CN116047617B - 一种井间地质特征识别方法及装置 - Google Patents
一种井间地质特征识别方法及装置 Download PDFInfo
- Publication number
- CN116047617B CN116047617B CN202310226760.3A CN202310226760A CN116047617B CN 116047617 B CN116047617 B CN 116047617B CN 202310226760 A CN202310226760 A CN 202310226760A CN 116047617 B CN116047617 B CN 116047617B
- Authority
- CN
- China
- Prior art keywords
- data
- background field
- multiple information
- wells
- plane
- 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
Images
Classifications
-
- 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/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/30—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electromagnetic waves
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种井间地质特征识别方法及装置,将井间电磁波数据表示为多重信息平面,场值分布具有明显的规律,在多重信息平面上开展数据处理和解释方法,能够充分利用场值平面的特征,在多重信息平面下对数据圆滑处理,考虑场值的变化趋势,能够有效去除噪声干扰,并保留异常的特征,提高了数据处理的效果;在多重信息平面下能够准确计算背景场,为相对衰减层析成像提供了基础;采用这种反演方法计算能够抑制虚假背景异常的产生,突出局部异常特征,从而提高解释的准确性,为进一步的井间地质特征识别提供了基础,使得识别结果更加准确。
Description
技术领域
本发明涉及地质特征识别技术领域,具体而言,涉及一种井间地质特征识别方法及装置。
背景技术
井间电磁波法是一种通过分析人工源电磁波在地下介质中的传播特性,研究地质体特征的地下物探方法。井间电磁波测量通过分布在两个钻孔中发射和接收电磁波,经过对数据处理解释后能够获得井间剖面的电磁参数分布特征。随着勘探深度和钻探成本的增加,井间电磁波测量提供的信息对于识别井间的地质特征尤为重要,因此在矿产资源勘查,工程与环境地质调查等方面具有广泛的应用。由于其特殊的观测方式,相邻观测剖面之间的场值具有明显变化的规律。井间电磁波法的一些解释方法,如射线交会法、层析成像等取得了良好的地质效果,然而没能充分利用场值的变化规律,同时还存在难以准确确定初始辐射场和天线方向因子的问题,导致反演结果的准确性和精确度降低,影响了地质解释的效果。
因此,如何提出一种准确性和精确度更高的井间地质特征识别方法,是目前亟待解决的问题。
发明内容
为了改善上述问题,本发明提供了一种井间地质特征识别方法及装置。
本发明实施例的第一方面,提供了一种井间地质特征识别方法,包括:
获取待识别地区井间电磁波的观测数据,建立多重信息平面;
在多重信息平面中对观测数据进行预处理;
在多重信息平面中进行背景场的计算,得到背景场曲线数据;
利用计算得到的背景场曲线数据和观测数据,进行相对衰减层析成像;
通过成像结果对井间地质特征进行识别。
可选地,所述建立多重信息平面的步骤,具体包括:
将不同频率、不同发射深度和不同接收深度的观测数据按照相应的三维坐标表示方式网格化得到多重信息空间,其中,x轴对应观测数据的频率,y轴对应观测数据的发射点深度,z轴对应观测数据的接收点深度;
沿着x轴切取信息空间,得到需要对应频率的信息平面。
可选地,所述在多重信息平面中对观测数据进行预处理的步骤,具体包括:
将观测数据的曲线异常按照多重信息空间的结构网格化,得到场值的平面数据;
对平面数据采用滑动平均法进行圆滑处理,得到圆滑后的平面网格数据;
最后在圆滑的平面数据上根据发射点深度和接收点深度进行插值得到圆滑后的曲线数据。
可选地,所述在多重信息平面中进行背景场的计算,得到背景场曲线数据的步骤,具体包括:
根据预先设置的梯度阈值,筛选出梯度值小于所述梯度阈值的数据点,作为计算背景场的数据点;
利用筛选出的数据点,采用趋势分析法得到背景场计算公式;
代入深度坐标得到背景场曲线数据。
可选地,所述利用计算得到的背景场曲线数据和观测数据,进行相对衰减层析成像的步骤,具体包括:
将所述背景场曲线数据和观测数据作为已知量,代入到相对衰减电磁波层析成像方法的线性方程组中;
对线性方程组求解,得到井间地下介质相对衰减的二维分布图像。
本发明实施例的第二方面,提供了一种井间地质特征识别装置,包括:
平面建立单元,用于获取待识别地区井间电磁波的观测数据,建立多重信息平面;
预处理单元,用于在多重信息平面中对观测数据进行预处理;
背景计算单元,用于在多重信息平面中进行背景场的计算,得到背景场曲线数据;
数据成像单元,用于利用计算得到的背景场曲线数据和观测数据,进行相对衰减层析成像;
特征识别单元,用于通过成像结果对井间地质特征进行识别。
可选地,所述平面建立单元,具体用于:
将不同频率、不同发射深度和不同接收深度的观测数据按照相应的三维坐标表示方式网格化得到多重信息空间,其中,x轴对应观测数据的频率,y轴对应观测数据的发射点深度,z轴对应观测数据的接收点深度;
沿着x轴切取信息空间,得到需要对应频率的信息平面。
可选地,所述预处理单元,具体用于:
将观测数据的曲线异常按照多重信息空间的结构网格化,得到场值的平面数据;
对平面数据采用滑动平均法进行圆滑处理,得到圆滑后的平面网格数据;
最后在圆滑的平面数据上根据发射点深度和接收点深度进行插值得到圆滑后的曲线数据。
可选地,所述背景计算单元,具体用于:
计算多重信息平面中平面数据的梯度值;
根据预先设置的梯度阈值,筛选出梯度值小于所述梯度阈值的数据点,作为计算背景场的数据点;
利用筛选出的数据点,采用趋势分析法得到背景场计算公式;
代入深度坐标得到背景场曲线数据。
可选地,所述数据成像单元,具体用于:
将所述背景场曲线数据和观测数据作为已知量,代入到相对衰减电磁波层析成像方法的线性方程组中;
对线性方程组求解,得到井间地下介质相对衰减的二维分布图像。
综上所述,本发明提供了一种井间地质特征识别方法及装置,将井间电磁波数据表示为多重信息平面,场值分布具有明显的规律,在多重信息平面上开展数据处理和解释方法,能够充分利用场值平面的特征,在多重信息平面下对数据圆滑处理,考虑场值的变化趋势,能够有效去除噪声干扰,并保留异常的特征,提高了数据处理的效果;在多重信息平面下能够准确计算背景场,为相对衰减层析成像提供了基础;采用这种反演方法计算能够抑制虚假背景异常的产生,突出局部异常特征,从而提高解释的准确性,为进一步的井间地质特征识别提供了基础,使得识别结果更加准确。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明实施例的井间地质特征识别方法的方法流程图;
图2为本发明实施例在多重信息平面计算出背景场数据的方法流程图;
图3为本发明实施例的井间地质特征识别方法的方法实施流程图;
图4为本发明实施例的井间地质特征识别装置的功能模块框图。
附图标记:
平面建立单元110;预处理单元120;背景计算单元130;数据成像单元140;特征识别单元150。
具体实施方式
井间电磁波测量通过在两个钻孔的不同深度点发射和接收电磁波,能够获得井间剖面的电磁参数分布。随着勘探深度和钻探成本的增加,井间电磁波测量提供的信息对于识别井间的地质特征尤为重要,因此在矿产资源勘查,工程与环境地质调查等方面具有广泛的应用。由于其特殊的观测方式,相邻观测剖面之间的场值具有明显变化的规律。井间电磁波法的一些解释方法,如射线交会法、层析成像等取得了良好的地质效果,然而没能充分利用场值的变化规律,同时还存在难以准确确定初始辐射场和天线方向因子的问题,导致反演结果的准确性和精确度降低,影响了地质特征识别的效果。
在以往的解释中,初始辐射场通常直接凭经验取值或者由“迭代线性拟合”得到,刘立振在论文“重建地下介质相对衰减的分布:新的成像方法”中,提出了“相对衰减层析成像”的方法,回避了估算辐射参数,但对于背景场的计算并没有详细阐述。
鉴于此,本发明设计者设计了一种井间地质特征识别方法及装置,充分了利用场值平面的特征,有效去除噪声数据,并保留异常的特征,提高了数据处理的效果;实现了一种在多重信息平面准确计算背景场的方法,为相对衰减层析成像提供了基础。相对传统经验取值或者“迭代线性拟合”方法,可以有效避免直接利用初始辐射常数和方向因子计算带来的误差;相对单独利用曲线数据计算背景场,能够充分利用井间电磁波的所有观测数据,获得准确的背景场;而相比双频透射电磁波电导率层析成像方法,不受电导率大小的限制,应用范围更为广泛。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
在本发明的描述中,需要说明的是,术语“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,或者是该发明产品使用时惯常摆放的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
在本发明的描述中,还需要说明的是,除非另有明确的规定和限定,术语“设置”、“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
请参阅图1,本实施例提供的一种井间地质特征识别方法,该方法包括:
步骤S101,获取待识别地区井间电磁波的观测数据,建立多重信息平面。
多重信息平面由实际工作时观测到的场值经过计算得到,按设计好的三维框架进行重排,形成立体的数据集。
具体地,将不同频率、不同发射深度和不同接收深度的观测数据按照相应的三维坐标表示方式网格化得到多重信息空间,其中,x轴对应观测数据的频率,y轴对应观测数据的发射点深度,z轴对应观测数据的接收点深度。
其中,坐标表示方式如表1所示,将离散的数据点(x,y,z,v)(对应于实际观测数据(f,z1,z2,v))网格化得到信息空间,空间信息包含了多种频率在不同发射点和观测点的数据,信息空间不同方向的切片图代表了不同模式下场值的平面图。
表1信息空间建立的坐标对应关系
通常采用单个频率的数据进行分析,可以沿着x轴切取信息空间,得到需要对应频率的信息平面。
步骤S102,在多重信息平面中对观测数据进行预处理。
在单个频率数据中,背景场和异常场的场值平面分布都具有规律变化的特征,不是孤立的异常点,这与随机干扰噪声是不同的,因此对场值的平面数据圆滑,可以达到保留场值特征同时去除噪声干扰的目的。
具体的做法是:
将观测数据的曲线异常按照多重信息空间的结构网格化,得到场值的平面数据;
对平面数据采用滑动平均法进行圆滑处理,得到圆滑后的平面网格数据;
最后在圆滑的平面数据上根据发射点深度和接收点深度进行插值得到圆滑后的曲线数据。
作为优选的实施方式,在此基础上,还可以进一步进行处理,补偿因距离衰减产生的影响。具体的处理方式如下式所示:
f(z1,z2)=Ds+20lgr
其中Ds为观测数据,预处理后的观测数据为f(z1,z2),发射点深度为z1,接收点深度为z2,发射点和接收点之间的距离为r。
步骤S103,在多重信息平面中进行背景场的计算,得到背景场曲线数据。
在多重信息平面数据中计算出背景场,是本发明区别于上述的“相对衰减层析成像”方法核心点之一。
AΔm=Δb
其中A是与射线长度有关的矩阵,上式与初始辐射场和方向因子无关,与Δb有关,因此关键在于如何求出Δb,显然,当Δm=0,则βb=β,说明在相对衰减为零的区域,实测场强即背景场。
可以看出,实现相对衰减层析成像的一个关键步骤是计算出背景场。在该方法中,对于背景的计算是在已有的观测数据曲线基础上,利用最小二乘多项式拟合或低通数字滤波等方法估算正常曲线其余部分,得到背景场的场值。本发明与之不同的是,在多重信息平面数据中计算出背景场,在这种表示模式下背景场和异常场的特征区别更为明显,因此更易准确计算出背景场,同时还能保留场值的规律变化特征。
本发明实施例中,在多重信息平面数据中计算出背景场的方法如图2所示,具体包括:
步骤S201,计算多重信息平面中平面数据的梯度值;
步骤S202,根据预先设置的梯度阈值,筛选出梯度值小于所述梯度阈值的数据点,作为计算背景场的数据点;
作为本发明实施例的优选实施方式,为了能够更好地拟合背景场,应该尽量使用背景场区域的数据点计算,为此将满足下式的点参与最小二乘拟合计算。
步骤S203,利用筛选出的数据点,采用趋势分析法得到背景场计算公式;
步骤S204,代入深度坐标得到背景场曲线数据。
本发明实施例所采用的趋势分析法,假设背景场G(z1,z2)可以用以下形式表示:
经过上述筛选后,数据点在深度坐标为(zi,zj)的观测值为f(zi,zj),则需要使多项式曲面和实际多重信息平面在最小二乘意义下得到最佳拟合,多项式系数aij应满足:
根据多元函数求极值法,令上式的偏导数为零,得到关于系数aij的线性方程组,求解后即得到多项式系数aij,再代入到背景场G(z1,z2)的公式中,便可计算各个点上的背景场。
在构建背景场公式时,p和q为多项式的阶数。其中,从背景场的特征分析,由于发射点深度和接收点深度对背景场的影响程度是一致的,因此设多项式的阶数p=q。根据背景场的复杂程度,选择合适的p和q值。作为优选的,可令p=q=2,当p=q=2时,背景场与均匀介质中的场值平面分布图相似,后续的成像结果中能更加突出显示异常区域。
通过上述过程得到背景场G(z1,z2)的表达式后,将实际发射点z1和接收点z2的深度代入到上述背景场计算方程的公式中,计算得到背景场曲线数据R(z1,z2)。
步骤S104,利用计算得到的背景场曲线数据和观测数据,进行相对衰减层析成像。
计算得到的背景场曲线数据后,代入相对层析成像的方程,即可以实现相对衰减层析成像的计算。具体的方法是:
将所述背景场曲线数据和观测数据作为已知量,代入到相对衰减电磁波层析成像方法的线性方程组中;
对线性方程组求解,得到井间地下介质相对衰减的二维分布图像。
需要主要的是,将所述背景场曲线数据和观测数据作为已知量,代入到相对衰减电磁波层析成像方法的线性方程组中并进行求解的过程,在论文“刘立振.重建地下介质相对衰减的分布:新的成像方法.地球物理学报,1990,33(01):98-105.”中有具体的说明,在此不再进行详细的说明。
步骤S105,通过成像结果对井间地质特征进行识别。
通过得到的成像结果,从图像上识别出井间地质特征中异常的区域。
本实施例提供的井间地质特征识别方法,其具体的实施过程如图3所示,首先建立多重信息平面,对观测数据预处理,消除因距离衰减产生的影响后,将处理后的数据在多重信息平面上网格化得到平面数据。然后计算平面数据的梯度,设置梯度阈值,利用梯度值小于阈值的数据点,采用趋势分析法得到背景场计算公式,代入深度坐标得到背景场曲线数据。最后通过背景场曲线数据和观测数据进行相对衰减层析成像。通过这样的方式,一方面在多重信息平面下对数据圆滑处理,实现有效去除噪声干扰,保留异常的特征;另一方面,在多重信息平面下采用趋势分析法准确计算背景场,将计算的背景场作为相对衰减层析成像的基础,然后采用相对衰减层析成像方法计算,抑制虚假背景异常的产生,突出局部异常特征,从而提高解释的准确性。
综上,本实施例提供的井间地质特征识别方法,将井间电磁波数据表示为多重信息平面,场值分布具有明显的规律,在多重信息平面上开展数据处理和解释方法,能够充分利用场值平面的特征,在多重信息平面下对数据圆滑处理,考虑场值的变化趋势,能够有效去除噪声干扰,并保留异常的特征,提高了数据处理的效果;在多重信息平面下能够准确计算背景场,为相对衰减层析成像提供了基础;采用这种反演方法计算能够抑制虚假背景异常的产生,突出局部异常特征,从而提高解释的准确性,为进一步的井间地质特征识别提供了基础,使得识别结果更加准确。
如图4所示,本发明实施提供的井间地质特征识别装置,所述装置包括:
平面建立单元110,用于获取待识别地区井间电磁波的观测数据,建立多重信息平面;
预处理单元120,用于在多重信息平面中对观测数据进行预处理;
背景计算单元130,用于在多重信息平面中进行背景场的计算,得到背景场曲线数据;
数据成像单元140,用于利用计算得到的背景场曲线数据和观测数据,进行相对衰减层析成像;
特征识别单元150,用于通过成像结果对井间地质特征进行识别。
作为本实施例的优选实施方式,所述平面建立单元110,具体用于:
将不同频率、不同发射深度和不同接收深度的观测数据按照相应的三维坐标表示方式网格化得到多重信息空间,其中,x轴对应观测数据的频率,y轴对应观测数据的发射点深度,z轴对应观测数据的接收点深度;
沿着x轴切取信息空间,得到需要对应频率的信息平面。
作为本实施例的优选实施方式,所述预处理单元120,具体用于:
将观测数据的曲线异常按照多重信息空间的结构网格化,得到场值的平面数据;
对平面数据采用滑动平均法进行圆滑处理,得到圆滑后的平面网格数据;
最后在圆滑的平面数据上根据发射点深度和接收点深度进行插值得到圆滑后的曲线数据。
作为本实施例的优选实施方式,所述背景计算单元130,具体用于:
计算多重信息平面中平面数据的梯度值;
根据预先设置的梯度阈值,筛选出梯度值小于所述梯度阈值的数据点,作为计算背景场的数据点;
利用筛选出的数据点,采用趋势分析法得到背景场计算公式;
代入深度坐标得到背景场曲线数据。
作为本实施例的优选实施方式,所述数据成像单元150,具体用于:
将所述背景场曲线数据和观测数据作为已知量,代入到相对衰减电磁波层析成像方法的线性方程组中;
对线性方程组求解,得到井间地下介质相对衰减的二维分布图像。
本发明实施例提供的井间地质特征识别装置,用于实现上述井间地质特征识别方法,因此具体实施方式与上述方法相同,在此不再赘述。
综上所述,本发明提供了一种井间地质特征识别方法及装置,将井间电磁波数据表示为多重信息平面,场值分布具有明显的规律,在多重信息平面上开展数据处理和解释方法,能够充分利用场值平面的特征,在多重信息平面下对数据圆滑处理,考虑场值的变化趋势,能够有效去除噪声干扰,并保留异常的特征,提高了数据处理的效果;在多重信息平面下能够准确计算背景场,为相对衰减层析成像提供了基础;采用这种反演方法计算能够抑制虚假背景异常的产生,突出局部异常特征,从而提高解释的准确性,为进一步的井间地质特征识别提供了基础,使得识别结果更加准确。
在本申请所公开的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和框图显示了根据本申请的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本申请各个实施例中的各功能模块可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
Claims (8)
1.一种井间地质特征识别方法,其特征在于,包括:
获取待识别地区井间电磁波的观测数据,建立多重信息平面;
在多重信息平面中对观测数据进行预处理;
在多重信息平面中进行背景场的计算,得到背景场曲线数据;
利用计算得到的背景场曲线数据和观测数据,进行相对衰减层析成像;
通过成像结果对井间地质特征进行识别;
所述在多重信息平面中进行背景场的计算,得到背景场曲线数据的步骤,具体包括:
计算多重信息平面中平面数据的梯度值;
根据预先设置的梯度阈值,筛选出梯度值小于所述梯度阈值的数据点,作为计算背景场的数据点;
利用筛选出的数据点,采用趋势分析法得到背景场计算公式;
代入深度坐标得到背景场曲线数据;
背景场G(z1,z2)表示为:
筛选后,数据点在深度坐标为(zi,zj)的观测值为f(zi,zj),使多项式曲面和实际多重信息平面在最小二乘意义下得到最佳拟合,多项式系数aij满足:
根据多元函数求极值法,令上式的偏导数为零,得到关于系数aij的线性方程组,求解后即得到多项式系数aij,再代入到背景场G(z1,z2)的公式中,便可计算各个点上的背景场;
在构建背景场公式时,p和q为多项式的阶数,且多项式的阶数p=q;
得到背景场G(z1,z2)的表达式后,将实际发射点z1和接收点z2的深度代入到上述背景场计算方程的公式中,计算得到背景场曲线数据R(z1,z2)。
2.根据权利要求1所述的井间地质特征识别方法,其特征在于,所述建立多重信息平面的步骤,具体包括:
将不同频率、不同发射深度和不同接收深度的观测数据按照相应的三维坐标表示方式网格化得到多重信息空间,其中,x轴对应观测数据的频率,y轴对应观测数据的发射点深度,z轴对应观测数据的接收点深度;
沿着x轴切取信息空间,得到需要对应频率的信息平面。
3.根据权利要求2所述的井间地质特征识别方法,其特征在于,所述在多重信息平面中对观测数据进行预处理的步骤,具体包括:
将观测数据的曲线异常按照多重信息空间的结构网格化,得到场值的平面数据;
对平面数据采用滑动平均法进行圆滑处理,得到圆滑后的平面网格数据;
最后在圆滑的平面数据上根据发射点深度和接收点深度进行插值得到圆滑后的曲线数据。
4.根据权利要求3所述的井间地质特征识别方法,其特征在于,所述利用计算得到的背景场曲线数据和观测数据,进行相对衰减层析成像的步骤,具体包括:
将所述背景场曲线数据和观测数据作为已知量,代入到相对衰减电磁波层析成像方法的线性方程组中;
对线性方程组求解,得到井间地下介质相对衰减的二维分布图像。
5.一种井间地质特征识别装置,其特征在于,包括:
平面建立单元,用于获取待识别地区井间电磁波的观测数据,建立多重信息平面;
预处理单元,用于在多重信息平面中对观测数据进行预处理;
背景计算单元,用于在多重信息平面中进行背景场的计算,得到背景场曲线数据;
数据成像单元,用于利用计算得到的背景场曲线数据和观测数据,进行相对衰减层析成像;
特征识别单元,用于通过成像结果对井间地质特征进行识别;
所述背景计算单元,具体用于:
计算多重信息平面中平面数据的梯度值;
根据预先设置的梯度阈值,筛选出梯度值小于所述梯度阈值的数据点,作为计算背景场的数据点;
利用筛选出的数据点,采用趋势分析法得到背景场计算公式;
代入深度坐标得到背景场曲线数据;
背景场G(z1,z2)表示为:
筛选后,数据点在深度坐标为(zi,zj)的观测值为f(zi,zj),使多项式曲面和实际多重信息平面在最小二乘意义下得到最佳拟合,多项式系数aij满足:
根据多元函数求极值法,令上式的偏导数为零,得到关于系数aij的线性方程组,求解后即得到多项式系数aij,再代入到背景场G(z1,z2)的公式中,便可计算各个点上的背景场;
在构建背景场公式时,p和q为多项式的阶数,且多项式的阶数p=q;
得到背景场G(z1,z2)的表达式后,将实际发射点z1和接收点z2的深度代入到上述背景场计算方程的公式中,计算得到背景场曲线数据R(z1,z2)。
6.根据权利要求5所述的井间地质特征识别装置,其特征在于,所述平面建立单元,具体用于:
将不同频率、不同发射深度和不同接收深度的观测数据按照相应的三维坐标表示方式网格化得到多重信息空间,其中,x轴对应观测数据的频率,y轴对应观测数据的发射点深度,z轴对应观测数据的接收点深度;
沿着x轴切取信息空间,得到需要对应频率的信息平面。
7.根据权利要求6中所述的井间地质特征识别装置,其特征在于,所述预处理单元,具体用于:
将观测数据的曲线异常按照多重信息空间的结构网格化,得到场值的平面数据;
对平面数据采用滑动平均法进行圆滑处理,得到圆滑后的平面网格数据;
最后在圆滑的平面数据上根据发射点深度和接收点深度进行插值得到圆滑后的曲线数据。
8.根据权利要求7所述的井间地质特征识别装置,其特征在于,所述数据成像单元,具体用于:
将所述背景场曲线数据和观测数据作为已知量,代入到相对衰减电磁波层析成像方法的线性方程组中;
对线性方程组求解,得到井间地下介质相对衰减的二维分布图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310226760.3A CN116047617B (zh) | 2023-03-10 | 2023-03-10 | 一种井间地质特征识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310226760.3A CN116047617B (zh) | 2023-03-10 | 2023-03-10 | 一种井间地质特征识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116047617A CN116047617A (zh) | 2023-05-02 |
CN116047617B true CN116047617B (zh) | 2023-06-27 |
Family
ID=86120352
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310226760.3A Active CN116047617B (zh) | 2023-03-10 | 2023-03-10 | 一种井间地质特征识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116047617B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2707419C1 (ru) * | 2019-02-28 | 2019-11-26 | Ярослав Викторович Савватеев | Способ георадиолокационного зондирования и устройство для его осуществления |
CN110794464A (zh) * | 2019-11-07 | 2020-02-14 | 中国地质科学院地球物理地球化学勘查研究所 | 一种去除时间域航空电磁数据背景场的方法 |
CN111830578A (zh) * | 2020-08-14 | 2020-10-27 | 中国地质科学院地球物理地球化学勘查研究所 | 一种地下管线精确探测的电磁波ct方法 |
CN112363236A (zh) * | 2020-10-15 | 2021-02-12 | 中国地质大学(武汉) | 一种基于pde的重力场数据等效源延拓与数据类型转换方法 |
CN115097534A (zh) * | 2022-07-06 | 2022-09-23 | 中国地质调查局水文地质环境地质调查中心 | 干热岩压裂改造空间形态识别方法和系统 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2084930C1 (ru) * | 1993-07-22 | 1997-07-20 | Товарищество с ограниченной ответственностью "Радионда Лтд" | Способ радиоволнового межскважинного просвечивания |
US6593746B2 (en) * | 2001-08-27 | 2003-07-15 | Larry G. Stolarczyk | Method and system for radio-imaging underground geologic structures |
CN100429531C (zh) * | 2006-01-20 | 2008-10-29 | 中国石油天然气集团公司 | 目标最小化的三维电磁快速反演方法 |
US8812237B2 (en) * | 2009-02-05 | 2014-08-19 | Schlumberger Technology Corporation | Deep-reading electromagnetic data acquisition method |
AU2014391098B2 (en) * | 2014-04-16 | 2017-05-11 | Halliburton Energy Services, Inc. | Time-lapse electromagnetic monitoring |
CN106291725B (zh) * | 2015-05-13 | 2018-11-30 | 核工业北京地质研究院 | 一种快速反演地下地质体空间位置的方法 |
CN109799540B (zh) * | 2018-12-25 | 2020-12-11 | 核工业北京地质研究院 | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 |
CN111983703B (zh) * | 2020-07-24 | 2023-07-25 | 中国石油天然气集团有限公司 | 井间电磁测量流体成像方法、系统及装置 |
CN113484920B (zh) * | 2021-08-17 | 2023-05-19 | 成都理工大学 | 一种频域电磁测深资料二维结构化反演方法 |
CN114398816A (zh) * | 2022-01-18 | 2022-04-26 | 科吉思石油技术咨询(北京)有限公司 | 全三维井周速度场反演的定量解释及三维可视化方法 |
CN114629503A (zh) * | 2022-03-29 | 2022-06-14 | 安徽惠洲地质安全研究院股份有限公司 | 物探类数据有损压缩方法、分析方法及压缩分析系统 |
-
2023
- 2023-03-10 CN CN202310226760.3A patent/CN116047617B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2707419C1 (ru) * | 2019-02-28 | 2019-11-26 | Ярослав Викторович Савватеев | Способ георадиолокационного зондирования и устройство для его осуществления |
CN110794464A (zh) * | 2019-11-07 | 2020-02-14 | 中国地质科学院地球物理地球化学勘查研究所 | 一种去除时间域航空电磁数据背景场的方法 |
CN111830578A (zh) * | 2020-08-14 | 2020-10-27 | 中国地质科学院地球物理地球化学勘查研究所 | 一种地下管线精确探测的电磁波ct方法 |
CN112363236A (zh) * | 2020-10-15 | 2021-02-12 | 中国地质大学(武汉) | 一种基于pde的重力场数据等效源延拓与数据类型转换方法 |
CN115097534A (zh) * | 2022-07-06 | 2022-09-23 | 中国地质调查局水文地质环境地质调查中心 | 干热岩压裂改造空间形态识别方法和系统 |
Non-Patent Citations (1)
Title |
---|
Imaging complex geology through challenging surface terrain – a case study from West China;Huiwen Xie等;《First Break》;第33卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116047617A (zh) | 2023-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6885943B2 (en) | Simultaneous resolution enhancement and dip correction of resistivity logs through nonlinear iterative deconvolution | |
US9046509B2 (en) | Method and system for estimating rock properties from rock samples using digital rock physics imaging | |
EP3111041B1 (en) | Automatic method for three-dimensional structural interpretation of borehole images acquired in high-angle and horizontal wells | |
US20160223704A1 (en) | Workflow adjustment methods and systems for logging operations | |
US11230922B2 (en) | Fracture interpretation with resistivity and sonic logs in biaxial anisotropic formations | |
US9105075B1 (en) | Enhancing seismic features using an optical filter array | |
EP0930519B1 (en) | Method and apparatus for evaluating the resistivity of invaded formations at high apparent dip angle | |
US20220341835A1 (en) | Methods and Systems for Determining Reservoir and Fracture Properties | |
US11307323B2 (en) | Methods and systems to analyze bed boundary detection | |
US8600115B2 (en) | Borehole image reconstruction using inversion and tool spatial sensitivity functions | |
Bitenc et al. | Evaluation of wavelet denoising methods for small-scale joint roughness estimation using terrestrial laser scanning | |
US9829591B1 (en) | Determining seismic stratigraphic features using a symmetry attribute | |
WO2022056147A1 (en) | Methods and apparatus for determining dielectric constant and resistivity with electromagnetic propagation measurements | |
EP3217354B1 (en) | Seismic modeling system providing seismic survey data frequency domain inpainting and related methods | |
CN112709567A (zh) | 用电磁测量值确定各向异性的地下性质 | |
US20220237891A1 (en) | Method and system for image-based reservoir property estimation using machine learning | |
WO2018183223A1 (en) | Methods and systems of determining parameters characterizing porous media from data gathered by a plurality of different tools | |
CN116047617B (zh) | 一种井间地质特征识别方法及装置 | |
US9869798B2 (en) | Method for estimating a logging tool response in a layered formation | |
Gryazin et al. | Numerical solution of a subsurface imaging inverse problem | |
US20220012853A1 (en) | Joint minimum entropy method for simultaneous processing and fusion of multi-physics data and images | |
Zhang et al. | An iterative algorithm for conductivity image reconstruction from crosswell EM measurements | |
CN112132372B (zh) | 用于神经网络的层测井信息无损处理方法和装置 | |
GB2495634A (en) | Resolution matched nonlinear resolution enhancement of well logs | |
CN116299701B (zh) | 断溶体储层识别方法、装置及电子设备 |
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 |