发明名称 用于地震勘探的照明度计算方法
摘要 针对现有技术中存在的用于地震勘探的照明度计算方法计算量大、计算效率低等问题,本发明提出了一种所述用于地震勘探的照明度计算方法,包括:通过布置多个检波点来检测各炮点的位置x<sub>s</sub>、炮点被激发后所产生的频率域波场的频率角速度ω、炮点被激发后的地表速度v;将多个球面波源转换成一个平面波源;对平面波入射角度θ在0°至180°范围内取多个值来对转换后的平面波源进行多角度入射,从而在多个方向计算用于表示照明度的S<sub>syn</sub>(ω,x<sub>s</sub>,p)。根据本发明示例性实施例的用于地震勘探的照明度计算方法,在正确地模拟球面波点震源的激发过程的同时,可显著降低计算量,减少计算时间,从而可应用于实际的野外地震勘探作业。
申请公布号 CN102608657B 申请公布日期 2014.12.17
申请号 CN201110424517.X 申请日期 2011.12.16
申请人 中国石油集团川庆钻探工程有限公司地球物理勘探公司 发明人 李亚林;何光明;胡善政;陈爱萍;敬龙江;龙资强;耿春;巫骏
分类号 G01V1/30(2006.01)I 主分类号 G01V1/30(2006.01)I
代理机构 北京铭硕知识产权代理有限公司 11286 代理人 王青芝
主权项 一种用于地震勘探的照明度计算方法,其特征在于,所述用于地震勘探的照明度计算方法包括以下步骤:通过布置多个检波点来检测各炮点的位置x<sub>s</sub>、炮点被激发后所产生的频率域波场的频率角速度ω、炮点被激发后的频率波场的地表速度v;根据以下公式将多个球面波源转换成一个平面波源:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>S</mi><mi>syn</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><msub><mi>x</mi><mi>s</mi></msub><mo>=</mo><msub><mi>s</mi><mi>start</mi></msub></mrow><msub><mi>s</mi><mi>end</mi></msub></munderover><msup><mi>e</mi><mrow><mo>-</mo><mi>i&omega;p</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>s</mi></msub><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></msup><mi>&psi;</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>x</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000524670690000011.GIF" wi="1185" he="224" /></maths>其中,S<sub>syn</sub>(ω,x<sub>0</sub>,p)为转换后的频率域波场的平面波震源,其值用于表示照明度,ψ(ω,x<sub>s</sub>)表示球面波点震源,x0为震源位置参考点,s<sub>start</sub>和s<sub>end</sub>分别表示炮点的开始位置和结束位置,p为射线参数,其中,平面波入射角度记为θ,射线参数p为地表速度v/cosθ,其中,<img file="FDA0000524670690000012.GIF" wi="680" he="223" />ω为炮点被激发后所产生的频率域波场的频率角速度,ω<sub>M</sub>为用户设定的特定频段内最大频率所对应的频率角速度;对平面波入射角度θ在0°至180°范围内取多个值来对转换后的平面波源进行多角度入射,从而在多个方向计算用于表示照明度的S<sub>syn</sub>(ω,x<sub>s</sub>,p),其中,S<sub>syn</sub>(ω,x<sub>s</sub>,p)=|S<sub>syn</sub>(ω,x<sub>0</sub>,p)|<sup>2</sup>。
地址 610213 四川省成都市双流县华阳镇华阳大道一段一号1-4