CN108005646A - 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 - Google Patents
基于随钻电磁波测井资料的地层各向异性电阻率提取方法 Download PDFInfo
- Publication number
- CN108005646A CN108005646A CN201711336574.6A CN201711336574A CN108005646A CN 108005646 A CN108005646 A CN 108005646A CN 201711336574 A CN201711336574 A CN 201711336574A CN 108005646 A CN108005646 A CN 108005646A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- electromagnetic wave
- inverting
- mover
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 23
- 241001074085 Scophthalmus aquosus Species 0.000 claims abstract description 19
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 230000003044 adaptive effect Effects 0.000 claims abstract description 7
- 238000012937 correction Methods 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims abstract 2
- 239000011435 rock Substances 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 10
- 230000004044 response Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 230000017105 transposition Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 230000007423 decrease Effects 0.000 claims 1
- 230000007613 environmental effect Effects 0.000 abstract description 2
- 238000005755 formation reaction Methods 0.000 description 14
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 102220350010 c.119C>A Human genes 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 102220206869 rs1057522105 Human genes 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/13—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling by electromagnetic energy, e.g. radio frequency
Landscapes
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于随钻电磁波测井资料的地层各向异性电阻率提取方法，包括步骤：s1.对输入的随钻电磁波测井曲线进行环境校正；s2.获取区块地质结构信息；s3.对测井数据滑动开窗处理；s4.建立单界面反演地层模型；s5.判断是否有上一滑动窗口的反演结果，若有则执行步骤s6，否则执行步骤s8；s6.将上一滑动窗口反演结果作为当前滑动窗口的反演初始值；s7.采用非线性最优化算法对随钻电磁波测井资料反演，若实测资料与模拟资料拟合差在容许精度误差允许之内，执行步骤s9，否则执行步骤s8；s8.采用自适应多初值选取方法，并结合非线性最优化算法对随钻电磁波测井资料反演；s9.对多个反演结果进行优化选取。本发明可为实时地质导向和储层评价提供准确电阻率信息。
Description
技术领域
本发明涉及石油勘探开发领域，属于电测井方法范畴，具体的说是涉及一种基于随钻电磁波测井资料的地层各向异性电阻率提取方法。
背景技术
随钻电磁波测井因其分辨率高、对各向异性电阻率敏感，在地质结构相对明确的油藏中得到广泛应用。但因受地层邻近地层界面、相对倾角及各向异性等影响，不同探测深度的视电阻率曲线分离严重，存在异常电阻率“犄角”，导致视电阻率值远大于地层真实值。因此，地层真实电阻率信息的准确、实时反演是地质导向与储层评价首要解决的关键问题之一。
目前，随钻电磁波测井反演主要集中在钻后储层精细评价方面，而关于随钻过程中实时电阻率曲线校正研究极少。从最优化的角度来看，随钻电磁波测井资料的实时处理主要存在以下三个问题：1)反演计算中采用的多层模型速度慢；2)多层模型反演参数多，Jacobian计算量大；3)存在多个局部极小值，且已知地层信息少，反演结果易困于局部极小值。
发明内容
本发明的目的在于提出一种基于随钻电磁波测井资料的地层各向异性电阻率提取方法，以期为实时地质导向和储层评价提供准确电阻率信息。
本发明为了实现上述目的，采用如下技术方案：
基于随钻电磁波测井资料的地层各向异性电阻率提取方法，包括如下步骤：
s1.对输入的随钻电磁波测井曲线进行井眼环境校正；
s2.获取区块的地质结构信息；
s3.对测井数据进行滑动开窗处理，将初始复杂地层模型转化为一系列开窗后的地层模型；
s4.对开窗后的地层模型进一步简化，将邻近层做等效处理；根据区块层厚，在随钻电磁波测井过程中将开窗后的地层模型简化为单界面地层模型，即单界面反演地层模型；
s5.判断是否有上一滑动窗口的反演结果，如果有，则执行步骤s6，否则执行步骤s8；
s6.将上一滑动窗口反演的结果作为当前滑动窗口的反演初始值；
s7.采用非线性最优化算法对随钻电磁波测井资料进行反演，如果实测资料与模拟资料拟合差在容许精度误差允许之内，则执行步骤s9，否则执行步骤s8；
s8.采用自适应多初值选取方法，并结合非线性最优化算法对随钻电磁波测井资料进行反演；
s9.对步骤s8中的多个反演结果进行优化选取，确定最优地层水平和垂直电阻率。
优选地，所述步骤s7和步骤s8中，非线性最优化算法具体为：
采用Taylor展开，对非线性反演问题进行线性处理，同时配合正则化的手段；
反演所用的代价函数C(m)为：
C(m)＝||W·(S(m)-d)||2+λ||m-mref||2+ξ||n||2 (1)
公式(1)中，d表示实测数据，即相位差和幅度比视电导率曲线，m为待反演的参数向量，||W·(S(m)-d)||2表示实测数据d与正演响应S(m)向量之差的L2范数；
||m-mref||2表示已知模型参考矢量mref与待反演的参数矢量m的差异；
ξ||n||2表示反演值与上一滑动窗口反演结果mpre的差异，n＝mk-mpre；
公式(1)中，λ和ξ表示正则化参数，W为测量数据的权重系数矩阵；
采用Gauss-Newton最优化方法对公式(1)进行求解，令公式(1)对m的导数为零，第k次迭代中，可得：
公式(2)中，Δm为反演模型的下降方向，为雅克比矩阵，T表示矩阵的转置，ν表示下降方向上的最优步长；正则化参数λ采用式(3)进行自动更新：
优选地，所述步骤s8中，自适应多初值方法具体为：
s8.1基于仪器响应对各个参数敏感性的不同与步骤s2给出的地质结构信息，分别确定各待反演参数初值选取的数量；各待反演参数的初值选取方式，参照步骤s8.2-s8.5；
s8.2根据地层结构信息与井眼轨迹，确定地层相对倾角范围θ1～θ2，选取初值个数为n1；
s8.3根据井眼轨迹和仪器的探测特性，确定仪器到临近地层界面距离范围H1～H2，选取初值个数为n2；
s8.4根据实测曲线的视电阻率值，定给定模型水平电阻率Rh，参数的初值个数为1；
s8.5根据目的层电阻率，建立相应高阻围岩模型和低阻围岩模型，选取初值个数为n3；
s8.6将步骤s8.2-s8.5选取的初值进行自由组合，构成n1*n2*n3个初始反演地层模型，对各个初始反演地层模型进行正则化Gauss-Newton反演迭代，获取n1*n2*n3个反演结果。
优选地，所述步骤s4中，近邻层做等效处理的具体过程为：
将围岩等效为各向同性地层，将目的层等效为各向异性地层。
优选地，所述步骤s4中，单界面地层模型中待反演的参数包括围岩电阻率、仪器到邻近界面的距离、当前层的各向异性电阻率以及仪器与地层的相对倾角。
本发明具有如下优点：
本发明摒弃了以往所采用的1D多层反演模型，考虑到随钻电磁波测井的探测深度较小，当地层层厚较大时，邻近围岩对响应的贡献很小甚至可以忽略，在随钻电磁波测井过程中可将地层简化为两层，即围岩等效为各向同性地层，目的层等效为各向异性地层，以简化求解参数的个数和反演的难度。传统三层模型简化后的模型只需要邻近围岩电阻率、仪器到界面的距离、当前层的各向异性电阻率及仪器与地层的相对倾角五个反演参数；计算一个点只需要进行n1*n2*n3次迭代(约60次)，1D反演算法每秒钟可以处理6个测量点，极大提高了反演速度，可以解决随钻过程中实时电阻率曲线校正问题。
附图说明
图1为本发明中基于随钻电磁波测井资料的地层各向异性电阻率提取方法的流程图；
图2为本发明中复杂地层模型示意图；
图3为本发明中开窗处理后的等效地层模型示意图；
图4为本发明中简化后的单界面反演地层模型示意图；
图5为本发明中五层地层模型对应的随钻电磁波测井响应图；
图6为本发明中五层地层模型对应的随钻电磁波反演结果图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明：
结合图1所示，本发明实施例给出了一种基于随钻电磁波测井资料的地层各向异性电阻率提取方法，其包括如下步骤：
s1.对输入的随钻电磁波测井曲线进行环境校正，以消除井眼对测井响应的影响。
s2.获取区块的地质结构信息，诸如层厚、电阻率及地层倾斜等变化范围。
s3.对测井数据进行滑动开窗处理，将初始的复杂地层模型转化为一系列开窗后的地层模型，如图2所示和图3所示，将原始复杂三维问题简化为一系列的一维问题。每个窗口均为无限延伸的1维水平层状各向异性地层，窗口长度取决于仪器长度，一般设置为1～2m。
s4.对开窗后的地层模型(多层模型)进一步简化，将邻近层做等效处理，即将围岩等效为各向同性地层，将目的层等效为各向异性地层，以简化求解参数的个数和反演的难度。
根据区块层厚，在随钻电磁波测井过程中将开窗后的地层模型简化成单界面地层模型，即单界面反演地层模型如图4所示。
单界面地层模型中待反演的参数包括围岩电阻率(Rs)、仪器到邻近界面的距离(H)、当前层的各向异性电阻率(Rh,Rv)及仪器与地层的相对倾角(θ)。
s5.判断是否有上一滑动窗口的反演结果，如果有，则执行步骤s6，否则执行步骤s8。
s6.将上一滑动窗口反演的结果作为当前滑动窗口的反演初始值。
考虑到地层性质横向非均质变化缓慢，如有上个窗口结果，可将上一滑动窗口反演的结果作为当前窗口的初始值，其优点是代价函数收敛快，计算效率高。
s7.采用非线性最优化算法对随钻电磁波测井资料进行反演，如果实测资料与模拟资料拟合差在容许精度误差允许之内，则执行步骤s9，否则执行步骤s8。
s8.因代价函数存在多个局部极小值，为保证反演算法获得全局最优解，采用自适应多初值选取方法，并结合非线性最优化算法对随钻电磁波测井资料进行反演。
s9.对步骤s8中的多个反演结果进行优化选取，确定最优的地层水平和垂直电阻率。
在步骤s7和步骤s8中，非线性最优化算法具体为：
采用Taylor展开，对非线性反演问题进行线性处理，同时配合正则化的手段，以简化计算量和提高计算速度；反演所用的代价函数C(m)为：
C(m)＝||W·(S(m)-d)||2+λ||m-mref||2+ξ||n||2 (1)
公式(1)中，d表示实测数据，即相位差和幅度比视电导率曲线，m为待反演的参数向量，||W·(S(m)-d)||2表示实测数据d与正演响应S(m)向量之差的L2范数；
||m-mref||2表示已知模型参考矢量mref与待反演的参数矢量m的差异；
ξ||n||2表示反演值与上一滑动窗口反演结果mpre的差异，n＝mk-mpre；
公式(1)中，λ和ξ表示正则化参数，W为测量数据的权重系数矩阵；
采用Gauss-Newton最优化方法对公式(1)进行求解，令公式(1)对m的导数为零，第k次迭代中可得：
公式(2)中，Δm为反演模型的下降方向，为雅克比矩阵，T表示矩阵的转置，ν表示下降方向上的最优步长；正则化参数λ的采用式(3)进行自动更新：
在步骤s8中，自适应多初值方法具体为：
s8.1基于仪器响应对各个参数敏感性的不同与步骤s2给出的地质结构信息，分别确定各待反演参数初值选取的数量；各待反演参数的初值选取方式，参照步骤s8.2-s8.5。
s8.2根据地层结构信息与井眼轨迹，确定地层相对倾角范围θ1～θ2，选取初值个数为n1，例如可以选取5个。
s8.3根据井眼轨迹和仪器的探测特性，确定仪器到临近地层界面距离范围H1～H2，选取初值个数n2，例如可以选取4个。
s8.4根据实测曲线的视电阻率的平均值，确定给定模型水平电阻率Rh，即该参数的初值个数为1。
s8.5根据目的层电阻率，建立相应高阻围岩模型和低阻围岩模型，选取初值个数n3，例如可以选取3个。
s8.6将步骤s8.2-s8.5选取的初值进行自由组合，构成n1*n2*n3个初始反演地层模型，对各个初始反演地层模型进行正则化Gauss-Newton反演迭代，获取n1*n2*n3个反演结果。
如图5和图6所示，下面以五层地层模型对应的随钻电磁波测井响应与对该五层地层模型进行反演后的结果进行对比分析，以验证本发明方法的有效性。
其中，图5和图6中，各个参数的表示含义如下：
A28H和P28H表示频率为2MHz，源距为28inch的幅度比和相位差电阻率；
A40H和P40H表示频率为2MHz，源距为40inch的幅度比和相位差电阻率；
Rh(Model)和Rv(Model)表示原始地层模型水平和垂直电阻率；
Rh(Inv)和Rv(Inv)表示反演得到的水平和垂直电阻率。
在图5中，仪器自上而下穿过地层，与地层法向的相对倾角为89度。实线为相位差和幅度比视电阻率曲线，虚线为原始地层模型水平和垂直电阻率。在图6中，实线为原始地层模型水平和垂直电阻率，虚点代表反演得到的水平和垂直电阻率。
实验证明，反演的电阻率与地层模型值一致性高，验证了反演算法的准确性和可靠性。
当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替代、明显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。
Claims (5)
1.基于随钻电磁波测井资料的地层各向异性电阻率提取方法，其特征在于，包括如下步骤：
s1.对输入的随钻电磁波测井曲线进行井眼环境校正；
s2.获取区块的地质结构信息；
s3.对测井数据进行滑动开窗处理，将初始复杂地层模型转化为一系列开窗后的地层模型；
s4.对开窗后的地层模型进一步简化，将邻近层做等效处理；根据区块层厚，在随钻电磁波测井过程中将开窗后的地层模型简化为单界面地层模型，即单界面反演地层模型；
s5.判断是否有上一滑动窗口的反演结果，如果有，则执行步骤s6，否则执行步骤s8；
s6.将上一滑动窗口反演的结果作为当前滑动窗口的反演初始值；
s7.采用非线性最优化算法对随钻电磁波测井资料进行反演，如果实测资料与模拟资料拟合差在容许精度误差允许之内，则执行步骤s9，否则执行步骤s8；
s8.采用自适应多初值选取方法，并结合非线性最优化算法对随钻电磁波测井资料进行反演；
s9.对步骤s8中的多个反演结果进行优化选取，确定最优的地层水平和垂直电阻率。
2.根据权利要求1所述的基于随钻电磁波测井资料的地层各向异性电阻率提取方法，其特征在于，所述步骤s7和步骤s8中，非线性最优化算法具体为：
采用Taylor展开，对非线性反演问题进行线性处理，同时配合正则化的手段；
反演所用的代价函数C(m)为：
C(m)＝||W·(S(m)-d)||2+λ||m-mref||2+ξ||n||2 (1)
公式(1)中，d表示实测数据，即相位差和幅度比视电导率曲线，m为待反演的参数向量，||W·(S(m)-d)||2表示实测数据d与正演响应S(m)向量之差的L2范数；
||m-mref||2表示已知模型参考矢量mref与待反演的参数矢量m的差异；
ξ||n||2表示反演值与上一滑动窗口反演结果mpre的差异，n＝mk-mpre；
λ和ξ表示正则化参数，W为测量数据的权重系数矩阵；
采用Gauss-Newton最优化方法对公式(1)进行求解，令公式(1)对m的导数为零，第k次迭代中可得：
<mrow>
<mi>&Delta;</mi>
<mi>m</mi>
<mo>=</mo>
<msup>
<mi>m</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>m</mi>
<mi>k</mi>
</msup>
<mo>=</mo>
<mo>-</mo>
<mi>v</mi>
<mfrac>
<mrow>
<msup>
<mover>
<mover>
<mi>J</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&OverBar;</mo>
</mover>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>m</mi>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msup>
<mi>W</mi>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>W</mi>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>(</mo>
<msup>
<mi>m</mi>
<mi>k</mi>
</msup>
<mo>)</mo>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>m</mi>
<mi>k</mi>
</msup>
<mo>-</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&xi;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>m</mi>
<mi>k</mi>
</msup>
<mo>-</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mover>
<mover>
<mi>J</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&OverBar;</mo>
</mover>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>m</mi>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msup>
<mi>W</mi>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>W</mi>
<mo>&CenterDot;</mo>
<mover>
<mover>
<mi>J</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<msup>
<mi>m</mi>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&lambda;</mi>
<mi>I</mi>
<mo>+</mo>
<mi>&xi;</mi>
<mi>I</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(2)中，Δm为反演模型的下降方向，为雅克比矩阵，T表示矩阵的转置，ν表示下降方向上的最优步长；正则化参数λ采用公式(3)进行自动更新：
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<mn>0.1</mn>
<mo>&times;</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>W</mi>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mover>
<mover>
<mi>J</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>m</mi>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>m</mi>
<mo>-</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
3.根据权利要求1所述的基于随钻电磁波测井资料的地层各向异性电阻率提取方法，其特征在于，所述步骤s8中，自适应多初值方法具体为：
s8.1基于仪器响应对各个参数敏感性的不同与步骤s2给出的地质结构信息，分别确定各待反演参数初值选取的数量；各待反演参数的初值选取方式，参照步骤s8.2-s8.5；
s8.2根据地层结构信息与井眼轨迹，确定地层相对倾角范围θ1～θ2，选取初值个数为n1；
s8.3根据井眼轨迹和仪器的探测特性，确定仪器到临近地层界面距离范围H1～H2，选取初值个数为n2；
s8.4根据实测曲线的视电阻率值，定给定模型水平电阻率Rh，参数的初值个数为1；
s8.5根据目的层电阻率，建立相应高阻围岩模型和低阻围岩模型，选取初值个数为n3；
s8.6将步骤s8.2-s8.5选取的初值进行自由组合，构成n1*n2*n3个初始反演地层模型，对各个初始反演地层模型进行正则化Gauss-Newton反演迭代，获取n1*n2*n3个反演结果。
4.根据权利要求1所述的基于随钻电磁波测井资料的地层各向异性电阻率提取方法，其特征在于，所述步骤s4中，近邻层做等效处理的具体过程为：
将围岩等效为各向同性地层，将目的层等效为各向异性地层。
5.根据权利要求1所述的基于随钻电磁波测井资料的地层各向异性电阻率提取方法，其特征在于，所述步骤s4中，单界面地层模型中待反演的参数包括围岩电阻率、仪器到邻近界面的距离、当前层的各向异性电阻率以及仪器与地层的相对倾角。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711336574.6A CN108005646B (zh) | 2017-12-14 | 2017-12-14 | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711336574.6A CN108005646B (zh) | 2017-12-14 | 2017-12-14 | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108005646A true CN108005646A (zh) | 2018-05-08 |
| CN108005646B CN108005646B (zh) | 2021-03-02 |
Family
ID=62058791
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711336574.6A Active CN108005646B (zh) | 2017-12-14 | 2017-12-14 | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108005646B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108952690A (zh) * | 2018-08-01 | 2018-12-07 | 中国石油大学(华东) | 基于随钻方位电磁波测井资料的地层界面实时提取方法 |
| CN109577963A (zh) * | 2018-10-19 | 2019-04-05 | 中国石油天然气股份有限公司 | 模拟各向异性地层阵列侧向测井响应的装置及其构建方法 |
| CN111980662A (zh) * | 2020-07-14 | 2020-11-24 | 中国石油大学(华东) | 一种斜井各向异性地层阵列侧向测井资料快速处理方法 |
| CN112100886A (zh) * | 2020-08-28 | 2020-12-18 | 中国石油天然气集团有限公司 | 随钻电磁波测井响应计算的完全匹配层构造方法及系统 |
| CN113868919A (zh) * | 2021-09-30 | 2021-12-31 | 西南石油大学 | 一种随钻电磁波测井3d模拟简化方法 |
| CN113887034A (zh) * | 2021-09-27 | 2022-01-04 | 西南石油大学 | 一种适用于2d地层模型的随钻电磁波测井跨维模拟方法 |
| CN116398126A (zh) * | 2023-06-07 | 2023-07-07 | 山东万洋石油科技有限公司 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2831644A1 (en) * | 2012-03-27 | 2015-02-04 | Services Pétroliers Schlumberger | Anisotropy processing in low angle wells |
| CN104775811A (zh) * | 2015-04-13 | 2015-07-15 | 中国海洋石油总公司 | 一种地层各向异性信息的提取及校正方法和系统 |
| WO2017078708A1 (en) * | 2015-11-04 | 2017-05-11 | Halliburton Energy Services, Inc. | Conductivity-depth transforms of electromagnetic telemetry signals |
| CN107002488A (zh) * | 2014-12-31 | 2017-08-01 | 哈利伯顿能源服务公司 | 使用各向异性电容率和电阻率的基于多分量信号的测量值而进行的地层测井 |
| CN107002489A (zh) * | 2015-01-07 | 2017-08-01 | 哈里伯顿能源服务公司 | 用于电阻率反演的函数地球模型参数化 |
| CN107191181A (zh) * | 2017-05-24 | 2017-09-22 | 中国石油大学(华东) | 一种基于电磁散射的井周界面探测方法 |
-
2017
- 2017-12-14 CN CN201711336574.6A patent/CN108005646B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2831644A1 (en) * | 2012-03-27 | 2015-02-04 | Services Pétroliers Schlumberger | Anisotropy processing in low angle wells |
| CN107002488A (zh) * | 2014-12-31 | 2017-08-01 | 哈利伯顿能源服务公司 | 使用各向异性电容率和电阻率的基于多分量信号的测量值而进行的地层测井 |
| CN107002489A (zh) * | 2015-01-07 | 2017-08-01 | 哈里伯顿能源服务公司 | 用于电阻率反演的函数地球模型参数化 |
| CN104775811A (zh) * | 2015-04-13 | 2015-07-15 | 中国海洋石油总公司 | 一种地层各向异性信息的提取及校正方法和系统 |
| WO2017078708A1 (en) * | 2015-11-04 | 2017-05-11 | Halliburton Energy Services, Inc. | Conductivity-depth transforms of electromagnetic telemetry signals |
| CN107191181A (zh) * | 2017-05-24 | 2017-09-22 | 中国石油大学(华东) | 一种基于电磁散射的井周界面探测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 夏培等: "水平井各向异性地层双侧向测井响应数值模拟", 《石油天然气学报（江汉石油学院学报）》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108952690A (zh) * | 2018-08-01 | 2018-12-07 | 中国石油大学(华东) | 基于随钻方位电磁波测井资料的地层界面实时提取方法 |
| CN108952690B (zh) * | 2018-08-01 | 2022-01-25 | 中国石油大学(华东) | 基于随钻方位电磁波测井资料的地层界面实时提取方法 |
| CN109577963A (zh) * | 2018-10-19 | 2019-04-05 | 中国石油天然气股份有限公司 | 模拟各向异性地层阵列侧向测井响应的装置及其构建方法 |
| CN111980662A (zh) * | 2020-07-14 | 2020-11-24 | 中国石油大学(华东) | 一种斜井各向异性地层阵列侧向测井资料快速处理方法 |
| CN111980662B (zh) * | 2020-07-14 | 2023-11-07 | 中国石油大学(华东) | 一种斜井各向异性地层阵列侧向测井资料快速处理方法 |
| CN112100886A (zh) * | 2020-08-28 | 2020-12-18 | 中国石油天然气集团有限公司 | 随钻电磁波测井响应计算的完全匹配层构造方法及系统 |
| CN113887034A (zh) * | 2021-09-27 | 2022-01-04 | 西南石油大学 | 一种适用于2d地层模型的随钻电磁波测井跨维模拟方法 |
| CN113868919A (zh) * | 2021-09-30 | 2021-12-31 | 西南石油大学 | 一种随钻电磁波测井3d模拟简化方法 |
| CN116398126A (zh) * | 2023-06-07 | 2023-07-07 | 山东万洋石油科技有限公司 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
| CN116398126B (zh) * | 2023-06-07 | 2023-08-25 | 山东万洋石油科技有限公司 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108005646B (zh) | 2021-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108005646A (zh) | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 | |
| US8364404B2 (en) | System and method for displaying data associated with subsurface reservoirs | |
| CN109753755B (zh) | 一种确定储层含水饱和度的方法 | |
| CN108345719B (zh) | 一种大斜度井/水平井随钻电磁波测井资料钻后反演方法 | |
| CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
| US20130204534A1 (en) | Method Of Estimating A Subterranean Formation Property | |
| CN108073765A (zh) | 一种水平井常规随钻测井地层界面识别与边界距反演方法 | |
| BR112018072718B1 (pt) | Método e aparelho para estimativa de um parâmetro invertido de uma formação subsuperficial | |
| CN106407574A (zh) | 一种多分量随钻方位电磁波仪器的快速正反演处理方法 | |
| EP2253971A1 (en) | Method of monitoring a hydrocarbon reservoir | |
| CN103883318B (zh) | 一种中低矿化度地层水储层的参数反演方法及装置 | |
| DE112016005364T5 (de) | Verfahren zum Auswählen eines Erdmodells von einer Vielzahl von Erdmodellen | |
| WO2017069650A1 (en) | Reservoir monitoring using galvanically excited transient electromagnetic fields | |
| CN107784159A (zh) | 一种储层电阻率各向异性系数的确定方法 | |
| Luo et al. | An ensemble-based framework for proactive geosteering | |
| CN110532507B (zh) | 一种提高水平井致密油储层钻遇率的方法 | |
| Kuchuk et al. | Determination of in situ two-phase flow properties through downhole fluid movement monitoring | |
| CN109885927B (zh) | 一种地层径向电阻率连续反演方法 | |
| CN108952690A (zh) | 基于随钻方位电磁波测井资料的地层界面实时提取方法 | |
| CN107526117B (zh) | 基于自动编码和超限学习联合网络的声波速度预测方法 | |
| CN110488367B (zh) | 一种基于阵列侧向测井资料的电阻率反演初值选取方法 | |
| CN105464650A (zh) | 一种随钻测井解释方法 | |
| CN113775326A (zh) | 可动水饱和度的评价方法、装置、电子设备及介质 | |
| CN113887034B (zh) | 一种适用于2d地层模型的随钻电磁波测井跨维模拟方法 | |
| US11060397B2 (en) | Disposing a carrier downhole in a wellbore to evaluate an earth formation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |