CN109343134A - A kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system - Google Patents
A kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system Download PDFInfo
- Publication number
- CN109343134A CN109343134A CN201811421881.9A CN201811421881A CN109343134A CN 109343134 A CN109343134 A CN 109343134A CN 201811421881 A CN201811421881 A CN 201811421881A CN 109343134 A CN109343134 A CN 109343134A
- Authority
- CN
- China
- Prior art keywords
- data
- depth
- standard deviation
- measuring point
- transient electromagnetic
- 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.)
- Pending
Links
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/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to a kind of data analysis interpretation method and system, belong to geological exploration in coal mine technical field, are specifically related to a kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system.Data Analysis Services of the invention are established on the basis of the apparent resistivity and depth data that initial data is calculated, and sense datum is needed not rely on.Using standard deviation standardized analysis method, the case where apparent resistivity isoline is in concentric circles in TA data are effectively had modified, the analyticity of data is improved.
Description
Technical field
The present invention relates to a kind of data analysis interpretation method and system, belong to geological exploration in coal mine technical field, specifically
It is related to a kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system.
Background technique
Transient electromagnetic method is a kind of geophysical exploration method that the electromagnetic field using artificial field source excitation is detected,
With electromagnetic coupling theory for basic detection principle.According to the difference of working space, aviation, ground, Mine transient electromagnetic can be divided into
Method.Mine Transient Electromagnetic Method be mainly used in underground exploration, wherein coal mine down-hole tunnel roof and floor, heading end forward probe and
Working face detection etc. is widely used.Its detection process mainly has the following steps: using the rectangular loop line of earth-free multiturn as
Emission source feeds steady current thereto, excites primary field in surrounding space；Electric current is sometime cut off, primary field disappears therewith
It loses, the drastic change of primary field excites induced current in surrounding medium and generates secondary field；This secondary field data is received using probe,
And the electrical property feature of surrounding medium is analyzed with this data.
To Transient Electromagnetic Method in Mine data carry out computation of apparent resistivity and when convert deeply, it is disconnected to obtain depth-apparent resistivity
Face.Because of working environment complexity, the reasons such as device influences or ambient field influences cause apparent resistivity pseudosection feature complicated and changeable, hold
Difficulty easily is caused to analysis interpretation work, to need to carry out conversion process to data convenient for recognizing apparent resistance section.
The Chinese invention patent of Publication No. CN103344999A discloses a kind of utilization ratio approach elimination transient electromagnetic
The method of ambient field in data, the method which uses are as follows: rule of thumb judge the direction that ambient field is distributed with survey line, it is right
Data press regular grid interpolation along X, Y-direction, and statistical average is sought to survey line in empirically determined direction, will be on survey line
The data and average value of every bit make ratio, repeat above step to all surveys line, then draw, judge mapping result and experience
Judge whether direction is consistent, judges whether that the ratio for needing to synthesize both direction is drawn again.The invention needs to establish to be sentenced in experience
It is higher to explanation personnel's skill requirement on the basis of disconnected, and it be easy to cause erroneous judgement.
The Chinese patent of Publication No. CN103278855 discloses a kind of elimination tunnel and landform to direct-current exploration apparent resistance
The method that rate influences, specifically includes and establishes model according to roadway shape and cutting depth, hypsography elevation etc.；By exploration area
Resistivity of the minimum value of layer resistivity as the homogeneous half space containing tunnel and with landform, by current electrode according to direct-current exploration
The type of device and cloth pole mode of use, place corresponding position in a model, and the every movement of current electrode is primary, carry out one time three
Dimensional Simulation of Air Motion obtains the potential value of ambient field, then interpolation obtains potential difference of the potential difference between receiving electrode as ambient field,
It is located at current potential formula in horizontal homogeneous half space from point electric current source, derives apparent resistivity formula, actual measurement potential difference is subtracted
As the potential difference of ambient field, the potential difference of no tunnel and the influence of topography is obtained, substitutes into the half space view of current source promise breaking underground
Resistivity formula obtains relative apparent resistivity.For the invention is mainly for direct-current exploration method, discussion is grounding requirement
Under, the influence in tunnel and landform, and problem with grounding is not present in transient electromagnetic exploration, in addition this method needs per treatment are done primary
Three-dimensional simulation, calculation amount are huge, it is difficult to meet the requirement of fast calculation analysis.Therefore, intelligence in the prior art is supervised
Control system improves, and is current technical problem in the urgent need to address to meet the needs of different application scene.
Summary of the invention
A brief summary of one or more aspects is given below to provide to the basic comprehension in terms of these.This general introduction is not
The extensive overview of all aspects contemplated, and be both not intended to identify critical or decisive element in all aspects also non-
Attempt to define the range in terms of any or all.Its unique purpose is to provide the one of one or more aspects in simplified form
A little concepts are with the sequence for more detailed description given later.
In the prior art, rule of thumb judge the direction that ambient field is distributed with survey line, rule mesh is pressed along X, Y-direction to data
Statistical average is sought to survey line in lattice interpolation, empirically determined direction, and the data of every bit on survey line and average value are made
Ratio repeats above step to all surveys line, then draws, judge whether mapping result is consistent with micro-judgment direction, judges
The ratio for whether needing to synthesize both direction is drawn again, which needs to establish on the basis of micro-judgment, to the personnel of explanation
Skill requirement is higher, and be easy to cause erroneous judgement.In response to this problem, the present invention proposes a kind of Transient Electromagnetic Method in Mine data analysis
Method, this method do interpolation calculation to depth-apparent resistivity data, seek even depth average value and standard deviation, are based on statistics
Theory is standardized amendment to data using average value and standard deviation, after standard deviation standardization again gridding and at
Map analysis makes full use of statistics regular, carries out processing analysis to data, disclose its internal relation, avoid artificial experience because
The deficiency of element.
The prior art needs to establish model, using the minimum value of exploration area formation resistivity as containing tunnel and with the equal of landform
The resistivity of even half space, every measurement once carry out a Three-dimensional simulation, obtain ambient field, and actual measurement potential difference is subtracted back
Scape field potential is poor, substitutes into half space apparent resistivity formula, obtains relative apparent resistivity.The prior art is mainly for direct-current exploration side
For method, discussion is the influence in tunnel and landform under grounding requirement, and problem with grounding is not present in transient electromagnetic exploration, in addition
This method is per treatment to need to do a Three-dimensional simulation, and calculation amount is huge, it is difficult to meet the requirement of fast calculation analysis.Needle
To this problem, the present invention proposes a kind of Transient Electromagnetic Method in Mine data analysing method, and this method is obtained based on preliminary treatment
Depth-apparent resistivity data is further processed, and has evaded the problem of Three-dimensional simulation calculates complicated and time consumption, is based on statistics
Theory, to data progress standard deviation standardization and at figure, realization highlights exception.
To solve the above problems, the scheme of the invention is:
A kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system, comprising:
Depth-apparent resistivity data of instrument acquisition is arranged successively by measuring point sequence；
Depth-the apparent resistivity data is pressed into regular grid interpolation；
Interpolation is by being weighted combination to depth-apparent resistivity data come what is carried out, and calculation formula is as follows:
In formula, C (x) is interpolation result, d (xi) it is original depth-resistivity data, λiFor interpolation weights, λiPass through following formula
It acquires:
Wherein parameter γijCalculation formula is as follows:
γ (h)=Var [d (x+h)-d (x)]/2
Equally distributed data point in investigative range can be obtained by interpolation, provides data point branch for subsequent anomaly extracting
Support.
Wherein, regular grid sizing grid on line direction is generally measuring point away from sizing grid on depth direction
About 2-10m, general optional 5m, the data file after exporting interpolation.It can be obtained in depth bounds and be uniformly distributed by interpolation
Data volume, be convenient for subsequent progress anomaly extracting explanation.
Even depth data C is extracted by certain depth intervalih；
By following formula, each group of even depth statistical average P is soughtih；
Wherein CihThe data of even depth h are corresponded to for measuring point i, n is measuring point number.After one group of depth-averaged value has been calculated, then count
Next group is calculated, until all depth datas calculate average value.
Seek each group of even depth data standard difference σih；
Each group of even depth data standard difference σihIt can be calculated with following formula:
In formula, n is measuring point number, CihThe data of even depth h, P are corresponded to for measuring point iihFor the flat of all measuring point even depth data
Mean value.
To the even depth data CihC ' is obtained after carrying out standard deviation standardizationih；
Reciprocity depth data carries out standard deviation standardization, and formula is as follows:
By the even depth data CihGained C ' after being standardized with standard deviationihReplacement；
By the even depth data CihWith C 'ihReplacement keeps depth column data constant, only marks apparent resistivity numerical value
Data after quasi- difference standardization are replaced, and obtain depth-standard deviation standardized data
Depth-C ' is drawn using mapping softwareihSectional view is for interpretive analysis.
Depth-C ' is drawn using mapping softwares such as SurferihSectional view carries out analytic solution to sectional view in conjunction with other data
It releases, solves corresponding geological problem.
As can be seen from the above description, the beneficial effects of the present invention are:
Present invention is generally directed to existing Mine transient electromagnetic Data processings, and apparent resistivity curve feature is changeable, directly solve
It releases and is difficult to the case where obtaining better effects, propose that a kind of standardize using standard deviation carries out Mine transient electromagnetic Data Management Analysis
Method, it is changeable to solve Mine transient electromagnetic exploration apparent resistivity pseudosection feature, explains the high problem of difficulty.Number of the invention
According to analysis processing establish on the basis of the apparent resistivity and depth data that initial data is calculated, need not rely on experience because
Element.Using standard deviation standardized analysis method, the feelings that apparent resistivity isoline in TA data is in concentric circles are effectively had modified
Condition improves the analyticity of data.
Detailed description of the invention
It is incorporated herein and the attached drawing for forming part of specification instantiates the embodiment of the present invention, and attached drawing and explanation
Book is further used for explaining the principle of the present invention together and one of ordinary skill in the art is enabled to make and use the disclosure.
Fig. 1 instantiates the flow chart of one of embodiment of the present invention Transient Electromagnetic Method in Mine data analysing method；
Fig. 2 instantiates tunnel and the signal of anomalous body relative position in the embodiment of the present invention；
Fig. 3 instantiates the apparent resistivity using conventional treatment method of the comparative example in the embodiment of the present invention into figure；
Fig. 4 is instantiated in the embodiment of the present invention and is used standard deviation standardization exception into figure；
Fig. 5 instantiate the apparent resistivity obtained using conventional treatment method in the comparative example in the embodiment of the present invention at
Figure；
Fig. 6 instantiates the standard deviation that uses in the embodiment of the present invention and is standardized into figure.
The embodiment of the present invention is described with reference to the accompanying drawings.
Specific embodiment
Embodiment
The present embodiment provides firstly a kind of Transient Electromagnetic Method in Mine data analysing method.It is detailed below in conjunction with Fig. 1
The Transient Electromagnetic Method in Mine data analysing method of the embodiment of the present invention is described.
In step slo, depth-apparent resistivity data is arranged successively by measuring point sequence；
In the present embodiment, depth-view electricity is obtained after computation of apparent resistivity and time and depth transfer on the basis of initial data
Resistance rate data, these data files may be disordered arrangements, for convenient for analysis, unification is successively arranged data by measuring point sequence
Column.
In step S20, depth-apparent resistivity data is pressed into regular grid interpolation；
To need depth-apparent resistivity data by regular grid interpolation, the grid on line direction convenient for comparing processing
Size is generally measuring point away from sizing grid is about 2-10m on depth direction, general optional 5m, the data file after exporting interpolation.
Specific interpolation the following steps are included:
Interpolation is by being weighted combination to depth-apparent resistivity data come what is carried out, and calculation formula is as follows:
In formula, C (x) is interpolation result, d (xi) it is original depth-resistivity data, λiFor interpolation weights, λiPass through following formula
It acquires:
Wherein parameter γijCalculation formula is as follows:
γ (h)=Var [d (x+h)-d (x)]/2
Equally distributed data point in investigative range can be obtained by interpolation, provides data point branch for subsequent anomaly extracting
It supports, in formula, d (x+h) is x+h depth-resistivity data, and d (x) is x depth-resistivity data, and Var is that variance calculates function.
In step s 30, even depth data C is extracted by certain depth intervalih；
In previous step interpolation file basis, even depth data, the optional 2- of depth interval are extracted by certain depth interval
10m is typically chosen 5m, this step depth interval is generally consistent with back interpolated depths direction size of mesh opening；Being consistent can
Interpolated data is directly utilized, guarantees the accuracy of data, secondly can guarantee that extracting data in depth bounds is uniformly distributed, after convenient
It is continuous to explain.
In step s 40, each group of even depth statistical average P is soughtih；
For all measuring point even depth data, its average value P is soughtih, Mean Value Formulas are as follows:
Wherein CihThe data of even depth h are corresponded to for measuring point i, n is measuring point number.After one group of depth-averaged value has been calculated, then count
Next group is calculated, until all depth datas calculate average value.
(5) each group of even depth data standard difference σ is soughtih；
Each group of even depth data standard difference σihIt can be calculated with following formula:
In formula, n is measuring point number, CihThe data of even depth h, P are corresponded to for measuring point iihFor the flat of all measuring point even depth data
Mean value.
(6) reciprocity depth data CihC ' is obtained after carrying out standard deviation standardizationih；
Reciprocity depth data carries out standard deviation standardization, and formula is as follows:
Symbol meaning is consistent with front in above formula.
(7) by even depth data CihWith C 'ihReplacement；
By even depth data CihWith C 'ihReplacement keeps depth column data constant, only by apparent resistivity numerical value standard deviation
Data after standardization are replaced, and obtain depth-standard deviation standardized data.
(8) depth-C ' is drawn using mapping softwareihSectional view is for interpretive analysis.
Depth-C ' is drawn using mapping softwares such as SurferihSectional view carries out analytic solution to sectional view in conjunction with other data
It releases, solves corresponding geological problem.
Fig. 2-4 is the standardized treatment effect of standard deviation carried out according to embodiments of the present invention.In Fig. 2-4, middle theory
It sets out, is verified using method of the method for numerical simulation to the present embodiment, to examine its effect.
In the present embodiment, it is verified using the positive algorithm of finite difference.
First from numerical simulation angle, geological model is established, is based on the positive algorithm of Three dimensional finite difference, carries out mine
Transient electromagnetic forward probe electromagnetic field responds tentative calculation.Model built as shown in Fig. 2, anomalous body be positioned at meeting head on front 50m just
Cube shape post non of low resistance body, geometric dimension are 30m × 30m × 30m.Shoulder-bed resistivity (SBR) is 100 Ω m, and abnormal body resistivity is 1
Ω·m。
According to Conventional processing steps, to analog result calculate apparent resistivity and at obtain apparent resistivity as shown in Figure 3 at
Figure.
Fig. 4 is then the standardization standard deviation method proposed using the present embodiment into figure.
Solid line boxes position is anomalous body range in model in Fig. 3-4, although from figure 3, it can be seen that conventional treatment regards
Resistivity is set to low-resistance response at figure abnormal posture, but low resistance region range is significantly greater than anomalous body range, and low-resistance region is in surround
The semicircular met head on, with anomalous body be meet head on front cube-shaped post non of low resistance body coincide it is poor.It standardizes and becomes through standard deviation in Fig. 4
After changing, anomalous body position coincide substantially in abnormal ranges and model.It can be seen that being construed to figure, standard deviation compared to conventional resistive rate
Standardization is abnormal preferably to improve Explanation Accuracy at figure.
Fig. 5-6 is the validity based on field data verifying interpretation of anomaly at figure Plays difference standardized method.
Fig. 5 and Fig. 6 is to be handled using conventional apparent resistivity at figure and standard deviation standardization exception analysis method respectively
Cheng Tu.
Fig. 5 is conventional apparent resistivity into figure, and apparent resistivity isoline is substantially in the concentric circles put centered on meeting head on, difficult
With delineation abnormality range.Fig. 6 is Abnormal Map after standard deviation standardization, is irised out altogether abnormal at 2.Position solid line boxes are helped on right side in figure
For known ponding goaf range, it can be seen that No. 2 exceptions and goaf range are coincide preferable in Fig. 6；Wherein to No. 1 in Fig. 6
Exception is drilled verifying, water outlet 5.8m3/h。
Comparison diagram 5 and Fig. 6 are explained it can be found that being coal seam within the scope of forward probe, actual resistivity should not be
Concentric circles distribution, is difficult to accurately explain front hydrological geological conditions with Fig. 5, and standard deviation standardizes interpretation of anomaly knot
Fruit and practical hydrological geological conditions coincide preferable.
As can be seen from the above description, the scheme of the present embodiment solves Mine transient electromagnetic exploration apparent resistivity pseudosection feature
It is changeable, explain the high problem of difficulty.This implementation Data Analysis Services establish the apparent resistivity being calculated in initial data and depth
Degree needs not rely on sense datum on the basis of.Using standard deviation standardized analysis method, effectively have modified in TA data
The case where apparent resistivity isoline is in concentric circles, improves the analyticity of data.
In the present embodiment, although the above method to be illustrated to and is described as a series of actions to simplify to explain, answer
It is appreciated and understood that, the order that these methods are not acted is limited, because according to one or more embodiments, some movements can be by not
Occur with order and/or with from depicted and described herein or not shown herein and describe but those skilled in the art can be with
Other movements understood concomitantly occur.
It is noted that " one embodiment ", " embodiment ", " example embodiment ", " some embodiments " etc. in specification
Reference instruction described embodiment may include a particular feature, structure, or characteristic, but each embodiment may not necessarily include
The a particular feature, structure, or characteristic.Moreover, such phrase is not necessarily referring to the same embodiment.In addition, ought retouch in conjunction with the embodiments
When stating a particular feature, structure, or characteristic, regardless of whether being expressly recited, such feature, structure are realized in conjunction with other embodiments
Or characteristic will be in the knowledge of those skilled in the art.
Offer is to make any person skilled in the art all and can make or use this public affairs to the previous description of the disclosure
It opens.The various modifications of the disclosure all will be apparent for a person skilled in the art, and as defined herein general
Suitable principle can be applied to other variants without departing from the spirit or scope of the disclosure.The disclosure is not intended to be limited as a result,
Due to example described herein and design, but should be awarded and principle disclosed herein and novel features phase one
The widest scope of cause.
Claims (10)
1. a kind of Transient Electromagnetic Method in Mine data analysis interpretation method characterized by comprising
Even depth data are extracted from depth-apparent resistivity data by certain depth interval；
Standard deviation standardization is carried out to the even depth data and data replace depth-apparent resistivity data with treated
In resistivity value to obtain depth-standard deviation standardized data；
Geologic interpretation is carried out based on the depth-standard deviation standardized data.
2. a kind of Transient Electromagnetic Method in Mine data analysis interpretation method according to claim 1, which is characterized in that packet
It includes:
Depth-apparent resistivity data is arranged successively by measuring point sequence and carries out interpolation by regular grid；
Even depth data are extracted by certain depth interval based on interpolation file.
3. a kind of Transient Electromagnetic Method in Mine data analysis interpretation method according to claim 2, which is characterized in that described
Depth interval is consistent with the longitudinal grid size in back interpolation.
4. a kind of Transient Electromagnetic Method in Mine data analysis interpretation method according to claim 1, which is characterized in that institute
Stating the progress standard deviation standardization of even depth data includes:
Measuring point even depth statistical average P is sought based on following formulaih；
Wherein, CihThe data of even depth h are corresponded to for measuring point i, n is measuring point number；
Each group of even depth data standard difference σ is sought based on following formulaih:
In formula, n is measuring point number, CihThe data of even depth h, P are corresponded to for measuring point iihFor the average value of all measuring point even depth data.
Based on following formula equity depth data CihCarry out standard deviation standardization:
C ' in formulaihFor the data of standard deviation standardization；By the even depth data CihWith C 'ihReplacement keeps depth columns
According to constant, the data after only standardizing apparent resistivity numerical value with standard deviation are replaced, and obtain depth-standard deviation normalized number
According to.
5. a kind of Transient Electromagnetic Method in Mine data analysis interpretation method according to claim 1, which is characterized in that utilize
Depth-C ' draws in the mapping softwares such as SurferihSectional view carries out analysis interpretation to sectional view in conjunction with other data, solves corresponding
Geological problem.
6. a kind of Transient Electromagnetic Method in Mine data analysis interpretation system characterized by comprising
Data extraction module extracts even depth data by certain depth interval from depth-apparent resistivity data；
Standard deviation computing module carries out standard deviation standardization to the even depth data and data is replaced deeply with treated
Resistivity value in degree-apparent resistivity data is to obtain depth-standard deviation standardized data；
Geologic interpretation module carries out geologic interpretation based on the depth-standard deviation standardized data.
7. a kind of Transient Electromagnetic Method in Mine data analysis interpretation system according to claim 6, which is characterized in that described
Data extraction module includes:
Depth-apparent resistivity data is arranged successively by measuring point sequence and carries out interpolation by regular grid by interpolating unit；
Extracting unit extracts even depth data by certain depth interval based on interpolation file.
8. a kind of Transient Electromagnetic Method in Mine data analysis interpretation system according to claim 7, which is characterized in that described
Depth interval is consistent with the longitudinal grid size in back interpolation.
9. a kind of Transient Electromagnetic Method in Mine data analysis interpretation system according to claim 6, which is characterized in that institute
Stating the progress standard deviation standardization of standard deviation computing module equity depth data includes:
Measuring point even depth statistical average P is sought based on following formulaih；
Wherein, CihThe data of even depth h are corresponded to for measuring point i, n is measuring point number；
Each group of even depth data standard difference σ is sought based on following formulaih:
In formula, n is measuring point number, CihThe data of even depth h, P are corresponded to for measuring point iihFor the average value of all measuring point even depth data.
Based on following formula equity depth data CihCarry out standard deviation standardization:
C ' in formulaihFor the data of standard deviation standardization；By the even depth data CihWith C 'ihReplacement keeps depth columns
According to constant, the data after only standardizing apparent resistivity numerical value with standard deviation are replaced, and obtain depth-standard deviation normalized number
According to.
10. a kind of Transient Electromagnetic Method in Mine data analysis interpretation system according to claim 6, which is characterized in that institute
It states geologic interpretation module and draws depth-C ' using mapping softwares such as SurferihSectional view carries out sectional view in conjunction with other data
Analysis interpretation solves corresponding geological problem.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811421881.9A CN109343134A (en) | 2018-11-27 | 2018-11-27 | A kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811421881.9A CN109343134A (en) | 2018-11-27 | 2018-11-27 | A kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109343134A true CN109343134A (en) | 2019-02-15 |
Family
ID=65318142
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811421881.9A Pending CN109343134A (en) | 2018-11-27 | 2018-11-27 | A kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109343134A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110133726A (en) * | 2019-04-25 | 2019-08-16 | 中铁二院工程集团有限责任公司 | Railway tunnel airborne electromagnetic prospecting arrangement of measuring-line method |
| CN111965721A (en) * | 2020-09-07 | 2020-11-20 | 中国煤炭地质总局地球物理勘探研究院 | Transient electromagnetic data multi-attribute fusion technology |
| CN113608267A (en) * | 2021-06-08 | 2021-11-05 | 华北科技学院(中国煤矿安全技术培训中心) | Double-circuit zero-flux transient electromagnetic detection device |
| CN114778948A (en) * | 2022-06-17 | 2022-07-22 | 中铁大桥科学研究院有限公司 | Method for monitoring resistivity of rock mass of flowing water tunnel and related equipment |
Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1528413A2 (en) * | 2003-10-27 | 2005-05-04 | Fugro N.V. | A method for the generation and application of anisotropic elastic parameters |
| CN1743872A (en) * | 2005-09-26 | 2006-03-08 | 大庆油田有限责任公司 | Multi-parameter dimension-reducing oil-gas-water-layer identifying method |
| CN102914797A (en) * | 2012-10-16 | 2013-02-06 | 中国石油天然气股份有限公司 | Method and device for obtaining formation anisotropy coefficient |
| CN103344999A (en) * | 2013-07-09 | 2013-10-09 | 中煤科工集团西安研究院 | Method for eliminating ambient fields in transient electromagnetic data through ratio method |
| CN104422969A (en) * | 2013-09-04 | 2015-03-18 | 中国石油天然气集团公司 | Method for reducing non-uniqueness of electromagnetic sounding inversion result |
| CN105116453A (en) * | 2015-08-14 | 2015-12-02 | 中国石油天然气股份有限公司 | Transient electromagnetic exploration method and device for natural gas hydrate in permafrost zone |
| CN105652325A (en) * | 2016-01-05 | 2016-06-08 | 吉林大学 | Exponential fit-adaptive Kalman-based ground-air electromagnetic data de-noising method |
| CN106324689A (en) * | 2016-06-24 | 2017-01-11 | 杭州迅美科技有限公司 | Resistivity anisotropy recognition method in horizontal well stratum environment |
| CN106556863A (en) * | 2015-09-29 | 2017-04-05 | 中国石油天然气股份有限公司 | Porosity prediction method based on Depth Domain prestack angle gathers |
| CN106814402A (en) * | 2016-12-22 | 2017-06-09 | 中石化石油工程技术服务有限公司 | Transient electromagnetic signal Prestack Noise Suppression Methods |
| US20170315199A1 (en) * | 2016-04-29 | 2017-11-02 | Schlumberger Technology Corporation | Use of echo shapes in nuclear magnetic resonance log data acquisition and quality control |
| CN107728231A (en) * | 2017-10-18 | 2018-02-23 | 科为联合(青岛)能源技术研究院有限公司 | One kind prediction nuclear magnetic resonance log T2 T2The method of distribution |
| CN107748394A (en) * | 2017-09-30 | 2018-03-02 | 中南大学 | A kind of two-parameter inversion algorithm for RMT data |
| CN107765319A (en) * | 2017-10-09 | 2018-03-06 | 中国矿业大学(北京) | A kind of Mine Transient Electromagnetic Method resistivity correction method |
| CN107807396A (en) * | 2017-09-15 | 2018-03-16 | 中国石油天然气股份有限公司 | Method and device for determining formation matrix resistivity |
| CN108122229A (en) * | 2017-12-25 | 2018-06-05 | 重庆大学 | A kind of anti-image quality evaluating method of full reference resistance |
| CN108873085A (en) * | 2018-08-02 | 2018-11-23 | 中国石油天然气集团有限公司 | A kind of sedimentary rock thin-layer electric magnetic data processing method and device |
-
2018
- 2018-11-27 CN CN201811421881.9A patent/CN109343134A/en active Pending
Patent Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1528413A2 (en) * | 2003-10-27 | 2005-05-04 | Fugro N.V. | A method for the generation and application of anisotropic elastic parameters |
| CN1743872A (en) * | 2005-09-26 | 2006-03-08 | 大庆油田有限责任公司 | Multi-parameter dimension-reducing oil-gas-water-layer identifying method |
| CN102914797A (en) * | 2012-10-16 | 2013-02-06 | 中国石油天然气股份有限公司 | Method and device for obtaining formation anisotropy coefficient |
| CN103344999A (en) * | 2013-07-09 | 2013-10-09 | 中煤科工集团西安研究院 | Method for eliminating ambient fields in transient electromagnetic data through ratio method |
| CN104422969A (en) * | 2013-09-04 | 2015-03-18 | 中国石油天然气集团公司 | Method for reducing non-uniqueness of electromagnetic sounding inversion result |
| CN105116453A (en) * | 2015-08-14 | 2015-12-02 | 中国石油天然气股份有限公司 | Transient electromagnetic exploration method and device for natural gas hydrate in permafrost zone |
| CN106556863A (en) * | 2015-09-29 | 2017-04-05 | 中国石油天然气股份有限公司 | Porosity prediction method based on Depth Domain prestack angle gathers |
| CN105652325A (en) * | 2016-01-05 | 2016-06-08 | 吉林大学 | Exponential fit-adaptive Kalman-based ground-air electromagnetic data de-noising method |
| US20170315199A1 (en) * | 2016-04-29 | 2017-11-02 | Schlumberger Technology Corporation | Use of echo shapes in nuclear magnetic resonance log data acquisition and quality control |
| CN106324689A (en) * | 2016-06-24 | 2017-01-11 | 杭州迅美科技有限公司 | Resistivity anisotropy recognition method in horizontal well stratum environment |
| CN106814402A (en) * | 2016-12-22 | 2017-06-09 | 中石化石油工程技术服务有限公司 | Transient electromagnetic signal Prestack Noise Suppression Methods |
| CN107807396A (en) * | 2017-09-15 | 2018-03-16 | 中国石油天然气股份有限公司 | Method and device for determining formation matrix resistivity |
| CN107748394A (en) * | 2017-09-30 | 2018-03-02 | 中南大学 | A kind of two-parameter inversion algorithm for RMT data |
| CN107765319A (en) * | 2017-10-09 | 2018-03-06 | 中国矿业大学(北京) | A kind of Mine Transient Electromagnetic Method resistivity correction method |
| CN107728231A (en) * | 2017-10-18 | 2018-02-23 | 科为联合(青岛)能源技术研究院有限公司 | One kind prediction nuclear magnetic resonance log T2 T2The method of distribution |
| CN108122229A (en) * | 2017-12-25 | 2018-06-05 | 重庆大学 | A kind of anti-image quality evaluating method of full reference resistance |
| CN108873085A (en) * | 2018-08-02 | 2018-11-23 | 中国石油天然气集团有限公司 | A kind of sedimentary rock thin-layer electric magnetic data processing method and device |
Non-Patent Citations (3)
| Title |
|---|
| 卓永鴻: ""標凖差舆標凖化"", 《百度文库HTTPS://WENKU.BAIDU.COM/VIEW/8C36BD797275A417866FB84AE45C3B3567ECDD05.HTML》 * |
| 李明星: ""塔里木盆地北缘侏罗系烧变岩富水性精细探测"", 《煤矿开采》 * |
| 陈研 等主编: "《资源与环境系统分析（第2版）》", 31 August 2017, 中国农业大学出版社 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110133726A (en) * | 2019-04-25 | 2019-08-16 | 中铁二院工程集团有限责任公司 | Railway tunnel airborne electromagnetic prospecting arrangement of measuring-line method |
| CN110133726B (en) * | 2019-04-25 | 2024-03-22 | 中铁二院工程集团有限责任公司 | Method for arranging exploration survey lines of railway tunnel aviation electromagnetic method |
| CN111965721A (en) * | 2020-09-07 | 2020-11-20 | 中国煤炭地质总局地球物理勘探研究院 | Transient electromagnetic data multi-attribute fusion technology |
| CN111965721B (en) * | 2020-09-07 | 2023-05-30 | 中国煤炭地质总局地球物理勘探研究院 | Transient electromagnetic data multi-attribute fusion technology |
| CN113608267A (en) * | 2021-06-08 | 2021-11-05 | 华北科技学院(中国煤矿安全技术培训中心) | Double-circuit zero-flux transient electromagnetic detection device |
| CN114778948A (en) * | 2022-06-17 | 2022-07-22 | 中铁大桥科学研究院有限公司 | Method for monitoring resistivity of rock mass of flowing water tunnel and related equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109343134A (en) | A kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system | |
| CN116665067B (en) | Ore finding target area optimization system and method based on graph neural network | |
| CN107884825B (en) | Uncertainty modeling method based on seismic multi-attribute | |
| CN109490985A (en) | A kind of geophysical reconnaissance system and method | |
| CN113051779B (en) | Numerical simulation method of three-dimensional direct-current resistivity method | |
| CN112949134A (en) | Earth-well transient electromagnetic inversion method based on non-structural finite element method | |
| CN108037490A (en) | Ground Penetrating Radar Linear Positioning Accuracy Measurement Methods and system | |
| CN106443803B (en) | Ocean controllable source electromagnetic response computational methods based on emitter actual measurement morphological data | |
| CN115238550B (en) | Land electric field numerical simulation calculation method for landslide rainfall of self-adaptive unstructured grid | |
| CN110333543A (en) | Post non of low resistance body explanation and imaging method and system based on reflection coefficient analysis | |
| CN110458129A (en) | Nonmetallic mine recognition methods based on depth convolutional neural networks | |
| Chen et al. | A total error-based multiquadric method for surface modeling of digital elevation models | |
| CN105785476A (en) | High density electrical resistivity imaging method | |
| CN105044782A (en) | Method for obtaining total organic carbon content of marine underground medium | |
| CN114114438A (en) | Quasi-three-dimensional inversion method for loop source ground-air transient electromagnetic data | |
| CN103344999B (en) | A kind of method utilizing ratio approach to eliminate the ambient field in transient electromagnetic data | |
| CN113406707A (en) | Magnetotelluric multi-scale and multi-time-period detection method | |
| Celaya et al. | Inversion of Time-Lapse Surface Gravity Data for Detection of 3-D CO 2 Plumes via Deep Learning | |
| Li et al. | Improved hybrid particle swarm optimizer with sine-cosine acceleration coefficients for transient electromagnetic inversion | |
| CN115586577A (en) | Full-time conversion method for fixed-source transient electromagnetic non-central-point observation data | |
| Parkan et al. | Individual tree segmentation in deciduous forests using geodesic voting | |
| Raith et al. | Visual Eddy Analysis of the Agulhas Current. | |
| Lebedev et al. | Comparative analysis of the North Atlantic surface circulation reproduced by three different methods | |
| CN105550442B (en) | Data processing and D integral pin-fin tube method based on the transformation of transient electrical magnetic moment | |
| CN114488327A (en) | Horizontal magnetic field and vertical magnetic field combined measurement method in well based on ground base point |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190215 |
|
| RJ01 | Rejection of invention patent application after publication |