CN107272067A - Method is found oil in a kind of optically-active - Google Patents
Method is found oil in a kind of optically-active Download PDFInfo
- Publication number
- CN107272067A CN107272067A CN201710413368.4A CN201710413368A CN107272067A CN 107272067 A CN107272067 A CN 107272067A CN 201710413368 A CN201710413368 A CN 201710413368A CN 107272067 A CN107272067 A CN 107272067A
- Authority
- CN
- China
- Prior art keywords
- frequency
- electric field
- data
- oil
- observation
- 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.)
- Granted
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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Method is found oil the present invention relates to a kind of optically-active, is comprised the following steps：S1：An observation system is provided, the observation system is made up of main frame, three magnetic probes, four electrodes and gps antenna；S2：The electromagnetic field data of acquisition time sequence, the sample frequency of observation data is 15Hz, and the Continuous Observation time is more than 4 hours, and the electromagnetic field data includes five track datas：Ex (t), Ey (t), Hx (t), Hy (t) and Hz (t)；S3：The electromagnetic field observation Ex (t) of time-domain, Ey (t), Hx (t), Hy (t) are changed into electromagnetic field Ex (f), Ey (f), Hx (f), the Hy (f) of frequency domain, by the way of segmentation conversion, the result that multistage is changed is overlapped acquisition average value；S4：Electric field intensity azimuth is calculated according to calculation formula θ (f)=(180/3.1415926) arctan [︱ Ey (f) ︱/︱ Ex (f) ︱]；S5：Draw the azimuthal frequency response curve of electric field intensity and explain oiliness.The optically-active oil prospecting method the step of provide can direct detection go out the position of oil-gas reservoir, reduce exploration cost.
Description
Technical field
The present invention relates to field of geophysical exploration, more particularly to a kind of optically-active oil prospecting method.
Background technology
In petroleum exploration field, applicable geophysical method has earthquake, gravity, magnetic force and electrical method method, seismic method
Based on the elastic difference of material, gravity based on density variation, magnetic force based on magnetic contrast, electrical method is then with underground
Based on the resistivity contrasts of medium.These methods are all to study oil gas field indirectly by analyzing geologic structure, and can not be straight
Connect the presence for judging oil gas field.The specific method of electrical method method include a variety of mutation, such as induced polarization method, complex resistivity method and when
Frequency electromagnetic method.Induced polarization method, complex resistivity method and time-frequency electromagnetic method are in addition to being capable of the resistivity of study of rocks, moreover it is possible to grind
Study carefully the polarizability parameter of rock.Polarizability is a parameter for reflecting effect of induced polarization, and the effect mainly reflects troilite,
The not substantive characteristics of oil-gas reservoir.Sum it up, up to the present, also none of these methods being capable of direct detection oil-gas reservoir.
The content of the invention
The purpose of the present invention is can to go out the position of oil-gas reservoir there is provided a kind of with direct detection for above-mentioned present situation, is dropped simultaneously
The optically-active oil prospecting method of low exploration cost.
The technical solution adopted by the present invention：Method is found oil in a kind of optically-active, is comprised the following steps：S1：One observation system is provided,
The observation system is made up of main frame, three magnetic probes, four electrodes and gps antenna, and four electrodes are used to observe two mutually
Vertical electric field component Ex, Ey, three magnetic probes are used to observing orthogonal three magnetic-field components Hx, Hy, Hz, wherein Ex with
Hx is parallel, and Ey is parallel with Hy；S2：The electromagnetic field data of acquisition time sequence, the sample frequency of observation data is 15Hz, continuous to see
The survey time is more than 4 hours, and the electromagnetic field data includes five track datas：Ex (t), Ey (t), Hx (t), Hy (t) and Hz (t),
Wherein E represents electric field, and H represents magnetic field, and x represents North and South direction reference axis, is northwards positive direction, y represents east-west direction reference axis,
It is positive direction eastwards, z represents reference axis vertically downward, and t represents the time；S3：Using Fast Fourier Transform (FFT), to time-domain
Electromagnetic field observation Ex (t), Ey (t), Hx (t), Hy (t) change into electromagnetic field Ex (f), Ey (f), Hx (f), the Hy of frequency domain
(f), by the way of segmentation conversion, the result that multistage is changed is overlapped acquisition average value；S4：After being superimposed in frequency domain
Electric field component Ex (f) and Ey (f), according to calculation formula θ (f)=(180/3.1415926) arctan [︱ Ey (f) ︱/︱ Ex
(f) ︱] to calculate θ (f) in electric field intensity azimuth, formula be the azimuth that frequency is f electric field intensitys；And S5：Sat using frequency to be horizontal
Mark, electric field intensity azimuth is ordinate, draws the azimuthal frequency response curve of electric field intensity.
The effect of the technical program is：The optically-active oil prospecting method the step of provide can direct detection go out the position of oil-gas reservoir
Put, reduce exploration cost.
Further, in step S1, the main frame is V5-2000 electromagnetism measurement apparatus.
Further, in step S2, the Continuous Observation time is 8 hours.
Further, in step S3, the quantity for the data point chosen is changed every time into N, if sample frequency is Fs, it is superimposed
Ex (f) afterwards is a data sequence, and the frequency that i-th of data (i is since 1) represents is：F=(i-1) * Fs/N；For handing over
For flowing signal, the modulus value divided by N/2 of i-th complex data are exactly the amplitude to signal that should be under frequency, are believed for direct current
For number, the modulus value divided by N of i-th complex data are exactly the amplitude to signal that should be under frequency, the frequency of the 1st data point
For 0, expression direct current signal, the frequency Fs/N of the 2nd data point, the frequency 2Fs/N of the 3rd, until N/2 data point.
Brief description of the drawings
Fig. 1 is optical activity schematic diagram, wherein, 1 represents Electromagnetic Wave Propagation direction, and 2 represent direction of an electric field, and 3 represent oil reservoir, 4
Direction of an electric field during without oil reservoir is represented, 5 represent direction of an electric field when having oil reservoir.
Fig. 2 is the observation system of data acquisition.
The electric field anglec of rotation-frequency curve (1 Hangzhou somewhere when Fig. 3 is without oil reservoir；2 Hebei Nanpi County somewhere).
Fig. 4 is the electric field anglec of rotation-frequency curve when having oil reservoir.
Embodiment
The principle and feature of the present invention are described below in conjunction with accompanying drawing, the given examples are served only to explain the present invention, and
It is non-to be used to limit the scope of the present invention.
Referring to Fig. 1, be the schematic diagram that method is found oil in a kind of optically-active for providing of the present invention, i.e. electromagnetic wave after oil reservoir,
The direction of electric field will rotate a certain angle.Wherein 1 is Electromagnetic Wave Propagation direction；2 be direction of an electric field；3 be oil reservoir；4 be oil-free
Direction of an electric field during Tibetan；5 be direction of an electric field when having oil reservoir.
Optically-active oil prospecting method comprises the following steps：
1. the collection of time series electric field data
An observation system is provided, referring to Fig. 2, the observation system was by main frame, three magnetic probes, four electrodes and GPS days
Line is constituted, and four electrodes are used to observe two orthogonal electric field component Ex, Ey, and three magnetic probes are mutually perpendicular to for observation
Three magnetic-field components Hx, Hy, Hz.Ex is parallel with Hx, and Ey is parallel with Hy.Specifically, the main frame is Canadian phoenix earth thing
The V5-2000 electromagnetism measurement apparatus of reason company, the electromagnetic field of the V5-2000 electromagnetism measurement apparatus recordable time sequences.
Operation is as follows：On 028 measuring point of the L01 surveys line in some work area, the sample frequency of observation data is 15Hz, i.e., every
Second records equally spaced 15 data, and the Continuous Observation time is more than 4 hours, can obtain three data file L01-
028A.TBL, L01-028A.TSH and L01-028A.TSL.L01-028A.TBL is gauge outfit file, have recorded the one of field construction
A little information.L01-028A.TSH is electromagnetic field of high frequency record data, and sample frequency is 2400Hz, and L01-028A.TSL is low-frequency electrical
Magnetic field data, sample frequency is 15Hz.
Electromagnetic field data includes five track datas：Ex (t), Ey (t), Hx (t), Hy (t), Hz (t), wherein E represent electric field, H
Magnetic field is represented, x represents North and South direction reference axis, is northwards positive direction, y represents east-west direction reference axis, is eastwards positive direction, z
Reference axis vertically downward is represented, t represents the time.Ex (t) means that the electric field observation in moment t North and South direction, according to this class
Push away.
In the present embodiment, the Continuous Observation time is 8 hours.
2. conversion of the electric field time-domain to frequency domain
Using Fast Fourier Transform (FFT), the electromagnetic field observation Ex (t) of time-domain, Ey (t), Hx (t), Hy (t) are converted
Into the electromagnetic field Ex (f) of frequency domain, Ey (f), Hx (f), Hy (f).By the way of segmentation conversion, the purpose of segmentation conversion is pressure
Random noise processed, improves the reliability of data, the result that multistage is changed is overlapped into acquisition average value.Conversion is chosen every time
The quantity of data point be N, if sample frequency is Fs, the Ex (f) after superposition is also a data sequence, multiple comprising 1024
Number, plural angle represents phase, and the frequency that i-th of data (i is since 1) represents is：F=(i-1) * Fs/N；I-th of plural number
The modulus value divided by N/2 of data are exactly the amplitude (being divided by N) to signal that should be under frequency for direct current signal.1st data
The frequency of point is 0, represents direct current signal, the frequency Fs/N of the 2nd data point, the frequency 2Fs/N of the 3rd, until N/2 data
Point, thereafter data without.
3. electric field intensity is azimuthal calculates
Electric field component Ex (f) and Ey (f) after being superimposed in frequency domain, calculate electric field intensity azimuth.Its public affairs calculated
Formula is that θ (f) is that frequency is f electric field intensitys in θ (f)=(180/3.1415926) arctan [︱ Ey (f) ︱/︱ Ex (f) ︱], formula
Azimuth.Namely with negating tangent, the azimuth that acquisition is represented with radian again after Ey (f) mould divided by Ex (f) mould.
4. Drawing of Curve and the interpretation and evaluation of oiliness
Using frequency as abscissa, electric field intensity azimuth is ordinate, draws the azimuthal frequency response of electric field intensity bent
Line, as shown in Figure 3 and Figure 4.In figure, 8Hz nearby curve to it is recessed be that because Schumann resonance, Schumann resonance is the earth and ionosphere
The inherent characteristic of the resonator of composition.Fig. 3 be Hangzhou somewhere and two, Hebei Nanpi County somewhere observation station the electric field anglec of rotation with
Frequency and the curve map changed, electric field data is in same time segment record.It can be seen that frequency 1. by 0.5Hz to
7Hz, being basically unchanged of angle of electric field represents do not have oil in the range of the certain depth of underground.Electric field when Fig. 4 is Xinjiang oil field
The anglec of rotation-frequency curve, from small to large, the angle of electric field tapers to 40 degree by 55 degree to frequency, and change is due to underground
Oil caused by, actually the observation station be above some oil field, it is theoretical to be coincide with actual, illustrate that method is found oil in optically-active
Validity.
In this way, the optically-active oil prospecting method the step of provide can direct detection go out the position of oil-gas reservoir, reduce and explore into
This.
The foregoing is only presently preferred embodiments of the present invention, be not intended to limit the invention, it is all the present invention spirit and
Within principle, any modification, equivalent substitution and improvements made etc. should be included in the scope of the protection.
Claims (4)
1. method is found oil in a kind of optically-active, it is characterised in that comprised the following steps：
S1：An observation system is provided, the observation system is made up of main frame, three magnetic probes, four electrodes and gps antenna, four
Individual electrode is used to observe two orthogonal electric field component Ex, Ey, and three magnetic probes are used to observe orthogonal three magnetic
Field component Hx, Hy, Hz, wherein Ex are parallel with Hx, and Ey is parallel with Hy；
S2：The electromagnetic field data of acquisition time sequence, the sample frequency of observation data is 15Hz, and the Continuous Observation time is more than 4
Hour, the electromagnetic field data includes five track datas：Ex (t), Ey (t), Hx (t), Hy (t) and Hz (t), wherein E represent electric field,
H represents magnetic field, and x represents North and South direction reference axis, is northwards positive direction, y represents east-west direction reference axis, is eastwards positive direction, z
Reference axis vertically downward is represented, t represents the time；
S3：Using Fast Fourier Transform (FFT), the electromagnetic field observation Ex (t) of time-domain, Ey (t), Hx (t), Hy (t) are changed into
Electromagnetic field Ex (f), Ey (f), Hx (f), the Hy (f) of frequency domain, by the way of segmentation conversion, the result that multistage is changed is carried out
Superposition obtains average value；
S4：Electric field component Ex (f) and Ey (f) after being superimposed in frequency domain, according to calculation formula θ (f)=(180/
3.1415926) it is that frequency is f electric field intensitys that arctan [︱ Ey (f) ︱/︱ Ex (f) ︱], which calculates θ (f) in electric field intensity azimuth, formula,
Azimuth；
S5：Using frequency as abscissa, electric field intensity azimuth is ordinate, draws the azimuthal frequency response of electric field intensity bent
Line, the curve contains the information that oil-gas reservoir whether there is.
2. method is found oil in optically-active as claimed in claim 1, it is characterised in that：In step S1, the main frame is V5-2000 electromagnetism
Measuring instrument.
3. method is found oil in optically-active as claimed in claim 1, it is characterised in that：In step S2, the Continuous Observation time is 8 small
When.
4. method is found oil in optically-active as claimed in claim 1, it is characterised in that：In step S3, the data point that conversion is chosen every time
Quantity be N, if sample frequency be Fs, the Ex (f) after superposition is a data sequence, i-th of data (i is since 1) table
The frequency shown is：F=(i-1) * Fs/N；For AC signal, the modulus value divided by N/2 of i-th complex data are exactly correspondence
The amplitude of signal under the frequency, for direct current signal, the modulus value divided by N of i-th complex data are exactly to should frequency
Under signal amplitude, the frequency of the 1st data point is 0, represents direct current signal, the frequency Fs/N of the 2nd data point, the 3rd
Frequency 2Fs/N, until N/2 data point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710413368.4A CN107272067B (en) | 2017-06-05 | 2017-06-05 | A kind of optically-active oil prospecting method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710413368.4A CN107272067B (en) | 2017-06-05 | 2017-06-05 | A kind of optically-active oil prospecting method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107272067A true CN107272067A (en) | 2017-10-20 |
| CN107272067B CN107272067B (en) | 2019-06-18 |
Family
ID=60065754
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710413368.4A Active CN107272067B (en) | 2017-06-05 | 2017-06-05 | A kind of optically-active oil prospecting method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107272067B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1580818A (en) * | 2003-08-01 | 2005-02-16 | 中国石油天然气集团公司 | Artificial source time frequency electro magnetic bathymetry |
| US20100212909A1 (en) * | 2009-02-20 | 2010-08-26 | Anatoly Baumstein | Method For Analyzing Multiple Geophysical Data Sets |
| WO2013104888A1 (en) * | 2012-01-10 | 2013-07-18 | The Secretary Of State For Defence | Signature identification |
| CN103809203A (en) * | 2012-11-06 | 2014-05-21 | 核工业北京地质研究院 | Outdoor detection working method resolving electromagnetic interference |
-
2017
- 2017-06-05 CN CN201710413368.4A patent/CN107272067B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1580818A (en) * | 2003-08-01 | 2005-02-16 | 中国石油天然气集团公司 | Artificial source time frequency electro magnetic bathymetry |
| US20100212909A1 (en) * | 2009-02-20 | 2010-08-26 | Anatoly Baumstein | Method For Analyzing Multiple Geophysical Data Sets |
| WO2013104888A1 (en) * | 2012-01-10 | 2013-07-18 | The Secretary Of State For Defence | Signature identification |
| CN103809203A (en) * | 2012-11-06 | 2014-05-21 | 核工业北京地质研究院 | Outdoor detection working method resolving electromagnetic interference |
Non-Patent Citations (1)
| Title |
|---|
| 何展翔: "圈定油气藏边界的井地电法研究", 《中国优秀博硕士学位论文全文数据库（博士）基础科学辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107272067B (en) | 2019-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Streich | Controlled-source electromagnetic approaches for hydrocarbon exploration and monitoring on land | |
| US7362102B2 (en) | Electromagnetic surveying for resistive or conductive bodies | |
| Wirianto et al. | A feasibility study of land CSEM reservoir monitoring in a complex 3-D model | |
| Cho et al. | Miniature array analysis of microtremors | |
| AU2006214069B2 (en) | System and method for using time-distance characteristics in acquisition, processing and imaging of t-CSEM data | |
| US8762062B2 (en) | Method for marine geoelectrical exploration with electrical current focusing | |
| CN103995301A (en) | Method and device for evaluating total organic carbon content in shale gas reservoir | |
| Guérin et al. | Geophysical characterisation of karstic networks–Application to the Ouysse system (Poumeyssen, France) | |
| CN114236621A (en) | Multi-technology cooperation based nondestructive detection method for potential hazards of dikes | |
| Pueyo Anchuela et al. | Geophysical characterization of buried active faults: the Concud Fault (Iberian Chain, NE Spain) | |
| Zhang et al. | A new algorithm for frequency-dependent shear-wave splitting parameters extraction | |
| Rahimi et al. | Application of microtremor horizontal to vertical spectra ratio (MHVSR) and multichannel analysis of surface wave (MASW) for shallow bedrock mapping for transportation projects | |
| CN105807326B (en) | The system and method that a kind of utilization sky wave carries out deep prospecting | |
| RU2733095C2 (en) | Method of three-dimensional objects searching by tm-polarization geoelectrics methods | |
| CN107272067B (en) | A kind of optically-active oil prospecting method | |
| Karamzadeh et al. | Small-aperture array as a tool to monitor fluid injection-and extraction-induced microseismicity: applications and recommendations | |
| Wilt et al. | Cross-borehole and surface-to-borehole electromagnetic induction for reservoir characterization | |
| RU2408036C1 (en) | Focused current marine geoelectric prospecting method | |
| Lahti | Audiomagnetotelluric (AMT) measurements: A new tool for mineral exploration and upper crustal research at the Geological Survey of Finland | |
| Thiel et al. | Eucla MT profile and modelling | |
| Amjadi et al. | An application of magnetotelluric data inversion in a stratovolcano region | |
| Saraev et al. | Multiparameter monitoring of electromagnetic earthquake precursors in the frequency range 0.1 Hz-1 MHz | |
| Summers | On the frequency response of induction anomalies | |
| RU2446417C2 (en) | Three-dimensional frequency-time electrical prospecting method (ftem-3d) | |
| Hickey | Application of Land-Based Controlled Source EM Method to Hydraulic Fracture Monitoring |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |