CN112327383B - 基于梯度张量特征值的位场数据地质构造位置圈定方法 - Google Patents
基于梯度张量特征值的位场数据地质构造位置圈定方法 Download PDFInfo
- Publication number
- CN112327383B CN112327383B CN202011336514.6A CN202011336514A CN112327383B CN 112327383 B CN112327383 B CN 112327383B CN 202011336514 A CN202011336514 A CN 202011336514A CN 112327383 B CN112327383 B CN 112327383B
- Authority
- CN
- China
- Prior art keywords
- geological structure
- gradient
- gradient tensor
- field data
- potential field
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000002159 abnormal effect Effects 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000000926 separation method Methods 0.000 claims abstract description 5
- 230000000694 effects Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于梯度张量特征值的位场数据地质构造位置圈定方法,包括如下步骤:S1:获取实测位场数据,采用基于向上延拓的异常剥离法实现区域场和局部场的异常分离,并采用空间域‑频率域导数计算方法和拉普拉斯方程计算梯度分量;S2:利用所述梯度分量建立梯度张量矩阵;S3:计算所述梯度张量矩阵的特征值;S4:利用所述特征值定义地质构造位置圈定方法。该基于梯度张量特征值的位场数据地质构造位置圈定方法,具有更高的水平分辨率,圈定出的断裂构造更加清晰,可有效地圈定出深浅地质构造对应的边界位置,提升了地质构造解释精度。
Description
技术领域
本发明涉及地球科学技术领域,特别提供了一种基于梯度张量特征值的位场数据地质构造位置圈定方法。
背景技术
位场数据相比其他地球物理测量数据具有更高的水平分辨率,在海洋、陆地区域地质构造划分中得到了广泛地应用,为区域构造位置圈定、断裂识别提供了一种重要的手段。基于位场数据建立高精度的地质构造位置圈定方法可更好地服务于海洋和陆地资源的勘探,进一步指导其它地球物理方法的精细勘探和解释。现在存在大量的文献和工作进行了地质构造位置圈定方法的研究,在一定假设条件下获得了较好的地质构造水平位置圈定效果,但在多地质体空间叠加异常的处理中会产生虚假边界位置、精度低等问题。针对大面积、高精度探测获取的位场数据,发展地质构造位置圈定方法是现在研究的重点方向,如何进一步提高地质构造位置圈定精度、分辨率以及可靠性,成为亟待解决的问题。
发明内容
鉴于此,本发明的目的在于针对上述现有技术的不足,提供一种基于梯度张量特征值的位场数据地质构造位置圈定方法,以提高地质构造位置圈定精度、分辨率以及可靠性。
本发明提供的技术方案是:基于梯度张量特征值的位场数据地质构造位置圈定方法,包括如下步骤:
S1:获取实测位场数据,采用基于向上延拓的异常剥离法实现区域场和局部场的异常分离,并采用空间域-频率域导数计算方法和拉普拉斯方程计算梯度分量Gxx、Gxy、Gxz、Gyx、Gyy、Gyz、Gzx、Gzy、Gzz;
S2:利用所述梯度分量建立梯度张量矩阵G,其中,所述梯度张量矩阵G如公式(1)所示:
S3:利用公式(2)和公式(3)计算所述梯度张量矩阵G的特征值λ1和λ2:
λ1=A+B (2)
I1=GxxGyy+GyyGzz+GxxGzz-Gxy 2-Gxz 2-Gyz 2;
I2=Gxx(GyyGzz-Gyz 2)+Gxy(GyzGxz-GxyGzz)+Gxz(GxyGyz-GxzGyy);
S4:利用所述特征值λ1和λ2定义地质构造位置圈定方法NL1和NL2,其中:
本发明提供的基于梯度张量特征值的位场数据地质构造位置圈定方法,基于梯度张量特征值建立归一化导数地质构造圈定方法,综合利用了多方向梯度分量对地质构造的空间位置刻画能力,定义了归一化导数计算方法可有效地均衡深浅地质构造引起的位场异常,在正负位场异常同时存在的条件下同样能够很好地进行地质构造空间位置的圈定。在实际数据处理过程中,采用定义的地质构造圈定方法可直接利用极大值位置圈定出构造边界位置,相比传统方法具有更高的水平分辨率和精度。
本发明提供的基于梯度张量特征值的位场数据地质构造位置圈定方法,将包含全部梯度分量的特征值进行构造位置圈定方法的定义,提高了构造位置圈定精度。本发明突破了传统地质构造位置圈定方法采用部分梯度分量进行地质构造位置圈定,定义了包含全部梯度分量的特征值的地质构造圈定方法,在多地质体叠加异常复杂条件下具有更高的空间分辨率和精度,避免了虚假边界位置的产生,提高了地质构造空间位置圈定的可靠性,为分析海洋和陆地构造演化规律和油气勘探远景区圈定提供技术支撑。
附图说明
下面结合附图及实施方式对本发明作进一步详细的说明:
图1为本发明提供的基于梯度张量特征值的位场数据地质构造位置圈定方法的流程图;
图2为实施例中吉林某地区实测位场数据异常图;
图3a为地质构造位置圈定方法NL1的地质构造位置圈定结果图;
图3b为地质构造位置圈定方法NL2的地质构造位置圈定结果图;
图3c为Thdr地质构造位置圈定结果图;
图3d为Theta地质构造位置圈定结果图。
具体实施方式
下面将结合具体的实施方案对本发明进行进一步的解释,但并不局限本发明。
如图1所示,本发明提供了一种基于梯度张量特征值的位场数据地质构造位置圈定方法,包括如下步骤:
S1:获取实测位场数据,采用基于向上延拓的异常剥离法实现区域场和局部场的异常分离,并采用空间域-频率域导数计算方法和拉普拉斯方程计算梯度分量Gxx、Gxy、Gxz、Gyx、Gyy、Gyz、Gzx、Gzy、Gzz;
S2:利用所述梯度分量建立梯度张量矩阵G,其中,所述梯度张量矩阵G如公式(1)所示:
S3:利用公式(2)和公式(3)计算所述梯度张量矩阵G的特征值λ1和λ2:
λ1=A+B (2)
I1=GxxGyy+GyyGzz+GxxGzz-Gxy 2-Gxz 2-Gyz 2;
I2=Gxx(GyyGzz-Gyz 2)+Gxy(GyzGxz-GxyGzz)+Gxz(GxyGyz-GxzGyy);
S4:利用所述特征值λ1和λ2定义地质构造位置圈定方法NL1和NL2,其中:
实施例:
以吉林某地区的实测位场异常数据为例,通过基于梯度张量特征值的位场数据地质构造位置圈定方法进行断裂构造空间展布的确定,其中,图2为该地区实测位场数据异常图。
具体步骤如下:
S1:获取该地区实测位场数据,采用基于向上延拓的异常剥离法实现区域场和局部场的异常分离,并采用空间域-频率域导数计算方法和拉普拉斯方程计算梯度分量Gxx、Gxy、Gxz、Gyx、Gyy、Gyz、Gzx、Gzy、Gzz;
S2:利用所述梯度分量建立梯度张量矩阵G,其中,所述梯度张量矩阵G如公式(1)所示:
S3:计算所述梯度张量矩阵G的特征值λ1和λ2:
λ1=A+B (2)
I1=GxxGyy+GyyGzz+GxxGzz-Gxy 2-Gxz 2-Gyz 2;
I2=Gxx(GyyGzz-Gyz 2)+Gxy(GyzGxz-GxyGzz)+Gxz(GxyGyz-GxzGyy);
S4:采用归一化导数构造位置圈定方法NL1进行实测位场数据的地质构造边界位置圈定,其中,参数k根据异常特征选择为0.1,图3a为该方法的应用效果,NL1表达式为:
采用归一化导数构造位置圈定方法NL2进行实测位场数据的地质构造边界位置圈定,其中,参数k根据异常特征选择为0.1,图3b为该方法的应用效果,NL2表达式为:
对比例:
为了分析本方法的地质构造圈定精度和分辨率,选择现有的Thdr和Theta图进行方法效果对比,其中,Thdr的表达式为:
Theta的表达式为:
图3c为Thdr图地质构造位置圈定结果,图3d为Theta图地质构造位置圈定结果。
经对比可以看出,相比传统的地质构造圈定方法,本发明提供的方法具有更高的水平分辨率,圈定出的断裂构造更加清晰,本发明提供的方法可有效地圈定出深浅地质构造对应的边界位置,提升了地质构造解释精度。上面结合附图对本发明的实施方式做了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (1)
1.基于梯度张量特征值的位场数据地质构造位置圈定方法,其特征在于,包括如下步骤:
S1:获取实测位场数据,采用基于向上延拓的异常剥离法实现区域场和局部场的异常分离,并采用空间域-频率域导数计算方法和拉普拉斯方程计算梯度分量Gxx、Gxy、Gxz、Gyx、Gyy、Gyz、Gzx、Gzy、Gzz;
S2:利用所述梯度分量建立梯度张量矩阵G,其中,所述梯度张量矩阵G如公式(1)所示:
S3:利用公式(2)和公式(3)计算所述梯度张量矩阵G的特征值λ1和λ2:
λ1=A+B (2)
I1=GxxGyy+GyyGzz+GxxGzz-Gxy 2-Gxz 2-Gyz 2;
I2=Gxx(GyyGzz-Gyz 2)+Gxy(GyzGxz-GxyGzz)+Gxz(GxyGyz-GxzGyy);
S4:利用所述特征值λ1和λ2定义地质构造位置圈定方法NL1和NL2,其中:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011336514.6A CN112327383B (zh) | 2020-11-25 | 2020-11-25 | 基于梯度张量特征值的位场数据地质构造位置圈定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011336514.6A CN112327383B (zh) | 2020-11-25 | 2020-11-25 | 基于梯度张量特征值的位场数据地质构造位置圈定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112327383A CN112327383A (zh) | 2021-02-05 |
CN112327383B true CN112327383B (zh) | 2021-12-31 |
Family
ID=74309541
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011336514.6A Expired - Fee Related CN112327383B (zh) | 2020-11-25 | 2020-11-25 | 基于梯度张量特征值的位场数据地质构造位置圈定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112327383B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108508490A (zh) * | 2018-03-07 | 2018-09-07 | 吉林大学 | 一种基于解析信号的磁张量梯度数据均衡边界识别方法 |
CN110967745A (zh) * | 2018-09-29 | 2020-04-07 | 中国石油化工股份有限公司 | 一种火成岩的深度域速度建模方法 |
CN111007571A (zh) * | 2019-11-28 | 2020-04-14 | 吉林大学 | 一种基于三维构造张量的航磁数据地质体边界识别方法 |
CN111123380A (zh) * | 2019-12-18 | 2020-05-08 | 山东大学 | 基于重磁梯度数据张量不变量的目标深度估计方法及*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008154679A1 (en) * | 2007-06-18 | 2008-12-24 | Commonwealth Scientific And Industrial Research Organisation | Method and apparatus for detection using magnetic gradient tensor |
-
2020
- 2020-11-25 CN CN202011336514.6A patent/CN112327383B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108508490A (zh) * | 2018-03-07 | 2018-09-07 | 吉林大学 | 一种基于解析信号的磁张量梯度数据均衡边界识别方法 |
CN110967745A (zh) * | 2018-09-29 | 2020-04-07 | 中国石油化工股份有限公司 | 一种火成岩的深度域速度建模方法 |
CN111007571A (zh) * | 2019-11-28 | 2020-04-14 | 吉林大学 | 一种基于三维构造张量的航磁数据地质体边界识别方法 |
CN111123380A (zh) * | 2019-12-18 | 2020-05-08 | 山东大学 | 基于重磁梯度数据张量不变量的目标深度估计方法及*** |
Non-Patent Citations (1)
Title |
---|
全张量重力梯度数据的综合分析与处理解释;袁园;《中国博士学位论文全文数据库 基础科学辑》;20150815;正文第37-42、72页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112327383A (zh) | 2021-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2573467C (en) | Computer-based method for while-drilling modeling and visualization of layered subterranean earth formations | |
CN111007571B (zh) | 一种基于三维构造张量的航磁数据地质体边界识别方法 | |
EP3105420B1 (en) | Log processing and fracture characterization in biaxially anisotropic formations | |
WO2016138874A1 (zh) | 位场构造格架自动提取方法 | |
JP2013527926A (ja) | 地球物理学的データセット中の異常検出のための窓型統計的分析 | |
OA11145A (en) | Method for automatic detection of planar heterogeneities crossing the stratification of an environment | |
CA2740636A1 (en) | Windowed statistical analysis for anomaly detection in geophysical datasets | |
CN111221048B (zh) | 基于跨孔电阻率ct多尺度反演的孤石边界识别与成像方法 | |
Shafiabadi et al. | Identification of reservoir fractures on FMI image logs using Canny and Sobel edge detection algorithms | |
US9354353B2 (en) | Dip correction in invaded zones | |
US20220237891A1 (en) | Method and system for image-based reservoir property estimation using machine learning | |
Srivardhan | Stratigraphic correlation of wells using discrete wavelet transform with fourier transform and multi-scale analysis | |
CN112327383B (zh) | 基于梯度张量特征值的位场数据地质构造位置圈定方法 | |
CN113466947A (zh) | 一种应用于超导瞬变电磁的坏点自动去除方法 | |
CN112180447B (zh) | 一种消除储层强反射屏蔽的方法及*** | |
WO2021163582A1 (en) | Systems and methods for generating continuous grain size logs from petrographic thin section images | |
EP3136132B1 (en) | Direct resistivity determination | |
Fan et al. | An elaborately designed virtual frame to level aeromagnetic data | |
Kassenaar | An application of principal components analysis to borehole geophysical data | |
NO343122B1 (no) | Analyse av geologiske strukturer basert på seismiske attributter | |
Al-Hadithi et al. | Using source parameter imaging technique to the aeromagnetic data to estimate the basement depth of Tharthar Lake and surrounding area in Central Iraq | |
Wang et al. | Wavelet Transform Based Morphological Matching Area Selection for Underwater Gravity Gradient-Aided Navigation | |
CN113552639B (zh) | 物探数据边界异常增强方法 | |
Fu et al. | Accelerated Geophysical Inversion for Airborne Transient Electromagnetic Data Using GPU | |
CN115797652A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211231 |
|
CF01 | Termination of patent right due to non-payment of annual fee |