发明名称 一种检测地层中裂缝发育程度的方法
摘要 本发明提供一种检测地层中裂缝发育程度的方法,所述方法利用横波在裂缝发育地层中传播时发生横波分裂现象,导致井下三分量检波器的两个水平分量接收到的横波初至波的波形随着地层中裂缝的发育程度变化而变化,通过计算两个水平分量上的横波初至波互相关误差系数来定量检测地层裂缝的发育程度,为降低勘探开发的风险,提高油气藏开采效率提供参考数据。
申请公布号 CN103645506A 申请公布日期 2014.03.19
申请号 CN201310553812.4 申请日期 2013.11.08
申请人 中国石油大学(北京) 发明人 李向阳;蔡志光;撒振宇;宋建军;陈海峰
分类号 G01V1/30(2006.01)I;G01V1/40(2006.01)I 主分类号 G01V1/30(2006.01)I
代理机构 北京汇泽知识产权代理有限公司 11228 代理人 张瑾
主权项 1.一种检测地层中裂缝发育程度的方法,所述方法包括以下步骤:(1)根据所欲观测的井所在的现场地形,按照如下要求设计三维三分量VSP观测方案:从井口所在位置O点指向可控震源所在位置S点的观测方位角度取值α∈{α<sub>j,j=1,2,......,M</sub>},α<sub>j</sub>∈[0,360),α的M个取值将圆周均匀分割为M等份,根据对观测精度的要求,M取值从集合{6、8、10、12、16、18、36}中选取一个合适的数值;井下三分量检波器所放置的位置G与井口位置O的距离为z<sub>0</sub>米,z<sub>0</sub>的取值保证井下三分量检波器的正常工作并尽量接近要观测的地层;可控震源所在位置S与井口位置O之间的距离为x<sub>0</sub>米,x<sub>0</sub>的取值满足当激发可控震源时,井中三分量检波器能够接收到较高信噪比的有效信号;设定井下三分量检波器的离散采样间隔dt与离散采样点数N,dt与N的取值保证记录到完整的可控震源引发的位置G处的质点振动;(2)对于α∈{α<sub>j,j=1,2,......,M</sub>}的每一个取值α<sub>j</sub>,将可控震源放置在与井口O距离为x<sub>0</sub>米方位角度为α<sub>j</sub>处;将井下三分量检波器放置于与井口位置O的距离为z<sub>0</sub>米的井下,并使井下三分量检波器的R分量指向的方位角度与α<sub>j</sub>保持一致;激发可控震源产生沿着α<sub>j</sub>方位角度偏振的横波,并同时以井下三分量检波器记录可控震源引发的井下三分量检波器所在位置G处的质点振动,当α取值为α=α<sub>j</sub>时井下三分量检波器两个水平分量R、T记录到的质点振动记为R<sub>j</sub>(t<sub>k</sub>),T<sub>j</sub>(t<sub>k</sub>),k=1,2,......,N;(3)对步骤(2)中得到的所有观测方位的数据记录{R<sub>j</sub>(t<sub>k</sub>),T<sub>j</sub>(t<sub>k</sub>)},j=1,2,......,M,k=1,2,......,N,编程计算直达波初至的到达时间{P<sub>j</sub>},j=1,2,......,M,编程的具体算法为:对于α=α<sub>j</sub>的每一个取值得到的井下三分量检波器的两个水平分量的数据记录R<sub>j</sub>(t<sub>k</sub>),T<sub>j</sub>(t<sub>k</sub>),k=1,2,......,N,扫描计算{r<sub>k</sub>=F<sub>k</sub>/E<sub>k</sub>,k=L+1,L+2,......,N-L}中的最大值所对应的k的取值,则得到P<sub>j</sub>=kdt;在公式{r<sub>k</sub>=F<sub>k</sub>/E<sub>k</sub>,k=L+1,L+2,......,N-L}中,<![CDATA[<math><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mrow><mo>(</mo><msubsup><mi>R</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>-</mo><mi>l</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>T</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>-</mo><mi>l</mi></mrow></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><![CDATA[<math><mrow><msub><mi>F</mi><mi>k</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mrow><mo>(</mo><msubsup><mi>R</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>+</mo><mi>l</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>T</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>+</mo><mi>l</mi></mrow></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>L的取值由可控震源激发的地震子波的延续时间t0决定,L=t0/dt,L<N;(4)根据{P<sub>j</sub>}的值,求取在所有观测方位角度观测的R分量上的横波初至与T分量上横波初至的互相关误差系数<![CDATA[<math><mrow><msub><mi>C</mi><mi>j</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mrow><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msub><mi>R</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>+</mo><mi>l</mi></mrow></msub><mo>)</mo></mrow><msub><mi>T</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>+</mo><mi>l</mi></mrow></msub><mo>)</mo></mrow><mo>|</mo></mrow><mrow><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msubsup><mi>R</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></msqrt><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msubsup><mi>T</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></msqrt></mrow></mfrac><mo>,</mo><mi>k</mi><mo>=</mo><msub><mi>P</mi><mi>j</mi></msub><mo>/</mo><mo>,</mo></mrow></math>]]></maths>L的取值同上所述,j=1,2,……,M;(5)计算所有C<sub>j</sub>,j=1,2,……,M的平均值<img file="FDA0000410992070000024.GIF" wi="371" he="149" />则sigma值的大小用于定量描述所观测地层中的裂缝发育程度,sigma的取值范围为sigma∈(0,1),当sigma值增大时,代表地层中的裂缝发育程度提高,当sigma值减小时代表地层中的裂缝发育程度降低;(6)流程结束。
地址 102249 北京市昌平区府学路18号