WO2021243967A1 - 一种三维电阻率层析成像方法及系统 - Google Patents
一种三维电阻率层析成像方法及系统 Download PDFInfo
- Publication number
- WO2021243967A1 WO2021243967A1 PCT/CN2020/129416 CN2020129416W WO2021243967A1 WO 2021243967 A1 WO2021243967 A1 WO 2021243967A1 CN 2020129416 W CN2020129416 W CN 2020129416W WO 2021243967 A1 WO2021243967 A1 WO 2021243967A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- dimensional
- data
- resistivity
- module
- coordinate
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 100
- 238000003325 tomography Methods 0.000 title claims abstract description 27
- 238000001514 detection method Methods 0.000 claims abstract description 41
- 230000004927 fusion Effects 0.000 claims abstract description 24
- 238000006243 chemical reaction Methods 0.000 claims abstract description 23
- 238000012847 principal component analysis method Methods 0.000 claims abstract description 8
- 238000004590 computer program Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000000513 principal component analysis Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000001174 ascending effect Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims 1
- 230000002547 anomalous effect Effects 0.000 abstract 1
- 238000010291 electrical method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000001154 acute effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000002247 constant time method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
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
- 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
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
- G01V2210/665—Subsurface modeling using geostatistical modeling
- G01V2210/6652—Kriging
Definitions
- the invention belongs to the field of geophysical prospecting, and in particular relates to a three-dimensional resistivity tomography imaging method and system.
- the reflection signal of the geological radar high-frequency electromagnetic wave carries the dielectric constant information of the formation, which has strong resolving power but limited detection depth; the high-density electrical method obtains a large amount of data in a single detection and responds well to high-resistance anomalies, but it is affected by the terrain.
- the ups and downs have a greater impact.
- the first aspect of the present invention provides a three-dimensional resistivity tomography method, which can compare resistivity data obtained by at least two geophysical prospecting methods to select resistivity data points with the same exploration area coordinates , And based on the principal component analysis method to perform data fusion on the data points obtained by at least two geophysical methods, and finally realize the formation of a three-dimensional model through three-dimensional coordinate conversion, improve the detection accuracy, and make the detection results have good intuitiveness and visibility sex.
- a three-dimensional electrical resistivity tomography method including:
- the resistivity data obtained after the fusion is converted into a three-dimensional coordinate to form a three-dimensional model.
- the second aspect of the present invention provides a three-dimensional electrical resistivity tomography system.
- a three-dimensional electrical resistivity tomography system including:
- a two-dimensional resistivity data acquisition module which is used to separately explore regions containing geological anomalies by adopting at least two exploration methods to obtain two-dimensional resistivity data of the corresponding detection plane;
- the same coordinate data point extraction module which is used to unify the coordinate system of the resistivity data points obtained by all exploration methods, and extract data points with the same coordinates;
- a data fusion module which is used to perform data fusion on the extracted resistivity data at the same position using principal component analysis
- the three-dimensional conversion module is used to perform three-dimensional coordinate conversion on the resistivity data obtained after fusion to form a three-dimensional model.
- the third aspect of the present invention provides a computer-readable storage medium.
- a computer-readable storage medium has a computer program stored thereon, and when the program is executed by a processor, the steps in the above-mentioned three-dimensional resistivity tomography method are realized.
- the fourth aspect of the present invention provides a computer device.
- a computer device includes a memory, a processor, and a computer program stored on the memory and capable of running on the processor.
- the processor implements the above-mentioned three-dimensional resistivity tomography method when the program is executed. A step of.
- the present invention compares the resistivity data obtained by at least two geophysical prospecting methods to select resistivity data points with the same exploration area coordinates, and performs data fusion on the data points obtained by the at least two geophysical prospecting methods based on the principal component analysis method Finally, the formation of a three-dimensional model is realized through three-dimensional coordinate conversion, which improves the detection accuracy, and at the same time makes the detection results have good intuitiveness and visibility.
- the three-dimensional coordinate conversion imaging method of the present invention has good visibility and can aggregate the cross-sectional data of multiple two-dimensional planes into a three-dimensional model, which can very intuitively reflect the reality of abnormal objects in the detection area. Circumstances, it is also convenient for later interpretation, analysis and guidance.
- Figure 1 is a flow chart of a three-dimensional electrical resistivity tomography method according to an embodiment of the present invention
- FIG. 2 is a schematic diagram of the structure of a three-dimensional electrical resistivity tomography system according to an embodiment of the present invention
- Fig. 3(a) is the coordinate conversion mode under the first line survey mode within the tunnel range of the embodiment of the present invention
- Fig. 3(b) is the coordinate conversion mode in the second way of surveying the line in the tunnel area according to the embodiment of the present invention
- Fig. 3(c) is the coordinate conversion mode in the third way of surveying the line in the tunnel range according to the embodiment of the present invention.
- Fig. 3(d) is the coordinate conversion mode under the fourth line survey mode within the tunnel range according to the embodiment of the present invention.
- the component analysis method performs data fusion on the extracted resistivity data at the same position; the resistivity data obtained after the fusion is converted into a three-dimensional coordinate to form a three-dimensional model.
- the geological data first analyze the geological data to roughly determine the location of geological anomalies such as water-filled caves and water-conducting faults, and then conduct geophysical exploration of the geological anomaly area, that is, cross-hole resistivity CT detection and high-density electrical detection.
- the two-dimensional resistivity cross-section distribution map of the cross-hole resistivity detection and the two-dimensional resistivity cross-section distribution map of the high-density electrical method detection are obtained by the geophysical inversion method.
- the coordinate system is unified, so that the coordinates of the data points obtained by the two methods for the same location in the area are the same.
- extract the data points with the same coordinates in the exploration area from the data obtained by the high-density electrical method and the data obtained by the cross-hole method, and use the method based on principal component analysis for data fusion.
- the resistivity data of multiple two-dimensional sections are obtained. Then, the obtained data points in the two-dimensional coordinate system are converted to the three-dimensional coordinate system, and the data points in multiple three-dimensional coordinates are collectively imaged through the Kriging interpolation method to form a three-dimensional model.
- the cross-hole resistivity can be improved through data fusion.
- the CT method has a poor effect near the hole, and can accurately locate the low-resistance abnormal body adjacent to the electrode, reducing the cross-hole resistance.
- the distribution range and number of high-resistance pseudo-abnormal bodies in the hole attachment; for the side-by-side abnormal bodies, through data fusion, the problem of poor cross-hole resistivity CT imaging in the horizontal direction and high-density electrical method in the vertical direction can be improved, so that The separation of the two geoelectrical anomalies makes the image more visible.
- the two complement each other and improve the image interpretation capability of the ERT method.
- the three-dimensional electrical resistivity tomography method process includes:
- the resistivity data points detected by one of the exploration methods Take the resistivity data points detected by one of the exploration methods as the benchmark, and transform the coordinates of the resistivity data points detected by the other exploration method. For example: Taking the resistivity data points obtained by the cross-hole method as a reference, the high-density method resistivity data points are transformed into coordinates. Because in the field detection, the position of the starting point of the survey line of cross-hole detection and high-density electrical detection may be different, the resistivity data points obtained by the two methods are located at the coordinate origin of the coordinate system, which leads to two methods The initial coordinates are different.
- the invention uses the data detected by the cross-hole method as a reference, and adjusts the position of the coordinate origin of the data point coordinate system of the high-density electrical resistivity to keep it consistent with the coordinate system where the cross-hole method is located.
- First carry out data centering that is, standardize the resistivity sample data obtained by the cross-hole method and the high-density electrical method to eliminate errors caused by different dimensions, self-variation or large differences in values;
- the sample points of the centralized data are projected onto the eigenvector base with the largest eigenvalue, and the result of the fused resistivity data can be considered as the comprehensive nature of the cross-hole detection data and the high-density detection data.
- R 1 is the position matrix
- R 2 is the data point matrix
- X, Y, Z are the final three-dimensional coordinates with O as the origin
- X 1 is the horizontal distance from the starting point of the survey line to the origin of the coordinates
- Y 1 is the distance coordinate from the starting point of the survey line
- the longitudinal depth of the origin, Z 1 is the vertical height from the starting point of the survey line to the origin of the coordinates
- X' is the horizontal length of the original data point
- Y' is the detection depth of the original data point.
- the direction of the survey line is the positive direction of the coordinate axis, Y'is usually negative, ⁇ is the angle between the tangent of the survey line along the XOZ plane and the Z axis (acute angle), and ⁇ is the angle between the horizontal direction and the longitudinal direction of the survey line (take Acute angle).
- the acquired resistivity data points in the three-dimensional rectangular coordinate system under multiple planes are kriging interpolation to form a three-dimensional model with a smooth curved surface.
- the three-dimensional electrical resistivity tomography system of this embodiment includes:
- a two-dimensional resistivity data acquisition module which is used to separately explore areas containing geological anomalies using at least two exploration methods to obtain two-dimensional resistivity data of the corresponding detection plane;
- geological analysis is performed to determine the geological anomaly area.
- geological and hydrological phenomena such as geological structure, topography, overburden, karst cave development, water-conducting faults, etc. in the area.
- the coordinate system of the resistivity data points obtained by one survey method is used as a reference, and the resistivity data points obtained by other survey methods are coordinate transformed.
- a data fusion module which is used to perform data fusion on the extracted resistivity data at the same position by using principal component analysis
- the data fusion module includes:
- Data centralization module which is used to standardize the resistivity data obtained by all exploration methods
- the covariance matrix solving module is used to solve the covariance matrix between the variables after centralization, and to measure whether the deviation change trends of the two variables are consistent;
- the eigenvalue and eigenvector solution module is used to find the eigenvalues and eigenvectors of the covariance, arrange the eigenvalues in ascending order, select the largest one among them, and find the corresponding eigenvector;
- the data projection module is used to project the sample points of the centralized data onto the eigenvector base with the largest eigenvalue to obtain the fused resistivity data result.
- a three-dimensional conversion module which is used to perform three-dimensional coordinate conversion on the resistivity data obtained after fusion to form a three-dimensional model.
- the converted three-dimensional coordinates form a three-dimensional model by kriging interpolation.
- the resistivity data obtained by at least two geophysical prospecting methods are compared to select resistivity data points with the same exploration area coordinates, and the data points obtained by the at least two geophysical prospecting methods are data fused based on the principal component analysis method, and finally Through the conversion of three-dimensional coordinates, the formation of three-dimensional models is realized, the detection accuracy is improved, and the detection results have good intuitiveness and visibility.
- This embodiment provides a computer-readable storage medium on which a computer program is stored, and when the program is executed by a processor, the steps in the three-dimensional resistivity tomography method described in the first embodiment are implemented.
- the resistivity data obtained by at least two geophysical prospecting methods are compared to select resistivity data points with the same exploration area coordinates, and the data points obtained by the at least two geophysical prospecting methods are data fused based on the principal component analysis method, and finally Through the conversion of three-dimensional coordinates, the formation of three-dimensional models is realized, the detection accuracy is improved, and the detection results have good intuitiveness and visibility.
- This embodiment provides a computer device, including a memory, a processor, and a computer program stored in the memory and capable of running on the processor.
- a computer program as described in the first embodiment is implemented. Steps in the three-dimensional electrical resistivity tomography method.
- the resistivity data obtained by at least two geophysical prospecting methods are compared to select resistivity data points with the same exploration area coordinates, and the data points obtained by the at least two geophysical prospecting methods are data fused based on the principal component analysis method, and finally Through the conversion of three-dimensional coordinates, the formation of three-dimensional models is realized, the detection accuracy is improved, and the detection results have good intuitiveness and visibility.
- the embodiments of the present invention can be provided as a method, a system, or a computer program product. Therefore, the present invention may adopt the form of a hardware embodiment, a software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may be in the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, optical storage, etc.) containing computer-usable program codes.
- These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing equipment to work in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture including the instruction device.
- the device implements the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
- These computer program instructions can also be loaded on a computer or other programmable data processing equipment, so that a series of operation steps are executed on the computer or other programmable equipment to produce computer-implemented processing, so as to execute on the computer or other programmable equipment.
- the instructions provide steps for implementing functions specified in a flow or multiple flows in the flowchart and/or a block or multiple blocks in the block diagram.
- the program can be stored in a computer readable storage medium. During execution, it may include the processes of the above-mentioned method embodiments.
- the storage medium may be a magnetic disk, an optical disc, a read-only memory (Read-Only Memory, ROM), or a random access memory (Random Access Memory, RAM), etc.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Claims (10)
- 一种三维电阻率层析成像方法,其特征在于,包括:采用至少两种勘探方式分别对含有地质异常体区域进行勘探,获得相应探测平面的二维电阻率数据;统一所有勘探方式获得的电阻率数据点的坐标系,提取坐标相同的数据点;利用主成分分析法对提取的相同位置的电阻率数据进行数据融合;将融合后得到的电阻率数据进行三维坐标转换形成三维模型。
- 如权利要求1所述的三维电阻率层析成像方法,其特征在于,利用主成分分析法对提取的相同位置的电阻率数据进行数据融合的过程为:数据中心化,即所有勘探方式获得的电阻率数据进行标准化;求解中心化后变量间的协方差矩阵,衡量两个变量的偏差变化趋势是否一致;求出协方差的特征值和特征向量,将特征值按照从小到大的顺序排列,选择其中最大的一个,求出其对应的特征向量;将中心化数据的样本点投影到特征值最大的特征向量基上,得到的融合后的电阻率数据结果。
- 如权利要求1所述的三维电阻率层析成像方法,其特征在于,以一种勘探方式获得的电阻率数据点的坐标系为基准,对其他勘探方式获得的电阻率数据点进行坐标转化。
- 如权利要求1所述的三维电阻率层析成像方法,其特征在于,转换后的三维坐标通过克里金插值法形成三维模型。
- 一种三维电阻率层析成像系统,其特征在于,包括:二维电阻率数据获取模块,其用于采用至少两种勘探方式分别对含有地质 异常体区域进行勘探,获得相应探测平面的二维电阻率数据;相同坐标数据点提取模块,其用于统一所有勘探方式获得的电阻率数据点的坐标系,提取坐标相同的数据点;数据融合模块,其用于利用主成分分析法对提取的相同位置的电阻率数据进行数据融合;三维转换模块,其用于将融合后得到的电阻率数据进行三维坐标转换形成三维模型。
- 如权利要求5所述的三维电阻率层析成像系统,其特征在于,所述数据融合模块,包括:数据中心化模块,其用于所有勘探方式获得的电阻率数据进行标准化;协方差矩阵求解模块,其用于求解中心化后变量间的协方差矩阵,衡量两个变量的偏差变化趋势是否一致;特征值和特征向量求解模块,其用于求出协方差的特征值和特征向量,将特征值按照从小到大的顺序排列,选择其中最大的一个,求出其对应的特征向量;数据投影模块,其用于将中心化数据的样本点投影到特征值最大的特征向量基上,得到的融合后的电阻率数据结果。
- 如权利要求5所述的三维电阻率层析成像系统,其特征在于,在所述相同坐标数据点提取模块中,以一种勘探方式获得的电阻率数据点的坐标系为基准,对其他勘探方式获得的电阻率数据点进行坐标转化。
- 如权利要求5所述的三维电阻率层析成像系统,其特征在于,在所述三维转换模块中,转换后的三维坐标通过克里金插值法形成三维模型。
- 一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-4中任一项所述的一种三维电阻率层析成像方法中的步骤。
- 一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-4中任一项所述的一种三维电阻率层析成像方法中的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17/431,051 US11960048B2 (en) | 2020-06-03 | 2020-11-17 | Three-dimensional electrical resistivity tomography method and system |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010495471.X | 2020-06-03 | ||
CN202010495471.XA CN111722292B (zh) | 2020-06-03 | 2020-06-03 | 一种三维电阻率层析成像方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2021243967A1 true WO2021243967A1 (zh) | 2021-12-09 |
Family
ID=72565871
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2020/129416 WO2021243967A1 (zh) | 2020-06-03 | 2020-11-17 | 一种三维电阻率层析成像方法及系统 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11960048B2 (zh) |
CN (1) | CN111722292B (zh) |
WO (1) | WO2021243967A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220090486A1 (en) * | 2020-09-24 | 2022-03-24 | Halliburton Energy Services, Inc. | Borehole localization relative to objects and subterranean formations |
CN114509822A (zh) * | 2022-01-20 | 2022-05-17 | 中铁二院工程集团有限责任公司 | 铁路隧道的地空电磁法阵列勘察方法及其测线布置方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108873074B (zh) * | 2018-04-18 | 2020-03-24 | 浙江大学 | 一种电极随机分布式高密度电阻率测量方法及勘探系统 |
CN111722292B (zh) | 2020-06-03 | 2021-06-01 | 山东大学 | 一种三维电阻率层析成像方法及系统 |
CN113552635B (zh) * | 2021-06-25 | 2023-03-03 | 山东大学 | 固源瞬变电磁数据融合隧道三维超前预报方法及系统 |
CN114397329A (zh) * | 2022-01-24 | 2022-04-26 | 中国矿业大学 | 一种煤层瓦斯含量原位动态探测系统及方法 |
CN117272213B (zh) * | 2023-11-21 | 2024-02-02 | 中南大学 | 地下污染物的地物化综合参数扫面方法、装置、设备及介质 |
CN118094438A (zh) * | 2024-04-18 | 2024-05-28 | 江西四联节能环保股份有限公司 | 一种基于大数据的配电电缆运行分析方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016018492A1 (en) * | 2014-07-30 | 2016-02-04 | Battelle Memorial Institute | Method of imaging the electrical conductivity distribution of a subsurface |
CN110361792A (zh) * | 2019-06-24 | 2019-10-22 | 山东大学 | 一种地球物理数据融合及成像方法、介质与设备 |
CN110823962A (zh) * | 2019-11-14 | 2020-02-21 | 山东大学 | 一种滑坡体的三维成像方法及系统 |
CN111123404A (zh) * | 2020-01-13 | 2020-05-08 | 安徽理工大学 | 一种地震和直流电法巷道超前探测的数据融合方法 |
CN111722292A (zh) * | 2020-06-03 | 2020-09-29 | 山东大学 | 一种三维电阻率层析成像方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6295512B1 (en) * | 1998-05-01 | 2001-09-25 | John Bryant | Subsurface mapping apparatus and method |
US10197702B2 (en) * | 2015-04-27 | 2019-02-05 | Pgs Geophysical As | Seismic guided inversion of electromagnetic survey data |
US11112516B2 (en) * | 2018-04-30 | 2021-09-07 | Schlumberger Technology Corporation | Data fusion technique to compute reservoir quality and completion quality by combining various log measurements |
-
2020
- 2020-06-03 CN CN202010495471.XA patent/CN111722292B/zh active Active
- 2020-11-17 US US17/431,051 patent/US11960048B2/en active Active
- 2020-11-17 WO PCT/CN2020/129416 patent/WO2021243967A1/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016018492A1 (en) * | 2014-07-30 | 2016-02-04 | Battelle Memorial Institute | Method of imaging the electrical conductivity distribution of a subsurface |
CN110361792A (zh) * | 2019-06-24 | 2019-10-22 | 山东大学 | 一种地球物理数据融合及成像方法、介质与设备 |
CN110823962A (zh) * | 2019-11-14 | 2020-02-21 | 山东大学 | 一种滑坡体的三维成像方法及系统 |
CN111123404A (zh) * | 2020-01-13 | 2020-05-08 | 安徽理工大学 | 一种地震和直流电法巷道超前探测的数据融合方法 |
CN111722292A (zh) * | 2020-06-03 | 2020-09-29 | 山东大学 | 一种三维电阻率层析成像方法及系统 |
Non-Patent Citations (1)
Title |
---|
YU, JIFENG ET AL.: " Chapter 8: Principal Components and Factor Analysis", GEOMATHEMATICS: METHODS AND APPLICATIONS , 28 February 2019 (2019-02-28), CN, pages 110 - 111, XP009532664, ISBN: 978-7-5646-3694-4 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220090486A1 (en) * | 2020-09-24 | 2022-03-24 | Halliburton Energy Services, Inc. | Borehole localization relative to objects and subterranean formations |
CN114509822A (zh) * | 2022-01-20 | 2022-05-17 | 中铁二院工程集团有限责任公司 | 铁路隧道的地空电磁法阵列勘察方法及其测线布置方法 |
CN114509822B (zh) * | 2022-01-20 | 2023-04-07 | 中铁二院工程集团有限责任公司 | 铁路隧道的地空电磁法阵列勘察方法及其测线布置方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111722292B (zh) | 2021-06-01 |
US11960048B2 (en) | 2024-04-16 |
CN111722292A (zh) | 2020-09-29 |
US20220308251A1 (en) | 2022-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021243967A1 (zh) | 一种三维电阻率层析成像方法及系统 | |
CN108345049B (zh) | 地下工程不良地质探测多方法约束反演与联合解释方法 | |
CN110823962B (zh) | 一种滑坡体的三维成像方法及系统 | |
Cao et al. | Automatic fracture detection based on Terrestrial Laser Scanning data: A new method and case study | |
CN106780730A (zh) | 三维地质模型的构建方法和装置 | |
US10385658B2 (en) | In-situ wellbore, core and cuttings information system | |
CN111221048B (zh) | 基于跨孔电阻率ct多尺度反演的孤石边界识别与成像方法 | |
CN111812730B (zh) | 一种用于滑坡探测的电阻率数据融合三维成像方法及系统 | |
US20230184983A1 (en) | Vector-resistivity-based real-time advanced detection method for water-bearing hazard body | |
CN111986321B (zh) | 三维地质建模方法、装置、设备及存储介质 | |
CN108873073B (zh) | 一种基于网络并行电法的三维跨孔电阻率层析成像方法 | |
US20140095078A1 (en) | Method and system for presenting seismic information | |
CN110073246B (zh) | 与质量控制有关的改进的方法 | |
BR112018017328B1 (pt) | Método e sistema para caracterização de uma formação subterrânea | |
CN112433253A (zh) | 一种动水软弱破碎围岩隧道三维探测方法 | |
CN105220694B (zh) | 一种工程抛石量估算方法 | |
CN113077546B (zh) | 一种面向矿产资源储量估算的搜索椭球体自动设置方法 | |
CN114820969B (zh) | 一种三维地质模型构建方法 | |
CN116484670A (zh) | 一种基于综合物探的三维地质建模与受力分析方法及系统 | |
CN115577616A (zh) | 基于深度学习的碳酸岩缝洞地震刻画方法及装置 | |
US20220090486A1 (en) | Borehole localization relative to objects and subterranean formations | |
US20210381365A1 (en) | Borehole localization relative to objects and subterranrean formations | |
CN105279305B (zh) | 一种地面三维激光扫描技术测站选取方法 | |
CN106951644A (zh) | 基于gocad的软基砂土液化范围确定方法 | |
CN112347901A (zh) | 一种基于三维激光扫描技术岩体块体分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20939241 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 20939241 Country of ref document: EP Kind code of ref document: A1 |
|
32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 22/06/2023) |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 20939241 Country of ref document: EP Kind code of ref document: A1 |