CN105089663B - 一种大斜度/水平井地层电阻率各向异性校正方法 - Google Patents
一种大斜度/水平井地层电阻率各向异性校正方法 Download PDFInfo
- Publication number
- CN105089663B CN105089663B CN201510463526.8A CN201510463526A CN105089663B CN 105089663 B CN105089663 B CN 105089663B CN 201510463526 A CN201510463526 A CN 201510463526A CN 105089663 B CN105089663 B CN 105089663B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- resistivity
- anisotropy
- mfrac
- 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.)
- Active
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 48
- 238000012937 correction Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000004458 analytical method Methods 0.000 claims abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 14
- 239000011435 rock Substances 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 6
- 239000008398 formation water Substances 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 2
- 238000005755 formation reaction Methods 0.000 claims 15
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004457 water analysis Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010224 classification analysis Methods 0.000 description 1
- 230000007547 defect Effects 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
- 230000005251 gamma ray Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种大斜度/水平井地层电阻率各向异性校正方法。其包括如下步骤,步骤一,根据电阻率宏观各向异性和微观各向异性的各项参数以及之间的相互关系,建立水平和垂直电阻率解释图版;步骤二,根据对水平和垂直电阻率解释图版的分析,采用迭代算法确定地层电阻率宏观和微观总的各向异性系数,从而进行地层电阻率各向异性校正,计算出代表目的层真实电阻率的水平电阻率。有效解决了图版法无法反应储层的非均质性,不适合测井资料的逐点解释,自动校正法因校正图版数量少,且比较离散,并非完全符合实际情况,具有很大的区域性等不足之处的问题,十分适合于测井资料的逐点解释,能够及时服务和满足生产上的快速解释,具有很好的应用前景。
Description
技术领域
本发明涉及大斜度/水平井测井解释评价技术领域,具体为一种大斜度/水平井地层电阻率各向异性校正方法。
背景技术
随着石油工业的迅速发展,勘探开发的油气藏越来越复杂,水平井技术以其投资回收率高和适用范围广的优点得到了全世界的青睐。在未来的钻井类型中,大部分将会是水平井。然而水平井不同于直井,水平井测井响应受其井眼附近地层的几何形状、测量方位、异常侵入剖面、重力引起的仪器偏心和各向异性等影响严重,在进行水平井地层评价之前,需对原始测井数据进行必要的校正处理。在诸多的影响因素中,电阻率测井受地层电性各向异性影响最大。
目前,电阻率各向异性校正主要有图版法和计算机自动校正法。校正图版根据理论计算或实验结果制作而成,图版中所有参数都是代表一个层的各采样点的平均特性,其校正结果只是把大斜度/水平井目的层电阻率平均抬高或降低一个相同的比例,与测井曲线特有的连续性和非均质性相违背,只具有理论上定性解释和提供参考的价值,无法运用到实际测井资料的逐点解释中。自动校正法是在图版法的基础上,对校正图版中的曲线进行合理采样读值拟合,用标准化公式来代替大量的图版曲线,并编制校正处理程序库实现计算机自动校正,但由于各向异性图版校正数量少,且比较离散,并非完全符合实际情况,具有很大的区域性特点,校正方法适应性差。
发明内容
针对现有技术中存在的问题,本发明提供一种能够确定地层电阻率各向异性系数,得出能够代表目的层真实电阻率的水平电阻率的大斜度/水平井地层电阻率各向异性校正方法。
本发明是通过以下技术方案来实现:
一种大斜度/水平井地层电阻率各向异性校正方法,其特征在于,包括如下步骤,
步骤一,根据电阻率宏观各向异性和微观各向异性的各项参数以及之间的相互关系,建立水平和垂直电阻率解释图版;
步骤二,根据对水平和垂直电阻率解释图版的分析,采用迭代算法确定地层电阻率宏观和微观总的各向异性系数,从而进行地层电阻率各向异性校正,计算出代表目的层真实电阻率的水平电阻率。
优选的,步骤一中建立水平和垂直电阻率解释图版时,包括如下步骤,
1.1在大斜度/水平井中,当水平井井眼轨迹离地层层面的距离小于测井仪器的探测深度或与地层层面相交,测井仪器探测的是砂/泥岩互层,测量值为砂岩和泥岩电阻率的综合信息,地层电阻率各向异性表现为宏观各向异性;此时,
地层水平电阻率Rh表示如下,
地层垂直电阻率Rv表示如下,
电阻率各向异性系数λ表示如下,
其中,Rsd为在砂/泥岩模型中砂岩电阻率,单位为(Ω.m);Rsh为泥岩电阻率,单位为(Ω.m);Hsd为砂岩累计厚度,单位为(m);Hsh为泥岩累计厚度,单位为(m);
1.2在大斜度/水平井中,当水平井井眼轨迹到地层层界面的距离超出测 井仪器的探测深度时,视地层为无限厚的均匀各向异性介质,地层电阻率各向异性表现为微观各向异性;则得到,
地层视电阻率Ra表示如下,
其中,为各向异性地层的平均电阻率,单位为(Ω.m);θ为井斜角单位为(°);
1.3根据步骤1.1和1.2宏观各向异性和微观各向异性系数与水平电阻率、垂直电阻率和地层视电阻率等的关系,建立地层的水平和垂直电阻率解释图版。
进一步,步骤二中采用迭代算法进行地层电阻率各向异性校正时,包括如下步骤,
2.1利用自然伽马测井曲线得到地层中的泥质体积含量Vsh如下,
其中,ΔGR为自然电位相对值;GR为自然伽马测井值,GRmin为纯砂岩段的自然伽马值,GRmax为纯泥岩段的自然伽马值;GCUR为希尔奇系数,老地层取2,第三系新地层取3.7;
2.2结合区域的经验值确定泥岩电阻率Rsh,或用岩石完全含水时电阻率Ro近似替代,其中,a、m为阿尔奇公式岩电参数,Rw为地层水电阻率,φ为孔隙度;
2.3给定各向异性系数λ初始值λ=λk和给定误差,结合井斜角θ、步骤2.1计算的泥质体积含量Vsh和步骤2.2计算的泥岩电阻率Rsh,根据公式 (1)~(4)得到中间变量的水平电阻率Rh ·,Rh ·=f(λ,θ,Vsh,Rsh);
2.4将各向异性系数λ、井斜角θ和步骤2.3得到的水平电阻率Rh ·,结合公式(4)得到中间变量的地层视电阻率Ra ·;
2.5利用步骤2.4计算的Ra ·与实测地层视电阻率Ra比较,计算ε,ε=‖Ra-Ra ·‖;若ε大于给定误差,则λ=λk+1并重复步骤2.3~2.5,直到ε满足不大于给定误差的要求,输出电阻率宏观和微观总的各向异性系数λ,并计算出代表目的层真实电阻率的水平电阻率Rh。
再进一步,阿尔奇公式岩电参数a和m由岩电实验数据分析确定。
再进一步,地层水电阻率Rw由水分析资料确定。
再进一步,孔隙度φ由声波时差测井曲线计算确定。
再进一步,岩石完全含水时电阻率通过利用声波时差计算孔隙度而确定。
再进一步,当取希尔奇系数GCUR时,老地层取2,第三系新地层取3.7。
与现有技术相比,本发明具有以下有益的技术效果:
本发明通过对电阻率宏观各向异性和微观各向异性的分类分析,建立了水平和垂直电阻率解释图版,结合井斜角、自然伽马和声波时差测井曲线以及泥岩电阻率,泥岩电阻率由区域的经验值确定或用利用声波时差计算孔隙度后得到的岩石完全含水时电阻率近似替代,采用迭代算法确定地层电阻率宏观和微观总的各向异性系数,从而进行地层电阻率各向异性校正,得到能够代表目的层真实电阻率的水平电阻率,进而更好地应用于大斜度/水平井测井解释评价中。有效解决了图版法无法反应储层的非均质性,不适合测井资料的逐点解释,自动校正法因校正图版数量少,且比较离散,并非完全符合实际情况,具有很大的区域性等不足之处的问题,十分适合于测井资料的逐点解释,能够及时服务和满足生产上的快速解释,具有很好的应用前景。
附图说明
图1为本发明实施例提供的一种大斜度/水平井地层电阻率各向异性校正方法的流程图。
图2为本发明实施例所述的水平和垂直电阻率解释图版。
图3为本发明实施例所述的大斜度/水平井地层电阻率各向异性校正处理解释成果图。
具体实施方式
下面结合具体的实施例对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。
本发明中,地层电阻率各向异性可分为宏观各向异性和微观各向异性,描述地层电阻率各向异性的参数主要有,水平电阻率Rh,垂直电阻率Rv以及电阻率各向异性系数λ。在大斜度井/水平井中,仪器测量的地层电阻率Ra是水平电阻率Rh和垂直电阻率Rv的综合反映。
当水平井井眼轨迹离地层层面较近也就是水平井井眼轨迹离地层层面的距离小于测井仪器的探测深度,或与地层层面相交,测井仪器探测的是砂/泥岩互层,测量值为砂岩和泥岩电阻率的综合信息,地层电阻率各向异性表现为宏观各向异性。设在砂/泥岩模型中砂岩电阻率为Rsd,泥岩电阻率为Rsh,砂岩累计厚度为Hsd,泥岩累计厚度为Hsh。此时地层的水平电阻率Rh可表示为:
地层垂直电阻率Rv可表示为:
电阻率各向异性系数λ:
当水平井井眼轨迹远离地层层界面时,也就是水平井井眼轨迹到地层层界面的距离超出测井仪器的探测深度时,假设地层为无限厚的均匀各向异性介质,地层电阻率各向异性表现为微观各向异性。设垂直层界面方向的电阻率为Rv,平行层界面方向的电阻率为Rh,可以推导出地层视电阻率Ra与水平电阻率Rh和垂直电阻率Rv的关系:
式中:为各向异性地层的平均电阻率,单位为(Ω.m);θ为井斜角,单位为(°)。
通过上述电阻率宏观各向异性和微观各向异性的理论研究,结合井斜角、自然伽马和声波时差测井曲线以及泥岩电阻率,泥岩电阻率由区域的经验值确定或用岩石完全含水时电阻率Ro近似替代,采用迭代算法确定地层电阻率宏观和微观总的各向异性系数λ,从而进行地层电阻率各向异性校正,计算出可以代表目的层真实电阻率的水平电阻率Rh。
在本实施例中,如图1所示,采用迭代算法确定地层电阻率宏观和微观总的各向异性系数λ,计算出可以代表目的层真实电阻率的水平电阻率Rh的具体实现步骤如下:
(1)利用自然伽马测井曲线计算地层中的泥质体积含量Vsh,其中,ΔGR为自然电位相对值,无量纲,计算值为小数;GR为自然伽马测井值,单位为(API);GRmin为纯砂岩段的自然伽马值,单位为(API);GRmax为纯泥岩段的自然伽马值,单位为(API);GCUR为希尔奇系数,老地层2,第三系新地层3.7;Vsh为地层中的泥质体积含量,小数。
(2)结合区域的认识确定泥岩电阻率Rsh,也可用岩石完全含水时电阻率Ro近似替代,其中a,m为阿尔奇公式岩电参数,无量纲,由岩电实验数据分析确定,Rw为地层水电阻率,单位为(Ω.m),由水分析资料确定;φ为孔隙度,小数,由声波时差测井曲线计算确定。
(3)对各向异性系数λ赋初始值,结合井斜角θ、步骤(1)计算的泥质体积含量Vsh和步骤(2)计算的泥岩电阻率Rsh,利用公式(1)~(4)计算出中间变量的水平电阻率Rh ·,Rh ·=f(λ,θ,Vsh,Rsh)。
(4)将各向异性系数λ、井斜角θ和步骤(3)计算的水平电阻率Rh ·代入公式(4)计算出中间变量的地层视电阻率Ra ·。
(5)利用步骤(4)计算的Ra ·与实测地层电阻率Ra比较,计算ε,ε=‖Ra-Ra ·‖。若ε大于给定的误差,则重新从步骤(3)开始执行,直到ε满足误差要求为止,输出各向异性系数λ,并计算出可以代表目的层真实电阻率的水平电阻率Rh。
通过电阻率宏观各向异性和微观各向异性的理论研究,建立了水平和垂直电阻率解释图版。如图2所示,在水平和垂直电阻率解释图版中,横轴是水平电阻率与泥岩电阻率的比值,纵轴是垂直电阻率与泥岩电阻率的比值,图中一系列的实线代表的是泥质体积含量,从右到左泥质体积含量Vsh从0逐渐增大到100%,一系列虚线代表的是砂岩电阻率与泥岩电阻率的比值,从下到上Rsd/Rsh的值从0.1逐渐增大到1000Ω.m。对图进行分析,泥质体积含量Vsh可以利用岩性测井曲线如自然伽马测井计算确定,即确定图中的一条实线。然后在该条实线中任意给定A点,就可以从图版中确定A点对应 的水平电阻率Rh和垂直电阻率Rv,再代入到公式(4)计算A点的地层视电阻率,并与水平井实测电阻率进行比较,如果两者差别较大,则重新调整A点,直到A点计算的地层视电阻率与水平井实测电阻率接近为止,最后输出各向异性系数λ,并计算出可以代表目的层真实电阻率的水平电阻率Rh。
实际资料处理过程中,地层电阻率各向异性校正通过编写程序实现。如图3所示,大斜度/水平井地层电阻率各向异性校正处理解释成果图,计算的各向异性系数主要分布在1~2之间。受各向异性的影响,校正后确定的水平电阻率比原始深侧向电阻率低,与自然伽马有很好的一致性,能很好的反映储层的含气性。结合阿尔奇公式采用相同的岩电参数和地层水电阻率,分别利用原始深侧向电阻率和校正后确定的水平电阻率计算了含水饱和度,图中第8道黑色曲线为深侧向电阻率计算的含水饱和度,全井段近似直线,4250m以上含水饱和度平均为50%,4250m以下含水饱和度平均为45%,全井段解释为气层、差气层和干层。该井在3865~3875m、4026~4036m、4250~4260m、4370~4380m、4480~4490m、4590~4600m和4720~4730m井段射孔,日产气5.44万方,水20.5方。利用深侧向电阻率计算的含水饱和度,与试气产水结论不吻合。各向异性校正后确定的水平电阻率计算的含水饱和度如图中第8道红色曲线,对应射孔井段的6、9、11、17、22、25和28号层的含水饱和度分别为59.02%、67.29%、53.61%、60.64%、51.47%、57.49%和60.56%,同时发现2、4和8号层计算的含水饱和度分别为58.82%、57.39%和59.01%,结合研究区气水识别界限,将2、4、6和8号层的原解释结论从气层改为气水同层。从射孔井段看,试气出水主要来源于6号层,表明修改后的解释结论与试气吻合,说明了各向异性校正后的水平电阻率可以代表目的层真实电阻率,能准确的反映储层的含气性特征,可 以应用于水平井测井地层评价。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种大斜度/水平井地层电阻率各向异性校正方法,其特征在于,包括如下步骤,
步骤一,根据电阻率宏观各向异性和微观各向异性的各项参数以及之间的相互关系,建立水平和垂直电阻率解释图版;
步骤二,根据对水平和垂直电阻率解释图版的分析,采用迭代算法确定地层电阻率宏观和微观总的各向异性系数,从而进行地层电阻率各向异性校正,计算出代表目的层真实电阻率的水平电阻率;
步骤一中建立水平和垂直电阻率解释图版时,包括如下步骤,
1.1在大斜度/水平井中,当水平井井眼轨迹离地层层面的距离小于测井仪器的探测深度或与地层层面相交,测井仪器探测的是砂/泥岩互层,测量值为砂岩和泥岩电阻率的综合信息,地层电阻率各向异性表现为宏观各向异性;此时,
地层水平电阻率Rh表示如下,
<mrow>
<msub>
<mi>R</mi>
<mi>h</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
地层垂直电阻率Rv表示如下,
<mrow>
<msub>
<mi>R</mi>
<mi>v</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
电阻率各向异性系数λ表示如下,
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<msqrt>
<mfrac>
<msub>
<mi>R</mi>
<mi>v</mi>
</msub>
<msub>
<mi>R</mi>
<mi>h</mi>
</msub>
</mfrac>
</msqrt>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<msubsup>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
<mn>2</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mfrac>
<msubsup>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
<mn>2</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,Rsd为在砂/泥岩模型中砂岩电阻率,单位为Ω.m;Rsh为泥岩电阻率,单位为Ω.m;Hsd为砂岩累计厚度,单位为m;Hsh为泥岩累计厚度,单位为m;
1.2在大斜度/水平井中,当水平井井眼轨迹到地层层界面的距离超出测井仪器的探测深度时,视地层为无限厚的均匀各向异性介质,地层电阻率各向异性表现为微观各向异性;则得到,
地层视电阻率Ra表示如下,
<mrow>
<msub>
<mi>R</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
</msqrt>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&lambda;R</mi>
<mi>h</mi>
</msub>
</mrow>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为各向异性地层的平均电阻率,单位为Ω.m;θ为井斜角单位为度;
1.3根据步骤1.1和1.2宏观各向异性和微观各向异性系数与水平电阻率、垂直电阻率和地层视电阻率的关系,建立地层的水平和垂直电阻率解释图版;
步骤二中采用迭代算法进行地层电阻率各向异性校正时,包括如下步骤,
2.1利用自然伽马测井曲线得到地层中的泥质体积含量Vsh如下,
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mi>G</mi>
<mi>C</mi>
<mi>U</mi>
<mi>R</mi>
<mo>&times;</mo>
<mi>&Delta;</mi>
<mi>G</mi>
<mi>R</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mi>G</mi>
<mi>C</mi>
<mi>U</mi>
<mi>R</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
其中,ΔGR为自然电位相对值;GR为自然伽马测井值,GRmin为纯砂岩段的自然伽马值,GRmax为纯泥岩段的自然伽马值;GCUR为希尔奇系数;
2.2结合区域的经验值确定泥岩电阻率Rsh,或用岩石完全含水时电阻率Ro近似替代,其中,a、m为阿尔奇公式岩电参数,Rw为地层水电阻率,φ为孔隙度;
2.3给定各向异性系数λ初始值λ=λk和给定误差,结合井斜角θ、步骤2.1计算的泥质体积含量Vsh和步骤2.2计算的泥岩电阻率Rsh,根据公式(1)~(4)得到中间变量的水平电阻率Rh ●,Rh ●=f(λ,θ,Vsh,Rsh);
2.4将各向异性系数λ、井斜角θ和步骤2.3得到的水平电阻率Rh ●,结合公式(4)得到中间变量的地层视电阻率Ra ●;
2.5利用步骤2.4计算的Ra ●与实测地层视电阻率Ra比较,计算ε,ε=||Ra-Ra ●||;若ε大于给定误差,则λ=λk+1并重复步骤2.3~2.5,直到ε满足不大于给定误差的要求,输出电阻率宏观和微观总的各向异性系数λ,并计算出代表目的层真实电阻率的水平电阻率Rh。
2.根据权利要求1所述的一种大斜度/水平井地层电阻率各向异性校正方法,其特征在于,阿尔奇公式岩电参数a和m由岩电实验数据分析确定。
3.根据权利要求1所述的一种大斜度/水平井地层电阻率各向异性校正方法,其特征在于,地层水电阻率Rw由水分析资料确定。
4.根据权利要求1所述的一种大斜度/水平井地层电阻率各向异性校正方法,其特征在于,孔隙度φ由声波时差测井曲线计算确定。
5.根据权利要求1所述的一种大斜度/水平井地层电阻率各向异性校正方法,其特征在于,岩石完全含水时电阻率通过利用声波时差计算孔隙度而确定。
6.根据权利要求1所述的一种大斜度/水平井地层电阻率各向异性校正方法,其特征在于,当取希尔奇系数GCUR时,老地层取2,第三系新地层取3.7。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510463526.8A CN105089663B (zh) | 2015-07-31 | 2015-07-31 | 一种大斜度/水平井地层电阻率各向异性校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510463526.8A CN105089663B (zh) | 2015-07-31 | 2015-07-31 | 一种大斜度/水平井地层电阻率各向异性校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105089663A CN105089663A (zh) | 2015-11-25 |
CN105089663B true CN105089663B (zh) | 2017-11-07 |
Family
ID=54570894
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510463526.8A Active CN105089663B (zh) | 2015-07-31 | 2015-07-31 | 一种大斜度/水平井地层电阻率各向异性校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105089663B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106522932B (zh) * | 2016-09-28 | 2019-07-09 | 中国石油天然气股份有限公司 | 获取水平井各向异性地层电阻率的方法及装置 |
CN106990448B (zh) * | 2017-04-05 | 2018-11-16 | 中国石油天然气股份有限公司 | 一种海相碎屑岩水平井隔夹层测井解释方法及装置 |
CN108804728B (zh) * | 2017-05-02 | 2021-07-20 | 中国石油化工股份有限公司 | 水平井地层储层分级分析方法及计算机可读存储介质 |
CN107784159B (zh) * | 2017-09-19 | 2021-09-28 | 中国石油天然气集团公司 | 一种储层电阻率各向异性系数的确定方法 |
CN107907913A (zh) * | 2017-11-10 | 2018-04-13 | 西南石油大学 | 一种利用成像、倾角资料校正地层电阻率的方法 |
CN107748393B (zh) * | 2017-11-10 | 2019-04-16 | 西南石油大学 | 一种基于数值模拟的地层倾角对电阻率影响的校正方法 |
CN110320562B (zh) * | 2018-03-30 | 2020-11-13 | 中国石油化工股份有限公司 | 页岩气储层水平井中声波时差的校正方法 |
CN109236284B (zh) * | 2018-09-06 | 2022-05-10 | 中国石油天然气集团有限公司 | 一种感应与侧向联测确定高陡地层真电阻率的方法 |
CN112727441A (zh) * | 2021-03-11 | 2021-04-30 | 大庆油田有限责任公司 | 一种适用于薄油层水平井油层对比方法 |
CN113107464B (zh) * | 2021-05-11 | 2024-05-07 | 中国石油天然气集团有限公司 | 一种水平井步进式水淹层识别测井方法 |
CN113356834B (zh) * | 2021-07-14 | 2023-04-18 | 中国石油天然气股份有限公司 | 井眼声波时差的校正方法、校正装置和电子设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101004454A (zh) * | 2006-01-20 | 2007-07-25 | 中国石油天然气集团公司 | 目标最小化的三维电磁快速反演方法 |
CN101775981A (zh) * | 2009-01-09 | 2010-07-14 | 中国石油天然气集团公司 | 一种确定地层真电阻率的方法 |
CN101899973A (zh) * | 2010-05-10 | 2010-12-01 | 中国石油天然气集团公司 | 一种地层水电阻率的测量方法及装置 |
CN102608666A (zh) * | 2012-03-23 | 2012-07-25 | 中煤科工集团西安研究院 | 瞬变电磁数据快速准确的深度反演方法 |
CN102635347A (zh) * | 2012-03-30 | 2012-08-15 | 中国电子科技集团公司第二十二研究所 | 一种定量将薄互地层等效成水平和垂直电阻率的方法 |
CN104775811A (zh) * | 2015-04-13 | 2015-07-15 | 中国海洋石油总公司 | 一种地层各向异性信息的提取及校正方法和*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8433518B2 (en) * | 2009-10-05 | 2013-04-30 | Schlumberger Technology Corporation | Multilevel workflow method to extract resistivity anisotropy data from 3D induction measurements |
-
2015
- 2015-07-31 CN CN201510463526.8A patent/CN105089663B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101004454A (zh) * | 2006-01-20 | 2007-07-25 | 中国石油天然气集团公司 | 目标最小化的三维电磁快速反演方法 |
CN101775981A (zh) * | 2009-01-09 | 2010-07-14 | 中国石油天然气集团公司 | 一种确定地层真电阻率的方法 |
CN101899973A (zh) * | 2010-05-10 | 2010-12-01 | 中国石油天然气集团公司 | 一种地层水电阻率的测量方法及装置 |
CN102608666A (zh) * | 2012-03-23 | 2012-07-25 | 中煤科工集团西安研究院 | 瞬变电磁数据快速准确的深度反演方法 |
CN102635347A (zh) * | 2012-03-30 | 2012-08-15 | 中国电子科技集团公司第二十二研究所 | 一种定量将薄互地层等效成水平和垂直电阻率的方法 |
CN104775811A (zh) * | 2015-04-13 | 2015-07-15 | 中国海洋石油总公司 | 一种地层各向异性信息的提取及校正方法和*** |
Non-Patent Citations (1)
Title |
---|
水平井地层电阻率各向异性校正方法研究;罗少成等;《测井技术》;20090430;126-129、147 * |
Also Published As
Publication number | Publication date |
---|---|
CN105089663A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105089663B (zh) | 一种大斜度/水平井地层电阻率各向异性校正方法 | |
CN108518212B (zh) | 一种计算页岩气藏复杂裂缝网络非稳态产量的方法 | |
Honarpour et al. | Characterization of critical fluid, rock, and rock-fluid properties-impact on reservoir performance of liquid-rich shales | |
CN103256046B (zh) | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 | |
Worthington | Net pay—what is it? What does it do? How do we quantify it? How do we use it? | |
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
CN110321595B (zh) | 一种测井提取静态品质系数的断层封闭性评价方法 | |
CN110320139B (zh) | 缝洞型储层裂缝孔隙度定量评价方法及*** | |
CN107829731B (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
CN110554441A (zh) | 一种基于钙质含量校正的致密砂岩含气饱和度计算方法 | |
Zou et al. | Porosity type analysis and permeability model for micro-trans-pores, meso-macro-pores and cleats of coal samples | |
CN105931125B (zh) | 一种致密油分段多簇体积压裂水平井产量预测方法 | |
CN110954944A (zh) | 一种断层圈闭含油高度地震预测方法 | |
CN108804728B (zh) | 水平井地层储层分级分析方法及计算机可读存储介质 | |
CN110927794A (zh) | 一种致密储层裂缝识别与孔隙度定量计算的方法 | |
CN104834003B (zh) | 一种非常规致密储层的相控压缩系数地震预测方法 | |
CN112698399B (zh) | 一种砂砾岩井测震联动约束高效储层定量预测方法与*** | |
CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
Yasin et al. | Study on brittleness templates for shale gas reservoirs-A case study of Longmaxi shale in Sichuan Basin, southern China | |
Zhang et al. | Postfracturing permeability prediction for CBM well with the analysis of fracturing pressure decline | |
Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
Farouk et al. | Integrating petrophysical and geomechanical rock properties for determination of fracability of an Iraqi Tight Oil Reservoir | |
Yang et al. | Study on the stimulation effectiveness evaluation of large-scale hydraulic fracturing simulation experiment based on optical scanning technology | |
CN110529106B (zh) | 一种利用测井资料确定煤层显微组分含量的方法 | |
CN112943229A (zh) | 一种储气库盖层突破压力连续预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |