CN107783198B - A kind of magnetotelluric inverting data-bias imaging method - Google Patents

A kind of magnetotelluric inverting data-bias imaging method Download PDF

Info

Publication number
CN107783198B
CN107783198B CN201710817823.7A CN201710817823A CN107783198B CN 107783198 B CN107783198 B CN 107783198B CN 201710817823 A CN201710817823 A CN 201710817823A CN 107783198 B CN107783198 B CN 107783198B
Authority
CN
China
Prior art keywords
data
inverting
magnetotelluric
wave
electromagnetic
Prior art date
Application number
CN201710817823.7A
Other languages
Chinese (zh)
Other versions
CN107783198A (en
Inventor
李竹强
林治模
于会臻
尹克敏
胡加山
王雨洁
郝志伟
王有涛
唐付良
陈学国
Original Assignee
中国石油化工股份有限公司
中国石油化工股份有限公司胜利油田分公司勘探开发研究院
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 中国石油化工股份有限公司, 中国石油化工股份有限公司胜利油田分公司勘探开发研究院 filed Critical 中国石油化工股份有限公司
Priority to CN201710817823.7A priority Critical patent/CN107783198B/en
Publication of CN107783198A publication Critical patent/CN107783198A/en
Application granted granted Critical
Publication of CN107783198B publication Critical patent/CN107783198B/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves

Abstract

The present invention relates to a kind of magnetotelluric inverting data-bias imaging methods of petroleum and natural gas exploration, and magnetotelluric inverting data are carried out migration imaging, obtain migrated section, improve resolution ratio.The technical solution adopted by the present invention is that: obtain the magnetotelluric inverting section of high quality;Carry out wave field separation and transformation;Extract apparent reflectance;Generate migration model;Human-computer interaction obtains high-resolution migrated section.The electromagnetic inversion migrated section that the present invention obtains can preferably improve the resolution ratio of electric layers, can carry out joint interpretation with earthquake, to improve structure interpretation precision.

Description

A kind of magnetotelluric inverting data-bias imaging method

Technical field

The present invention relates to the processing of geophysical information in oil exploration and geological syntheses interpretation technique fields, and in particular to A kind of magnetotelluric inverting data-bias imaging method.

Background technique

Magneto-electrotelluric survey method in geophysical exploration usually is used to study one kind of earth interior electrical structure Important geophysics means can detect each point Electrical distribution from tens meters of ground to the hundreds of kms of earth mantle, to identification basin The geologic objectives such as ground structure, basin mountain contact relation, fracture play an important role.

Large size overlapping oil-gas bearing basin is after more phase tectonic movements and develops, and construction is complicated, and early stage rift valley or rift rank The strong magmatic activity of more phases generally occurs for section, saves numerous igneous lithofacies.Therefore in certain igneous rocks area, thick-layer Igneous rock has stronger shielding and absorption to seismic signal, so that seismic data igneous rock inside story is in blank or mixed and disorderly anti-more It penetrates, description and the target for constraining igneous rock are preferred.

The natural Magnetotelluric signal of MT method acquisition has not in oil exploration by shield layer with high resistivity, to high-conductivity layer point Strong, the lateral resolution capability of ability is distinguished compared with strong, exploration expenditure is low, the advantages such as easy for construction, but due to its bulk effect, longitudinal resolution Ability weakens rapidly with the increase of depth, and conventional inversion result resolution ratio is difficult to and the material matches such as earthquake, Bu Nengman Foot carries out the requirement of weight magnetoelectricity shake joint exploration research.

Currently used raising MT resolution of inversion approach is more: high precision collecting, high Precision Processing, joint inversion etc., But resolution ratio still cannot be satisfactory, still there is room for promotion.

Summary of the invention

Magnetotelluric is improved the purpose of the present invention is to provide a kind of magnetotelluric inverting data-bias imaging method to differentiate Rate;This method can effectively improve the longitudinal resolution of magnetotelluric inverting section, can be provided with earthquake, drilling well etc. by migration imaging Material carries out comprehensive study, improves accuracy of identification and reliability, reduces exploration risk.

In order to achieve the above objectives, the present invention adopts the following technical scheme that, a kind of magnetotelluric inverting data-bias imaging Its step of method includes:

1) magnetotelluric high Precision Processing technology is utilized, high quality two dimensional inversion data section is obtained;

2) according to polarization mode, wave field separation is carried out, the mask data of different polarization modes is obtained;

3) by isolated electromagnetic data, wave field transformation is carried out, apparent resistivity and impedance phase data are changed into wave field;

