发明名称 基于地震记录的震源激发信号一致性评价方法
摘要 本发明涉及一种基于地震记录的震源激发信号一致性评价方法,首先提取震源地震记录中近偏移距数据中的直达波,确定不同道数据直达波信号的到时差,然后进行时差校正,将多道直达波信号进行叠加,再通过分析不同记录中的直达波信号的相似系数、相对时延以及能量差异,从而实现对不同震源的激发效果的一致性评价。解决了震源激发信号一致性的检测问题。避免了因各震源点处震源激发信号不一致而引起的地下构造探测结果较大误差的出现。不仅适用于可控震源也适用于炸药震源的一致性评价,不仅能用于单可控震源、组合震源,还可用于相控震源现场质量监控及震源维护,为地震数据处理过程对震源进行一致补偿提供必要的参考信息。
申请公布号 CN104536046A 申请公布日期 2015.04.22
申请号 CN201410827020.6 申请日期 2014.12.26
申请人 吉林大学 发明人 姜弢;贾海青;马瑞琪;徐学纯;林君
分类号 G01V1/30(2006.01)I;G01V1/36(2006.01)I 主分类号 G01V1/30(2006.01)I
代理机构 长春吉大专利代理有限责任公司 22201 代理人 王立文
主权项 一种基于地震记录的震源激发信号一致性评价的方法,包括以下步骤:a、输入采集系统一致的两个地震记录,记为U<sub>1</sub>(t,x)、U<sub>2</sub>(t,x),其中,t为时间,x为正整数,表示地震记录中各道数据的道号。记U<sub>i</sub>(t,x)=(u<sub>i</sub>(t,1),u<sub>i</sub>(t,2),…,u<sub>i</sub>(t,n)),i=1,2,u<sub>i</sub>(t,n)表示第i个地震记录第n道数据;b、在两炮地震记录中分别截取包含清晰直达波信号的近偏移距道地震记录<img file="FDA0000645454570000011.GIF" wi="364" he="72" />对应的时间范围为0≤t≤t<sub>0</sub>,对应的地震道范围为x<sub>min</sub>≤x≤x<sub>max</sub>,即<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>U</mi><mi>i</mi><mi>z</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mo>{</mo><msub><mi>U</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>|</mo><mn>0</mn><mo>&le;</mo><mi>t</mi><mo>&le;</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>x</mi><mi>min</mi></msub><mo>&le;</mo><mi>x</mi><mo>&le;</mo><msub><mi>x</mi><mi>max</mi></msub><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA0000645454570000012.GIF" wi="928" he="92" /></maths>c、记<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>U</mi><mi>i</mi><mi>z</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msubsup><mi>u</mi><mi>i</mi><mi>z</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><msub><mi>x</mi><mi>min</mi></msub><mo>)</mo></mrow><mo>,</mo><msubsup><mi>u</mi><mi>i</mi><mi>z</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><msub><mi>x</mi><mrow><mi>min</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>u</mi><mi>i</mi><mi>z</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><msub><mi>x</mi><mi>max</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000645454570000013.GIF" wi="978" he="75" /></maths>设d为<img file="FDA0000645454570000014.GIF" wi="160" he="76" />中最小偏移道道号,令<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>S</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>u</mi><mi>i</mi><mi>z</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>d</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000645454570000015.GIF" wi="334" he="65" /></maths>d、以S<sub>i</sub>(t)为参考道,对<img file="FDA0000645454570000016.GIF" wi="172" he="71" />的所有道数据作互相关运算,公式为<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><msub><mi>S</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CircleTimes;</mo><msubsup><mi>u</mi><mi>i</mi><mi>z</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000645454570000017.GIF" wi="453" he="78" /></maths>k=x<sub>min</sub>,x<sub>min</sub>+1,…,x<sub>max</sub>;e、记R<sub>i,k</sub>最大值对应时刻为τ<sub>i,k</sub>,将<img file="FDA0000645454570000018.GIF" wi="166" he="75" />时移‑τ<sub>i,k</sub>得到<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>u</mi><mi>i</mi><mi>&tau;</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close='' separators=''><mtable><mtr><mtd><msubsup><mi>u</mi><mi>i</mi><mi>z</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>+</mo><msub><mi>&tau;</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mn>0</mn><mo>&le;</mo><mi>t</mi><mo>&le;</mo><msub><mi>t</mi><mn>0</mn></msub><mo>-</mo><msub><mi>&tau;</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi></mrow></msub></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo></mtd><mtd><msub><mi>t</mi><mn>0</mn></msub><mo>-</mo><msub><mi>&tau;</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>&lt;</mo><mi>t</mi><mo>&lt;</mo><msub><mi>t</mi><mn>0</mn></msub></mtd></mtr></mtable><mo>,</mo></mfenced></mrow>]]></math><img file="FDA0000645454570000019.GIF" wi="852" he="180" /></maths>k=x<sub>min</sub>,x<sub>min</sub>+1,…,x<sub>max</sub>,输出对齐直达波同相轴后的记录<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>U</mi><mi>i</mi><mi>&tau;</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msubsup><mi>u</mi><mi>i</mi><mi>&tau;</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><msub><mi>x</mi><mi>min</mi></msub><mo>)</mo></mrow><mo>,</mo><msubsup><mi>u</mi><mi>i</mi><mi>&tau;</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><msub><mi>x</mi><mrow><mi>min</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>u</mi><mi>i</mi><mi>&tau;</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><msub><mi>x</mi><mi>max</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00006454545700000110.GIF" wi="977" he="70" /></maths>f、叠加<img file="FDA00006454545700000111.GIF" wi="175" he="77" />的所有道数据,计算得到对应的直达波信号<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>p</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><msub><mi>x</mi><mi>min</mi></msub></mrow><msub><mi>x</mi><mi>max</mi></msub></munderover><msubsup><mi>u</mi><mi>i</mi><mi>&tau;</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00006454545700000112.GIF" wi="433" he="146" /></maths>g、对p<sub>1</sub>(t)和p<sub>2</sub>(t)作互相关计算,<img file="FDA00006454545700000113.GIF" wi="397" he="75" />记R<sub>p</sub>的最大值对应的时刻为t<sub>p</sub>,则两地震记录中震源相对时延Δt=t<sub>p</sub>;h、求取p<sub>1</sub>(t)和p<sub>2</sub>(t)的相似系数,公式为<img file="FDA00006454545700000114.GIF" wi="686" he="167" />&lt;·&gt;表示内积运算;i、定义p<sub>1</sub>(t)和p<sub>2</sub>(t)的能量差异系数为ω=(W<sub>2</sub>‑W<sub>1</sub>)/W<sub>1</sub>*100%,其中,<img file="FDA00006454545700000115.GIF" wi="673" he="152" />j表示在一定采样率的情况下得到的离散时间序列;j、根据地震勘探中对震源的实际要求,确定基于地震记录的震源激发信号一致性判定准则为:|Δt|<0.14ms、γ≥0.9、|ω|≤5%,若以上三个条件同时满足,则判定两个地震记录的震源激发信号一致性强,反之,则两个地震记录的震源激发信号一致性较弱。
地址 130012 吉林省长春市前进大街2699号
您可能感兴趣的专利