二维几何法上、下底界面的几何法.docxVIP

二维几何法上、下底界面的几何法

地震资料的解释是地震勘探过程的最终环节，也是地震地质专家的一项重要任务。时间变化方法是地震数据解释的重要因素。正确的时间-深度转换方法是提高地震解释精度和地质效率的必要条件。

以往将水平叠加地震时间剖面解释为地震地质剖面,在一小段剖面中将新生界、基岩到煤层视为二层均匀介质,且取新生界平均速度为vQ=1800m/s左右,基岩到煤层平均速度为vP=3600m/s左右.按分层作图法:在水平叠加时间剖面上的某相邻桩号中,分别读取两桩号处新生界底界面反射时间tQ,煤层底板反射时间tm,再以桩号为圆心,以tQ×vQ/2为半径画弧,这两段弧的公切线之切点即为两桩号在新生界底界的反射点;再以这两个反射点为圆心,以(tm-tQ)×vP/2为半径画弧,这两段弧的公切线之切点即为两桩号在煤层底板上的反射点.这个做法不符合射线原理,因为在煤层底板上反射点的射线不能反射回原发射桩号处,这是一个假反射点;同时,手工作图慢且不准.根据我队以前的科研成果,假设煤层倾角αm≤45°,试用新生界底界倾角αQ=-25°～+25°作为入射角θ给上一个终止条件迭代扫描.这个方法程序冗长、运算量大、速度慢,且当煤层倾角较陡时,会形成死循环或程序运行不稳定.还有假设新生界底界与煤层倾向一致或相反,从而给出在新生界底界的入射角θ和煤层倾角αm的特例公式,均无法正确适应各种情况.因此,作者提出一种二层介质的时-深转换模型.

1原理

1.1新生界底界面倾角公式

如图1,AB为基准面的剖面线上的两桩号,已知A桩号到新生界底界的时间为tQA,到煤层底板的时间为tmA;B桩号到新生界底界的时间为tQB,到煤层底板的时间为tmB;新生界的速度为vQ,基岩到煤层的速度为vP;平面波在这二层均匀介质中传播,射线互相平行.

如图1,新生界底界倾向与煤层倾向一致且上倾,只是煤层倾角更陡的情况.C,D分别为过A,B两桩号的射线在新生界底界面上的反射点,C′,D′分别为过A,B两桩号的射线在新生界底界面上的透射点,E,F分别为过A,B两桩号的射线在煤层底板上的反射点.过B作AC的垂线,垂足为A′,过B作AC′的垂线,垂足为A″,显然A′C与BD用的时间相等,则sinαQ=AA′AB=(tQB?tQA)×vQ/2AB,αQsinαQ=AA′AB=(tQB-tQA)×vQ/2AB,αQ即为新生界底界倾角.射线与等时面正交,则A″,B在等时面上,E,F也在等时面上,故A″E与BF等时.则sinω=(tmB?tmA)×vQ/2AB?BHsinω=(tmB-tmA)×vQ/2AB?BΗ为铅垂线,据两锐角的两对应边垂直,则两锐角相等,得ω=∠ABA″=∠HBD′和αQ=∠ABA′=∠HBD以及γ=∠QFE,又∠DBD′=θ(入射角),故θ=ω-αQ.过F作水平线GF,作CD的平行线QF,αm=∠GFE(煤层倾角),αQ=∠GFQ(新生界底界之倾角),则αm=αQ+γ.

总之用前面公式所求得ω,αQ均为负值,ω的绝对值大于αQ的绝对值;而几何证明过程中ω,αQ均视为正角,则θ角为正角,γ也为正角.因而可将证明结果用如下形式表示

θ=αQ?ω,(1)αm=αQ?γ.(2)θ=αQ-ω,(1)αm=αQ-γ.(2)

对于煤层倾向与新生界底界倾向的其它关系,C′E⊥EF和D′F⊥EF的充要条件是:斜率互为负倒数即k1k2+1=0(k1,k2分别为两直线C′E和EF之斜率),据这一原则、用通式(1),(2)、模型实例(图2).Qbasic语言编程验算结果如表1.

1.2ce东方的射线路径a-c-e,tqa为唯一的介质速度,tqa为6.

图2为当两介质层倾角不同时反射波射线与界面的关系模型实例.A,B坐标分别为A(300,50),B(400,50).上、下层介质速度的二分之一分别为:vQ=900,vm=1800.上层上倾时:tQA=80,tQB=35,tmA=160,tmB=160-j×5(j=-2,0,13)上层下倾时:tQA=35,tQB=80,tmA=80,tmB=100+j×5(j=-2,2,10).

从图1可以看出,C′E⊥EF即A到E的射线路径A-C′-E是唯一的,故所用数学模型正确.

2新生界底界面

在已知A,B桩号(AB),A桩号到新生界底界的时间为tQA,到煤层底板的时间为tmA;B桩号到新生界底界的时间为tQB,到煤层底板的时间为tmB;新生界的速度为vQ,基岩到煤层的速度为vm,视为平面波在均匀二层介质中传播.

sinαQ=AA′AB=(tQB?tQA)×vQ/2AB?sinω=(tmB?