4) wave field converted out is handled, inverting obtains apparent emission rate;

5) inverting data and apparent reflectance are subjected to joint inversion, generate migration model;

6) migration imaging is carried out to migration model, obtains high-resolution electromagnetic inversion migrated section.

Further, the step 1) comprises the concrete steps that: (1) carrying out accurate pre-processing, the fine pre- place to MT data Reason includes flying spot and static shift correction;(2) carry out TE TM high-precision two-dimensional quantitative inversion;(3) two dimension joint high-precision is carried out Inverting obtains the electrical inverting section of high quality two dimension;I.e. first with drilling well, the electrically data such as measurement of appearing, research area is built Vertical region Electrical Models;Then carry out magnetotelluric acquired original data and carry out flying spot correction and static shift correction;It is finally right Pretreatment data carries out one-dimensional, two-dimentional joint inversion and obtains inverting data section.

Further, the step 2) comprises the concrete steps that: selection polarization mode first carries out wave field separation respectively, obtains TE Polarity electromagnetic reflectance data, TM polarization mode ELECTROMAGNETIC REFLECTION data, TEM polarity electromagnetic reflectance data;Then the polarization side TE is obtained The apparent reflectance data of formula, the apparent reflectance data of TM polarization mode, the apparent reflectance data of TEM polarization mode.

Further, the step 3) comprises the concrete steps that: carry out wave field separation for the polarization mode in step 2 respectively, It needs to reasonably select the parameters such as regularization factors, time sampling points, time sampling interval, obtains TE polarity electromagnetic transformation wave, TM polarity electromagnetic converts wave, TE and TM polarity electromagnetic converts wave, using electromagnetic wave slow wave propagation equation:

Using the fast wave propagation equations of electromagnetic wave:

Propagation equation is moved using electromagnetic field all-wave:

Further, the step 4) comprises the concrete steps that: first the wave field in step 3 being edited, denoised and is mentioned High-resolution Processing;And deconvolution is carried out respectively, propose approximate apparent reflectance;Then approximate apparent reflectance is carried out non-linear anti- Drill acquisition more accurately apparent emission rate.

Further, the step 5) comprises the concrete steps that: obtaining high-precision firstly the need of two-dimentional continuous media inverting is carried out Inverting data volume;Then preference pattern depth is consistent with inverting depth parameter, selects reasonable sampling step length, according to electromagnetic vector Concerning With Fast-slow Waves wave theory obtains Pseudo velocity model, as migration model.

Further, the step 6) comprises the concrete steps that: first migration model data .dat is converted into .sgy data, It needs to select suitable wavelet, correction reflectivity, obtains stable migration model;The migration model .sgy number that will finally obtain According to progress migration imaging obtains high-resolution inverting migrated section.

Inverting electromagnetism section is converted to Pseudo velocity section migration imaging by the present invention, realizes conventional magnetotelluric inverting number According to high-resolution offset imaging method, the inversion imaging precision of MT data is improved, can be widely applied to oil-gas exploration Middle Fine structural interpretation and joint interpretation of gravel.

Detailed description of the invention

Fig. 1 is a kind of flow diagram of magnetotelluric inverting data-bias imaging method of the present invention.

Fig. 2 is MT data inversion process figure in the embodiment of the present invention.

Fig. 3 is MQ-E048 magnetotelluric inverting sectional view in the embodiment of the present invention.

Fig. 4 is that MQ-E048 extracts schematic diagram to certain road electromagnetic wave data apparent resistivity in the embodiment of the present invention.

Fig. 5 is MQ-E048 magnetotelluric inverting data-bias imaging effect figure in the embodiment of the present invention.

Specific embodiment

Detailed description for the present invention and technology contents, cooperation Detailed description of the invention it is as follows, however attached drawing only provides reference with Purposes of discussion is not intended to limit the present invention.

Fig. 1 is a kind of flow diagram of magnetotelluric inverting data-bias imaging method of the present invention, as shown in Figure 1, the party Method the following steps are included:

The first step establishes work area Electrical Models in conjunction with first branch apparent resistivity information using drilling well, the electrical information such as appear, Carry out firsthand information pretreatment.In this case by taking the work area MQ as an example, it is contemplated that the influence of power plant, landform needs to carry out fine pre- place Reason, eliminate flying spot, static displacement, then carry out TE TM polarization mode inverting work, obtain the inverting sectional view 3 of high quality, Specific flow chart is as shown in Figure 2.

