CN109407160A - A method of D integral pin-fin tube is carried out based on WEM vector finite element method - Google Patents
A method of D integral pin-fin tube is carried out based on WEM vector finite element method Download PDFInfo
- Publication number
- CN109407160A CN109407160A CN201811554392.0A CN201811554392A CN109407160A CN 109407160 A CN109407160 A CN 109407160A CN 201811554392 A CN201811554392 A CN 201811554392A CN 109407160 A CN109407160 A CN 109407160A
- Authority
- CN
- China
- Prior art keywords
- wem
- field
- fin tube
- integral pin
- finite element
- 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 invention discloses a kind of methods for carrying out D integral pin-fin tube based on WEM vector finite element method, belong to deep resource detection technology field, D integral pin-fin tube method are as follows: on the basis of considering that ionosphere effect and displacement current act on, it is decomposed into the path that primary field and secondary field solve respectively using by electromagnetism resultant field, overcomes the singularity of field source；In fast Hankel transform, using 2040 high-density sampling filter factors, retains enough high frequency oscillation information, improve arithmetic accuracy.The present invention can be improved Neoproterozoic Progress of Residual Basins forecast of distribution Three-dimensional simulation in ultra-deep and provide scientific basis for petroleum resources assessment in the depth bounds.
Description
Technical field
The invention belongs to deep resource detection technology fields, for solving three-dimensional geology of deep part stratum, the fine detection of construction
Numerical simulation and data interpretation, specifically, be related to it is a kind of based on WEM vector finite element method carry out D integral pin-fin tube method.
Background technique
Artificial source's Extremely Low Frequency Electromagnetic method (WEM) is up to several hundred amperes of electricity to subsurface injection using the emission source for being about 100km
Stream, and then penetrate field source in the inscribed transmitting-receiving of thousands of miles range and generate alternating electromagnetism field signal, area can be carried out in China
The three-dimensional Sounding Research On Deep and resource exploration of property.Ball medium is complicated three-dimensional structure, and collected data need to carry out
3-d inversion and explanation, and consideration emission source D integral pin-fin tube existing at present is modeled as considering ionosphere and cannot achieve high band
It calculates, can not also overcome the singularity of field source, arithmetic accuracy is very limited.
In view of this present invention is specifically proposed.
Summary of the invention
The technical problem to be solved in the present invention is that overcoming the deficiencies of the prior art and provide a kind of limited based on WEM vector
The method that first method carries out D integral pin-fin tube can improve Neoproterozoic Progress of Residual Basins forecast of distribution essence in ultra-deep by Three-dimensional simulation
It spends and provides scientific basis for petroleum resources assessment in the depth bounds.
In order to solve the above technical problems, the present invention is using the basic conception of technical solution:
The present invention realizes the D integral pin-fin tube of WEM method using Vector Finite-Element Method, firstly, the surprise in order to overcome field source
The opposite sex, present invention employs electromagnetism resultant field is decomposed into the path that primary field and secondary field solve respectively.In the solution of primary field
In, it is calculated using modified open source software Dipole1D, considers ionosphere effect and displacement current effect.In quick Hankel
In transformation, using 2040 high-density sampling filter factors, retains enough high frequency oscillation information, improve arithmetic accuracy.It solves
Complete primary field EPAfterwards, then secondary field E is calculatedS.In research process, taking time-harmonic factor is e-iωt, ignore the variation of dielectric constant e, magnetic
Conductance is taken as μ0.In view of field source is finite length grounded source, if dispatch from foreign news agency flux is JS.It can according to Maxwell equation group
To obtain double vorticity equations of electric field:
According to principle of stacking, electric field resultant field is decomposed into the form of primary field and secondary field sum
E=EP+ES (2)
Double vorticity equations of same available electric field primary field:
σ in formulaPFor the conductivity for calculating a field model (generally homogeneous half space or layer-cake model), σ
For the conductivity of true model.
Formula (1), which is subtracted formula (3), to be obtained
Wherein, σa=σ-σPIt is the difference of model conductivity and background conductance rate to be calculated, referred to as residual conduction rate.Always
Singularity Fields source item in field equation is by nonsingular function σaEPIt is replaced.Since secondary field is often weaker, and boundary condition is apart from different
When often body is remote enough, the secondary field that anomalous body generates is on boundary it is considered that being equal to zero.Therefore, by aerial top boundary, underground
The secondary electrical field E of bottom boundaries and four lateral boundariesSZero is taken as the boundary condition in Three-dimensional simulation.
Because seamed edge base divergence is 0, do not need to carry out Divergence correction in solution procedure, thus it is limited using trivector herein
Member solves above-mentioned boundary value problem.With weight function dot product boundary value problem is various respectively integrate again after be added obtain Weighted Residual, enable
Weighted Residual is equal to zero and obtains weak form.According to the golden method of gal the Liao Dynasty, weighting function is basic function.Since basic function is only in unit
It is inside not zero, element integral can be written as
Above equation shares 8 × Nelement(unit number) is a, and unknown number only has Nedge(seamed edge number) is a redundancy side
Journey group.It is combined in global seamed edge by local seamed edge, the total number of equation group is reduced to NedgeIt is a, become with determining solution
Equation group.Last overall situation matrix equation can be written as
KIJ、PIJ、MIJ、WIJFor global coefficient matrix.The secondary field value ES of boundary is forced in formula (6).Finally, formula
(6) only ESJ is unknown quantity in, and secondary field value ESJ can be obtained by solving above-mentioned equation group, is added that electric field can be obtained is total with primary field
?.Then according to Maxwell equation group, the magnetic field H of point position at earth's surface is acquired, the impedance at measuring point can be calculated, into
And obtain apparent resistivity and phase information.
The present invention is decomposed into one using by electromagnetism resultant field on the basis of considering that ionosphere effect and displacement current act on
The path that secondary field and secondary field solve respectively overcomes the singularity of field source；It is high using 2040 points in fast Hankel transform
Density sampling filter coefficient retains enough high frequency oscillation information, improves arithmetic accuracy.Have to realize WEM trivector
First method D integral pin-fin tube is limited, lays a solid foundation for subsequent WEM electromagnetic method 3-d inversion, Neoproterozoic Progress of Residual Basins in ultra-deep can be improved
Forecast of distribution Three-dimensional simulation and in the depth bounds petroleum resources assessment scientific basis is provided.
After adopting the above technical scheme, the present invention has the advantages that compared with prior art.
Very low frequencies Vector finite element method D integral pin-fin tube method proposed by the present invention, ionosphere is considered when simulating to emission source
It influences and displacement current acts on, realize the D integral pin-fin tube simulation of range in WEM emission system transmitting covering thousands of miles；Fast
2040 high-density sampling filter factors are used in fast Hankel transformation, improve high band electromagnetic field arithmetic accuracy, are met practical
The needs of simulation.
A specific embodiment of the invention is described in further detail with reference to the accompanying drawing.
Detailed description of the invention
The a part of attached drawing as the application, for providing further understanding of the invention, of the invention is schematic
Examples and descriptions thereof are used to explain the present invention, but does not constitute an undue limitation on the present invention.Obviously, the accompanying drawings in the following description
Only some embodiments to those skilled in the art without creative efforts, can be with
Other accompanying drawings can also be obtained according to these attached drawings.In the accompanying drawings:
Fig. 1 is three-dimensional WEM low-resistance exception body Model；
Fig. 2 is the apparent resistivity curve at coordinate origin；
Phase curve at Fig. 3 coordinate origin；
The apparent resistivity forward modeling pseudosection map of Fig. 4 low-resistance anomalous body；
The apparent resistivity forward modeling slice map of Fig. 5 low-resistance anomalous body.
It should be noted that these attached drawings and verbal description are not intended to the design model limiting the invention in any way
It encloses, but illustrates idea of the invention by referring to specific embodiments for those skilled in the art.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, the technical solution in embodiment is clearly and completely described, the following examples are intended to illustrate the invention, but
It is not intended to limit the scope of the invention.
Low-resistance exception body Model as shown in Figure 1 is established, the geometric dimension of low-resistance prism is 750m × 750m × 250m,
Upper top surface is 250m far from ground, and defining is coordinate origin at anomalous body center.There are two at anomalous body center left 600km
Orthogonal injection antenna, length are 100km, emission current 100A.The height in ionosphere is 100m, and the resistivity in ionosphere is close
It is seemingly 10-4 Ω m.The resistivity of air layer is set as 10-14 Ω m, and ball medium resistivity is set as 100 Ω m, and low-resistance is different
The resistivity of normal body is set as 10 Ω m.Frequency range used in forward modeling is 1~300Hz, totally 20 frequency points.Not according to two
Equidirectional polarized electric and magnetic fields carry out being converted to impedance tensor Zxx、Zxy、Zyx、Zyy, further find out apparent resistivity and phase
Position
Fig. 2 is the apparent resistivity map of (low-resistance anomalous body center) at coordinate origin under XY mode and YX mode, Fig. 3 respectively
It is the phase diagram under XY mode and YX mode at coordinate origin respectively.Fig. 4 illustrate by anomalous body geometric center along x-axis and
Apparent resistivity forward modeling pseudosection map on upper two sections along the y-axis direction.Fig. 5 be when frequency be 1.35,4.48,14.91,
When 49.53,164.57Hz, the spatial distribution of low-resistance anomalous body in the horizontal direction.In forward modeling pseudosection map and forward modeling slice map
In on, the presence of low-resistance anomalous body can be clearly visible, illustrate that WEM method is sensitiveer to low-resistance anomalous body.But post non of low resistance body
Size can not judge with depth, need further inverting that can just learn.In high band, since electromagnetic wave skin depth is small, depending on
Resistivity is close to the resistivity of country rock, and unobvious to abnormal precursor reactant, in middle low frequency, the influence of anomalous body is shown.
The above is only presently preferred embodiments of the present invention, is not intended to limit the present invention in any form, though
So the present invention has been disclosed as a preferred embodiment, and however, it is not intended to limit the invention, any technology people for being familiar with this patent
Member without departing from the scope of the present invention, when the technology contents using above-mentioned prompt make it is a little change or be modified to
The equivalent embodiment of equivalent variations, but anything that does not depart from the technical scheme of the invention content, it is right according to the technical essence of the invention
Any simple modification, equivalent change and modification made by above embodiments, in the range of still falling within the present invention program.
Claims (3)
1. a kind of method for carrying out D integral pin-fin tube based on WEM vector finite element method, which is characterized in that sent out using the signal of WEM method
The electromagnetic wave for penetrating source generation upwardly propagates to form sky wave, and sky wave propagation is to ionosphere, then by ionospheric reflection to underground medium, most
Ground is traveled to by underground medium afterwards；By this method to the Three-dimensional simulation of trial zone, help to understand the three-dimensional for surveying area
The WEM electromagnetic signal that geologic structure generates, achievable WEM electromagnetic signal observes data inversion on this basis, final to obtain oil
Gas ore field deep complex three-dimensional electrical structure.
2. a kind of method for carrying out D integral pin-fin tube based on WEM vector finite element method according to claim 1, feature exist
In considering ionosphere effect and displacement current when simulating to emission source and act on, it is thousands of to realize the transmitting covering of WEM emission system
The D integral pin-fin tube simulation of range, can be obtained each electromagnetic field information and apparent resistivity and impedance phase in kilometer.
3. a kind of method for carrying out D integral pin-fin tube based on WEM vector finite element method according to claim 1, feature exist
In 2040 high-density sampling filter factors are used in the method in fast Hankel transform, retain the shake of enough high frequencies
Information is swung, high band electromagnetic field arithmetic accuracy is improved, meets the needs of practical WEM field source signal imitation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811554392.0A CN109407160A (en) | 2018-12-19 | 2018-12-19 | A method of D integral pin-fin tube is carried out based on WEM vector finite element method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811554392.0A CN109407160A (en) | 2018-12-19 | 2018-12-19 | A method of D integral pin-fin tube is carried out based on WEM vector finite element method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109407160A true CN109407160A (en) | 2019-03-01 |
Family
ID=65460970
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811554392.0A Pending CN109407160A (en) | 2018-12-19 | 2018-12-19 | A method of D integral pin-fin tube is carried out based on WEM vector finite element method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109407160A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110346834A (en) * | 2019-07-22 | 2019-10-18 | 中国科学院地球化学研究所 | The forward modeling method of three-dimensional frequency domain controllable source electromagnetism, system |
| CN113505516A (en) * | 2021-02-25 | 2021-10-15 | 中国科学院地质与地球物理研究所 | Boundary fault interception method and device for three-dimensional forward modeling of low-frequency magnetotelluric method |
| US11378713B2 (en) * | 2019-07-25 | 2022-07-05 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method for collecting and processing tensor artificial-source electromagnetic signal data and device thereof |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109031432A (en) * | 2018-04-09 | 2018-12-18 | 中国科学院地质与地球物理研究所 | A kind of very low frequencies and magnetotelluric union measuring method |
-
2018
- 2018-12-19 CN CN201811554392.0A patent/CN109407160A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109031432A (en) * | 2018-04-09 | 2018-12-18 | 中国科学院地质与地球物理研究所 | A kind of very low frequencies and magnetotelluric union measuring method |
Non-Patent Citations (2)
| Title |
|---|
| 曹萌: "人工源极低频电磁法三维数据空间OCCAM反演研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
| 李勇 等: "基于二次场电导率分块连续变化的三维可控源电磁有限元数值模拟", 《地球物理学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110346834A (en) * | 2019-07-22 | 2019-10-18 | 中国科学院地球化学研究所 | The forward modeling method of three-dimensional frequency domain controllable source electromagnetism, system |
| CN110346834B (en) * | 2019-07-22 | 2020-11-17 | 中国科学院地球化学研究所 | Forward modeling method and system for three-dimensional frequency domain controllable source electromagnetism |
| US11378713B2 (en) * | 2019-07-25 | 2022-07-05 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method for collecting and processing tensor artificial-source electromagnetic signal data and device thereof |
| CN113505516A (en) * | 2021-02-25 | 2021-10-15 | 中国科学院地质与地球物理研究所 | Boundary fault interception method and device for three-dimensional forward modeling of low-frequency magnetotelluric method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | 3D vector finite-element electromagnetic forward modeling for large loop sources using a total-field algorithm and unstructured tetrahedral grids | |
| Chave et al. | The magnetotelluric method: Theory and practice | |
| Grayver et al. | 3D inversion and resolution analysis of land-based CSEM data from the Ketzin CO 2 storage formation | |
| Pedersen et al. | The gradient tensor of potential field anomalies: Some implications on data collection and data processing of maps | |
| Chang-Chun et al. | Review on airborne electromagnetic inverse theory and applications | |
| McMillan et al. | 3D parametric hybrid inversion of time-domain airborne electromagnetic data | |
| Liu et al. | 2D sequential inversion of total magnitude and total magnetic anomaly data affected by remanent magnetization | |
| Pedersen et al. | Groundwater exploration using combined controlled-source and radiomagnetotelluric techniques | |
| Cevallos et al. | Application of curvatures to airborne gravity gradient data in oil exploration | |
| Liang et al. | A new inversion method based on distorted born iterative method for grounded electrical source airborne transient electromagnetics | |
| Zhdanov et al. | Large-scale 3D inversion of marine magnetotelluric data: Case study from the Gemini prospect, Gulf of Mexico | |
| Xu et al. | How did the Dabie Orogen collapse? Insights from 3‐D magnetotelluric imaging of profile data | |
| CN109407160A (en) | A method of D integral pin-fin tube is carried out based on WEM vector finite element method | |
| Shockley | Ground vehicle navigation using magnetic field variation | |
| Nguyen et al. | Comparing large-scale 3D Gauss–Newton and BFGS CSEM inversions | |
| Dubey et al. | Computation of the gravity field and its gradient: Some applications | |
| Bang et al. | Imaging subsurface orebodies with airborne electromagnetic data using a recurrent neural network | |
| CN105891895B (en) | A kind of system and method determining sky wave propagation characteristic | |
| US9020205B2 (en) | Methods of multinary inversion for imaging objects with discrete physical properties | |
| Beiki et al. | Leveling HEM and aeromagnetic data using differential polynomial fitting | |
| Adrian et al. | Exploration of a copper ore deposit in elbistan/turkey using 2D inversion of the time-domain induced polarization data by using unstructured mesh | |
| Pirttijärvi et al. | Drone-based electromagnetic survey system for geophysical applications | |
| Teklesenbet | Multidimensional inversion of MT data from Alid Geothermal area, Eritrea. Comparison with geological structures and identification of a geothermal reservoir | |
| Wen et al. | Seismic structure-constrained inversion of CSAMT data for detecting karst caves | |
| Wang et al. | Application of the three-component transient electromagnetic method in gold exploration: a case study of the Heicigou gold deposit in east Kunlun |
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: 20190301 |
|
| RJ01 | Rejection of invention patent application after publication |