CN107272067B - A kind of optically-active oil prospecting method - Google Patents

A kind of optically-active oil prospecting method Download PDF

Info

Publication number
CN107272067B
CN107272067B CN201710413368.4A CN201710413368A CN107272067B CN 107272067 B CN107272067 B CN 107272067B CN 201710413368 A CN201710413368 A CN 201710413368A CN 107272067 B CN107272067 B CN 107272067B
Authority
CN
China
Prior art keywords
frequency
electric field
data
optically
oil
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
Application number
CN201710413368.4A
Other languages
Chinese (zh)
Other versions
CN107272067A (en
Inventor
陈紫静
胡文宝
陈清礼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yangtze University
Original Assignee
Yangtze University
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 Yangtze University filed Critical Yangtze University
Priority to CN201710413368.4A priority Critical patent/CN107272067B/en
Publication of CN107272067A publication Critical patent/CN107272067A/en
Application granted granted Critical
Publication of CN107272067B publication Critical patent/CN107272067B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric 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

The present invention relates to a kind of optically-actives to find oil method, comprising the following steps: S1: providing an observation system, the observation system is made of host, three magnetic probes, four electrodes and GPS antenna;S2: the electromagnetic field data of acquisition time sequence, the sample frequency for observing data is 15Hz, and the continuous observation time is greater than 4 hours, 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 converted to 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 converted 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: it draws the azimuthal frequency response curve of electric field intensity and explains oiliness.The step of optically-active oil prospecting method provides can direct detection go out the position of oil-gas reservoir, reduce exploration cost.

Description

A kind of optically-active oil prospecting method
Technical field
The present invention relates to field of geophysical exploration more particularly to a kind of optically-active to find oil method.
Background technique
In petroleum exploration field, applicable geophysical method has earthquake, gravity, magnetic force and electrical method method, seismic method Based on the elastic difference of substance, 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 indirectly to study oil gas field by analysis geologic structure, and can not be straight Connect the presence for judging oil gas field.The specific method of electrical method method includes a variety of mutation, as induced polarization method, complex resistivity method and when Frequency electromagnetic method.Induced polarization method, complex resistivity method and time-frequency electromagnetic method are other than capableing 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.To sum up, up to the present, there are no a kind of methods being capable of direct detection oil-gas reservoir.
Summary of the invention
The purpose of the present invention is being directed to above-mentioned status, a kind of position that can go out oil-gas reservoir with direct detection is provided, is dropped simultaneously Method is found oil in the optically-active of low exploration cost.
A kind of the technical solution adopted by the present invention: optically-active oil prospecting method, comprising the following steps: S1: providing an observation system, The observation system is made of host, three magnetic probes, four electrodes and GPS antenna, and four electrodes are for observing two mutually Vertical electric field component Ex, Ey, three magnetic probes for 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 for observing data is 15Hz, continuous to see The time is surveyed greater than 4 hours, the electromagnetic field data includes five track datas: Ex (t), Ey (t), Hx (t), Hy (t) and Hz (t), Wherein E indicates that electric field, H indicate magnetic field, and it is northwards positive direction that x, which indicates North and South direction reference axis, and y indicates east-west direction reference axis, It is eastwards positive direction, z indicates reference axis vertically downward, and t represents the time;S3: Fast Fourier Transform (FFT) is utilized, to time-domain Electromagnetic field observation Ex (t), Ey (t), Hx (t), Hy (t) are converted to 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 converted is overlapped acquisition average value;S4: by frequency domain after being superimposed Electric field component Ex (f) and Ey (f), according to calculation formula θ (f)=(180/3.1415926) arctan [︱ Ey (f) ︱/︱ Ex (f) ︱] electric field intensity azimuth is calculated, θ (f) is the azimuth that frequency is f electric field intensity in formula;And S5: with frequency for horizontal seat Mark, electric field intensity azimuth are ordinate, draw the azimuthal frequency response curve of electric field intensity.
The effect of the technical program is: the step of optically-active oil prospecting method provides can direct detection go out the position of oil-gas reservoir It sets, reduces exploration cost.
Further, in step S1, the host is V5-2000 electromagnetism measurement apparatus.
Further, in step S2, the continuous observation time is 8 hours.
Further, it in step S3, converts the quantity for the data point chosen every time as N, if sample frequency is Fs, is superimposed it Ex (f) afterwards is a data sequence, the frequency that i-th of data (i is since 1) indicates are as follows: f=(i-1) * Fs/N;For handing over It flows for signal, the modulus value of i-th of complex data is exactly the amplitude of signal under the corresponding frequency divided by N/2, and direct current is believed For number, the modulus value of i-th of complex data is exactly the amplitude of signal under the corresponding frequency, the frequency of the 1st data point divided by N It is 0, indicates direct current signal, the frequency Fs/N of the 2nd data point, the 3rd frequency 2Fs/N, until N/2 data point.
Detailed description of the invention
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, and 4 Direction of an electric field when 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.
Electric field rotates angle-frequency curve (1 Hangzhou somewhere when Fig. 3 is no oil reservoir;2 Hebei Nanpi County somewhere).
Fig. 4 is that electric field rotates angle-frequency curve when having oil reservoir.
Specific embodiment
The principle and features of the present invention will be described below with reference to the accompanying drawings, and the given examples are served only to explain the present invention, and It is non-to be used to limit the scope of the invention.
Referring to Fig. 1, find oil the schematic diagram of method for a kind of optically-active provided by the 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 when hiding;5 be direction of an electric field when having oil reservoir.
The optically-active find oil method the following steps are included:
1. the acquisition of time series electric field data
An observation system is provided, referring to Fig. 2, the observation system was by host, three magnetic probes, four electrodes and GPS days Line composition, four electrodes are mutually perpendicular to for observing two orthogonal electric field component Ex, Ey, three magnetic probes for observing Three magnetic-field components Hx, Hy, Hz.Ex is parallel with Hx, and Ey is parallel with Hy.Specifically, which is Canadian phoenix earth object The V5-2000 electromagnetism measurement apparatus of reason company, the electromagnetic field of the V5-2000 electromagnetism measurement apparatus recordable time sequence.
Operate as follows: on 028 measuring point of the L01 survey line in some work area, the sample frequency for observing data is 15Hz, i.e., often Second records equally spaced 15 data, and the continuous observation time is greater than 4 hours, and three data file L01- can be obtained 028A.TBL, L01-028A.TSH and L01-028A.TSL.L01-028A.TBL is gauge outfit file, has recorded the one of field construction A little information.L01-028A.TSH is electromagnetic field of high frequency record data, and sample frequency 2400Hz, L01-028A.TSL are low-frequency electricals Magnetic field data, sample frequency 15Hz.
Electromagnetic field data includes five track datas: Ex (t), Ey (t), and Hx (t), Hy (t), Hz (t), wherein E indicates electric field, H Indicate magnetic field, it is northwards positive direction that x, which indicates North and South direction reference axis, and it is eastwards positive direction that y, which indicates east-west direction reference axis, z Indicate reference axis vertically downward, t represents the time.Ex (t) means that the electric field observation in moment t North and South direction, according to this class It pushes away.
In this present embodiment, the continuous observation time is 8 hours.
2. electric field time-domain is to the conversion of frequency domain
Using Fast Fourier Transform (FFT), the electromagnetic field observation Ex (t) of time-domain, Ey (t), Hx (t), Hy (t) are converted At the electromagnetic field Ex (f) of frequency domain, Ey (f), Hx (f), Hy (f).By the way of segmentation conversion, the purpose for being segmented conversion is pressure Random noise processed improves the reliability of data, and the result that multistage is converted is overlapped acquisition average value.Conversion is chosen every time The quantity of data point be N, if sample frequency is Fs, Ex (f) after being superimposed is also a data sequence, multiple comprising 1024 Number, plural angle indicate phase, the frequency that i-th of data (i is since 1) indicates are as follows: f=(i-1) * Fs/N;I-th of plural number The modulus value of data is exactly the amplitude (being divided by N for direct current signal) of signal under the corresponding frequency divided by N/2.1st data The frequency of point is 0, indicates direct current signal, the frequency Fs/N of the 2nd data point, the 3rd frequency 2Fs/N, until N/2 data Point, data do not have to thereafter.
3. the azimuthal calculating of electric field intensity
By the electric field component Ex (f) and Ey (f) in frequency domain after being superimposed, electric field intensity azimuth is calculated.Its public affairs calculated Formula is θ (f)=(180/3.1415926) arctan [︱ Ey (f) ︱/︱ Ex (f) ︱], and it is f electric field intensity that θ (f), which is frequency, in formula Azimuth.Namely with the mould of Ey (f) divided by negating tangent, the azimuth that acquisition is indicated with radian after the mould of Ex (f) again.
4. the interpretation and evaluation of Drawing of Curve and oiliness
Using frequency as abscissa, electric field intensity azimuth is ordinate, and it is bent to draw the azimuthal frequency response of electric field intensity Line, as shown in Figure 3 and Figure 4.In figure, 8Hz nearby curve to it is recessed be because of Schumann resonance, Schumann resonance is the earth and ionosphere The inherent characteristic of the resonant cavity of composition.Fig. 3 be Hangzhou somewhere and two, Hebei Nanpi County somewhere observation station electric field rotation angle with Frequency and the curve graph changed, electric field data are 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 indicate within the scope of the certain depth of underground without petroleum.Electric field when Fig. 4 is Xinjiang oil field Angle-frequency curve is rotated, from small to large, the angle of electric field tapers to 40 degree by 55 degree to frequency, and variation is due to underground Petroleum caused by, actually the observation station be above some oil field, it is theoretical to coincide with practical, illustrate that method is found oil in optically-active Validity.
In this way, optically-active oil prospecting method the step of providing can direct detection go out the position of oil-gas reservoir, reduce exploration at This.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all in spirit of the invention and Within principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.

Claims (4)

1. method is found oil in a kind of optically-active, which comprises the following steps:
S1: providing an observation system, and the observation system is made of host, three magnetic probes, four electrodes and GPS antenna, and four A electrode is for observing two orthogonal electric field component Ex, Ey, and three magnetic probes are for observing orthogonal three magnetic Field component Hx, Hy, Hz, wherein Ex is parallel with Hx, and Ey is parallel with Hy;
S2: the electromagnetic field data of acquisition time sequence, the sample frequency for observing data is 15Hz, and the continuous observation time is greater than 4 Hour, the electromagnetic field data includes five track datas: Ex (t), Ey (t), Hx (t), Hy (t) and Hz (t), and wherein E indicates electric field, H indicates magnetic field, and it is northwards positive direction that x, which indicates North and South direction reference axis, and it is eastwards positive direction that y, which indicates east-west direction reference axis, z Indicate reference axis vertically downward, t represents the time;
S3: Fast Fourier Transform (FFT) is utilized, the electromagnetic field observation Ex (t) of time-domain, Ey (t), Hx (t), Hy (t) are converted to Electromagnetic field Ex (f), Ey (f), Hx (f), the Hy (f) of frequency domain are carried out the result of multistage conversion by the way of segmentation conversion Superposition obtains average value;
S4: by the electric field component Ex (f) and Ey (f) in frequency domain after being superimposed, according to calculation formula θ (f)=(180/ 3.1415926) arctan [︱ Ey (f) ︱/︱ Ex (f) ︱] calculates electric field intensity azimuth, and θ (f) is that the electric field that frequency is f is sweared in formula Measure azimuth;
S5: using frequency as abscissa, electric field intensity azimuth is ordinate, and it is bent to draw the azimuthal frequency response of electric field intensity Line, the curve contain the information that oil-gas reservoir whether there is.
2. method is found oil in optically-active as described in claim 1, it is characterised in that: in step S1, the host is V5-2000 electromagnetism Measuring instrument.
3. method is found oil in optically-active as described 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 described in claim 1, it is characterised in that: in step S3, convert the data point chosen every time Quantity be N, if sample frequency is Fs, the Ex (f) after being superimposed is a data sequence, the frequency that i-th of data indicates are as follows: F=(i-1) * Fs/N, i is since 1;For AC signal, the modulus value of i-th of complex data is exactly corresponding be somebody's turn to do divided by N/2 The amplitude of signal under frequency, for direct current signal, the modulus value of i-th of complex data is exactly under the corresponding frequency divided by N Signal amplitude, the frequency of the 1st data point is 0, indicate direct current signal, the frequency Fs/N of the 2nd data point, the 3rd Frequency 2Fs/N, until N/2 data point.
CN201710413368.4A 2017-06-05 2017-06-05 A kind of optically-active oil prospecting method Active CN107272067B (en)

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 CN107272067A (en) 2017-10-20
CN107272067B true 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)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1239922C (en) * 2003-08-01 2006-02-01 中国石油天然气集团公司 Artificial source time frequency electro magnetic bathymetry
US8352190B2 (en) * 2009-02-20 2013-01-08 Exxonmobil Upstream Research Company Method for analyzing multiple geophysical data sets
GB201200313D0 (en) * 2012-01-10 2012-02-22 Secr Defence Signature identification
CN103809203B (en) * 2012-11-06 2017-02-22 核工业北京地质研究院 Outdoor detection working method resolving electromagnetic interference

