CN116818842B - 油井地层电导率信息的获取方法 - Google Patents
油井地层电导率信息的获取方法 Download PDFInfo
- Publication number
- CN116818842B CN116818842B CN202311098855.8A CN202311098855A CN116818842B CN 116818842 B CN116818842 B CN 116818842B CN 202311098855 A CN202311098855 A CN 202311098855A CN 116818842 B CN116818842 B CN 116818842B
- Authority
- CN
- China
- Prior art keywords
- axial
- discrete
- conductivity
- electric field
- measuring points
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000003129 oil well Substances 0.000 title claims abstract description 35
- 230000005684 electric field Effects 0.000 claims abstract description 93
- 239000002184 metal Substances 0.000 claims abstract description 23
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 15
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 39
- 238000005755 formation reaction Methods 0.000 claims description 14
- 230000035945 sensitivity Effects 0.000 claims description 11
- 238000005259 measurement Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 3
- 239000004020 conductor Substances 0.000 claims 2
- 230000035699 permeability Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/041—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
本发明提供了一种地层电导率信息的获取方法,包括:S1、对金属套管内轴向的不同测点电磁场进行测量,得到个测点的轴向实测电场;计算接地导线源在所有离散圆环中心处激发的轴向初始电场;S2、建立离散圆环内轴向电场之间相互关系→形成第二类弗雷德方程组→得到个离散圆环内轴向电场;S3、得到套管内个测点轴向预测电场;S4、建立数据拟合项,结合光滑约束项建立反演目标函数;S5、计算数据拟合误差并进行判断:若数据拟合误差小于设定阈值,则输出层状模型电导率数据作为油井地层电导率信息,结束本次的油井地层电导率信息的获取;若数据拟合误差大于设定阈值则对电导率值进行更新,再基于更新后的电导率值返回S2。
Description
技术领域
本发明涉及电磁法油气勘探技术领域,尤其涉及一种油井地层电导率信息的获取方法。
背景技术
测井技术是一种地球物理手段,通过使用各种传感器记录数据,获取石油、煤及金属矿体开采过程中井下岩层的物理参数,从而服务与石油与矿产开采。
传统的电磁法测井将发射源与传感器全部置于金属套管内,电磁信号受到油气开采井金属套管严重屏蔽,如何从微弱电磁信号数据获取地层电导率信息一直是油气测井工程难题。
研究发现,有限元等理论计算方法可以获取金属套管对电磁测井数据的影响,从而去除套管对电磁场测井干扰,但计算效率低下无法满足快速获取油井地层电导率信息的需求。
地表接地长导线源具有大功率、信号衰减慢的特点,能够有效加强井内电磁传感器信号强度;同时,油井地层结构相对简单,基于层状模型对油井地层结构进行近似,使用矩量法计算层状模型背景下套管内电磁测井信号具有高效性。
利用地表接地导线源信号强、油井地层结构简单、矩量法高效的特点,提出一种获取油井地层电导率信息的套管内电磁测井方法及其数据解释方法以解决现有开发油气井电磁测井信号屏蔽与数据解释困难。
发明内容
本发明目的在于提供一种获取油井地层电导率信息的套管内电磁测井方法及其数据解释方法,旨在解决现有油田电磁测井装置信号微弱及数据解释困难的问题。
本发明提供了一种油井地层电导率信息的获取方法,包括以下步骤:
步骤一、基于布设于地表的接地导线源以及置于金属套管内的电磁传感器,对金属套管内的不同测点的轴向电磁场进行测量,得到个测点的轴向实测电场/>;
设定油井地层为层水平均匀层状模型并给出各水平均匀地层的顶部厚度与深度以及将每层水平均匀层状模型赋予电导率/>,将金属套管剖分为/>个离散圆环,将接地导线源离散设置为/>个电偶极子源,使用解析法计算接地导线源在所有离散圆环中心处激发的轴向初始电场/>;
其中:取大于等于1的自然数;/>取大于等于1的自然数;/>取大于等于1的自然数,且每个离散圆环均不跨越相邻两个地层;/>取大于等于1的自然数;
步骤二、基于离散圆环中心处激发的轴向初始电场建立离散圆环内轴向电场之间相互关系;
基于离散圆环内轴向电场之间相互关系形成第二类弗雷德方程组;
基于第二类弗雷德方程组求解方程组得到个离散圆环内轴向电场/>;
步骤三、基于个离散圆环内轴向电场/>插值得到套管内/>个测点轴向预测电场/>;
步骤四、基于个测点的轴向实测电场/>和套管内/>个测点轴向预测电场建立数据拟合项,结合光滑约束项建立反演目标函数/>;
步骤五、基于个测点的轴向实测电场/>和套管内/>个测点轴向预测电场计算数据拟合误差/>,并对数据拟合误差/>进行判断:若数据拟合误差小于设定阈值,则输出层状模型电导率数据作为油井地层电导率信息,结束本次的油井地层电导率信息的获取;若数据拟合误差/>大于设定阈值则进入步骤六;
步骤六、逐一改变每层水平均匀层状模型赋予电导率,计算并存贮因每层水平均匀层状模型赋予电导率/>改变而导致各测点轴向预测电场产生扰动的扰动值,形成测点轴向预测电场对每层水平均匀层状模型赋予电导率/>的灵敏度矩阵;
步骤七、对反演目标函数使用泰勒级数展开,忽略2次以上高次项,使用灵敏度矩阵分别表达反演目标函数/>的梯度矩阵/>以及反演目标函数/>的海森矩阵/>;
步骤八、基于反演目标函数的海森矩阵/>和反演目标函数/>的梯度矩阵/>,得到层状模型电导率更新量/>,从而将步骤一中每层水平均匀层状模型赋予电导率/>进行更新,得到更新后的电导率/>;
步骤九、基于更新后的电导率返回步骤二。
可选的,得到个测点的中心轴向实测电场/>的具体过程如下:
在地表布设接地导线源并通以固定频率电流,将电磁传感器置于金属套管内,并沿金属套管的延伸方向进行移动同时对金属套管内不同测点的轴向电磁场进行测量,得到个测点的中心轴向实测电场/>。
可选的,所述步骤二中建立离散圆环内轴向电场之间相互关系的具体过程如下:
忽略离散圆环的中心点径向电场;
对离散圆环中心处激发的轴向初始电场使用格林函数建立离散圆环内轴向电场之间相互关系。
可选的,所述步骤二中第二类弗雷德方程组的表达式如下:
;
其中,为/>数组,代表/>个离散圆环的轴向电场;/>为/>数组,代表个离散圆环的轴向初始电场;/>为/>数组,代表各离散圆环与其所在地层电导率差异;/>为/>矩阵,代表/>个离散圆环轴向电场之间的相互作用。
可选的,所述步骤二中所述的矩阵各元素具体为:
,
,
;
当不同单元之间相互作用时:
当 />;
当同一单元内部作用时,对单元进行二次剖分,得到:
当/>;
其中,为离散圆环的内半径,/>为离散圆环的外半径,/>为作用的离散圆环的编号,/>为被作用的离散圆环的编号,/>为/>号离散圆环中心的纵坐标,/>为/>号离散圆环中心的纵坐标,/>为/>号离散圆环所在地层的电导率,/>为积分变量,/>为单位虚数,/>为发射电流频率,/>为真空磁导率,/>为离散圆环的长度,/>为一阶贝塞尔函数,/>为波数。
可选的,所述步骤三中的由离散圆环电场插值得到所有测点预测电场的表达式如下:
;
其中,为/>个套管内测点轴向预测电场;/>为/>矩,其由三次样条插值函数决定。
可选的,所述步骤四中反演目标函数的表达式如下:
;
其中,为/>矩阵,代表/>个套管内测点轴向实测电场;/>为正则化参数,用于平衡数据拟合项/>与模型光滑项/>;/>是/>阶的差分矩阵;/>为/>阶矩阵,代表/>个地层电导率参数。
可选的,所述步骤五中的数据拟合误差的表达式如下:
。
可选的,所述步骤七中的反演目标函数的海森矩阵/>的表达式为:
;
所述梯度矩阵的表达式为:
;
其中,为灵敏度矩阵的转秩。
与现有技术相比,本发明具有以下有益效果:
本发明提出的技术方案能在油井金属套管内开展电磁法测井工作,所提装置能够从金属套管内测量反映地层电导率信息的电磁场信号,所采用矩量法对层状模型下套管内轴向电磁场进行快速数值模拟,所使用反演手段能够从所测量测井电磁场信号中获取地层电导率信息。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明实施例中一种油井地层电导率信息的获取方法的整体流程示意图;
图2(a)是本发明实施例中测井装置形式示意图;
图2(b)是本发明实施例中测井装置理论模型离散示意图;
图3是本发明实施例中单个离散圆环几何参数示意图;
图4(a)是本发明实施例中基于三次样条插值获取的第一种层状模型下套管内轴向预测电场的曲线;
图4(b)是本发明实施例中基于三次样条插值获取的第二种层状模型下套管内轴向预测电场的曲线;
图5是本发明实施例中反演过程迭代次数与数据拟合差曲线图;
图6是本发明实施例中理论油井地层模型与反演解释结果对比图。
具体实施方式
以下是本发明的具体实施例并结合附图,对本发明的技术方案作进一步的描述,但本发明并不限于这些实施例。
除非另有定义,本文所使用的所有的技术和科学术语与本发明的技术领域技术人员通常理解的含义相同。本文中在本发明说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。
本实施例:
本实施案例中首先对油井测井装置形式进行描述,其中发射源位于地表,电磁传感器位于金属套管内部;使用接地导线源在油井地层内产生激发电流,将电磁传感器沿油气开发井移动,从而实现对套管内电磁场进行测量;将油井地层近似为水平均匀层状模型,在给定各地层厚度后对每个地层给定初始电导率;将层状介质中套管剖分为离散圆环,使用解析计算层状模型背景下,接地导线在套管内形成的初始电流,使用矩量法实现套管内电场正演;以套管内实测电场为数据拟合目标,基于套管内电场正演计算各地层电导率对套管内电场的灵敏函数,使用线性化反演理论,迭代更改层状模型电导率信息,最终获取油井各地层电导率。
具体的,本发明所提供的一种油井地层电导率信息的获取方法,具体是指一种加强油井的电磁测井信号的装置形式与快速数据解释方法;参见图1所示,包括以下步骤:
步骤一、参见图2(a)所示,在地表布设接地导线源并通以固定频率电流,将电磁传感器置于金属套管内,并沿金属套管的高度方向进行移动同时对金属套管内轴向的不同测点电磁场进行测量,得到(/>取大于等于1的自然数)个测点的轴向实测电场/>;
参见图2(b)所示,设定油井地层为(/>取大于等于1的自然数)层水平均匀层状模型并给出各水平均匀地层的顶部厚度与深度以及将每层水平均匀层状模型赋予电导率值/>;将金属套管剖分为/>(/>取大于等于1的自然数,且每个离散圆环均不跨越相邻两个地层)个离散圆环;将接地导线源离散设置为/>(/>取大于等于1的自然数)个电偶极子源(具体的,本实施例中优选将接地导线源离散设置为4个电偶极子源),使用解析法(具体设置为格林函数/>)计算接地导线源在所有离散圆环中心处激发的轴向初始电场;
;
其中,为/>数组,代表地表发射源在/>个离散圆环的轴向初始电场;/>为矩阵,代表地表发射导线所离散的4个电偶极子各段偶极矩大小,/>为/>的格林函数矩阵,代表发射源离散得到的4个单位电偶极子源在/>个离散圆环产生的轴向初始电场。
步骤二、参见图3所示,忽略离散圆环内径向电场,认为圆环内轴向电场处处相等,对离散圆环中心处激发的轴向初始电场使用格林函数建立离散圆环内轴向电场之间相互关系,使用矩量法形成第二类弗雷德方程组(只针对离散圆环的轴向电场建立),使用直接解法求解方程组得到离散圆环内轴向电场;
进一步地,所述第二类弗雷德方程组的表达式如下:
;
其中,为/>数组,代表/>个离散圆环的轴向初始电场;/>为/>数组,代表各离散圆环与其所在地层电导率差异;/>为/>矩阵,代表/>个离散圆环轴向电场之间的相互作用。
进一步地,矩阵各元素具体为:
,
,
;
当不同单元之间相互作用时:
当/>;
进一步地,当同一单元内部作用时,对单元进行二次剖分,得到:
当/>;
其中,为离散圆环的内半径,/>为离散圆环的外半径,/>为作用的离散圆环的编号,/>为被作用的离散圆环的编号,/>为/>号离散圆环中心的纵坐标,/>为/>号离散圆环中心的纵坐标,/>为/>号离散圆环所在地层的电导率,/>为积分变量,/>为单位虚数,/>为发射电流频率,/>为真空磁导率,/>为离散圆环的长度,/>为一阶贝塞尔函数,/>为波数。
步骤三、使用三次样条函数将步骤二中所获取的个离散圆环内轴向电场插值得到套管内/>个测点轴向预测电场(参见图4(a)和图4(b)展示了基于三次样条插值获取的不同电导率层状模型下套管内轴向预测电场的曲线);
进一步地,由离散圆环电场插值得到所有测点预测电场的表达式如下:
;
其中,为/>个套管内测点轴向预测电场;/>为/>矩,其由三次样条插值函数决定。
步骤四、基于个测点的轴向实测电场/>和套管内/>个测点的轴向预测电场/>建立数据拟合项,结合光滑约束项建立反演目标函数/>:
;
其中,为/>矩阵,代表/>个套管内测点轴向实测电场;/>为正则化参数,用于平衡数据拟合项/>与模型光滑项/>;/>是/>阶的差分矩阵;/>为/>阶矩阵,代表/>个地层电导率参数;
其中,为避免使用格林函数对同一圆环内部作用带来的奇异积分项,选择将圆环进一步剖分为二等份。
步骤五、基于个测点的中心轴向实测电场/>和套管内/>个测点轴向预测电场/>计算数据拟合误差/>:
;
并对数据拟合误差进行判断:若数据拟合误差/>小于设定阈值,则输出层状模型电导率数据作为油井地层电导率信息,结束本次的油井地层电导率信息的获取;若数据拟合误差/>大于设定阈值则进入步骤六。
步骤六、逐一改变每层水平均匀层状模型赋予电导率值,计算并存贮因每层水平均匀层状模型赋予电导率值/>改变而导致各测点轴向预测电场产生扰动的扰动值,形成测点轴向预测电场对每层水平均匀层状模型赋予电导率值/>的灵敏度矩阵/>;
进一步地,灵敏度矩阵的主要形式为:
。
其中:为1至/>之间的测点编号(如图2(b)所示),/>为1至/>之间的地层编号(如图2(b)所示)。
步骤七、对反演目标函数使用泰勒级数展开,忽略2次以上高次项,使用灵敏度矩阵分别表达反演目标函数/>的梯度矩阵/>和反演目标函数/>的海森矩阵/>;
具体的,反演目标函数的海森矩阵/>的表达式为:
;
其中,为灵敏度矩阵的转秩。
反演目标函数的梯度矩阵/>表达式为:
。
步骤八、使用反演目标函数的海森矩阵除以反演目标函数的梯度矩阵,得到层状模型电导率更新量:
;
从而将步骤一中每层水平均匀层状模型赋予电导率值进行更新,得到更新后的电导率值/>:
。
具体的,参见图5展示了模型更迭过程中不同电导率层状模型所预测测点轴向电场数据与实际观测轴向电场数据的拟合差变化曲线。
步骤九、基于更新后的电导率值返回步骤二。
具体地,参见图6展示了最终反演获取的地层电导率结果与实际地层模型电导率对比。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种油井地层电导率信息的获取方法,其特征在于,包括以下步骤:
步骤一、基于布设于地表的接地导线源以及置于金属套管内的电磁传感器,对金属套管内轴向的不同测点电磁场进行测量,得到M个测点的轴向实测电场Eobs;
设定油井地层为K层水平均匀层状模型并给出各水平均匀地层的顶部厚度与深度以及将每层水平均匀层状模型赋予电导率σ,将金属套管剖分为N个离散圆环,将接地导线源离散设置为P个电偶极子源,使用解析法计算接地导线源在所有离散圆环中心处激发的轴向初始电场E0;
其中:M取大于等于1的自然数;K取大于等于1的自然数;N取大于等于1的自然数,且每个离散圆环均不跨越相邻两个地层;P取大于等于1的自然数;
步骤二、基于离散圆环中心处激发的轴向初始电场建立离散圆环内轴向电场之间相互关系;
基于离散圆环内轴向电场之间相互关系形成第二类弗雷德方程组;
基于第二类弗雷德方程组求解方程组得到N个离散圆环内轴向电场E,其中,第二类弗雷德方程组的表达式如下:
E=[σα⊙E]Γ+E0;
其中,E为1×N数组,代表N个离散圆环的轴向电场;其中,E0为1×N数组,代表N个离散圆环的轴向初始电场;σα为1×N数组,代表各离散圆环与其所在地层电导率差异;Γ为N×N矩阵,代表N个离散圆环轴向电场之间的相互作用;
步骤三、基于N个离散圆环内轴向电场E插值得到套管内M个测点轴向预测电场Epre;
步骤四、基于M个测点的轴向实测电场Eobs和套管内M个测点轴向预测电场Epre建立数据拟合项,结合光滑约束项建立反演目标函数Φ;
步骤五、基于M个测点的轴向实测电场Eobs和套管内M个测点轴向预测电场Epre计算数据拟合误差RMSE,并对数据拟合误差RMSE进行判断:若数据拟合误差RMSE小于设定阈值,则输出层状模型电导率数据作为油井地层电导率信息,结束本次的油井地层电导率信息的获取;若数据拟合误差RMSE大于设定阈值则进入步骤六;
步骤六、逐一改变每层水平均匀层状模型赋予电导率σ,计算并存贮因每层水平均匀层状模型赋予电导率σ改变而导致各测点轴向预测电场产生扰动的扰动值,形成测点轴向预测电场对每层水平均匀层状模型赋予电导率σ的灵敏度矩阵;
步骤七、对反演目标函数Φ使用泰勒级数展开,忽略2次以上高次项,使用灵敏度矩阵分别表达反演目标函数Φ的梯度矩阵g以及反演目标函数Φ的海森矩阵H,其中,反演目标函数Φ的表达式如下:
其中,Eobs为1×M矩阵,代表M个套管内测点轴向实测电场;β为正则化参数,用于平衡数据拟合项[Epre-Eobs]2与模型光滑项[σW]2;W是K×K阶的差分矩阵;σ为1×K阶矩阵,代表K个地层电导率参数;
步骤八、基于反演目标函数Φ的海森矩阵H和反演目标函数Φ的梯度矩阵g,得到层状模型电导率更新量δ,从而将步骤一中每层水平均匀层状模型赋予电导率σ进行更新,得到更新后的电导率σ1;
步骤九、基于更新后的电导率σ1返回步骤二。
2.根据权利要求1所述的油井地层电导率信息的获取方法,其特征在于,得到M个测点的轴向实测电场Eobs的具体过程如下:
在地表布设接地导线源并通以固定频率电流,将电磁传感器置于金属套管内,并沿金属套管的延伸方向进行移动同时对金属套管内轴向的不同测点电磁场进行测量,得到M个测点的轴向实测电场Eobs。
3.根据权利要求1所述的油井地层电导率信息的获取方法,其特征在于,所述步骤二中建立离散圆环内轴向电场之间相互关系的具体过程如下:
忽略离散圆环内径向电场;
对离散圆环中心处激发的轴向初始电场使用格林函数建立离散圆环内轴向电场之间相互关系。
4.根据权利要求1所述的油井地层电导率信息的获取方法,其特征在于,所述步骤二中所述的Γ矩阵各元素具体为:
当不同单元之间相互作用时:
当j≠n;
进一步地,当同一单元内部作用时,对单元进行二次剖分,得到:
当j=n;
其中,a为离散圆环的内半径,b为离散圆环的外半径,n为作用的离散圆环的编号,j为被作用的离散圆环的编号,Zn为n号离散圆环中心的纵坐标,Zj为j号离散圆环中心的纵坐标,σ*为j号离散圆环所在地层的电导率,λ为积分变量,i为单位虚数,ω为发射电流频率,μ为真空磁导率,dZ为离散圆环的长度,J1为一阶贝塞尔函数,u为波数,dλ为被积分变量λ的微分符号。
5.根据权利要求1所述的油井地层电导率信息的获取方法,其特征在于,所述步骤三中的由离散圆环轴向电场插值得到所有测点预测电场的表达式如下:
Epre=CE;
其中,Epre为1×M个套管内测点的轴向预测电场;C为N×M矩,其由三次样条插值函数决定。
6.根据权利要求1所述的油井地层电导率信息的获取方法,其特征在于,所述步骤五中的数据拟合误差RMSE的表达式如下:
其中,m为1至M之间的测点编号。
7.根据权利要求1所述的油井地层电导率信息的获取方法,其特征在于,所述步骤七中的反演目标函数Φ的海森矩阵H的表达式为:
H=JTJT+βWW;
所述梯度矩阵g的表达式为:
g=JT[Epre-Eobs]+βWWσ;其中,JT为灵敏度矩阵的转秩。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311098855.8A CN116818842B (zh) | 2023-08-30 | 2023-08-30 | 油井地层电导率信息的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311098855.8A CN116818842B (zh) | 2023-08-30 | 2023-08-30 | 油井地层电导率信息的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116818842A CN116818842A (zh) | 2023-09-29 |
CN116818842B true CN116818842B (zh) | 2023-12-05 |
Family
ID=88115310
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311098855.8A Active CN116818842B (zh) | 2023-08-30 | 2023-08-30 | 油井地层电导率信息的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116818842B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5210691A (en) * | 1990-05-08 | 1993-05-11 | Schlumberger Technology Corporation | Method and apparatus for producing a more accurate resistivity log from data recorded by an induction sonde in a borehole |
US5283520A (en) * | 1991-04-04 | 1994-02-01 | Martin Philip W | Method of determining thickness of magnetic pipe by measuring the time it takes the pipe to reach magnetic saturation |
CN101004454A (zh) * | 2006-01-20 | 2007-07-25 | 中国石油天然气集团公司 | 目标最小化的三维电磁快速反演方法 |
CN105863614A (zh) * | 2015-02-10 | 2016-08-17 | 中国石油集团长城钻探工程有限公司 | 一种三维感应测井数据实时处理方法 |
CN109001823A (zh) * | 2018-04-04 | 2018-12-14 | 杭州迅美科技有限公司 | 一种电磁大地透镜探测方法和探测装置 |
CN113655542A (zh) * | 2021-08-25 | 2021-11-16 | 中国地质调查局水文地质环境地质调查中心 | 一种基于地球物理的干热岩开发阶段储层信息获取方法 |
CN115629428A (zh) * | 2022-09-23 | 2023-01-20 | 中国海洋石油集团有限公司 | 一种泥浆侵入下电阻率及核测井联合反演储层参数的方法 |
CN115795994A (zh) * | 2022-09-29 | 2023-03-14 | 西安石油大学 | 基于Unet卷积神经网络的随钻方位电磁波测井数据反演方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7049821B2 (en) * | 2003-05-29 | 2006-05-23 | Schlumberger Technology Corporation | Determination of borehole geometry inside cased wells with crosswell electromagnetics |
US8726987B2 (en) * | 2010-10-05 | 2014-05-20 | Baker Hughes Incorporated | Formation sensing and evaluation drill |
US10928543B2 (en) * | 2014-07-31 | 2021-02-23 | Halliburton Energy Services, Inc. | Electrical log processing near casing point |
US10901111B2 (en) * | 2018-06-25 | 2021-01-26 | Halliburton Energy Services, Inc. | Adaptive workflows for artifact identification in electromagnetic pipe inspection |
-
2023
- 2023-08-30 CN CN202311098855.8A patent/CN116818842B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5210691A (en) * | 1990-05-08 | 1993-05-11 | Schlumberger Technology Corporation | Method and apparatus for producing a more accurate resistivity log from data recorded by an induction sonde in a borehole |
US5283520A (en) * | 1991-04-04 | 1994-02-01 | Martin Philip W | Method of determining thickness of magnetic pipe by measuring the time it takes the pipe to reach magnetic saturation |
CN101004454A (zh) * | 2006-01-20 | 2007-07-25 | 中国石油天然气集团公司 | 目标最小化的三维电磁快速反演方法 |
CN105863614A (zh) * | 2015-02-10 | 2016-08-17 | 中国石油集团长城钻探工程有限公司 | 一种三维感应测井数据实时处理方法 |
CN109001823A (zh) * | 2018-04-04 | 2018-12-14 | 杭州迅美科技有限公司 | 一种电磁大地透镜探测方法和探测装置 |
CN113655542A (zh) * | 2021-08-25 | 2021-11-16 | 中国地质调查局水文地质环境地质调查中心 | 一种基于地球物理的干热岩开发阶段储层信息获取方法 |
CN115629428A (zh) * | 2022-09-23 | 2023-01-20 | 中国海洋石油集团有限公司 | 一种泥浆侵入下电阻率及核测井联合反演储层参数的方法 |
CN115795994A (zh) * | 2022-09-29 | 2023-03-14 | 西安石油大学 | 基于Unet卷积神经网络的随钻方位电磁波测井数据反演方法 |
Non-Patent Citations (2)
Title |
---|
Three dimensional interpretations of single-well electromagnetic data for geothermal applications;Hung=Wen Tseng等;《Lawrence Berkeley National Laboratory》;第1-7页 * |
金属套管对井间电磁测井的影响;栗建军等;《石油物探》;第42卷(第2期);第256-260页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116818842A (zh) | 2023-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rücker et al. | pyGIMLi: An open-source library for modelling and inversion in geophysics | |
CN104537714B (zh) | 磁共振与瞬变电磁空间约束联合反演方法 | |
AU2011248989A1 (en) | Artifact reduction in method of iterative inversion of geophysical data | |
Liu et al. | Deep learning audio magnetotellurics inversion using residual-based deep convolution neural network | |
CN113221393A (zh) | 一种基于非结构有限元法的三维大地电磁各向异性反演方法 | |
CN107305600A (zh) | 最小二乘法电阻率三维近似反演技术 | |
CN109799540B (zh) | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 | |
CN110286416B (zh) | 一种基于物性函数的快速二维密度反演方法 | |
CN104992029B (zh) | 一种多尺度非均匀月壤层内离散随机介质建模方法 | |
Peng et al. | Targeted location of microseismic events based on a 3D heterogeneous velocity model in underground mining | |
Pardo et al. | Simulations of 3D DC borehole resistivity measurements with a goal-oriented hp finite-element method. Part II: through-casing resistivity instruments | |
Göktürkler | A hybrid approach for tomographic inversion of crosshole seismic first-arrival times | |
Ho | 3-D inversion of borehole-to-surface electrical data using a back-propagation neural network | |
CN113917560B (zh) | 一种三维重磁电震多参数协同反演方法 | |
CN110119586B (zh) | 轴向电导率各向异性瞬变电磁三分量三维fdtd正演方法 | |
CN115201927A (zh) | 基于交叉梯度约束的重磁三维联合反演的方法 | |
Lei et al. | A parallel conformal symplectic Euler algorithm for GPR numerical simulation on dispersive media | |
CN116818842B (zh) | 油井地层电导率信息的获取方法 | |
Yang et al. | 3-D DC resistivity modelling with arbitrary long electrode sources using finite element method on unstructured grids | |
Zhou et al. | Determination of pore size distribution in tight gas sandstones based on Bayesian regularization neural network with MICP, NMR and petrophysical logs | |
CN112666612B (zh) | 基于禁忌搜索的大地电磁二维反演方法 | |
CN112327374B (zh) | Gpu探地雷达复杂介质dgtd正演方法 | |
Fischer et al. | A cellular automata‐based deterministic inversion algorithm for the characterization of linear structural heterogeneities | |
Tadeu et al. | 3D scattering of waves by a cylindrical irregular cavity of infinite length in a homogeneous elastic medium | |
WO2016187238A1 (en) | Auto-validating earth interpretation and modeling system |
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 |