CN104375195B - Many source multi-component three-dimensional joint inversion methods of time-frequency electromagnetism - Google Patents
Many source multi-component three-dimensional joint inversion methods of time-frequency electromagnetism Download PDFInfo
- Publication number
- CN104375195B CN104375195B CN201310356385.0A CN201310356385A CN104375195B CN 104375195 B CN104375195 B CN 104375195B CN 201310356385 A CN201310356385 A CN 201310356385A CN 104375195 B CN104375195 B CN 104375195B
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- inversion
- source
- resistivity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention is a kind of many source multi-component three-dimensional joint inversion methods of the controllable source time-frequency electromagnetism of geophysical numerical and inverting.Time-frequency electromagnetism 3-d inversion initial resistivity model is determined according to known resistivity log data and seismic prospecting data, calculate the minimum and maximum coordinate in horizontal direction x, y direction, determine 3-d inversion Inverse range in the horizontal direction, select the vertical direction size of mesh opening of many source multi-component 3-d inversions of time-frequency electromagnetism, time-frequency Electromagnetic Launching source primary field in an initial model case is calculated, the Green tensor between each hexahedron of underground is calculated;The derivative of each time-frequency Electromagnetic Sources is calculated, the minimization of object function is made using the iterative algorithm of conjugate gradient, after iteration, complete many field source multi -components time-frequency electromagnetism 3-d inversions.Measured data of the present invention has obtained the resistivity distributed in three dimensions of work area underground medium after processing, meet the explanation needs of construction, tomography, trap.
Description
Technical field
The invention belongs to the electromagnetic prospecting data processing field of geophysical exploration, is a kind of geophysical numerical
Many source multi-component three-dimensional joint inversion methods with the controllable source time-frequency electromagnetism of inverting.
Background technology
A bathymetry of building in 90 years last century is used for oil reservoir prospecting and achieves extraordinary application effect.In the method
On the basis of, the characteristics of binding time domain transient electromagnetic sounding and FEM Sounding, employ time-frequency electromagnetic method within 2005.This
The characteristics of method of kind has TDEM and frequency domain electromagnetic methods.
Time-frequency electromagnetic method measures two components of Ex and Hz, with time-domain transient electromagnetic and frequency-domain sounding feature.In recent years,
The degree that become increasingly complex with exploration targets body, explores is more and more finer, explores and is transformed into three-dimensional (3D) by two-dimentional (2D), surveys
Line is become a plurality of survey line of webbed laying by a survey line.Traditional processing method is anti-with time-domain and frequency domain one-dimensional (1D)
Drill, based on qualitatively amplitude and phase parameter are asked for, these processing methods all can not carry out three-dimensional data process, netted laying
Time-frequency electromagnetic data has only carried out the inverting of three-dimensional, and target storage, the spatial of trap and scope could be understood more
Clear, this is extremely important to improving application effect of the time-frequency electromagnetic method in oil-gas exploration, and therefore, research time-frequency electromagnetism is three-dimensional anti-
Drill and be just highly desirable to.
When netted Hydrographic General Line is carried out, not all survey line shares a transmitting to time-frequency electromagnetic method, actually
It is multiple transmittings.Because after source position is fixed, the maximum offset in survey line range transmission source is substantially stationary, skew
Away from maximum offset is exceeded, the receiver on survey line is the signal that cannot collect transmitting.Therefore netted survey line is certainly existed
Multiple emission sources, when the data of all of measuring point do 3-d inversion on to netted survey line, actually more than one, field source is more
The 3 D electromagnetic inversion problem of component.
For the multicomponent 3-d inversion problem of the multi-source of time-frequency electromagnetic method, phase was not all carried out both at home and abroad at present
The research of pass, related document do not have yet.Time-frequency electromagnetic method belongs to controllable source electromagnetic method, and in controllable source electromagnetic method, treatment technology is sent out
The most fast just number ocean controllable source electromagnetic method of exhibition, mesh mainly be one-dimensional (1D), two-dimentional (2D) and three-dimensional (3D) inversion procedure
Technology, also relates to the problem of multiple emission sources, but mainly processes electric component, and magnetic component is for calculating magnetotelluric regarding
Resistivity.And the time-frequency electromagnetic method of the Field survey layout of marine electromagnetic and land is differed.
Time-frequency electromagnetism mainly uses one-dimensional processing method at present.As exploration targets becomes increasingly complex, go to carve with one-dimensional inversion
It is to meet required precision to draw 3 D complex objective body.One field source, the 3-d inversion of several surveys line there is also same asking
Topic, inversion result can only be the part in work area, it is impossible to once be finally inversed by the distribution of whole work area subsurface resistivity.And we
Actual observation mode is exactly multiple emission sources, multiple surveys line, multiple components, so only adopting multi-source, many surveys line, multi -components
3-d inversion method, could once all invertings be out by the distribution of whole work area underground resistance.
Due to there are multiple transmitting field sources, the length for launching field source is different, so the observation data phase of different transmitting field sources
Difference is very big, and just shown as in refutation process playing a part of how the data of different field sources is very crucial.Time-frequency electromagnetism
Need to calculate Jacobi Jacobian matrix in a large number in the 3-d inversion of method, if using traditional calculus of finite differences, that is, just drill twice
The variable quantity divided by resistivity is differed, this inversion algorithm calculating speed is very slow, be difficult to realize fast imaging.Quick calculating three
The Jacobi Jacobian matrix of dimension inverting is one of key of many source multi-component 3-d inversions of time-frequency electromagnetic method.How these skills are solved
Art problem, does not find relevant clue in the paper that has delivered.
Content of the invention
It is an object of the present invention to provide a kind of oil-gas possibility with satisfaction explanation tomography, construction, base rolling shape, trap is commented
Many source multi-component three-dimensional joint inversion methods of the time-frequency electromagnetism that valency needs.
The present invention is realized by following steps:According to distribution and the component type of actual measurement time-frequency Electromagnetic Launching frequency,
Select the data for participating in inverting;
The described scope for participating in the multicomponent tranmitting frequency of time-frequency electromagnetism multi-source in inverting is 0.01-100Hz.
The described multicomponent component of time-frequency electromagnetism multi-source in inverting of participating in is the electric field component Ex parallel with field source and hangs down
Straight magnetic-field component Hz.
2) the initial 1D electricity of time-frequency electromagnetism 3-d inversion is determined according to known resistivity log data and seismic prospecting data
Resistance rate model, thickness degree and resistivity value;
The deep lateral resistivity data of described Electric Log Data employing, the reflection line-ups of foundation seismic profile,
The height relation of electric logging data, determines that thickness degree and the number of plies of 1D background models, 1D background resistivity values size pass through electrical measurement
Well data determine, while ensureing total longitudinal conductance curve weight of total longitudinal conductance curve and 1D background resistivity models of log data
Close.
3) according to the coordinate in horizontal direction x, the y direction of all measuring points in ground, calculate the maximum in horizontal direction x, y direction
And min coordinates, determine 3-d inversion Inverse range in the horizontal direction, obtain along with 2000m according to the depth of objective body
The depth of 3-d inversion;With x, the Inverse range in tri- directions of y, z obtains 3-d inversion x, y divided by respective subdivision size of mesh opening,
The subdivision grid number in tri- directions of z;
Horizontal direction x of the 3-d inversion, the Inverse range of y include all of measuring point, and size of mesh opening is 100m.
Depth 2000m of the maximum inverting depth of vertical direction z of the 3-d inversion more than objective body, vertical direction z
By gradually increasing, the size of first grid is 50m to the size of mesh opening in direction, and the size of last grid is 500m.
4) size of the vertical direction z grid of many source multi-component 3-d inversions of correct time-frequency electromagnetism is selected;
Described selection is the principle according to only one of which resistivity value in a grid, the vertical z directions of 3-d inversion
Any one grid can not simultaneously be located at two stratum of 1D background layers in;
5) primary field of the time-frequency Electromagnetic Launching source in 1D resistivity initial models is calculated;
6) Green tensor of each hexahedron of underground to ground receiving point is calculated, calculates the lattice between each hexahedron of underground
Woods tensor;
7) the Frechet derivatives of each time-frequency Electromagnetic Sources, described calculating time-frequency electromagnetism Frechet Jacobian matrixs are calculated
Calculated by analytic formula:
In formula, EaIt is horizontal component of electric field Ex, HaIt is vertical magnetic field, σ is the derivative of resistivity,WithIt is electric field and magnetic field
Green tensor, EbBackground electric field.rjIt is the distance to measuring point, r is the distance to subdivision grid,It is relevant with Green tensor
Coefficient matrix, g is and ambient field EbThe coefficient matrix relevant with Green tensor.
8) according to regularization inversion principle calculating target function, described is calculated as follows:
In formulaIt is the error of fitting of the data in i-th time-frequency Electromagnetic Launching source;It is
The weight function in i time-frequency Electromagnetic Launching source, S (m) are to stabilize function;α is regularization parameter;AiSon is just being calculated;diWhen being i
During frequency Electromagnetic Launching source, the electric field E of time-frequency electromagnetic observation pointxWith magnetic field HzComponent;
9) the minimization of object function is made using the iterative algorithm of conjugate gradient;
When described the minimization of object function is calculated, in first time iterative process is minimized, regularization parameter value is
0.5, α span is 0.05-5.Automatically adjust the size of α in an iterative process, such as fit error increase, increase α;If
Error of fitting reduces, and α reduces.
Time-frequency electromagnetic data weighting matrix is calculated using equation below in described iterative algorithm
In formula:It is the error of fitting of the initial model in i-th time-frequency Electromagnetic Launching source;It is n-th time-frequency electromagnetism
The error of fitting of the initial model of emission source;Max is to calculate max function;Sum is summing function.
10) repeat step 9), when the error of fitting of the object function of the time-frequency electromagnetic inversion of multiple field sources reaches setting
After error criterion or iterations have exceeded the maximum iteration time for setting, many field source multi -components time-frequency electromagnetism are completed three-dimensional anti-
Drill.
The present invention synthesizes to model and measured data has done experiment process, the data inversion result of model synthesis and theoretical mould
Type is essentially the same, illustrates of the invention effective;The resistivity distributed in three dimensions of work area underground medium has been obtained after measured data process, full
Foot construction, tomography, the explanation of trap need.
Description of the drawings
Fig. 1 time-frequency electromagnetism electric fields and magnetic field amplitude curve, (left) Ex, (right) Hz.
Fig. 2 logs and total longitudinal conductance curve, (left) log, (right) longitudinal conductance curve.
The one-dimensional background resistivity models of Fig. 3.
Certain work area survey line of Fig. 4 and point position figure, before (left) rotating and translating, after (right) rotating and translating.
The 3D of many source multi-component 3D joint inversion results of the time-frequency electromagnetism of Fig. 5 work areas measured data shows.
Horizontal direction resistivity flat distribution map of Fig. 6 depth in 6000m.
The resistivity section figure of Fig. 7 surveys line 1.
Specific embodiment
The present invention is described in detail below in conjunction with accompanying drawing.
1) according to distribution and the component type of actual measurement time-frequency Electromagnetic Launching frequency, the data for participating in inverting are selected;When
The scope of the multicomponent tranmitting frequency of frequency electromagnetism multi-source is 0.01-100Hz.Fig. 1 is time-frequency electromagnetism electric field and magnetic field amplitude curve,
(left) Ex, (right) Hz.
2) the initial 1D electricity of time-frequency electromagnetism 3-d inversion is determined according to known resistivity log data and seismic prospecting data
Resistance rate model, thickness degree and resistivity value.Fig. 2 logs and total longitudinal conductance curve, (left) log (right) are indulged
Conductance plots, Fig. 2 (left) are according to Resistivity log Changing Pattern, set up the resistivity value and thickness degree of 1D backgrounds.
Fig. 2 (right) is the total longitudinal conductance curve comparison of the longitudinal conductance curve with the 1D background models that sets up of electric logging data, it is ensured that make
Therebetween it is fitted.Fig. 3 is one-dimensional background resistivity model.
3) according to the coordinate in horizontal direction x, the y direction of all measuring points in ground, calculate the maximum in horizontal direction x, y direction
And min coordinates, determine 3-d inversion Inverse range in the horizontal direction, obtain along with 2000m according to the depth of objective body
The depth of 3-d inversion;With x, the Inverse range in tri- directions of y, z obtains 3-d inversion x, y divided by respective subdivision size of mesh opening,
The subdivision grid number in tri- directions of z.
Fig. 4 is survey line and point position figure, before (left) rotating and translating, after (right) rotating and translating.Through rotation
After translation, the Inverse range in final x directions is -5000m 12600m, and the Inverse range in y directions is 0m 12400m, z
The Inverse range in direction is 0m 10000m.The size of the grid cell in x, y and z direction is all 200m, so 3-d inversion
Grid numerical digit 89 × 63 × 51.
4) size (parameter) of the vertical direction z grid of many source multi-component 3-d inversions of correct time-frequency electromagnetism is selected.
5) primary field of the time-frequency Electromagnetic Launching source in 1D resistivity initial models is calculated.
6) Green tensor of each hexahedron of underground to ground receiving point is calculated, calculates the lattice between each hexahedron of underground
Woods tensor.
7) the Frechet derivatives of each time-frequency Electromagnetic Sources are calculated.The computing formula of Frechet Jacobian matrixs is shown in formula
(1)、(2).
8) according to regularization inversion principle calculating target function, described computing formula is shown in formula (3).
9) the minimization of object function is made using the iterative algorithm of conjugate gradient, canonical in first time iterative process is minimized
The span for changing parameter alpha value for 0.5, α is 0.05-5.Automatically adjust the size of α in an iterative process, such as fit error increasing
Plus, increase α;If error of fitting reduces, α reduces.
10) repeat step 9), when the error of fitting of the object function of the time-frequency electromagnetic inversion of multiple field sources reaches setting
After error criterion or iterations have exceeded the maximum iteration time for setting, many field source multi -components time-frequency electromagnetism are completed three-dimensional anti-
Drill.Fig. 5 is that the 3D of many source multi-component 3D joint inversion results of time-frequency electromagnetism of certain work area measured data shows that Fig. 6 is through 20
After secondary iterative inversion, the 3 D resistivity distribution of work area underground medium has been obtained.Fig. 6 is horizontal direction electricity of the depth in 6000m
Resistance rate flat distribution map, Fig. 7 are the resistivity section figures of survey line 1.Basic electricity be can be seen that from 3-d inversion result above
Property layer be finally inversed by, with log data coincide.Base rolling shape is clearly finally inversed by by 3-d inversion.
Buried hill buried depth of top surface delimited for accurate, the migration imaging and Explanation Accuracy of earthquake improved.
Claims (10)
1. many source multi-components three-dimensional joint inversion methods of a kind of time-frequency electromagnetism, feature is through the following steps that realize:
1) according to distribution and the component type of actual measurement time-frequency Electromagnetic Launching frequency, the data for participating in inverting are selected;
2) the initial 1D resistivity of time-frequency electromagnetism 3-d inversion is determined according to known resistivity log data and seismic prospecting data
Model, thickness degree and resistivity value;
3) according to the coordinate in horizontal direction x, the y direction of all measuring points in ground, calculate the maximum of horizontal direction x, y direction and most
Little coordinate, determines 3-d inversion Inverse range in the horizontal direction, obtains three-dimensional according to the depth of objective body along with 2000m
The depth of inverting;With x, the Inverse range in tri- directions of y, z obtains 3-d inversion x, y, z tri- divided by respective subdivision size of mesh opening
The subdivision grid number in individual direction;
4) size of the vertical direction z grid of many source multi-component 3-d inversions of correct time-frequency electromagnetism is selected;
5) primary field of the time-frequency Electromagnetic Launching source in 1D resistivity initial models is calculated;
6) Green tensor of each hexahedron of underground to ground receiving point is calculated, calculates the Green between each hexahedron of underground
Amount;
7) the Frechet derivatives of each time-frequency Electromagnetic Sources are calculated;
8) according to regularization inversion principle calculating target function, described is calculated as follows:
In formulaIt is the error of fitting of the data in i-th time-frequency Electromagnetic Launching source;It is i-th
The weight function in time-frequency Electromagnetic Launching source, S (m) are to stabilize function;α is regularization parameter;AiSon is just being calculated;diIt is i time-frequency electricity
During magnetic emission source, the electric field E of time-frequency electromagnetic observation pointxWith magnetic field HzComponent;
9) the minimization of object function is made using the iterative algorithm of conjugate gradient;
10) repeat step 9), when the error of fitting of the object function of the time-frequency electromagnetic inversion of multiple field sources reaches the error of setting
After standard or iterations have exceeded the maximum iteration time for setting, many field source multi -components time-frequency electromagnetism 3-d inversions are completed.
2. method according to claim 1, feature is step 1) described in participation inverting in multicomponent of time-frequency electromagnetism multi-source
The scope of radio frequency rate is 0.01-100Hz.
3. method according to claim 1, feature is step 1) described in participation inverting in multicomponent point of time-frequency electromagnetism multi-source
Amount is the electric field component Ex parallel with field source and perpendicular magnetic component Hz.
4. method according to claim 1, feature is step 2) described in Electric Log Data using deep lateral resistivity
Data, according to the reflection line-ups of seismic profile, the height relation of electric logging data, determine the thickness degree and layer of 1D background models
Number, 1D background resistivity values size are determined by electric logging data, while ensureing the total longitudinal conductance curve and 1D back of the body of log data
Total longitudinal conductance curve co-insides of scape resistivity models.
5. method according to claim 1, feature is step 3) horizontal direction x of the 3-d inversion, the Inverse range of y includes
All of measuring point, size of mesh opening are 100m.
6. method according to claim 1, feature is step 3) the maximum inverting depth of vertical direction z of the 3-d inversion is big
In depth 2000m of objective body, the size of mesh opening in vertical direction z direction gradually increases, and the size of first grid is 50m, most
The size of a grid is 500m afterwards.
7. method according to claim 1, feature is step 4) described in selection be according to only one of which resistance in a grid
The principle of rate value, any one grid in the vertical z directions of 3-d inversion can not be located in two stratum of 1D background layers simultaneously.
8. method according to claim 1, feature is step 7) described in calculating time-frequency electromagnetism Frechet Jacobian matrixs pass through with
Lower analytic formula is calculated:
In formula, EaIt is horizontal component of electric field Ex, HaIt is vertical magnetic field, σ is the derivative of resistivity,WithIt is the Green in electric field and magnetic field
Tensor, EbBackground electric field, rjIt is the distance to measuring point, r is the distance to subdivision grid,It is the coefficient relevant with Green tensor
Matrix, g is and ambient field EbThe coefficient matrix relevant with Green tensor.
9. method according to claim 1, feature is step 9) described in the minimization of object function when calculating, minimizing first
In secondary iterative process, regularization parameter value is 0.05-5 for the span of 0.5, α;α big is automatically adjusted in an iterative process
Little, if fitting error increase, increase α;If error of fitting reduces, α reduces.
10. method according to claim 1, feature is step 9) described in iterative algorithm in time-frequency electricity is calculated using equation below
Magnetic data weighting matrix
In formula:It is the error of fitting of the initial model in i-th time-frequency Electromagnetic Launching source;It is that n-th time-frequency electromagnetism is sent out
Penetrate the error of fitting of the initial model in source;Max is to calculate max function;Sum is summing function.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310356385.0A CN104375195B (en) | 2013-08-15 | 2013-08-15 | Many source multi-component three-dimensional joint inversion methods of time-frequency electromagnetism |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310356385.0A CN104375195B (en) | 2013-08-15 | 2013-08-15 | Many source multi-component three-dimensional joint inversion methods of time-frequency electromagnetism |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104375195A CN104375195A (en) | 2015-02-25 |
CN104375195B true CN104375195B (en) | 2017-03-15 |
Family
ID=52554223
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310356385.0A Active CN104375195B (en) | 2013-08-15 | 2013-08-15 | Many source multi-component three-dimensional joint inversion methods of time-frequency electromagnetism |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104375195B (en) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107305600A (en) * | 2016-04-21 | 2017-10-31 | 新疆维吾尔自治区煤炭科学研究所 | Least square method resistivity three-dimensional approximate inversion technology |
CN106019394B (en) * | 2016-04-27 | 2019-04-05 | 中国地质科学院矿产资源研究所 | Three-dimensional parallel inversion method for nonlinear conjugate gradient of ocean magnetotelluric field |
CN106556876B (en) * | 2016-11-11 | 2018-05-15 | 山东大学 | A kind of three-dimensional NMR prestack inversion method based on the excitation of multifrequency off resonance |
CN108119764B (en) * | 2017-12-26 | 2020-02-07 | 东莞理工学院 | Time reversal self-adaptive grid gas transmission pipeline leakage detection method |
CN108304618B (en) * | 2018-01-05 | 2019-05-24 | 台州创兴环保科技有限公司 | A kind of gravimetric data and magnetotelluric data joint inversion method |
CN108345049B (en) * | 2018-02-12 | 2020-05-26 | 山东大学 | Underground engineering unfavorable geological detection multi-method constraint inversion and joint interpretation method |
CN108802834B (en) * | 2018-02-13 | 2020-12-22 | 中国科学院电子学研究所 | Underground target identification method based on joint inversion |
CN109541703A (en) * | 2018-11-30 | 2019-03-29 | 西南科技大学 | A kind of quantitative preferred method of magnetotelluric sounding curve |
CN109557601A (en) * | 2019-01-22 | 2019-04-02 | 青岛海洋地质研究所 | Reservoir parameter inversion method is combined in one-dimensional ocean controllable source electromagnetism and earthquake |
CN112415628B (en) * | 2019-08-20 | 2022-07-26 | 核工业二一六大队 | Sandstone-type uranium ore sampling plane distribution diagram forming method |
CN113447990B (en) * | 2020-03-24 | 2024-04-02 | 中国石油化工股份有限公司 | Method and device for observing abnormal electrical property of well site |
CN111983703B (en) * | 2020-07-24 | 2023-07-25 | 中国石油天然气集团有限公司 | Method, system and device for imaging electromagnetic measurement fluid between wells |
CN111856597B (en) * | 2020-08-05 | 2023-03-21 | 中国海洋大学 | Towed marine electromagnetic formation resistivity and receiving station position joint inversion method |
CN111880235B (en) * | 2020-08-05 | 2023-03-28 | 中国海洋大学 | Ocean electromagnetic formation anisotropic resistivity and emission source posture joint inversion method |
CN111856596B (en) * | 2020-08-05 | 2023-03-28 | 中国海洋大学 | Layered medium resistivity anisotropy ocean controllable source electromagnetic rapid inversion method |
CN111983689B (en) * | 2020-08-27 | 2023-06-02 | 青海省第四地质勘查院 | Near-source electromagnetic seismoelectric joint GR inversion method |
CN112379449B (en) * | 2020-10-30 | 2023-05-26 | 中国石油天然气集团有限公司 | Processing method and device for electromagnetic data of controllable source |
CN112596108B (en) * | 2020-11-24 | 2022-08-23 | 中国地质科学院地球物理地球化学勘查研究所 | AMT (automated mechanical Transmission) profile detection method, device and equipment |
CN114764151B (en) * | 2021-01-13 | 2023-06-23 | 中国石油化工股份有限公司 | Magnetotelluric frequency division chromatography inversion method |
CN114859421B (en) * | 2021-02-03 | 2024-05-31 | 中国科学院声学研究所 | Underwater buried target identification method based on multi-parameter simultaneous inversion |
CN113391362B (en) * | 2021-08-13 | 2021-10-29 | 成都理工大学 | Magnetotelluric profile three-dimensional structured inversion method based on corridor data constraint |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102798898A (en) * | 2012-08-20 | 2012-11-28 | 中国地质科学院矿产资源研究所 | Three-dimensional inversion method for nonlinear conjugate gradient of magnetotelluric field |
CN102798897A (en) * | 2012-08-20 | 2012-11-28 | 中国地质科学院矿产资源研究所 | Pit-well ground magnetotelluric field nonlinear conjugate gradient two-dimensional inversion method |
CN103576205A (en) * | 2013-11-17 | 2014-02-12 | 成都理工大学 | Land and air transient electromagnetism exploring method based on combined magnetism source technology |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2415511B (en) * | 2004-06-26 | 2008-09-24 | Statoil Asa | Processing electromagnetic data |
US20090150124A1 (en) * | 2007-12-07 | 2009-06-11 | Schlumberger Technology Corporation | Model based workflow for interpreting deep-reading electromagnetic data |
-
2013
- 2013-08-15 CN CN201310356385.0A patent/CN104375195B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102798898A (en) * | 2012-08-20 | 2012-11-28 | 中国地质科学院矿产资源研究所 | Three-dimensional inversion method for nonlinear conjugate gradient of magnetotelluric field |
CN102798897A (en) * | 2012-08-20 | 2012-11-28 | 中国地质科学院矿产资源研究所 | Pit-well ground magnetotelluric field nonlinear conjugate gradient two-dimensional inversion method |
CN103576205A (en) * | 2013-11-17 | 2014-02-12 | 成都理工大学 | Land and air transient electromagnetism exploring method based on combined magnetism source technology |
Also Published As
Publication number | Publication date |
---|---|
CN104375195A (en) | 2015-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104375195B (en) | Many source multi-component three-dimensional joint inversion methods of time-frequency electromagnetism | |
Grayver et al. | 3D inversion and resolution analysis of land-based CSEM data from the Ketzin CO 2 storage formation | |
CN101438176B (en) | Time lapse analysis with electromagnetic data | |
CN100429531C (en) | 3D electromagnetic fast inversion method of minimized target | |
Berge et al. | Electrical resistivity tomography investigations of multilayered archaeological settlements: part I–modelling | |
CN104280782B (en) | One-dimensional joint inversion method for time-frequency electromagnetic data and magnetotelluric data | |
CN108680964A (en) | A kind of normalization weight magnetoelectricity shake joint inversion method based on structural constraint | |
Yin et al. | 3D time-domain airborne EM forward modeling with topography | |
CN106170605A (en) | Proppant in induced breakage is positioned and the system and method for imaging | |
CN102707323A (en) | Controllable source audio-frequency magnetic field sounding method for geological exploration | |
Nguyen et al. | Comparing large-scale 3D Gauss–Newton and BFGS CSEM inversions | |
Park et al. | CSEM sensitivity study for Sleipner CO2-injection monitoring | |
Triantafilis et al. | Modeling the electrical conductivity of hydrogeological strata using joint-inversion of loop-loop electromagnetic data | |
CN107305600A (en) | Least square method resistivity three-dimensional approximate inversion technology | |
Ekinci et al. | Geophysical imaging survey in the south necropolis at the ancient city of Parion (Kemer-Biga), Northwestern Anatolia, Turkey: preliminary results | |
Papadopoulos et al. | Electrical resistivity tomography for the modelling of cultural deposits and geomophological landscapes at Neolithic sites: a case study from Southeastern Hungary | |
CN109164033A (en) | A method of the acquisition unsaturated soil waterpower conductance based on resistivity method | |
CN104422969A (en) | Method for reducing non-uniqueness of electromagnetic sounding inversion result | |
Pan et al. | Estimating field-scale soil water dynamics at a heterogeneous site using multi-channel GPR | |
Al‐Saadi et al. | Very‐high‐resolution electrical resistivity imaging of buried foundations of a Roman villa near Nonnweiler, Germany | |
Yu et al. | A combining regularization strategy for the inversion of airborne time-domain electromagnetic data | |
Yilmaz et al. | Preliminary results of an integrated archaeo‐geophysical survey on the basis of ancient finds unearthed by an illegal excavation at Kılıç Ören site (Isparta, Turkey) | |
Feng et al. | Contrast between 2D inversion and 3D inversion based on 2D high-density resistivity data | |
Xiong | 2.5 D forward for the transient electromagnetic response of a block linear resistivity distribution | |
Cosentino et al. | Study and monitoring of salt water intrusion in the coastal area between Mazara del Vallo and Marsala (South-Western Sicily) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |