(12)发明专利申请
(10)申请公布号 CN 107630697 A(43)申请公布日 2018.01.26
(21)申请号 201710884216.2(22)申请日 2017.09.26
(71)申请人 长江大学
地址 430100 湖北省武汉市蔡甸区蔡甸街
大学路111号(72)发明人 许巍 朱小露 刘迪仁 (74)专利代理机构 武汉开元知识产权代理有限
公司 42104
代理人 冯超 陈家安(51)Int.Cl.
E21B 49/00(2006.01)G06F 17/50(2006.01)
权利要求书2页 说明书6页 附图6页
CN 107630697 A(54)发明名称
基于随钻电磁波电阻率测井的地层电阻率联合反演方法(57)摘要
本发明公开了一种基于随钻电磁波电阻率测井的地层电阻率联合反演方法，该方法首先获取随钻电磁波电阻率测井幅度比和相位差数据；再根据幅度比和相位差数据进行加权求和，构建联合约束的反演目标函数；最后基于联合约束的反演目标函数，利用粒子群算法反演得到地层电阻率。本发明所提供的地层电阻率反演方法得到的地层电阻率曲线兼顾了传统幅度比电阻率曲线较好的径向探测深度及相位差电阻率曲线较好的纵向分辨率的优势，更加充分的利用了随钻电磁波电阻率测井信息，而反演工作量较传统反演方法缩减了一半。
CN 107630697 A
权 利 要 求 书
1/2页
1.一种基于随钻电磁波电阻率测井的地层电阻率联合反演方法，其特征在于：包括以下步骤：
1)获取随钻电磁波电阻率测井幅度比和相位差数据；2)根据幅度比和相位差数据进行加权求和，构建联合约束的反演目标函数E；3)基于联合约束的反演目标函数，利用粒子群算法反演得到地层电阻率。
2.根据权利要求1所述基于随钻电磁波电阻率测井的地层电阻率联合反演方法，其特征在于：所述步骤1)中，幅度比和相位差数据为实际测量数据和正演模拟数据，其中，幅度比和相位差数据包括正演模型幅度比信号、实测幅度比信号、正演模型相位差信号，实测相位差信号。
3.根据权利要求1所述基于随钻电磁波电阻率测井的地层电阻率联合反演方法，其特征在于：所述步骤2)中，联合约束的反演目标函数E为：
E＝||EATTS-EATTM||2+||△φS-△φM||2·cor2其中，EATTS为正演模型幅度比信号，EATTM为实测幅度比信号，E△φ为相位差反演目标函数，△φS为正演模型相位差信号，△φM为实测相位差信号，
cor为幅度比和相位差信号加权系数，表达式为：
4.根据权利要求1所述基于随钻电磁波电阻率测井的地层电阻率联合反演方法，其特征在于：所述步骤3)中，利用粒子群算法反演得到地层电阻率的步骤：
步骤S1、获取反演目标函数E，设定粒子数量Num、粒子分布范围F、最大迭代次数itermax及目标函数阈值ε；
步骤S2、获取给定范围内均匀分布的粒子初始位置及随机分布的粒子初始运动速度；步骤S3、计算初始粒子群个体极值P0best及全局极值N0best，初始模型参数X0＝N0best；步骤S4、根据粒子运动速度及位置更新表达式，获得更新的粒子运动速度及更新位置坐标；
步骤S5、以步骤4)获得的粒子运动速度及更新位置，确定当前位置，粒子群个体极值Pbest及全局极值Nbest，当前模型参数X＝Nbest；
步骤S6、将初始模型参数X0和当前模型参数X带入目标函数中，分别得到反演目标函数E的模型参数初始值E(X0)和当前值E(X)；
步骤S7、求取目标函数的模型参数初始值和当前值的差值△E＝E(X0)-E(X)，即为判断目标函数的模型参数初始值和当前值的差值是否符合预设规则，当目标函数的当前值是否小于或等于目标函数阈值ε，若是，则执行步骤S8，否则，返回执行步骤S4重新获取更新的粒子运动速度及位置；
步骤S8、将当前模型参数X赋予初始模型参数X0，然后执行步骤S9；步骤S9、判断迭代次数iter是否与所述最大迭代次数itermax相等，若是相等，则执行步骤S10，否则，返回步骤S2；
步骤S10、输出初始模型参数X0，作为反演得到的参数值。
5.根据权利要求4所述基于随钻电磁波电阻率测井的地层电阻率联合反演方法，其特
2
CN 107630697 A
权 利 要 求 书
2/2页
征在于：所述步骤S4)中，
粒子运动速度及位置更新表达式分别为：
其中，和分别为迭代更新前后第i个粒子运动速度；和分别为迭代更新前
后第i个粒子位置；c1＝2和c2＝2是加速常量，分别用于调节向全局最优粒子和个体最优粒子方向运动的最大步长；ω是权重系数，它根据下式从ωmax到ωmin线性减小：
其中，itermax是最大迭代次数，iter为当前迭代次数，ωmax＝1ωmin＝0。
3
CN 107630697 A
说 明 书
1/6页
基于随钻电磁波电阻率测井的地层电阻率联合反演方法
技术领域
[0001]本发明涉及油气藏开发技术领域，具体地指一种基于随钻电磁波电阻率测井的地层电阻率联合反演方法，其是利用随钻电磁波电阻率测井仪器所测量的相位差信号和幅度比信号联合反演地层电阻率的一种方法。
背景技术
[0002]作为主要的随钻测井技术之一，随钻电磁波电阻率测井在随钻地质导向及储层含油饱和度评价方面发挥了至关重要的作用。相对于其它直接测量地层电阻率的测井仪器，随钻电磁波电阻率测井实际测量记录的为幅度比和相位差两个信号，因此在进行地层电阻率评价时需要将幅度比和相位差信号反演成地层电阻率。[0003]目前，利用随钻电磁波电阻率测井资料反演地层电阻率时，通常采用最小二乘及模拟退火等算法。其中，最小二乘算法属于局部寻优算法，反演效果很大程度上依赖于模型参数的初始精度，目标函数容易陷入局部极值，而模拟退火算法虽然能够进行全局寻优，但依然没有完全摆脱初始模型参数的影响。由于不同地区井眼地层电阻率变化较大，在进行地层电阻率反演时，很多情况下无法准确获取较合适的初始模型参数。[0004]除了反演算法，反演目标函数的约束方式也对反演结果有较大影响。利用随钻电磁波电阻率测井资料反演地层电阻率时，通常在反演目标函数中对幅度比或相位差信号单独进行约束，从而反演得到幅度比电阻率或相位差电阻率。理论研究及实际应用结果表明，幅度比电阻率探测深度较深，但纵向分辨率较差且仅在低阻地层具有较好精度，而相位差电阻率纵向分辨率较好，但径向深度又相对较浅。
[0005]因此单独利用幅度比电阻率和相位差电阻率均无法较好反应地层真实电阻率信息，而组合利用幅度比电阻率和相位差电阻率时又存在纵向分辨率和径向探测深度不匹配的问题，这些问题都严重束缚了随钻电磁波电阻率测井的应用效果。发明内容
[0006]本发明现有技术存在的缺陷，提供了一种基于随钻电磁波电阻率测井的地层电阻率联合反演方法，该反演方法利用随钻电磁波电阻率测井幅度比和相位差信号联合约束反演目标函数，该方法具有全局寻优特性，不需要初始值或先验信息。[0007]为实现上述目的，本发明提供的一种基于随钻电磁波电阻率测井的地层电阻率联合反演方法，包括以下步骤：
[0008]1)获取随钻电磁波电阻率测井幅度比和相位差数据；[0009]2)根据幅度比和相位差数据进行加权求和，构建联合约束的反演目标函数E；其中，由于随钻电磁波测井幅度比和相位差信号幅值往往相差一个数量级以上，因此如果在目标函数中直接联合约束反演，幅值较大的数据会占主导作用，从而达不到联合约束反演的目的。
[0010]3)基于联合约束的反演目标函数，利用粒子群算法反演得到地层电阻率。
4
CN 107630697 A[0011]
说 明 书
2/6页
进一步地看，所述步骤1)中，幅度比和相位差数据为实际测量数据和正演模拟数
据，其中，幅度比和相位差数据包括正演模型幅度比信号、实测幅度比信号、正演模型相位差信号和实测相位差信号。[0012]再进一步地，所述步骤2)中，联合约束的反演目标函数为：[0013]E＝||EATTS-EATTM||2+||ΔφS-ΔφM||2·cor2[0014]其中，EATTS为正演模型幅度比信号，EATTM为实测幅度比信号，ΔφS为正演模型相位差信号，ΔφM为实测相位差信号，
[0015]cor为幅度比和相位差信号加权系数，表达式为：
[0016]
再进一步地，所述步骤3)中，利用粒子群算法反演得到地层电阻率的步骤：[0018]步骤S1、获取反演目标函数E，设定粒子数量Num、粒子分布范围F、最大迭代次数itermax及目标函数阈值ε；[0019]步骤S2、获取给定范围内均匀分布的粒子初始位置及随机分布的粒子初始运动速度；
[0020]步骤S3、计算初始粒子群个体极值P0best及全局极值N0best，初始模型参数X0＝N0best；
[0021]步骤S4、根据粒子运动速度及位置更新表达式，获得更新的粒子运动速度及更新位置坐标；
[0022]步骤S5、以步骤4)获得的粒子运动速度及更新位置，确定当前位置，粒子群个体极值Pbest及全局极值Nbest，当前模型参数X＝Nbest；[0023]步骤S6、将初始模型参数X0和当前模型参数X带入反演目标函数E中，分别得到目标函数E的模型参数初始值E(X0)和当前值E(X)；[0024]步骤S7、求取目标函数的模型参数初始值和当前值的差值ΔE＝E(X0)-E(X)，即为判断目标函数的模型参数初始值和当前值的差值是否符合预设规则，当目标函数的当前值是否小于或等于目标函数阈值ε，若是，则执行步骤S8，否则，返回执行步骤S4重新获取更新的粒子运动速度及位置；[0025]步骤S8、将当前模型参数X赋值给初始模型参数X0，然后执行步骤S9；[0026]步骤S9、判断迭代次数iter是否与所述最大迭代次数itermax相等，若是相等，则执行步骤S10，否则，返回步骤S2；[0027]步骤S10、输出初始模型参数X0，作为反演得到的地层电阻率值。[0028]利用粒子群算法对获取的幅度比和相位差数据进行多次联合约束反演，将反演得到的地层电阻率平均值作为最终的反演结果。[0029]粒子群算法(Particle Swarm Optimization，PSO)是由Kennedy和Eberhar(1995年)受鸟群觅食过程中个体和群体间的信息传递的启发提出的。当鸟群在特定区域内搜索唯一的食物时，所有的鸟都不知道食物在那里，而找到食物的最优策略就是搜寻目前离食物最近的鸟的周围区域。
[0030]基于上述鸟群觅食过程，粒子群算法将每个优化问题的解都看作是所搜寻空间中
5
[0017]
CN 107630697 A
说 明 书
3/6页
的一只鸟，称为“粒子”。每个粒子还有一个速度决定它们每一次迭代所飞翔的方向和距离，而在每一次迭代中，粒子通过跟踪个体极值和全局极值来来更新自己，最终获得使目标函数全局最小的最优解。[0031]再进一步地，所述步骤S4)中，粒子运动速度及位置更新表达式分别为：
[0032][0033][0034]
其中，分别为迭代更新前后第i个粒子运动速度；和分别为迭代更
新前后第i个粒子位置；c1＝2和c2＝2是加速常量，分别用于调节向全局最优粒子和个体最优粒子方向运动的最大步长；ω是权重系数，它根据下式从ωmax到ωmin线性减小：
[0035]
其中，itermax是最大迭代次数，iter为当前迭代次数，ωmax＝1，ωmin＝0。[0037]本发明的有益效果在于：
[0038]1)在本发明实施例的技术方案中，所提供的地层电阻率反演方法根据[0039]随钻电磁波电阻率测井仪器测量数据，建立了幅度比和相位差加权联合约束的反演目标函数，并通过粒子群算法反演得到地层电阻率。
[0040]2)与传统的利用幅度比和相位差单独约束的反演方法相比，本发明所提供的地层电阻率反演方法得到的地层电阻率曲线兼顾了传统幅度比电阻率曲线较好的径向探测深度及相位差电阻率曲线较好的纵向分辨率的优势，更加充分的利用了随钻电磁波电阻率测井信息，而反演工作量较传统反演方法缩减了一半。
[0041]3)粒子群算法与通常所采用的最小二乘法或模拟退火算法不同，它不仅搜索局部最优值，还搜索全局最优值，因而可以跳出局部极值。另外，不需要地层电阻率先验信息，也不依赖于模型参数的初始精度，只需要设置较为合理的所搜范围，即可反演得到较高精度的地层电阻率。
附图说明
[0042]图1为本发明实施实例所提供的地层电阻率反演流程示意图；[0043]图2为本发明实施实例所提供的粒子群算法的流程示意图；
[0044]图3为本发明实施实例所提供的随钻电磁波测井线圈系结构示意图；[0045]图4a为模型1中500kHz发射频率0.8m源距线圈系径向探测特性示意图；[0046]图4b为模型1中500kHz发射频率1m源距线圈系径向探测特性示意图；[0047]图4c为模型1中2MHz发射频率0.8m源距线圈系径向探测特性示意图；[0048]图4d为模型1中2MHz发射频率1m源距线圈系径向探测特性示意图；[0049]图5a为模型2中500kHz发射频率0.8m源距线圈系纵向分辨率示意图；[0050]图5b为模型2中500kHz发射频率1m源距线圈系纵向分辨率示意图；[0051]图5c为模型2中2MHz发射频率0.8m源距线圈系纵向分辨率示意图；[0052]图5d为模型2中2MHz发射频率1m源距线圈系纵向分辨率示意图。
[0036]
6
CN 107630697 A
说 明 书
4/6页
具体实施方式
[0053]为了更好地解释本发明，以下结合具体实施例进一步阐明本发明的主要内容，但本发明的内容不仅仅局限于以下实施例。[0054]实施例1
[0055]模型1为侵入地层模型，从仪器轴向沿径向方向将地层划分为侵入带和原状地层，为了计算不同侵入状况下随钻电磁波测井仪器响应，模型侵入带电阻率Rin和原状地层电阻率Rt分别为2.5Ω.m和10Ω.m，侵入深度计算范围为0.1至6m。为了反演得到侵入地层模型电阻率值，如图1～2所示：具体地层电阻率反演步骤如下：
[0056]S1)将上述侵入地层模型中随钻电磁波电阻率测井仪器响应值作为实际测量数据，而将均匀地层模型中仪器响应值作为正演模型数据。[0057]S2)根据幅度比和相位差数据进行加权求和，构建联合约束的反演目标函数E；[0058]E＝||EATTS-EATTM||2+||ΔφS-ΔφM||2·cor2[0059]其中，EATTS为正演模型幅度比信号，EATTM为实测幅度比信号，EΔφ为相位差反演目标函数，ΔφS为正演模型相位差信号，ΔφM为实测相位差信号，[0060]cor为幅度比和相位差信号加权系数，表达式为：
[0061]
S3)基于联合约束的反演目标函数，利用粒子群算法反演得到地层电阻率；其具体步骤如下：
[0063]步骤S301、获取反演目标函数E，设定粒子数量Num、粒子分布范围F、最大迭代次数itermax及目标函数阈值ε；[0064]步骤S302、获取给定范围内均匀分布的粒子初始位置及随机分布的粒子初始运动速度；
[0065]步骤S303、计算初始粒子群个体极值P0best及全局极值N0best，初始模型参数X0＝N0best；
[0066]步骤S304、根据粒子运动速度及位置更新表达式，获得更新的粒子运动速度及更新位置坐标；粒子运动速度及位置更新表达式分别为：
[0062]
[0067][0068]
其中，分别为迭代更新前后第i个粒子运动速度；和分别为迭代更
新前后第i个粒子位置；c1＝2和c2＝2是加速常量，分别用于调节向全局最优粒子和个体最优粒子方向运动的最大步长；ω是权重系数，它根据下式从ωmax到ωmin线性减小：
[0069][0070][0071]
其中，itermax是最大迭代次数，iter为当前迭代次数，ωmax＝1，ωmin＝0；
步骤S305、以步骤4)获得的粒子运动速度及更新位置，确定当前位置，粒子群个体
7
CN 107630697 A
说 明 书
5/6页
极值Pbest及全局极值Nbest，当前模型参数X＝Nbest；[0072]步骤S306、将初始模型参数X0和当前模型参数X带入目标函数中，分别得到反演目标函数E的模型参数初始值E(X0)和当前值E(X)；[0073]步骤S307、求取目标函数的模型参数初始值和当前值的差值ΔE＝E(X0)-E(X)，即为判断目标函数的模型参数初始值和当前值的差值是否符合预设规则，当目标函数的当前值是否小于或等于目标函数阈值ε，若是，则执行步骤S308，否则，返回执行步骤S304重新获取更新的粒子运动速度及位置；[0074]步骤S308、将当前模型参数X赋予初始模型参数X0，然后执行步骤S309；[0075]步骤S309、判断迭代次数iter是否与所述最大迭代次数itermax相等，若是相等，则执行步骤S310，否则，返回步骤S2；[0076]步骤S310、输出初始模型参数X0，作为反演得到的地层电阻率值Ra。[0077]根据径向积分几何因子Gr定义，Gr等于0.5时所对应的径向侵入深度值即为随钻电磁波电阻率测井仪器径向探测深度，其计算表达式如下：
[0078]
式中，Ra为随钻电磁波测井反演所得的地层电阻率值，将传统反演方法得到的幅度比电阻率和相位差电阻率及本发明中联合反演方法得到联合电阻率分别带入上式，便可计算得到不同的径向积分几何因子曲线。如图4所示，针对不同的线圈系结构及发射频率，利用径向积分几何因子Gr等于0.5所对应的径向探测深度值，发现线圈源距越大且发射频率越低时，仪器探测深度越大；在固定的线圈系源距及发射频率下，幅度比电阻率对应的径向探测深度最大，联合反演电阻率对应的径向探测深度次之，而相位差电阻率对应的径向探测深度最浅。结果表明，本发明中联合反演方法得到的视电阻率曲线所探测的径向深度居于传统反演方法所获得的幅度比电阻率和相位差电阻率径向探测深度之间，具有较好的径向探测特性。[0080]实施例2：纵向分辨率分析
[0081]本实施例的方法与实施例1的方法基本相同，不同之处在于：[0082]模型2为纵向分层的地层模型，背景地层电阻率和夹层电阻率分别为10Ω.m和1Ω.m，模型中背景地层与夹层交替出现，地层厚度范围为0.2至2m，根据传统反演方法得到的幅度比和相位差电阻率曲线及本发明中联合反演方法得到的视电阻率曲线如图5所示。如图所示，按照传统曲线半幅点划分地层厚度的方法，发现线圈源距越小且发射频率越高时，视电导率曲线纵向分辨率越好；在固定的线圈系源距及发射频率下，传统的相位差反演得到的视电导率曲线纵向分辨率最好，联合反演得到的视电导率曲线纵向分辨率与相位差反演得到的视电导率曲线纵向分辨率十分接近，而传统的幅度比反演得到的视电导率曲线纵向分辨率则明显较前两种方法差些。[0083]结果表明，本实施例中所提出的联合反演方法所获得的地层视电导率曲线兼顾了传统幅度比电阻率径向探测深度和相位差电阻率纵向分辨率的优势，能够比传统的单一视电阻率曲线更好的评价地层。
[0084]其它未详细说明的部分均为现有技术。尽管上述实施例对本发明做出了详尽的描述，但它仅仅是本发明一部分实施例，而不是全部实施例，人们还可以根据本实施例在不经
8
[0079]
CN 107630697 A
说 明 书
6/6页
创造性前提下获得其他实施例，这些实施例都属于本发明保护范围。
9
CN 107630697 A
说 明 书 附 图
1/6页
图1
10
CN 107630697 A
说 明 书 附 图
2/6页
图2
图3
11
CN 107630697 A
说 明 书 附 图
3/6页
图4a
图4b
12
CN 107630697 A
说 明 书 附 图
4/6页
图4c
图4d
13
CN 107630697 A
说 明 书 附 图
5/6页
图5a
图5b
14
CN 107630697 A
说 明 书 附 图
6/6页
图5c
图5d
15
因篇幅问题不能全部显示，请点此查看更多更全内容