Also Published As

Publication number Publication date
CN107272067A (en) 2017-10-20

Similar Documents

Publication Publication Date Title
US7362102B2 (en) Electromagnetic surveying for resistive or conductive bodies
US7483792B2 (en) Electromagnetic surveying for hydrocarbon reservoirs
CA2741011C (en) Method for determining electromagnetic survey sensor orientation
AU2006214069B2 (en) System and method for using time-distance characteristics in acquisition, processing and imaging of t-CSEM data
AU2013380971B2 (en) Sensors for passive electroseismic and seismoelectric surveying
Cho et al. Miniature array analysis of microtremors
BRPI0214678B1 (en) electromagnetic inspection method to inspect an area previously identified as potentially containing an underwater hydrocarbon reservoir
US10024995B2 (en) System and method for elevated source to borehole electromagnetic survey
CN103995301A (en) Method and device for evaluating total organic carbon content in shale gas reservoir
WO2010123696A2 (en) Marine source to borehole electromagnetic mapping of sub-bottom electrical resistivity
Guérin et al. Geophysical characterisation of karstic networks–Application to the Ouysse system (Poumeyssen, France)
Luzi et al. Site effects observed in the Norcia intermountain basin (Central Italy) exploiting a 20-year monitoring
WO2010129556A2 (en) Monitoring reservoirs using array based controlled source electromagnetic methods
Alumbaugh et al. A numerical sensitivity study of three dimensional imaging from a single borehole
CN107272067B (en) A kind of optically-active oil prospecting method
Zhang et al. A new algorithm for frequency-dependent shear-wave splitting parameters extraction
CN105807326B (en) The system and method that a kind of utilization sky wave carries out deep prospecting
Salem Multi-deconvolution analysis of potential field data
Karamzadeh et al. Small-aperture array as a tool to monitor fluid injection-and extraction-induced microseismicity: applications and recommendations
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
Shao‐Rui et al. Research on the Application of Borehole Seismometer Azimuth Detection Technology
Amjadi et al. An application of magnetotelluric data inversion in a stratovolcano region
Weidle et al. Identification of response and timing issues at permanent European broadband stations from automated data analysis

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