CN109471179A - A method of new processing CSAMT data - Google Patents
A method of new processing CSAMT data Download PDFInfo
- Publication number
- CN109471179A CN109471179A CN201811316566.XA CN201811316566A CN109471179A CN 109471179 A CN109471179 A CN 109471179A CN 201811316566 A CN201811316566 A CN 201811316566A CN 109471179 A CN109471179 A CN 109471179A
- Authority
- CN
- China
- Prior art keywords
- csamt
- data
- model
- electric field
- electric
- 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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
Abstract
A kind of method of new processing CSAMT data does not change the working method technology of existing CSAMT, it does not increase input, only from the data that existing CSAMT is acquired, extract electric field component amplitude data, the algorithm that a kind of pair of electric field component amplitude data directly carries out one-dimensional inversion is developed, the processing means of existing CSAMT data are enriched.
Description
Technical field
The present invention relates to mineral resources field, specially a kind of method of new processing CSAMT data.
Background technique
CSAMT method is the abbreviation of controlled-source audiomagnetotellurics method, and this method is by D.W.Strangway and Myron
A kind of geophysical exploration method that Goldstein is proposed the 1970s.This method uses earth lead or earth-free
Loop line calculates Ka Ni Asian TV Station resistivity in the electromagnetic field component that wave region measurement is mutually orthogonal as field source.CSAMT method user
Workshop source, all has than AMT the shortcomings that overcoming the randomness of AMT method field source on signal strength, working efficiency and operating accuracy
Very big raising.This method so far, is widely used, and obtain from proposition in mineral products, the hydrology and environment exploration
Effect is reconnoitred well.But CSAMT method also has the shortcomings that it is intrinsic: it uses artificial field source, but still according to the side of AMT
Method explains, it is easy to enter transition region and near region, be unsatisfactory for wave area condition, bring constructional error.
In order to solve the problems, such as CSAMT transition region and near region, Duroux proposes MELOS method, and this method only measures one
Electric component or a magnetic component calculate apparent resistivity, breach the limitation in far field, observation area is expanded to the area Liao Zhong.China
He Jishan academician propose wide area electromagnetic method (2010), it is believed that the ground of acquisition can't be reduced by only measuring an electric or magnetic component
Lower electrical property information, is based on this, he has developed wide area electromagnetic method, and has developed specific instrument and equipment well simultaneously.Yan Shu, Xue Guo
By force, Qiu Weizhong et al. also proposed CSAMT simple component data interpretation method (2017)；Therefore, using the electric field magnitude of CSAMT this
A but component carries out processing explanation, is feasible.
A kind of new data processing method based on existing instrument and method technology proposed by the invention, for existing instrument
The CSAMT electric field component amplitude data of device equipment directly carries out the algorithm of one-dimensional inversion, can enrich existing CSAMT data
Processing means.
Summary of the invention
Due to the geologic prospect research and production unit of most domestic, CSAMT preliminry basic research, this hair can be carried out at present
The bright working method technology for not changing existing CSAMT, does not increase input, and only from the data that existing CSAMT is acquired, extracts
Electric field component amplitude data develops the method and technology that a kind of pair of electric field component amplitude data directly carries out one-dimensional inversion, enriches existing
There are the processing means of CSAMT data；Technical problem solved by the invention is to provide a kind of side of new processing CSAMT data
Method, to solve the problems mentioned in the above background technology.
In the mathematic(al) representation of the electric and magnetic fields component of CSAMT, the expression formula of electric field magnitude be include underground be situated between
The resistivity information of matter, therefore in theory, direct inversion obtains the resistance of underground medium from electric field magnitude information
Rate information is feasible.
1, the one-dimensional layered medium CSAMT electric field component mathematic(al) representation of Electric Dipole is determined:
In formula:Expression and the equidirectional horizontal electric field component in source,For the electric dipole moment of dipole source,,For
Supply current,For eelctric dipole length,For permeability of free space,For angular frequency,,,For the first layer resistivity,For the radius vector at Electric Dipole midpoint to receiving point and the folder in Electric Dipole direction
Angle,For transmitting-receiving away from,WithRespectively zero and first order Bessel function,WithRespectively contact each layer
Electrical function.
2, the catalogue scalar functions constrained based on smooth model are determined:
In formula:For catalogue scalar functions,For model level roughness,It is fixed for model roughness across the lay matrix
Justice is the first derivative of model parameter,For regularization factors,It is poor for desired fitting,For data residual error,For
Model modification amount,It is Ex to conductivity first-order partial derivative matrix,For data weighting matrix.
3, model parameter is constrained:
Layering subdivision is carried out to model, it is assumed that model there are M layers, and the CSAMT simple component electric field magnitude data amount check measured is N
It is a, for the nonnegativity for efficiently reducing multi-solution He ensuring parameter, take:
In formula:To survey Ex data,For forward modeling response data,For modelThe conductivity of the prediction of layer
Value,It isThe current conductivity value of layer,It isThe standard deviation of a observation data.
4, partial derivative matrix is sought:
The method that the solution of partial derivative matrix uses difference coefficient, it may be assumed that
5, minimum are asked to catalogue scalar functions, and inversion result are obtained to initial model iterative inversion of modifying:
It obtains:
Above-mentioned equation group is solved, is obtained, inversion result is obtained to initial model iterative inversion of modifying.
6, application process
Program or software are formed by computer programming using above-mentioned algorithm, theoretical model and real data are carried out one-dimensional
The test of Inversion Calculation has obtained preferable application effect, and the algorithm for showing that the present invention is studied is feasible.
(1) theoretical model calculates
Two theoretical stratified models are devised, their electrical parameter is as shown in table 1.
1 theoretical model electrical parameter table of table
Inverting is carried out to above-mentioned two theoretical model using inversion algorithm of the invention, attached drawing 1 and attached drawing 2 are respectively G type and H-type
The inversion result of model.In figure, color solid line is theoretical model, and solid black lines are the result that inverting obtains.
From the point of view of the inversion result of two theoretical models, the electrical parameter of model is truly reduced substantially, shows this hair
It is bright it is researched and developed based on electric field magnitude one-dimensional inversion algorithm be correct and feasible.
(2) practical application
1. using the V8 instrument of existing Canadian phoenix company, fully according to conventional CSAMT working method in Anhui Province
Fanchang County is worked.In actual work, using scalar equipment, point is received and dispatched away from being 100m away from for 9km, working frequency
For 9600-1Hz.
2. only being extracted from initial data after obtaining CSAMT data using above-mentioned V8 instrument and scalar equipment
Electric field magnitude data carry out one-dimensional inversion to electric field magnitude data, obtain resistivity number using the step 1-5 of aforementioned present invention
According to
3. the resistivity data that inverting is obtained, according to actual electric wire away from being required with scale bar, formed using computer software
Isogram is electric field magnitude resistivity profile figure, sees attached drawing 3.
4. combining actual geological condition and other data, geologic interpretation is carried out to electric field magnitude resistivity profile figure, is drawn
Divide geological structure, sees attached drawing 4.
Compared with prior art, the beneficial effects of the present invention are: the present invention does not change the working method skill of existing CSAMT
Art is not increased input, and only from the data that existing CSAMT is acquired, extracts electric field component amplitude data, develops one kind to electricity
Field component amplitude data directly carries out the algorithm of one-dimensional inversion, enriches the processing means of existing CSAMT data.
Detailed description of the invention
Fig. 1 is G model inversion result figure of the present invention.
Fig. 2 is H-type model inversion result figure of the present invention.
Fig. 3 is electric field magnitude resistivity profile figure of the present invention.
Fig. 4 is section of synthetic interpretation figure of the present invention.
Specific embodiment
It is carried out below with reference to technical effect of the embodiment and attached drawing to design of the invention, specific structure and generation clear
Chu is fully described by, to be completely understood by the purpose of the present invention, feature and effect, it is clear that described embodiment is this hair
Bright a part of the embodiment, rather than whole embodiments, based on the embodiment of the present invention, those skilled in the art are not being paid
Other embodiments obtained, belong to the scope of protection of the invention under the premise of creative work, in addition, designing in patent
All connection/the connection relationships arrived not singly refer to that building directly connects, and addition can be passed through according to specific implementation situation by referring to
Or reduce all couple auxiliaries, the coupling structure that Lai Zucheng more has, in order to make implement technical means of the invention, creation characteristic,
Reach purpose to be easy to understand with effect, following further describes the present invention in conjunction with specific drawings.
As shown, the expression formula of electric field magnitude is to include in the mathematic(al) representation of the electric and magnetic fields component of CSAMT
There is the resistivity information of underground medium, therefore in theory, direct inversion obtains underground Jie from electric field magnitude information
The resistivity information of matter is feasible.
1, the one-dimensional layered medium CSAMT electric field component mathematic(al) representation of Electric Dipole is determined:
In formula:Expression and the equidirectional horizontal electric field component in source,For the electric dipole moment of dipole source,,For
Supply current,For eelctric dipole length,For permeability of free space,For angular frequency,,,For the first layer resistivity,For the radius vector at Electric Dipole midpoint to receiving point and the folder in Electric Dipole direction
Angle,For transmitting-receiving away from,WithRespectively zero and first order Bessel function,WithRespectively contact each layer
Electrical function.
2, the catalogue scalar functions constrained based on smooth model are determined:
In formula:For catalogue scalar functions,For model level roughness,It is fixed for model roughness across the lay matrix
Justice is the first derivative of model parameter,For regularization factors,It is poor for desired fitting,For data residual error,For
Model modification amount,It is Ex to conductivity first-order partial derivative matrix,For data weighting matrix.
3, model parameter is constrained:
Layering subdivision is carried out to model, it is assumed that model there are M layers, and the CSAMT simple component electric field magnitude data amount check measured is N
It is a, for the nonnegativity for efficiently reducing multi-solution He ensuring parameter, take:
In formula:To survey Ex data,For forward modeling response data,For modelThe conductivity of the prediction of layer
Value,It isThe current conductivity value of layer,It isThe standard deviation of a observation data.
4, partial derivative matrix is sought:
The method that the solution of partial derivative matrix uses difference coefficient, it may be assumed that
5, minimum are asked to catalogue scalar functions, and inversion result are obtained to initial model iterative inversion of modifying:
It obtains:
Above-mentioned equation group is solved, is obtained, inversion result is obtained to initial model iterative inversion of modifying.
6, application process
Program or software are formed by computer programming using above-mentioned algorithm, theoretical model and real data are carried out one-dimensional
The test of Inversion Calculation has obtained preferable application effect, and the algorithm for showing that the present invention is studied is feasible.
(1) theoretical model calculates
Two theoretical stratified models are devised, their electrical parameter is as shown in table 1.
1 theoretical model electrical parameter table of table
Inverting is carried out to above-mentioned two theoretical model using inversion algorithm of the invention, attached drawing 1 and attached drawing 2 are respectively G type and H-type
The inversion result of model.In figure, color solid line is theoretical model, and solid black lines are the result that inverting obtains.
From the point of view of the inversion result of two theoretical models, the electrical parameter of model is truly reduced substantially, shows this hair
It is bright it is researched and developed based on electric field magnitude one-dimensional inversion algorithm be correct and feasible.
(2) practical application
1. using the V8 instrument of existing Canadian phoenix company, fully according to conventional CSAMT working method in Anhui Province
Fanchang County is worked.In actual work, using scalar equipment, point is received and dispatched away from being 100m away from for 9km, working frequency
For 9600-1Hz.
2. only being extracted from initial data after obtaining CSAMT data using above-mentioned V8 instrument and scalar equipment
Electric field magnitude data carry out one-dimensional inversion to electric field magnitude data, obtain resistivity number using the step 1-5 of aforementioned present invention
According to
3. the resistivity data that inverting is obtained, according to actual electric wire away from being required with scale bar, formed using computer software
Isogram is electric field magnitude resistivity profile figure, sees attached drawing 3.
4. combining actual geological condition and other data, geologic interpretation is carried out to electric field magnitude resistivity profile figure, is drawn
Divide geological structure, sees attached drawing 4.
Compared with prior art, the beneficial effects of the present invention are: the present invention does not change the working method skill of existing CSAMT
Art is not increased input, and only from the data that existing CSAMT is acquired, extracts electric field component amplitude data, develops one kind to electricity
Field component amplitude data directly carries out the algorithm of one-dimensional inversion, enriches the processing means of existing CSAMT data.
The above shows and describes the basic principles and main features of the present invention and the advantages of the present invention, the technology of the industry
Personnel are it should be appreciated that the present invention is not limited to the above embodiments, and the above embodiments and description only describe this
The principle of invention, without departing from the spirit and scope of the present invention, various changes and improvements may be made to the invention, these changes
Change and improvement all fall within the protetion scope of the claimed invention, claimed range of the invention by appended claims and
Its equivalent thereof.
Claims (5)
1. a kind of method of new processing CSAMT data, it is characterised in that: include: the number in the electric and magnetic fields component of CSAMT
It learns in expression formula, the expression formula of electric field magnitude is the resistivity information for including underground medium: determining the one-dimensional stratiform of Electric Dipole
Medium CSAMT electric field component mathematic(al) representation:
In formula:Expression and the equidirectional horizontal electric field component in source,For the electric dipole moment of dipole source,,For power supply
Electric current,For eelctric dipole length,For permeability of free space,For angular frequency,,,
For the first layer resistivity,For the radius vector at Electric Dipole midpoint to receiving point and the angle in Electric Dipole direction,For transmitting-receiving away from,
WithRespectively zero and first order Bessel function,WithRespectively contact the function of each layer electrical property.
2. a kind of method of new processing CSAMT data according to claim 1 is based on modulus of smoothness it is characterized by: determining
The catalogue scalar functions of type constraint:
In formula:For catalogue scalar functions,For model level roughness,For model roughness across the lay matrix, definition
For the first derivative of model parameter,For regularization factors,It is poor for desired fitting,For data residual error,For model
Modification amount,It is Ex to conductivity first-order partial derivative matrix,For data weighting matrix.
3. a kind of method of new processing CSAMT data according to claim 1 is it is characterized by: carry out model parameter
Constraint:
Layering subdivision is carried out to model, it is assumed that model there are M layers, and the CSAMT simple component electric field magnitude data amount check measured is N
It is a, for the nonnegativity for efficiently reducing multi-solution He ensuring parameter, take:
In formula:To survey Ex data,For forward modeling response data,For modelThe conductivity value of the prediction of layer,
It isThe current conductivity value of layer,It isThe standard deviation of a observation data.
4. a kind of method of new processing CSAMT data according to claim 1 is it is characterized by: seek partial derivative matrix:
The method that the solution of partial derivative matrix uses difference coefficient, it may be assumed that
。
5. a kind of method of new processing CSAMT data according to claim 1 is it is characterized by: seek catalogue scalar functions
It is minimum, and inversion result is obtained to initial model iterative inversion of modifying:
It obtains:
Above-mentioned equation group is solved, is obtained, inversion result is obtained to initial model iterative inversion of modifying.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811316566.XA CN109471179A (en) | 2018-11-07 | 2018-11-07 | A method of new processing CSAMT data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811316566.XA CN109471179A (en) | 2018-11-07 | 2018-11-07 | A method of new processing CSAMT data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109471179A true CN109471179A (en) | 2019-03-15 |
Family
ID=65671871
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811316566.XA Pending CN109471179A (en) | 2018-11-07 | 2018-11-07 | A method of new processing CSAMT data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109471179A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110348568A (en) * | 2019-07-16 | 2019-10-18 | 山东科技大学 | A kind of deep Mined-Out Areas method suitable for strong electromagnetic area |
| CN113960674A (en) * | 2021-10-14 | 2022-01-21 | 湖北省水文地质工程地质勘察院有限公司 | Wide-area electromagnetic method two-dimensional inversion method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4616184A (en) * | 1984-06-27 | 1986-10-07 | The United States Of America As Represented By The United States Department Of Energy | CSAMT method for determining depth and shape of a sub-surface conductive object |
| CN103630139A (en) * | 2013-12-17 | 2014-03-12 | 哈尔滨工程大学 | Underwater vehicle all-attitude determination method based on magnetic gradient tensor measurement |
| CN106021859A (en) * | 2016-05-09 | 2016-10-12 | 吉林大学 | A controlled source audio-frequency magnetotellurics one-dimensional inversion method using an improved genetic algorithm |
-
2018
- 2018-11-07 CN CN201811316566.XA patent/CN109471179A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4616184A (en) * | 1984-06-27 | 1986-10-07 | The United States Of America As Represented By The United States Department Of Energy | CSAMT method for determining depth and shape of a sub-surface conductive object |
| CN103630139A (en) * | 2013-12-17 | 2014-03-12 | 哈尔滨工程大学 | Underwater vehicle all-attitude determination method based on magnetic gradient tensor measurement |
| CN106021859A (en) * | 2016-05-09 | 2016-10-12 | 吉林大学 | A controlled source audio-frequency magnetotellurics one-dimensional inversion method using an improved genetic algorithm |
Non-Patent Citations (3)
| Title |
|---|
| 汤井田 等: "基于频点CSAMT一维最小构造反演", 《物探化探计算技术》 * |
| 袁伟: "CSAMT一维正反演与场源效应研究", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》 * |
| 鲁杏 等: "激发极化测深数据反演方法及应用", 《地质科技情报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110348568A (en) * | 2019-07-16 | 2019-10-18 | 山东科技大学 | A kind of deep Mined-Out Areas method suitable for strong electromagnetic area |
| CN110348568B (en) * | 2019-07-16 | 2021-02-23 | 山东科技大学 | Deep goaf detection method suitable for strong electromagnetic interference area |
| CN113960674A (en) * | 2021-10-14 | 2022-01-21 | 湖北省水文地质工程地质勘察院有限公司 | Wide-area electromagnetic method two-dimensional inversion method |
| CN113960674B (en) * | 2021-10-14 | 2023-11-21 | 湖北省水文地质工程地质勘察院有限公司 | Wide-area electromagnetic method two-dimensional inversion method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110058317A (en) | Aviation transient electromagnetic data and aviation magnetotelluric data joint inversion method | |
| CN104597519B (en) | A kind of deep gold mine staged method of prospecting | |
| CN106019394B (en) | Three-dimensional parallel inversion method for nonlinear conjugate gradient of ocean magnetotelluric field | |
| Haber et al. | Modeling electromagnetic fields in the presence of casing | |
| CN105204073B (en) | A kind of tensor apparent conductivity measuring method | |
| WO2010104907A2 (en) | Electromagnetic survey using naturally occurring electromagnetic fields as a source | |
| CN109100808A (en) | A kind of horizontal magnetic polarization field detection method of multi-thread source transient electromagnetic | |
| Boteler | The evolution of Québec earth models used to model geomagnetically induced currents | |
| CN106443803B (en) | Ocean controllable source electromagnetic response computational methods based on emitter actual measurement morphological data | |
| CN103809204A (en) | Field audio magnetotelluric data acquisition method | |
| CN109471179A (en) | A method of new processing CSAMT data | |
| CN104123455B (en) | Three-dimensional dip inversion method for nonlinear conjugate gradient of magnetotelluric field | |
| CN109917466A (en) | A kind of method of measuring resistivity based on electromagnetic field vertical component | |
| Lin et al. | Three-dimensional conjugate gradient inversion of magnetotelluric full information data | |
| CN109188542A (en) | A kind of the remote of wave area correlation detection refers to magnetotelluric impedance computation method | |
| Kai-Feng et al. | Inversion of time-domain airborne EM data with IP effect based on Pearson correlation constraints | |
| Zhang et al. | 3D inversion of large-scale frequency-domain airborne electromagnetic data using unstructured local mesh | |
| CN110231661A (en) | The corresponding method of rock sample and ground observation measurement parameters of electrical induction | |
| CN116050046A (en) | Magnetotelluric fuzzy constraint inversion method based on cluster analysis | |
| Xue et al. | Full waveform inversion of transient electromagnetic data in the time domain | |
| Liu et al. | A study on directional resistivity logging-while-drilling based on self-adaptive hp-FEM | |
| Persova et al. | Numerical scheme for modelling the electromagnetic field in airborne electromagnetic survey taking into account follow currents in transmitter loop | |
| CN103135140A (en) | Computing method of center loop line transverse electric and magnetic field (TEM) whole period true resistivity without fringe effect | |
| Wang et al. | 3D finite-difference modeling algorithm and anomaly features of ZTEM | |
| Yang et al. | 3D stereotomography applied to the deep-sea data acquired in the South China Sea, Part II: The real case study |
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: 20190315 |
|
| RJ01 | Rejection of invention patent application after publication |