主权项 |
一种定向地震波延时参数优化设计方法,其特征在于,包括以下步骤:a、根据测区地质、地球物理资料,建立探区地下介质模型;b、震源点和观测排列按常规地震勘探方法设计,设震源点个数为N,编号顺序分别为1,2,...,N,其中N≥2,各震源点保持等间距d,检波器道数为n,道间距为D,震源函数为f(t),建立并求解弹性波方程,得到对应N个震源点的N个单炮记录;t为时间变量;c、求取目标反射层上覆介质的均方根速度v;d、设地震波主波束方向与法线夹角为θ,<img file="FDA0000934609230000011.GIF" wi="261" he="103" />其中延时时间τ为定向地震波形成方法的延时时间参数,再根据三角函数性质‑1≤sinθ≤1,则有<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mo>-</mo><mfrac><mi>d</mi><mi>v</mi></mfrac><mo>≤</mo><mi>τ</mi><mo>≤</mo><mfrac><mi>d</mi><mi>v</mi></mfrac></mrow>]]></math><img file="FDA0000934609230000012.GIF" wi="244" he="118" /></maths> (不等式1)其中延时参数τ取负值时代表信号超前,求解不等式1,得到不等式的解集t<sub>0</sub>≤τ≤t<sub>1</sub>;e、在[t<sub>0</sub>,t<sub>1</sub>]内,采用线性插值等间隔选取m个延时时间参数,记为τ<sub>1</sub>,τ<sub>2</sub>,...,τ<sub>i</sub>,...,τ<sub>m</sub>,其中τ<sub>1</sub>=t<sub>0</sub>,τ<sub>m</sub>=t<sub>1</sub>,m≥20,针对每个延时参数τ<sub>i</sub>(i∈[1,m]),将N个单炮记录按震源点编号分别延时0,2τ<sub>i</sub>,3τ<sub>i</sub>,……,(N‑1)τ<sub>i</sub>,然后叠加,得到对应延时参数τ<sub>i</sub>的定向地震记录H<sub>i</sub>;f、以H<sub>i</sub>中勘探目标反射波同相轴为中心,选取等长的反射波时窗,计算定向地震记录H<sub>i</sub>内各道目标反射信号能量E<sub>i,1</sub>,E<sub>i,2</sub>,...,E<sub>i,j</sub>,...,E<sub>i,n</sub>,其中j∈[1,n],根据公式<img file="FDA0000934609230000013.GIF" wi="291" he="207" />计算接收阵列内各道反射波能量密度p<sub>i,j</sub>;g、根据<img file="FDA0000934609230000014.GIF" wi="234" he="134" />计算H<sub>i</sub>内接收阵列接收的目标反射波总能量h<sub>i</sub>;根据[(τ<sub>1</sub>,h<sub>1</sub>),(τ<sub>2</sub>,h<sub>2</sub>),...,(τ<sub>i</sub>,h<sub>i</sub>),...,(τ<sub>m</sub>,h<sub>m</sub>)]绘制反射波能量随延时时间参数变化曲线,并从中找出前5个最大值点,设该五个最大值点对应的延时参数为:τ<sub>1</sub>',τ<sub>2</sub>',τ<sub>3</sub>',τ<sub>4</sub>',τ<sub>5</sub>';h、计算对应τ<sub>1</sub>',τ<sub>2</sub>',τ<sub>3</sub>',τ<sub>4</sub>',τ<sub>5</sub>'的5个定向地震记录的反射波能量密度方差,最小方差对应的延时时间即为定向地震波的最优延时时间参数。 |