Second step carries out wave field separation, obtains the mask data of different polarization modes according to polarization mode: firstly the need of To the TE obtained in the first step the data filtering of TM polarity electromagnetic, polarization data interpolation, polarization data restore impedance phase, one-dimensional Two-dimentional continuous media inverting data are identified, according to Regional Geologic Statistical error TE TM polarization mode;Secondly, carrying out Wave field separation obtains into the pole TE polarity electromagnetic reflectance data ud_te.dat, TM polarity electromagnetic reflectance data ud_tm.dat, TEM Change apparent reflectance data reflect_te.dat, the TM polarization mode of ELECTROMAGNETIC REFLECTION data ud_tem.dat, TE polarization mode The apparent reflectance data reflect_tem.dat of apparent reflectance data reflect_tm.dat, TEM polarization mode. ud_ Te.dat, ud_tm.dat and ud_tem.dat prepare for " wave field transformation " and " extracting apparent reflectance " operation.

Isolated electromagnetic data is carried out wave field transformation, apparent resistivity and impedance phase data is changed into fluctuation by third step .It is comprised the concrete steps that: firstly, setting transformation parameter, wherein polarization mode selects joint inversion, and regularization factors are set as 2* 10-7, the time is using number 5000, Temporal sampling 2ms;Then to second step obtain TE TM polarity electromagnetic reflectance data, into Traveling wave field transformation obtains TE polarization wave field transformation wave file f _ w_te.dat, TM polarization wave field transformation wave file f _ w_tm.dat And TEM polarization wave field transformation wave file f _ w_tem.dat, it prepares for the operation of " human-computer interaction promotion resolution ratio " module.

4th step is handled the wave field converted out, inverting obtains more accurately its specific steps of apparent emission rate It is: since the wave field resolution ratio converted out is very low, so to propose approximate apparent reflectance by deconvolution, then by non-thread Property inverting obtain more accurately apparent reflectance.By selection polarization mode and other parameters, corresponding data is inputted, is extracted Apparent reflectance operation, generates new reflect_te.dat, reflect_tm.dat, reflect_tem.dat file as defeated File out, to prepare in next step.

Inverting data and apparent reflectance are carried out joint inversion, generate migration model by the 5th step.It is comprised the concrete steps that: first First, formation speed model parameter is set, model depth must be consistent with front inverting depth profiled, and 10 meters of depth-sampling step-length, Seek the corresponding rate pattern of section, V_model_te.dat, V_model_tm.dat, V_model_Jem.dat etc..Then According to earthquake data format, SEGY format is converted to rate pattern, SEGY file consists of two parts, file header, data volume. Data volume is made of data together, and each track data body is made of trace header and sampled data two parts again, by f_w_ Te.dat is converted to f_w_te.sgy file.Then it modifies reflectivity and FORWARD AND INVERSE PROBLEMS to a certain track data of electromagnetic field data Work, Fig. 4 select 120 track datas, wavelet frequency 5hz, and complementary operation adjusts reflection coefficient, synthesizes new reflection coefficient and road Data, and then carry out the work to form migration model to complete section face.

6th step carries out migration imaging to migration model, obtains high-resolution electromagnetic inversion migrated section.It is specifically walked Suddenly it is: migration model data will be obtained, be loaded into earthquake correlation and show software, carries out migration imaging and show, to obtain high score The electromagnetic inversion migration imaging section of resolution.

Fig. 5 is the MQ-E048 section electromagnetism result map schematic diagram obtained using offset imaging method, is come from imaging results It sees, hence it is evident that than the better longitudinal resolution of former inverting sectional view 3 and lateral continuity, it was demonstrated that the reliability of the technology.

Above-described specific embodiment has carried out further the purpose of the present invention, technical scheme and beneficial effects Explanation, it should be understood that being not intended to limit the present invention the foregoing is merely a specific embodiment of the invention Protection scope, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in Within protection scope of the present invention.

Claims (6)

