用matlab程序实现沙泥岩指数预测.PDF

2016年第6期 科学管理 用matlab程序实现沙泥岩指数预测 姚姣凤　郑冲 河南省煤田地质局二队　河南　洛阳　471000 摘要：地震波的传播速度与地层的岩性、孔隙度、埋藏深度等因素密切相关。通过速度谱资料提取出层速度，利用 matlab编程，可以实现砂泥岩质量分数的自动预测。实际应用表明，此方法预测的砂泥岩质量分数结果与钻井资料一致。 关键词：砂泥岩体积物理模型　砂泥岩质量分数　matlab程序实现　岩性量板 matlab 1　理论基础 clear all 471000 file_fid=fopen（2.txt，r）； 砂、泥岩体积物理模型 为了找出岩石的整体速度与其中砂泥岩成分之间的关 TXT_FILE_HEAD=textscan（file_fid，%s %s %s，1）； matlab depth=TXT_FILE_HEAD{1}； 系，需要建立砂泥岩体积物理模型。沿地震波传播方向截 matlab 取一块边长为Z的立方体砂泥岩块，砂泥岩以互层形式分 index=TXT_FILE_HEAD{2}； 1 布其中，见图1a。若将岩块中的砂岩和泥岩分别集中，可 vstr=TXT_FILE_HEAD{3}； 用其等效体积来表示，见图1b。因此，若地震波以速度V ori_data=fscanf（file_fid，%f %f %f，[3 inf]）； Z 通过岩块，则可以认为，地震波通过Z的时间t等于通过等 ori_data=ori_data； a b v Z t fclose（file_fid）； 效砂泥岩段的时间。由以上理论可得到下式： t= Z /V= Zs/Vs+ Zh/Vh （1） [raw_number cloum_number]=size（ori_data） t= Z /V= Zs/Vs+ Zh/Vh 1 式中：V为地震波通过等效砂泥岩体的速度，m/s； niyan=[]； v m/s Vs Vh m/s shayan=[]； Vs、Vh分别为纯砂岩和纯泥岩的波速度，m/s。 Zs/Z= Ps for i=1：raw_number 令：Zs/Z= Ps 砂岩含量则： Zh/Z = 1-Ps 1/V = Ps/Vs + (1-Ps)/Vh 2 if ori_data（i，2）==1||ori_data（i，2）==4 Zh/Z = 1-Ps →1/V = Ps/Vs +（1-Ps）/Vh （2） 根据此式就可计算砂岩百分含量