发明名称 一种检波器组合参数计算方法
摘要 本发明涉及一种检波器组合参数计算方法。可根据表层调查的速度、试验接收子波的主频率、观测系统排列的长度、最浅目的层的深度等信息快速计算出最优的检波器组合方式。实现该方法首先需要建立各个组合参数与输出响应间的函数关系,该输出响应会因信号输入角度不同而不同,因此可以将有效信号区间的归一化响应累加值与噪音区间的响应累加值相除,作为目标函数,求取组合参数。具体包含三部分:1、建立子波主频、表层速度、入射角度、组合个数、检波器空间坐标、检波器加权系数与组合输出响应间的函数关系;2、建立优化目标函数;3、求取最优组合参数。本专利克服了当前检波器组合参数设计上的某些不足,寻找到了一种快速的组合参数定量化计算方法,同时具有良好地可操作性,可以为野外施工提供方案指导。
申请公布号 CN102830424B 申请公布日期 2016.02.17
申请号 CN201210328063.0 申请日期 2012.09.07
申请人 西南石油大学 发明人 刘福烈;徐峰;刘文革;敬龙江
分类号 G01V1/00(2006.01)I;G01V1/20(2006.01)I 主分类号 G01V1/00(2006.01)I
代理机构 代理人
主权项 一种检波器组合参数确定方法,依次包括以下步骤:①建立各种检波器组合参数与组合输出响应间的函数关系;②建立参数优化的目标函数;③对组合参数进行优化计算的程序实现;其特征在于步骤①中各种检波器组合参数与组合输出响应间有函数关系,表达为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&Phi;</mi><mo>=</mo><mi>F</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><msqrt><mrow><msup><mrow><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>c</mi><mi>i</mi></msub><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>&pi;f&Delta;t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>c</mi><mi>i</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>&pi;f&Delta;t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow>]]></math><img file="FDA0000800365990000011.GIF" wi="1589" he="217" /></maths>其中,F(f)为雷克子波频率域表达式,f<sub>0</sub>为其主频<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>2</mn><msqrt><mi>&pi;</mi></msqrt></mfrac><mfrac><msup><mi>f</mi><mn>2</mn></msup><mrow><msup><msub><mi>f</mi><mn>0</mn></msub><mn>3</mn></msup></mrow></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mi>f</mi><mn>2</mn></msup><mrow><msup><msub><mi>f</mi><mn>0</mn></msub><mn>2</mn></msup></mrow></mfrac></mrow></msup></mrow>]]></math><img file="FDA0000800365990000012.GIF" wi="652" he="263" /></maths>N为检波器组合个数,C<sub>i</sub>为检波器加权系数,△t<sub>i</sub>为各检波器到目标观测平面的时间差,表达为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&Delta;t</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>d</mi><mi>i</mi></msub><mi>v</mi></mfrac><mo>=</mo><mfrac><mrow><mo>|</mo><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>px</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>px</mi><mi>i</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>py</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>py</mi><mi>i</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>z</mi><mi>i</mi></msub><mo>-</mo><msub><mi>pz</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>pz</mi><mi>i</mi></msub><mo>|</mo></mrow><mrow><mi>v</mi><mi>r</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000800365990000013.GIF" wi="1694" he="198" /></maths>其中,d<sub>i</sub>为各检波器到目标观测平面的距离,v为检波器所在位置的表层速度,r为组合中心到观测平面的距离,px<sub>i</sub>,py<sub>i</sub>,pz<sub>i</sub>为过组合中心的观测平面法线在XYZ三个坐标轴上的投影长度,表达为:px<sub>i</sub>=rsinθcosβpy<sub>i</sub>=rsinθsinβpz<sub>i</sub>=rcosθ其中,θ为子波入射倾角,取值范围‑90°~90°,β为子波入射水平角度,取值范围0°~360°。
地址 610500 四川省成都市新都区新都大道8号