1. a kind of magnetotelluric inverting data-bias imaging method, which comprises the steps of:
1) magnetotelluric high Precision Processing technology is utilized, high quality two dimensional inversion data section is obtained;
2) according to polarization mode, wave field separation is carried out, the mask data of different polarization modes is obtained;
3) by isolated electromagnetic data, wave field transformation is carried out, apparent resistivity and impedance phase data are changed into wave field;
4) wave field converted out is handled, inverting obtains apparent reflectance;
5) inverting data and apparent reflectance are subjected to joint inversion, generate migration model;
6) migration model is handled, migration imaging, obtains high-resolution electromagnetic inversion migrated section;
It is further described in the step 3) are as follows: carry out wave field separation for the polarization mode in step 2 respectively, need rationally to select Regularization factors, time sampling points, time sampling interval parameter are selected, obtains TE polarity electromagnetic transformation wave, TM polarity electromagnetic becomes Wave is changed, TE and TM polarity electromagnetic converts wave, and equation is as follows:
Using electromagnetic wave slow wave propagation equation:
Using the fast wave propagation equations of electromagnetic wave:
Propagation equation is moved using electromagnetic field all-wave:
2. magnetotelluric inverting data-bias imaging method according to claim 1, which is characterized in that in the step 1) It is further broken into:
2.1) accurate pre-processing is carried out to MT data, the accurate pre-processing includes flying spot and static shift correction;
2.2) carry out TE TM high-precision two-dimensional quantitative inversion;
2.3) two dimension joint high-precision inverting is carried out, the electrical inverting section of high quality two dimension is obtained.
3. magnetotelluric inverting data-bias imaging method according to claim 1, which is characterized in that in the step 2 It is further broken into:
3.1) selection polarization mode carries out wave field separation respectively, and it is anti-to obtain TE polarity electromagnetic reflectance data, TM polarization mode electromagnetism Penetrate data and TEM polarity electromagnetic reflectance data;
3.2) selection polarization mode carries out wave field separation respectively, obtains apparent reflectance data, the TM polarization mode of TE polarization mode Apparent reflectance data and TEM polarization mode apparent reflectance data.
4. magnetotelluric inverting data-bias imaging method according to claim 1, which is characterized in that in the step 4) It is further broken into:
4.1) resolution processes are edited, denoised and is improved to the wave field in step 3, and carry out deconvolution respectively, proposed Approximate apparent reflectance;
4.2) approximate apparent reflectance is subjected to non-linear inversion and obtains more accurately apparent reflectance.
5. magnetotelluric inverting data-bias imaging method according to claim 1, which is characterized in that in the step 5) It is further broken into:
5.1) it needs to carry out two-dimentional continuous media inverting and obtains high-precision inverting data volume;
5.2) preference pattern depth is consistent with inverting depth parameter, selects reasonable sampling step length, according to electromagnetic vector Concerning With Fast-slow Waves wave Dynamic theory, obtains Pseudo velocity model, as migration model.
6. a kind of magnetotelluric inverting data-bias imaging method according to claim 1, which is characterized in that the step 6) it is further broken into:
6.1) migration model data .dat is converted into .sgy data, needs to select suitable wavelet, correction reflectivity, obtains steady Fixed migration model;
6.2) the migration model .sgy data that will be obtained carry out migration imaging, obtain high-resolution inverting migrated section.
CN201710817823.7A 2017-09-12 2017-09-12 A kind of magnetotelluric inverting data-bias imaging method CN107783198B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710817823.7A CN107783198B (en) 2017-09-12 2017-09-12 A kind of magnetotelluric inverting data-bias imaging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710817823.7A CN107783198B (en) 2017-09-12 2017-09-12 A kind of magnetotelluric inverting data-bias imaging method

Publications (2)

Publication Number Publication Date
CN107783198A CN107783198A (en) 2018-03-09
CN107783198B true CN107783198B (en) 2019-11-08

Family

ID=61438268

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710817823.7A CN107783198B (en) 2017-09-12 2017-09-12 A kind of magnetotelluric inverting data-bias imaging method

Country Status (1)

Country Link
CN (1) CN107783198B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010022296A1 (en) * 2008-08-20 2010-02-25 Lockheed Martin Corporation System and method to measure and track fluid movement in a reservoir using electromagnetic transmission
CN105116456A (en) * 2015-07-23 2015-12-02 吉林大学 Combined imaging method of structure and polarized attribute of underground object

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010022296A1 (en) * 2008-08-20 2010-02-25 Lockheed Martin Corporation System and method to measure and track fluid movement in a reservoir using electromagnetic transmission
CN105116456A (en) * 2015-07-23 2015-12-02 吉林大学 Combined imaging method of structure and polarized attribute of underground object

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Damage imaging of reinforced concrete structures using electromagnetic migration algorithm;Saeed Nojavan,et al.;《International Journal of Solids and Structures》;20061231;第5886-5908页 *
MT法电磁相位移偏移研究;李吉松,等;《石油地球物理勘探》;19940430;第189-200页 *
Underground imaging by frequency-domain electromagnetic migration;Michael S. Zhdanov,et al.;《GEOPHYSICS》;19961231;第666-682页 *
大地电磁偏移成像技术的研究进展及其解决复杂地质问题的能力;于鹏,等;《地震地质》;20010630;第239-242页 *
大地电磁全频段偏移成像结合反演技术;于鹏,等;《同济大学学报(自然科学版)》;20050430;第540-544页 *
大地电磁场成像方法综述与新进展;于鹏,等;《地球物理学进展》;20030331;第53-58页 *
有限差分法大地电磁多参数偏移成像;于鹏,等;《地球物理学报》;20010731;第552-562页 *

