CN103645506A - 一种检测地层中裂缝发育程度的方法 - Google Patents
一种检测地层中裂缝发育程度的方法 Download PDFInfo
- Publication number
- CN103645506A CN103645506A CN201310553812.4A CN201310553812A CN103645506A CN 103645506 A CN103645506 A CN 103645506A CN 201310553812 A CN201310553812 A CN 201310553812A CN 103645506 A CN103645506 A CN 103645506A
- Authority
- CN
- China
- Prior art keywords
- value
- component
- hole
- sigma
- vibroseis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000011161 development Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000015572 biosynthetic process Effects 0.000 title abstract description 16
- 239000002245 particle Substances 0.000 claims description 11
- 241000269793 Cryothenia peninsulae Species 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 5
- 230000010287 polarization Effects 0.000 claims description 5
- 235000007164 Oryza sativa Nutrition 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 235000009566 rice Nutrition 0.000 claims description 3
- 240000007594 Oryza sativa Species 0.000 claims 1
- 238000000605 extraction Methods 0.000 abstract 1
- 208000010392 Bone Fractures Diseases 0.000 description 24
- 206010017076 Fracture Diseases 0.000 description 24
- 238000010586 diagram Methods 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 125000001183 hydrocarbyl group Chemical group 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种检测地层中裂缝发育程度的方法,所述方法利用横波在裂缝发育地层中传播时发生横波***现象,导致井下三分量检波器的两个水平分量接收到的横波初至波的波形随着地层中裂缝的发育程度变化而变化,通过计算两个水平分量上的横波初至波互相关误差系数来定量检测地层裂缝的发育程度,为降低勘探开发的风险,提高油气藏开采效率提供参考数据。
Description
技术领域
本发明涉及地震勘探技术领域,是一种通过对三维三分量VSP中井下三分量检波器的两个水平分量记录的地震波场的横波初至能量的比值进行计算,从而对地层中裂缝发育程度进行定量检测的分析方法。本发明主要应用于地震勘探领域,为裂缝型油气藏勘探与开发提供地层裂缝发育程度信息,为降低勘探开发的风险,提高油气藏开采效率提供参考数据。
背景技术
在三维三分量VSP地震勘探数据中,当横波可控震源在井旁激发时,震源所产生的横波在地层中传播,穿透地层到达井中的三分量检波器所在位置,并在该点引起质点振动,该点的振动被放置于井中的三分量检波器记录下来。如果地层中裂缝发育,可控震源所激发的横波在穿透地层传播时会发生横波***现象,即横波会***为沿着裂缝走向偏振的快横波和垂直裂缝走向偏振的慢横波。通过计算可控震源激发、井中三分量检波器接收的三维三分量VSP地震数据中井下三分量检波器的水平分量的横波初至波互相关误差系数的大小,可以定量检测地层中裂缝发育的程度。
现有技术手段是通过地层的地质露头进行观测、或者对钻井中提取的岩心样本在实验室中观测、或者对测井资料进行分析来获得地层裂缝发育程度的相关信息。由于上述方法资料采集的局限性,只能够获得地层裂缝的局部信息。
现有技术手段也可以通过宽方位角纵波地震勘探获得纵波地震传播的方位各向异性属性来检测地层裂缝的走向方位角度信息,但是由于地震纵波对地层裂缝的发育情况不够敏感,因此通过纵波地震方位各向异性属性分析来对储层裂缝的发育特征只能够做到定性的或者是误差较大的半定量描述。
总之,现有的技术手段由于受到技术方法本身的限制,对地层裂缝发育程度的检测能力不能满足油田的生产需要。
发明内容
本发明通过计算三维三分量VSP地震勘探中井下三分量检波器水平分量的横波初至的互相关误差系数来定量检测地层裂缝发育程度的分析方法与流程。
当地层中发育有裂缝时,横波在地层中传播时将发生横波***现象,横波***导致井下三分量检波器的两个水平分量上接收到的横波初至波的能量不同,地层中裂缝发育程度越高,井下三分量检波器的两个水平分量上接收到的横波初至波的能量大小就越接近。通过计算井下三分量检波器的两个水平分量上所记录的地震波场中横波初至波互相关误差系数,可以定量检测地层中裂缝的发育程度。
图1是当横波穿过呈平行排列的垂直裂缝介质时发生横波***的原理示意图。当由横波可控震源激发的入射横波S穿过裂缝发育的地层时,发生了横波***,入射横波S分解为快横波S1和慢横波S2,其中快横波S1的偏振方向平行于裂缝走向,慢横波S2的偏振方向垂直于裂缝走向。
为了便于分析横波***现象,我们将三分量检波器的径向方向记录的地震波场记为R(t)、切向分量记录的地震波场记为T(t),三分量检波器的两个水平分量同时记录下了快横波S1(t)和慢横波S2(t)在观测点引起的质点振动,三分量检波器的两个水平分量R(t)与T(t)可以写成如下形式:
上式中,NR(t)与NT(t)分别表示三分量检波器的径向分量与横向分量记录的环境噪声;θ是地层中裂缝走向方位角度β与三分量检波器的径向分量指向的方位角度α的夹角,见图2中所示。
本发明的技术方案为:
一种检测地层中裂缝发育程度的方法,所述方法包括以下步骤:
(1)根据所欲观测的井所在的现场地形,按照如下要求设计三维三分量VSP观测方案:从井口所在位置O点指向可控震源所在位置S点的观测方位角度取值α∈{αj,j=1,2,......,M},αj∈[0,360),α的M个取值将圆周均匀分割为M等份,根据对观测精度的要求,M取值从集合{6、8、10、12、16、18、36}中选取一个合适的数值;井下三分量检波器所放置的位置G与井口位置O的距离为z0米,z0的取值保证井下三分量检波器的正常工作并尽量接近要观测的地层;可控震源所在位置S与井口位置O之间的距离为x0米,x0的取值满足当激发可控震源时,井中三分量检波器能够接收到较高信噪比的有效信号;设定井下三分量检波器的离散采样间隔dt与离散采样点数N,dt与N的取值保证记录到完整的可控震源引发的位置G处的质点振动;
(2)对于α∈{αj,j=1,2,......,M}的每一个取值αj,将可控震源放置在与井口O距离为x0米方位角度为αj处;将井下三分量检波器放置于与井口位置O的距离为z0米的井下,并使井下三分量检波器的R分量指向的方位角度与αj保持一致;激发可控震源产生沿着αj方位角度偏振的横波,并同时以井下三分量检波器记录可控震源引发的井下三分量检波器所在位置G处的质点振动,当α取值为α=αj时井下三分量检波器两个水平分量R、T记录到的质点振动记为Rj(tk),Tj(tk),k=1,2,......,N;
(3)对步骤(2)中得到的所有观测方位的数据记录{Rj(tk),Tj(tk)},j=1,2,......,M,k=1,2,......,N,编程计算直达波初至的到达时间{Pj},j=1,2,......,M,
编程的具体算法为:
对于α=αj的每一个取值得到的井下三分量检波器的两个水平分量的数据记录Rj(tk),Tj(tk),k=1,2,......,N,扫描计算{rk=Fk/Ek,k=L+1,L+2,......,N-L}中的最大值所对应的k的取值,则得到Pj=kdt;在公式{rk=Fk/Ek,k=L+1,L+2,......,N-L}中, L的取值由可控震源激发的地震子波的延续时间t0决定,L=t0/dt,L<N;
(4)根据{Pj}的值,求取在所有观测方位角度观测的R分量上的横波初至与T分量上横波初至的互相关误差系数 L的取值同上所述,j=1,2,……,M;
(5)计算所有Cj,j=1,2,……,M的平均值则sigma值的大小用于定量描述所观测地层中的裂缝发育程度,sigma的取值范围为sigma∈(0,1),当sigma值增大时,代表地层中的裂缝发育程度提高,当sigma值减小时代表地层中的裂缝发育程度降低;
(6)流程结束。
本发明的技术效果:利用横波在裂缝发育地层中传播时发生横波***现象,导致井下三分量检波器的两个水平分量接收到的横波初至波能量比值随着地层中裂缝的发育程度变化而变化,通过计算两个水平分量上的横波初至波互相关误差系数来定量检测地层裂缝的发育程度。
附图说明
图1为横波穿过呈平行排列的垂直裂缝介质时发生横波***的示意图。
图2为三分量检波器径向分量指向的方位角度α、裂缝走向β与θ的关系图。
图3为三维三分量VSP地震观测方法示意图。
图4为α取值(M=16)及其与井下检波器水平分量指向方位关系示意图。
图5与图6分别是理论合成的三维三分量VSP观测数据的水平分量R与T图。模型设计的裂缝方位为NE60°,地层速度100m/s,x=500m,z=500m。观测方位角度α的取值为:NE10°,NE20°,……,NE360°。井下三分量检波器的方向采用右手坐标系,径向分量与观测方位角度α保持一致,切向分量垂直于观测方位角度α,位于径向分量顺时针旋转90°的方位。合成数据地震子波采用的是40Hz的零相位雷克子波。
图7和图8是对图5和图6中的合成地震记录水平分量进行横波初至时间拾取的结果图。
图9是对图7和图8中的数据计算每个观测方位RT分量互相关误差系数数值图,从图9中可知所有观测方位RT分量互相关误差系数的平均值为sigma=0.62,表明地层中裂缝的发育程度为0.62。
具体实施方式
一种检测地层中裂缝发育程度的方法,所述方法包括以下步骤:
(1)根据所欲观测的井所在的现场地形,按照如下要求设计三维三分量VSP观测方案:参考图2所示,从井口所在位置O点指向可控震源(炮点)所在位置S点的观测方位角度取值α∈{αj,j=1,2,......,M},αj∈[0,360),α的M个取值将圆周均匀分割为M等份,根据对观测精度的要求,M取值可从集合{6、8、10、12、16、18、36}中选取一个合适的数值;井下三分量检波器所放置的位置G与井口位置O的距离为z0米,z0的取值应能够保证井下三分量检波器的正常工作并尽量接近要观测的地层;可控震源所在位置S与井口位置O之间的距离为x0米,x0的取值应满足当激发可控震源时,井中三分量检波器能够接收到较高信噪比的有效信号;设定井下三分量检波器的离散采样间隔dt与离散采样点数N,dt与N的取值应能保证记录到完整的可控震源引发的G点处的质点振动。
(2)对于α∈{αj,j=1,2,......,M}的每一个取值αj,将可控震源放置在与井口O距离为x0米方位角度为αj处;将井下三分量检波器放置于与井口位置O的距离为z0米的井下,并使井下三分量检波器的R分量指向的方位角度与αj保持一致(参考图3所示);激发可控震源产生沿着αj方位角度偏振的横波,并同时以井下三分量检波器记录可控震源引发的井下三分量检波器所在位置G处的质点振动。当α取值为α=αj时井下三分量检波器两个水平分量R、T记录到的质点振动记为Rj(tk),Tj(tk),k=1,2,......,N。
(3)对(2)中得到的所有观测方位的数据记录{Rj(tk),Tj(tk)},j=1,2,......,M,k=1,2,......,N,编程计算直达波初至的到达时间{Pj},j=1,2,......,M。
编程的具体算法为:
对于α=αj的每一个取值得到的井下三分量检波器的两个水平分量的数据记录Rj(tk),Tj(tk),k=1,2,......,N,扫描计算{rk=Fk/Ek,k=L+1,L+2,......,N-L}中的最大值所对应的k的取值,则得到Pj=kdt;在公式{rk=Fk/Ek,k=L+1,L+2,......,N-L}中, L的取值由可控震源激发的地震子波的延续时间t0决定,L=t0/dt,L<N;
(4)根据{Pj}的值,求取在所有观测方位角度观测的R分量上的横波初至与T分量上横波初至的互相关误差系数 L的取值同上所述,j=1,2,……,M;
(5)计算所有Cj,j=1,2,……,M的平均值则sigma值的大小可用于定量描述所观测地层中的裂缝发育程度,sigma的取值范围为sigma∈(0,1),当sigma值增大时,代表地层中的裂缝发育程度提高,当sigma值减小时降低;
(6)流程结束。
本发明利用横波在裂缝发育地层中传播时发生横波***现象,导致井下三分量检波器的两个水平分量接收到的横波初至波能量比值随着地层中裂缝的发育程度变化而变化,通过计算两个水平分量上的横波初至波互相关误差系数来定量检测地层裂缝的发育程度。
Claims (1)
1.一种检测地层中裂缝发育程度的方法,所述方法包括以下步骤:
(1)根据所欲观测的井所在的现场地形,按照如下要求设计三维三分量VSP观测方案:从井口所在位置O点指向可控震源所在位置S点的观测方位角度取值α∈{αj,j=1,2,......,M},αj∈[0,360),α的M个取值将圆周均匀分割为M等份,根据对观测精度的要求,M取值从集合{6、8、10、12、16、18、36}中选取一个合适的数值;井下三分量检波器所放置的位置G与井口位置O的距离为z0米,z0的取值保证井下三分量检波器的正常工作并尽量接近要观测的地层;可控震源所在位置S与井口位置O之间的距离为x0米,x0的取值满足当激发可控震源时,井中三分量检波器能够接收到较高信噪比的有效信号;设定井下三分量检波器的离散采样间隔dt与离散采样点数N,dt与N的取值保证记录到完整的可控震源引发的位置G处的质点振动;
(2)对于α∈{αj,j=1,2,......,M}的每一个取值αj,将可控震源放置在与井口O距离为x0米方位角度为αj处;将井下三分量检波器放置于与井口位置O的距离为z0米的井下,并使井下三分量检波器的R分量指向的方位角度与αj保持一致;激发可控震源产生沿着αj方位角度偏振的横波,并同时以井下三分量检波器记录可控震源引发的井下三分量检波器所在位置G处的质点振动,当α取值为α=αj时井下三分量检波器两个水平分量R、T记录到的质点振动记为Rj(tk),Tj(tk),k=1,2,......,N;
(3)对步骤(2)中得到的所有观测方位的数据记录{Rj(tk),Tj(tk)},j=1,2,......,M,k=1,2,......,N,编程计算直达波初至的到达时间{Pj},j=1,2,......,M,
编程的具体算法为:
对于α=αj的每一个取值得到的井下三分量检波器的两个水平分量的数据记录Rj(tk),Tj(tk),k=1,2,......,N,扫描计算{rk=Fk/Ek,k=L+1,L+2,......,N-L}中的最大值所对应的k的取值,则得到Pj=kdt;在公式{rk=Fk/Ek,k=L+1,L+2,......,N-L}中, L的取值由可控震源激发的地震子波的延续时间t0决定,L=t0/dt,L<N;
(4)根据{Pj}的值,求取在所有观测方位角度观测的R分量上的横波初至与T分量上横波初至的互相关误差系数 L的取值同上所述,j=1,2,……,M;
(5)计算所有Cj,j=1,2,……,M的平均值则sigma值的大小用于定量描述所观测地层中的裂缝发育程度,sigma的取值范围为sigma∈(0,1),当sigma值增大时,代表地层中的裂缝发育程度提高,当sigma值减小时代表地层中的裂缝发育程度降低;
(6)流程结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310553812.4A CN103645506B (zh) | 2013-11-08 | 2013-11-08 | 一种检测地层中裂缝发育程度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310553812.4A CN103645506B (zh) | 2013-11-08 | 2013-11-08 | 一种检测地层中裂缝发育程度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103645506A true CN103645506A (zh) | 2014-03-19 |
CN103645506B CN103645506B (zh) | 2015-09-02 |
Family
ID=50250754
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310553812.4A Expired - Fee Related CN103645506B (zh) | 2013-11-08 | 2013-11-08 | 一种检测地层中裂缝发育程度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103645506B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104297786A (zh) * | 2014-10-17 | 2015-01-21 | 中国石油天然气集团公司 | 检测地层裂缝走向方位角度的方法和装置 |
CN105572739A (zh) * | 2016-02-19 | 2016-05-11 | 中国石油集团川庆钻探工程有限公司 | 碳酸盐岩孔洞裂缝性发育特征判断方法 |
CN105572732A (zh) * | 2014-10-08 | 2016-05-11 | 中国石油化工股份有限公司 | 一种属性变化率逐次晋级的裂缝发育带检测方法 |
CN107992690A (zh) * | 2017-12-06 | 2018-05-04 | 西南石油大学 | 一种诱导应力干扰下多裂缝扩展形态均衡程度的评价方法 |
CN108277728A (zh) * | 2018-01-04 | 2018-07-13 | 中国恩菲工程技术有限公司 | 道路面层病害检测方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102053266A (zh) * | 2009-11-09 | 2011-05-11 | 中国石油化工股份有限公司 | 地下裂缝预测方法 |
CN102096097A (zh) * | 2009-12-14 | 2011-06-15 | 武侯区巅峰机电科技研发中心 | 快慢横波***的裂缝方位角测试方法 |
CN102879800A (zh) * | 2011-07-15 | 2013-01-16 | 中国石油天然气集团公司 | 一种横波***裂缝检测的方法 |
-
2013
- 2013-11-08 CN CN201310553812.4A patent/CN103645506B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102053266A (zh) * | 2009-11-09 | 2011-05-11 | 中国石油化工股份有限公司 | 地下裂缝预测方法 |
CN102096097A (zh) * | 2009-12-14 | 2011-06-15 | 武侯区巅峰机电科技研发中心 | 快慢横波***的裂缝方位角测试方法 |
CN102879800A (zh) * | 2011-07-15 | 2013-01-16 | 中国石油天然气集团公司 | 一种横波***裂缝检测的方法 |
Non-Patent Citations (2)
Title |
---|
张学映: "多层裂缝介质转换波横波***分析技术及应用", 《新疆石油地质》 * |
徐天吉: "横波***裂缝检测方法在川西深层气藏中的应用", 《地球物理学进展》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572732A (zh) * | 2014-10-08 | 2016-05-11 | 中国石油化工股份有限公司 | 一种属性变化率逐次晋级的裂缝发育带检测方法 |
CN104297786A (zh) * | 2014-10-17 | 2015-01-21 | 中国石油天然气集团公司 | 检测地层裂缝走向方位角度的方法和装置 |
CN104297786B (zh) * | 2014-10-17 | 2017-06-13 | 中国石油天然气集团公司 | 检测地层裂缝走向方位角度的方法和装置 |
CN105572739A (zh) * | 2016-02-19 | 2016-05-11 | 中国石油集团川庆钻探工程有限公司 | 碳酸盐岩孔洞裂缝性发育特征判断方法 |
CN107992690A (zh) * | 2017-12-06 | 2018-05-04 | 西南石油大学 | 一种诱导应力干扰下多裂缝扩展形态均衡程度的评价方法 |
CN108277728A (zh) * | 2018-01-04 | 2018-07-13 | 中国恩菲工程技术有限公司 | 道路面层病害检测方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN103645506B (zh) | 2015-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103645505B (zh) | 一种定量检测地层中裂缝发育程度的方法 | |
Maxwell | Microseismic imaging of hydraulic fracturing: Improved engineering of unconventional shale reservoirs | |
CN102879800B (zh) | 一种横波***裂缝检测的方法 | |
CN102053266B (zh) | 地下裂缝预测方法 | |
CN102455439B (zh) | 基于克希霍夫积分法的绕射波场分离方法 | |
CN103645506B (zh) | 一种检测地层中裂缝发育程度的方法 | |
CN102053261A (zh) | 一种地震数据处理方法 | |
RU2539745C1 (ru) | Способ сейсмического мониторинга в процесса разработки месторождений углеводородов на акваториях | |
CN102053259A (zh) | 一种油气藏识别方法 | |
CN102053263A (zh) | 调查表层结构的方法 | |
CN116378648B (zh) | 一种基于随钻声波前视的近钻头地层探测方法及装置 | |
CA2961168A1 (en) | Integrating vertical seismic profile data for microseismic anisotropy velocity analysis | |
CN107179555B (zh) | 随钻地震钻头震源侧帮地质构造探测方法 | |
CN103869366A (zh) | 一种确定裂隙裂缝走向的方法及装置 | |
EP2917857B1 (en) | Fracture characterization from refraction travel time data | |
CN104484573A (zh) | 一种地层的刚性系数的计算新方法 | |
CN104297786B (zh) | 检测地层裂缝走向方位角度的方法和装置 | |
CN103645508B (zh) | 一种检测地层裂缝走向方位角度的方法 | |
CN103645510B (zh) | 一种定量检测地层裂缝走向方位角度的方法 | |
CN110780345A (zh) | 一种隧道超前地震探测地震数据的三维速度分析方法 | |
Hu et al. | Quantitative comparative analysis of monitoring effect of microseismic fracturing borehole and surface | |
Sinha et al. | Hydraulic fracture characterization using cased-hole sonic data | |
CN112946753B (zh) | 近地表地层结构的分析***及方法 | |
Howard | Geothermal resource and reservoir investigations of US Bureau of Reclamation leaseholds at East Mesa, Imperial Valley, California | |
RU2339978C1 (ru) | Способ скважинной сейсморазведки |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150902 |
|
CF01 | Termination of patent right due to non-payment of annual fee |