Also Published As

Publication number Publication date
CN107783198A (en) 2018-03-09

Similar Documents

Publication Publication Date Title
Conyers Ground-penetrating radar for archaeology
Dix Seismic velocities from surface measurements
Jongmans et al. Geophysical investigation of landslides: a review
Nabighian The analytic signal of two-dimensional magnetic bodies with polygonal cross-section: its properties and use for automated anomaly interpretation
Gawthorpe et al. Ground penetrating radar: application to sandbody geometry and heterogeneity studies
Smith et al. Gridding with continuous curvature splines in tension
Selim et al. Shallow seismic refraction, two-dimensional electrical resistivity imaging, and ground penetrating radar for imaging the ancient monuments at the western shore of Old Luxor city, Egypt
Jongmans et al. Geophysical investigation of a large landslide in glaciolacustrine clays in the Trièves area (French Alps)
Nabighian et al. The historical development of the magnetic method in exploration
Grasmueck et al. Full-resolution 3D GPR imaging
Jongmans et al. EURO-SEISTEST: determination of the geological structure of the Volvi basin and validation of the basin response
Colella et al. High-resolution imaging of the High Agri Valley Basin (Southern Italy) with electrical resistivity tomography
Soupios et al. Use of engineering geophysics to investigate a site for a building foundation
CN104977618B (en) A kind of method evaluated shale gas reservoir and find dessert area
Forte et al. Integrated seismic tomography and ground-penetrating radar (GPR) for the high-resolution study of burial mounds (tumuli)
Dabas et al. Simultaneous use of electrostatic quadrupole and GPR in urban context: investigation of the basement of the Cathedral of Girona (Catalunya, Spain)
Boeniger et al. Improving the interpretability of 3D GPR data using target–specific attributes: application to tomb detection
Wilson et al. From outcrop to flow simulation: Constructing discrete fracture models from a LIDAR survey
Wisén et al. Laterally and mutually constrained inversion of surface wave seismic data and resistivity data
Stanley et al. Regional magnetotelluric surveys in hydrocarbon exploration, Parana Basin, Brazil
Longde et al. Petroleum exploration and development practices of sedimentary basins in China and research progress of sedimentology
Yordkayhun et al. 3D seismic traveltime tomography imaging of the shallow subsurface at the CO2SINK project site, Ketzin, GermanyCO2SINK 3D seismic traveltime tomography
Guérin et al. Geophysical characterisation of karstic networks–Application to the Ouysse system (Poumeyssen, France)
Vallee et al. Metalliferous mining geophysics—State of the art after a decade in the new millennium
Porsani et al. Ground-penetrating radar profiles over multiple steel tanks: Artifact removal through effective data processing

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
CB03 Change of inventor or designer information

Inventor after: Li Zhuqiang

Inventor after: Chen Xueguo

Inventor after: Lin Zhimo

Inventor after: Yu Huizhen

Inventor after: Yin Kemin

Inventor after: Hu Jiashan

Inventor after: Wang Yujie

Inventor after: Hao Zhiwei

Inventor after: Wang Youtao

Inventor after: Tang Fuliang

Inventor before: Li Zhuqiang

Inventor before: Wei Min

Inventor before: Guo Tao

Inventor before: Feng Guozhi

Inventor before: Wang Yujie

Inventor before: Shi Xiupeng

Inventor before: Zhang Jianhua

Inventor before: Qian Huanju

Inventor before: Yang Guojie

Inventor before: Wang Yuelei

Inventor before: Li Qin

Inventor before: Chen Xueguo

Inventor before: Wang Shuhua

Inventor before: Xiang Peng

Inventor before: Lin Huixi

Inventor before: Tan Shaoquan

Inventor before: Lin Zhimo

Inventor before: Xu Youde

Inventor before: Yu Huizhen

Inventor before: Wang Youtao

CB03 Change of inventor or designer information
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant