CN105137495B - A kind of method and system of oil and gas detection - Google Patents
A kind of method and system of oil and gas detection Download PDFInfo
- Publication number
- CN105137495B CN105137495B CN201510500584.3A CN201510500584A CN105137495B CN 105137495 B CN105137495 B CN 105137495B CN 201510500584 A CN201510500584 A CN 201510500584A CN 105137495 B CN105137495 B CN 105137495B
- Authority
- CN
- China
- Prior art keywords
- mrow
- curvature
- msup
- polarizability
- resistivity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
An embodiment of the present invention provides a kind of gas-oil detecting method and system, the described method includes：Obtain the resistivity curve and polarizability curve of observation station in formation at target locations；Calculate the resistivity curvature and polarizability curvature of the depth point under the observation station；Calculate the resistivity curvature of the depth point and the related coefficient of polarizability curvature；The related coefficient isopleth sectional drawing of the formation at target locations is drawn, identifies that the isopleth region that correlation coefficient value is more than first threshold in the isopleth sectional drawing is hydrocarbon-zone.Calculating and correlation analysis based on resistivity curvature in electromagnetic survey and polarizability curvature, effectively reject false anomaly, and extraction reflects whether the effective information of oily.Effective hydrocarbon-zone is identified according to correlation analysis, enhances the ability for differentiating effective anomaly, improves the accuracy of identification of hydrocarbon-zone.
Description
Technical field
The invention belongs to geophysics petroleum exploration domain, more particularly to a kind of method and system of oil and gas detection.
Background technology
In petroleum exploration domain, the either exploration of subtle pool or structural deposit, gas-oil detecting method is always
Praised highly and paid attention to by oil-gas exploration expert.In seismic prospecting, the seismic amplitude of hydrocarbon-bearing formation, frequency, speed etc. are often utilized
The exception that seismic properties occur, identifies or judges oil-gas reservoir.Than " bright spot " technology as everyone knows, which is applied to always
Identification of hydrocarbon target, but seismic properties research can also not accurately identify every time, oil gas drilling success rate does not obtain still
Significantly improve, mainly also reside in the wave impedance and filling oil gas and no filling wave impedance difference unobvious of oil gas and water.Ground
Geochemistry method is called unique directly oil prospecting method, and this method directly measures hydrocarbon component content, but due to simply on ground
Table gathers, and away from too far away, effect is often not so good as people's will for oil-gas reservoir with deep under ground.
Since rock electricity is close with Hydrocarbon Relationship, the resistivity contrasts of oil gas and water can reach 1000 times, fill oil gas
With no filling resistivity difference clearly, therefore, electrical method accounts for critical role in oil-gas recognition.As identification and evaluation oil-containing
The most basic parameter of gas reservoir, resistivity have been a concern in electric logging technology, either well logging or ground electromagnetic
The resistivity of measurement can reflect the oily situation of reservoir targets to a certain extent.
During the application is realized, inventor has found that at least there are the following problems in the prior art：Contain fire in stratum
The non-hydrocarbon formations of high resistivity as rock stratum are mixed in oil-gas Layer, if directly evaluating reservoir targets by resistivity value
Oily situation, then obtained abnormal resistivity value is likely to be to be produced by the non-hydrocarbon formations of high resistivity, this will lead
Cause the error of evaluation result.
The content of the invention
In the prior art by the method for resistivity value anomaly evaluation oil-gas reservoir, there are larger error, for above-mentioned feelings
Condition, the present invention propose a kind of method and system of oil and gas detection, and the method and system of the oil and gas detection are specifically realized in
's：
A kind of gas-oil detecting method, the described method includes：
Obtain the resistivity curve and polarizability curve of observation station in formation at target locations；
According to the resistivity curve and polarizability curve, depth point in the formation at target locations under observation station is calculated
Resistivity curvature and polarizability curvature；
Calculate the resistivity curvature of the depth point under the observation station and the related coefficient of polarizability curvature；
According to the related coefficient of the resistivity curvature of the depth point under the observation station and polarizability curvature, described in drafting
The related coefficient isopleth sectional drawing of formation at target locations；Identify that correlation coefficient value is more than first threshold in the isopleth sectional drawing
Isopleth region is hydrocarbon-zone.
Optionally, in an embodiment of the present invention, the electricity for calculating the depth point in the formation at target locations under observation station
The calculation formula of resistance rate curvature includes：
Wherein, kρFor the resistivity curvature of the depth point, ρ " is two of the resistivity curve at the depth point
Order derivative, ρ ' are first derivative of the resistivity curve at the depth point.
Optionally, in an embodiment of the present invention, the pole for calculating the depth point in the formation at target locations under observation station
The calculation formula of rate curvature includes：
Wherein, kηFor the polarizability curvature of the depth point, η " is two of the polarizability curve at the depth point
Order derivative, η ' are first derivative of the polarizability curve at the depth point.
Optionally, in an embodiment of the present invention, the resistivity curvature for calculating the depth point under the observation station with
And the calculation formula of the related coefficient of polarizability curvature includes：
Wherein, riThe phase relation of the resistivity curvature and polarizability curvature that are the observation station at i-th of depth point
Number, n are the number of the depth point of the observation station,The resistivity curvature for being the observation station at i-th of depth point,
The polarizability curvature for being the observation station at i-th of depth point.
Optionally, in an embodiment of the present invention, it is described to obtain the resistivity curve of observation station and polarization in formation at target locations
Rate curve, including：
Obtain the electromagnetic field amplitude of observation station and the observation data of electromagnetic field phase in formation at target locations；
The observation data inversion is obtained in formation at target locations using electromagnetic inversion method the resistivity curve of observation station and
Polarizability curve.
Optionally, in an embodiment of the present invention, the electromagnetic inversion method includes the Aukma method of inversion.
A kind of oil and gas detection system, including：
Acquiring unit, for obtaining the resistivity curve and polarizability curve of observation station in formation at target locations；
Curvature estimation unit, for according to the resistivity curve and polarizability curve, calculating in the formation at target locations
The resistivity curvature and polarizability curvature of depth point under observation station；
Related coefficient computing unit, the resistivity curvature and polarizability for calculating the depth point under the observation station are bent
The related coefficient of rate；
Sectional drawing drawing unit, for the resistivity curvature and polarizability curvature according to the depth point under the observation station
Related coefficient, draw the related coefficient isopleth sectional drawing of the formation at target locations, identify the correlation of the isopleth sectional drawing
The isopleth region that factor v is more than first threshold is hydrocarbon-zone.
Optionally, in an embodiment of the present invention, the curvature estimation unit includes：The resistivity curvature of the depth point
Calculation formula include：
Wherein, kρFor the resistivity curvature of the depth point, ρ " is two of the resistivity curve at the depth point
Order derivative, ρ ' are first derivative of the resistivity curve at the depth point.
Optionally, in an embodiment of the present invention, the curvature estimation unit includes：The polarizability curvature of the depth point
Calculation formula include：
Wherein, kηFor the polarizability curvature of the depth point, η " is two of the polarizability curve at the depth point
Order derivative, η ' are first derivative of the polarizability curve at the depth point.
Optionally, in an embodiment of the present invention, the related coefficient computing unit includes：Depth under the observation station
The resistivity curvature of point and the calculation formula of the related coefficient of polarizability curvature include：
Wherein, riThe phase relation of the resistivity curvature and polarizability curvature that are the observation station at i-th of depth point
Number, n are the number of the depth point of the observation station,The resistivity curvature for being the observation station at i-th of depth point,
The polarizability curvature for being the observation station at i-th of depth point.
Optionally, in an embodiment of the present invention, the acquiring unit includes：
Subelement is obtained, for obtaining the electromagnetic field amplitude of observation station in formation at target locations and the observation number of electromagnetic field phase
According to；
Inverting unit, for obtaining the electricity of observation station in formation at target locations to the observation data inversion using electromagnetic inversion method
Resistance rate curve and polarizability curve.
Optionally, in an embodiment of the present invention, the electromagnetic inversion method includes the Aukma method of inversion.
The above-mentioned technical proposal of introduction of the embodiment of the present invention has the advantages that：It is bent based on resistivity in electromagnetic survey
Calculating and correlation analysis of the rate with polarizability curvature, effectively reject false anomaly, and extraction reflects whether the effective information of oily.
Effective hydrocarbon-zone is identified according to correlation analysis, the ability for differentiating effective anomaly is enhanced, improves oil-gas bearing area
The accuracy of identification in domain.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is attached drawing needed in technology description to be briefly described, it should be apparent that, drawings in the following description are the present invention
Some embodiments, for those of ordinary skill in the art, without creative efforts, can also basis
These attached drawings obtain other attached drawings.
Fig. 1 is the flow chart for the gas-oil detecting method that one embodiment of the invention provides；
Fig. 2 is the method flow diagram of resistivity curve and polarizability curve that observation station is obtained in one embodiment of the invention；
Fig. 3 is observation station and the structure diagram of depth point in the earth formation in one embodiment of the invention；
Fig. 4 is the related coefficient isopleth sectional drawing in certain work area in one embodiment of the invention；
Fig. 5 is the structure diagram for the oil and gas detection system that one embodiment of the invention provides；
Fig. 6 is the structure diagram for the acquiring unit that one embodiment of the invention provides.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, the technical solution in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is
Part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
All other embodiments obtained without making creative work, belong to the scope of protection of the invention.
As two basic parameters of identification and evaluation oil-bearing reservoir, resistivity and polarizability all to a certain extent can be with
Reflect the oily situation of target reservoir.Resistivity method is based on the conductivity difference of rock and ore, by observing and studying
Stratum oily situation is analyzed in the distribution for the underground current field manually established.Polarizability is different from resistivity, it is not material
Intrinsic physical quantity.Polarizability method is to feed excitation current field to rock ore bed, and a series of charge movement will occur in rock ore deposit
Process, therefore polarization electromotive force is produced, stratum oily situation can be analyzed by observing and study this polarization electromotive force.Usual feelings
Under condition, in oil-bearing reservoir, resistivity and polarizability can produce a degree of exception, and in non-oil and gas reservoir, electricity
Resistance rate and polarizability will not be abnormal there is a situation where all occurring, but oil-bearing reservoir these extremely be not it is particularly evident, it is necessary to
It can just be extracted by some methods.In view of this, gas-oil detecting method proposed by the present invention combines resistivity and pole
The curvature correlation properties of rate, by analyzing both correlations, judge the oily situation of target reservoir.
The method flow diagram of gas-oil detecting method proposed by the present invention is as shown in Figure 1, specific method step is as follows：
Step S1：Obtain the resistivity curve and polarizability curve of observation station in formation at target locations.
The resistivity curve or polarizability curve of formation at target locations can generally be obtained by electromagnetic inversion method.Fig. 2 is to obtain
The resistivity curve of formation at target locations and the method flow diagram of polarizability curve are taken, can specifically include following two steps：
Step S101：Obtain the observation data of electromagnetic field amplitude and electromagnetic field phase in formation at target locations.
, it is necessary to obtain the electromagnetic field amplitude and electromagnetic field phase of observation station described in formation at target locations before electromagnetic inversion is carried out
The observation data of position.To obtain resistivity observation data instance, formation at target locations the earth can be specifically hit to using two supply terminals
Power supply, so as to form underground artificial electric field.The survey line of Resistivity testing instrument is passed through underground afterwards, two electrodes is measured respectively and leads to
Enter the potential difference of each depth point in underground, calculate and record the electromagnetic field amplitude and electromagnetic field phase value of each Frequency point on survey line.
Choose the multiple observation stations of formation at target locations to measure, record, form the electromagnetic field amplitude and electromagnetic field phase of formation at target locations
Observe data.The observation station can be along the distance between a line of observation, number and observation station of observation station according to reality
Border demand is chosen.
Step S102：The resistivity of observation station in formation at target locations is obtained to the observation data inversion using electromagnetic inversion method
Curve and polarizability curve.
The observation data inversion is obtained in formation at target locations using electromagnetic inversion method the resistivity curve of observation station and
Polarizability curve.The electromagnetic inversion method can include one-dimensional, the two-dimentional and 3 D electromagnetic method of inversion.The one-dimensional electromagnetic inversion
Method can include the Bostick methods of inversion, the gradient method of inversion, the Gauss-Newton method of inversion, least square method and simulated annealing.
The two-dimensional electromagnetic method of inversion can include Aukma (hereinafter referred to as OCCAM) method of inversion, the fast relaxation method of inversion (RRI methods) with
And Nonlinear Conjugate Gradient Methods (NLCG methods).The one-dimensional electromagnetic inversion method can include OCCAM methods, the fast relaxation method of inversion
(RRI methods), Nonlinear Conjugate Gradient Methods (NLCG methods) and the artificial neural network method of inversion (ANN methods).Due to the method for the present invention
Resistivity-polarizability related coefficient isopleth sectional drawing of depth-dimension of survey line length two can be subsequently obtained, therefore at this
The two-dimensional electromagnetic method of inversion can be used in inventive embodiments.
The OCCAM methods of inversion can be fitted all observation data, and computing is stablized, and step is simple, fast convergence rate, because
This embodiment of the present invention carries out inverting using the OCCAM methods of inversion to the resistivity observation data and polarizability observation data.Institute
State the OCCAM methods of inversion use Lagrange multiplier come balance model is smooth and the object function of data fitting degree for：
Wherein, μ-1For Lagrange multiplier, d observes data vector for resistivity or polarizability, and F is forward operator, | | Rm |
| it is the roughness of model, R is roughness matrix, | | Wd-WF [m] | |2For the second order norm of standard, observation data d and mould are represented
The fitting difference X of type forward response F [m]2,For X2Desired value, W be n × n data normalization diagonal matrix.
Inverting target is to obtain smooth enough resistivity curve or polarizability curve a model, can specifically be led to
Multiple interative computation is crossed to realize, the iteration expression formula of kth time is：
mk+1=[μ RTR+(WJk)TWJk]-1(WJk)TWdk (2)
In an iteration computing, first have to calculate the partial derivative matrix J of this iterationk, then give a series of μ
Value, obtains iterative model m (μ), using the forward response F of Finite Element computation model by Cholesky breakdowns (2)
[m], it is poor using the opposite fitting of formula (3) computation model, according to the optimal μ values of the definite current iteration of fitting difference and model.
Step S2：According to the resistivity curve and polarizability curve of the observation station, calculate and seen in the formation at target locations
The resistivity curvature and polarizability curvature of depth point under measuring point.
Step S102 invertings obtain the resistivity curve and polarizability curve of each observation station, according to the observation station
Resistivity curve and polarizability curve, calculate resistivity curvature and the pole of depth point in the formation at target locations under observation station
Rate curvature.Fig. 3 is the structure diagram of observation station and depth point in the earth formation, as shown in the figure, observation station A, B along
Line of observation direction, the line of observation can be located at ground.The depth point is the Inversion Calculation point under the observation station along depth direction,
As shown in figure 3, point A1, A2, A3 and A4 are the Inversion Calculation depth points of observation station A, point B1, B2, B3 and B4 are with observation station B
Inversion Calculation depth point.The calculation formula of the resistivity curvature of the depth point is：
Wherein, kρFor the resistivity curvature of the depth point, ρ " is two of the resistivity curve at the depth point
Order derivative, ρ ' are first derivative of the resistivity curve at the depth point.
The calculation formula of the polarizability curvature of the depth point is：
Wherein, kηFor the polarizability curvature of the depth point, η " is two of the polarizability curve at the depth point
Order derivative, η ' are first derivative of the polarizability curve at the depth point.
In step s 2, reflect whether to contain by the calculating of resistivity curvature in electromagnetic survey and polarizability curvature, extraction
The effective information of oil gas
Step S3：Calculate the resistivity curvature of the depth point under the observation station and the related coefficient of polarizability curvature.
There are corresponding resistivity curvature and polarizability curvature in each depth point under formation at target locations observation station.
Calculate the resistivity curvature of the depth point under the observation station and the related coefficient of polarizability curvature.Depth under the observation station
Spending the calculation formula of the resistivity curvature of point and the related coefficient of polarizability curvature includes：
Wherein, riThe phase relation of the resistivity curvature and polarizability curvature that are the observation station at i-th of depth point
Number, n are the number of the depth point of the observation station,The resistivity curvature for being the observation station at i-th of depth point,
The polarizability curvature for being the observation station at i-th of depth point.
Step S3 is analyzed both correlations, is had based on the electromagnetism rate curvature that step S2 is calculated and polarizability curvature
Effect rejects false anomaly, enhances the ability for differentiating effective anomaly, improves the accuracy of identification of hydrocarbon-zone.
Step S4：According to the related coefficient of the resistivity curvature of the depth point under the observation station and polarizability curvature,
The related coefficient isopleth sectional drawing of the objective area is drawn, when the related coefficient isopleth numerical value is more than first threshold
When, the isopleth region is hydrocarbon-zone.
The related coefficient isopleth sectional drawing is along the multiple depth point resistivity curvature of depth direction according to each observation station
And the result of the related coefficient of polarizability curvature is drawn and is formed.The related coefficient isopleth sectional drawing can reflect along survey line
Direction (direction of the observation station in straight line) and the related coefficient of depth direction resistivity curvature and polarizability curvature
Change.The abscissa of the related coefficient isopleth sectional drawing can be survey line length, and ordinate can be depth.Ordinary circumstance
Under, hydrocarbon-zone has the characteristics that resistivity curvature anomalies are high and polarizability curvature anomalies are high, therefore both are abnormal at the same time high
When, both related coefficients are also of a relatively high.Resistivity curvature and polarizability curvature without hydrocarbon-zone will not go out at the same time
It is now abnormal, therefore both related coefficients are smaller.In conclusion the related coefficient and mesh of resistivity curvature and polarizability curvature
The probability of mark stratum oily is directly proportional, and related coefficient is bigger, and the probability of formation at target locations oily is bigger.Described related coefficient etc.
When value line number value is more than first threshold, the isopleth region is hydrocarbon-zone.The first threshold is according to known prospect pit
Data is configured.
Fig. 4 is the related coefficient isopleth sectional drawing in certain work area, in a length of 6 km of survey line, depth in the stratum of 5 kms,
The distribution of related coefficient isopleth is as shown in the figure.May determine that depth is 3 kms, and it is oil gas storage at 3.5 kms that survey line length, which is,
Layer center.According to the data of known prospect pit, it is 0.2 to set first threshold, then isopleth region of the related coefficient more than 0.2 is
Hydrocarbon-zone.Hydrocarbon-zone in Fig. 4 is 0.2 area encompassed of isopleth.According to above-mentioned hydrocarbon-zone really
It is fixed, it can be determined that W1 is dry-well, and W2 is oil-producing well.
It is the oil and gas detection system to be described below with the corresponding oil and gas detection system of above-described embodiment method, Fig. 5
Structure diagram, as shown in the figure, the system 50 includes：Acquiring unit 51, curvature estimation unit 52, related coefficient computing unit
53, sectional drawing drawing unit 54, wherein,
Acquiring unit 51, for obtaining the resistivity curve and polarizability curve of observation station in formation at target locations.
Curvature estimation unit 52, for according to the resistivity curve and polarizability curve, calculating the formation at target locations
The resistivity curvature and polarizability curvature of depth point under middle observation station.
Related coefficient computing unit 53, for calculating the resistivity curvature and polarizability of the depth point under the observation station
The related coefficient of curvature.
Sectional drawing drawing unit 54, it is bent for the resistivity curvature and polarizability according to the depth point under the observation station
The related coefficient of rate, draws the related coefficient isopleth sectional drawing of the objective area, when the related coefficient isopleth numerical value
During more than first threshold, the isopleth region is hydrocarbon-zone.
Fig. 6 is the structure diagram of the acquiring unit 51, as shown in the figure, the acquiring unit 51 further includes：Obtain son
Unit 61, inverting unit 62, wherein,
Subelement 61 is obtained, for obtaining electromagnetic field amplitude and the observation of electromagnetic field phase of observation station in formation at target locations
Data.
Inverting unit 62, for obtaining observation station in formation at target locations to the observation data inversion using electromagnetic inversion method
Resistivity curve and polarizability curve.
The above-mentioned technical proposal of introduction of the embodiment of the present invention has the advantages that：It is bent based on resistivity in electromagnetic survey
Calculating and correlation analysis of the rate with polarizability curvature, effectively reject false anomaly, and extraction reflects whether the effective information of oily.
Effective hydrocarbon-zone is identified according to correlation analysis, the ability for differentiating effective anomaly is enhanced, improves oil-gas bearing area
The accuracy of identification in domain.
Those skilled in the art will also be appreciated that the various illustrative components, blocks that the embodiment of the present invention is listed
(illustrative logical block), unit, and step can pass through the knot of electronic hardware, computer software, or both
Conjunction is realized.To clearly show that the replaceability of hardware and software (interchangeability), above-mentioned various explanations
Property component (illustrative components), unit and step universally describe their function.Such work(
Can be that specific application and the design requirement of whole system are depended on to realize by hardware or software.Those skilled in the art
Various methods can be used to realize the function, but this realization is understood not to for every kind of specific application
Beyond the scope of protection of the embodiment of the present invention.
Various illustrative logical blocks described in the embodiment of the present invention, or unit can by general processor,
Digital signal processor, application-specific integrated circuit (ASIC), field programmable gate array (FPGA) or other programmable logic devices,
Described function is realized or operated in discrete gate or transistor logic, discrete hardware components, or the design of any of the above described combination.
General processor can be microprocessor, alternatively, the general processor can also be any traditional processor, controller,
Microcontroller or state machine.Processor can also be realized by the combination of computing device, such as digital signal processor and micro-
Processor, multi-microprocessor, one or more microprocessors combine a digital signal processor core, or any other like
Configuration realize.
The step of method or algorithm described in the embodiment of the present invention can be directly embedded into hardware, processor perform it is soft
Part module or the combination of both.Software module can be stored in RAM memory, flash memory, ROM memory, EPROM storages
Other any form of storaging mediums in device, eeprom memory, register, hard disk, moveable magnetic disc, CD-ROM or this area
In.Exemplarily, storaging medium can be connected with processor, so that processor can read information from storaging medium, and
Can be to storaging medium stored and written information.Alternatively, storaging medium can also be integrated into processor.Processor and storaging medium can
To be arranged in ASIC, ASIC can be arranged in user terminal.Alternatively, processor and storaging medium can also be arranged at use
In different components in the terminal of family.
In one or more exemplary designs, the described above-mentioned function of the embodiment of the present invention can be in hardware, soft
Part, firmware or any combination of this three are realized.If realized in software, these functions can store and computer-readable
On medium, or with one or more instruction or code form be transmitted on the medium of computer-readable.Computer readable medium includes electricity
Brain storaging medium and easy to allow computer program to be transferred to other local telecommunication medias from a place.Storaging medium can be with
It is that any general or special computer can be with the useable medium of access.For example, such computer readable media can include but
It is not limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage devices, or other
What can be used for carrying or store with instruct or data structure and it is other can be by general or special computer or general or specially treated
The medium of the program code of device reading form.In addition, any connection can be properly termed computer readable medium, example
Such as, if software is to pass through a coaxial cable, fiber optic cables, double from a web-site, server or other remote resources
Twisted wire, Digital Subscriber Line (DSL) or with defined in being also contained in of the wireless way for transmitting such as infrared, wireless and microwave
In computer readable medium.The disk (disk) and disk (disc) include Zip disk, radium-shine disk, CD, DVD, floppy disk
And Blu-ray Disc, disk is usually with magnetic duplication data, and disk usually carries out optical reproduction data with laser.Combinations of the above
It can also be included in computer readable medium.
The above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications also should
It is considered as protection scope of the present invention.
Claims (12)
- A kind of 1. gas-oil detecting method, it is characterised in that the described method includes：Obtain the resistivity curve and polarizability curve of observation station in formation at target locations；According to the resistivity curve and polarizability curve, the resistance of the depth point in the formation at target locations under observation station is calculated Rate curvature and polarizability curvature；Calculate the resistivity curvature of the depth point under the observation station and the related coefficient of polarizability curvature；According to the related coefficient of the resistivity curvature of the depth point under the observation station and polarizability curvature, the target is drawn The related coefficient isopleth sectional drawing on stratum；Identify that correlation coefficient value is more than the equivalence of first threshold in the isopleth sectional drawing Line region is hydrocarbon-zone.
- 2. gas-oil detecting method according to claim 1, it is characterised in that described to calculate observation station in the formation at target locations Under the calculation formula of resistivity curvature of depth point include：<mrow> <msub> <mi>k</mi> <mi>&rho;</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <msup> <mi>&rho;</mi> <mrow> <mo>&prime;</mo> <mo>&prime;</mo> </mrow> </msup> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mi>&rho;</mi> <mrow> <mo>&prime;</mo> <mn>2</mn> </mrow> </msup> <mo>)</mo> </mrow> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> </msup> </mfrac> <mo>)</mo> </mrow> </mrow>Wherein, kρFor the resistivity curvature of the depth point, ρ " leads for second order of the resistivity curve at the depth point Number, ρ ' is first derivative of the resistivity curve at the depth point.
- 3. gas-oil detecting method according to claim 1, it is characterised in that described to calculate observation station in the formation at target locations Under the calculation formula of polarizability curvature of depth point include：<mrow> <msub> <mi>k</mi> <mi>&eta;</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <msup> <mi>&eta;</mi> <mrow> <mo>&prime;</mo> <mo>&prime;</mo> </mrow> </msup> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mi>&eta;</mi> <mrow> <mo>&prime;</mo> <mn>2</mn> </mrow> </msup> <mo>)</mo> </mrow> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> </msup> </mfrac> <mo>)</mo> </mrow> </mrow>Wherein, kηFor the polarizability curvature of the depth point, η " leads for second order of the polarizability curve at the depth point Number, η ' is first derivative of the polarizability curve at the depth point.
- 4. gas-oil detecting method according to claim 1, it is characterised in that the depth point calculated under the observation station Resistivity curvature and the calculation formula of related coefficient of polarizability curvature include：<mrow> <msup> <mi>r</mi> <mi>i</mi> </msup> <mo>=</mo> <mfrac> <mrow> <mi>n</mi> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&rho;</mi> <mi>i</mi> </msubsup> <msubsup> <mi>k</mi> <mi>&eta;</mi> <mi>i</mi> </msubsup> <mo>-</mo> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&rho;</mi> <mi>i</mi> </msubsup> <mo>&CenterDot;</mo> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&eta;</mi> <mi>i</mi> </msubsup> </mrow> <mrow> <msqrt> <mrow> <mi>n</mi> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>k</mi> <mi>&rho;</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&rho;</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>&CenterDot;</mo> <msqrt> <mrow> <mi>n</mi> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>k</mi> <mi>&eta;</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&eta;</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> </mfrac> </mrow>Wherein, riThe related coefficient of the resistivity curvature for being the observation station at i-th of depth point and polarizability curvature, n are The number of the depth point of the observation station,The resistivity curvature for being the observation station at i-th of depth point,To be described Polarizability curvature of the observation station at i-th of depth point.
- 5. gas-oil detecting method according to claim 1, it is characterised in that the electricity for obtaining observation station in formation at target locations Resistance rate curve and polarizability curve, including：Obtain the electromagnetic field amplitude of observation station and the observation data of electromagnetic field phase in formation at target locations；The resistivity curve of observation station and polarization in formation at target locations are obtained to the observation data inversion using electromagnetic inversion method Rate curve.
- 6. gas-oil detecting method according to claim 5, it is characterised in that the electromagnetic inversion method includes Aukma inverting Method.
- A kind of 7. oil and gas detection system, it is characterised in that including：Acquiring unit, for obtaining the resistivity curve and polarizability curve of observation station in formation at target locations；Curvature estimation unit, for according to the resistivity curve and polarizability curve, calculating and being observed in the formation at target locations The resistivity curvature and polarizability curvature of depth point under point；Related coefficient computing unit, for calculating the resistivity curvature and polarizability curvature of the depth point under the observation station Related coefficient；Sectional drawing drawing unit, for the resistivity curvature and the phase of polarizability curvature according to the depth point under the observation station Relation number, draws the related coefficient isopleth sectional drawing of the formation at target locations, identifies the related coefficient of the isopleth sectional drawing The isopleth region that numerical value is more than first threshold is hydrocarbon-zone.
- 8. oil and gas detection system according to claim 7, it is characterised in that the curvature estimation unit includes：The depth Spending the calculation formula of the resistivity curvature of point includes：<mrow> <msub> <mi>k</mi> <mi>&rho;</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <msup> <mi>&rho;</mi> <mrow> <mo>&prime;</mo> <mo>&prime;</mo> </mrow> </msup> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mi>&rho;</mi> <mrow> <mo>&prime;</mo> <mn>2</mn> </mrow> </msup> <mo>)</mo> </mrow> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> </msup> </mfrac> <mo>)</mo> </mrow> </mrow>Wherein, kρFor the resistivity curvature of the depth point, ρ " leads for second order of the resistivity curve at the depth point Number, ρ ' is first derivative of the resistivity curve at the depth point.
- 9. oil and gas detection system according to claim 7, it is characterised in that the curvature estimation unit includes：The depth Spending the calculation formula of the polarizability curvature of point includes：<mrow> <msub> <mi>k</mi> <mi>&eta;</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <msup> <mi>&eta;</mi> <mrow> <mo>&prime;</mo> <mo>&prime;</mo> </mrow> </msup> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mi>&eta;</mi> <mrow> <mo>&prime;</mo> <mn>2</mn> </mrow> </msup> <mo>)</mo> </mrow> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> </msup> </mfrac> <mo>)</mo> </mrow> </mrow>Wherein, kηFor the polarizability curvature of the depth point, η " leads for second order of the polarizability curve at the depth point Number, η ' is first derivative of the polarizability curve at the depth point.
- 10. oil and gas detection system according to claim 7, it is characterised in that the related coefficient computing unit includes：Institute Stating the calculation formula of the resistivity curvature of the depth point under observation station and the related coefficient of polarizability curvature includes：<mrow> <msup> <mi>r</mi> <mi>i</mi> </msup> <mo>=</mo> <mfrac> <mrow> <mi>n</mi> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&rho;</mi> <mi>i</mi> </msubsup> <msubsup> <mi>k</mi> <mi>&eta;</mi> <mi>i</mi> </msubsup> <mo>-</mo> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&rho;</mi> <mi>i</mi> </msubsup> <mo>&CenterDot;</mo> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&eta;</mi> <mi>i</mi> </msubsup> </mrow> <mrow> <msqrt> <mrow> <mi>n</mi> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>k</mi> <mi>&rho;</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&rho;</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>&CenterDot;</mo> <msqrt> <mrow> <mi>n</mi> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>k</mi> <mi>&eta;</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mi>i</mi> <mi>n</mi> </munderover> <msubsup> <mi>k</mi> <mi>&eta;</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> </mfrac> </mrow>Wherein, riThe related coefficient of the resistivity curvature for being the observation station at i-th of depth point and polarizability curvature, n are The number of the depth point of the observation station,The resistivity curvature for being the observation station at i-th of depth point,To be described Polarizability curvature of the observation station at i-th of depth point.
- 11. oil and gas detection system according to claim 7, it is characterised in that the acquiring unit includes：Subelement is obtained, for obtaining the electromagnetic field amplitude of observation station in formation at target locations and the observation data of electromagnetic field phase；Inverting unit, for obtaining the resistivity of observation station in formation at target locations to the observation data inversion using electromagnetic inversion method Curve and polarizability curve.
- 12. oil and gas detection system according to claim 11, it is characterised in that it is anti-that the electromagnetic inversion method includes Aukma Drill method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510500584.3A CN105137495B (en) | 2015-08-14 | 2015-08-14 | A kind of method and system of oil and gas detection |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510500584.3A CN105137495B (en) | 2015-08-14 | 2015-08-14 | A kind of method and system of oil and gas detection |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105137495A CN105137495A (en) | 2015-12-09 |
| CN105137495B true CN105137495B (en) | 2018-05-04 |
Family
ID=54722901
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510500584.3A Active CN105137495B (en) | 2015-08-14 | 2015-08-14 | A kind of method and system of oil and gas detection |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105137495B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108828680B (en) * | 2018-07-25 | 2020-02-14 | 中国石油天然气集团有限公司 | Method and device for determining formation resistivity and polarizability |
| CN108828681B (en) * | 2018-08-01 | 2020-01-07 | 中国石油天然气集团有限公司 | Method and device for determining formation resistivity and polarizability |
| CN112083505B (en) * | 2019-12-04 | 2024-02-23 | 安徽省勘查技术院(安徽省地质矿产勘查局能源勘查中心) | Comprehensive physical logging method and system based on laser emission system |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4329647A (en) * | 1981-06-04 | 1982-05-11 | Petroleum Physics Corporation | Method for determining distance and direction from an open well to a cased well using resistivity and directional survey data |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8237446B2 (en) * | 2008-08-27 | 2012-08-07 | Baker Hughes Incorporated | High resolution resistivity earth imager |
| US8812237B2 (en) * | 2009-02-05 | 2014-08-19 | Schlumberger Technology Corporation | Deep-reading electromagnetic data acquisition method |
| CN102053281B (en) * | 2009-11-10 | 2013-01-23 | 中国石油化工集团公司 | Method for realizing oil-gas detection by applying long-offset distance transient electromagnetic array method |
| CN102012525B (en) * | 2010-12-21 | 2013-06-12 | 吉林大学 | Distribution type multi-parameter deep electromagnetic section imaging system and measuring method thereof |
| CN102565866A (en) * | 2012-02-08 | 2012-07-11 | 蔡运胜 | Geophysical prospecting two-dimensional electrical sounding data chromatography inversion processing technology |
| CN104614779B (en) * | 2015-01-16 | 2017-07-21 | 中国矿业大学 | A kind of multi-parameter electromagnetic method dynamic monitoring system and its method |
-
2015
- 2015-08-14 CN CN201510500584.3A patent/CN105137495B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4329647A (en) * | 1981-06-04 | 1982-05-11 | Petroleum Physics Corporation | Method for determining distance and direction from an open well to a cased well using resistivity and directional survey data |
Non-Patent Citations (2)
| Title |
|---|
| Oil detection within the Middle Jurassic sediments in the south of Western Siberia;Yuriy Tcimbaluk etl.;《SEG Denver 2014 Annual Meeting》;20141231;第1619-1623页 * |
| 大地电磁激电效应油气检测试验;罗卫锋 等;《石油地球物理勘探》;20111231;第46卷(第6期);第978-983页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105137495A (en) | 2015-12-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA3140364C (en) | Systems and methods for electromagnetic waterfront surveillance in a vicinity of an oil well | |
| He et al. | Transition probability‐based stochastic geological modeling using airborne geophysical data and borehole data | |
| CN105190366B (en) | The continuous capillary pressure curve of surface lower stratum is determined using saturation degree and nuclear magnetic resonance log data | |
| Hansen et al. | Accounting for imperfect forward modeling in geophysical inverse problems—exemplified for crosshole tomography | |
| CN104863574B (en) | A kind of Fluid Identification Method suitable for tight sandstone reservoir | |
| Morton et al. | Integrated interpretation for pressure transient tests in discretely fractured reservoirs | |
| CN105137495B (en) | A kind of method and system of oil and gas detection | |
| US8554482B2 (en) | Monitoring reservoirs using array based controlled source electromagnetic methods | |
| Li et al. | Electrical imaging of hydraulic fracturing fluid using steel-cased wells and a deep-learning method | |
| CN107861162B (en) | Microelectrode logging data-based natural crack identification method and system | |
| Nam et al. | Assessment of Delaware and Groningen effects on dual-laterolog measurements with a self-adaptive hp finite-element method | |
| Yang et al. | 1D subsurface electromagnetic fields excited by energized steel casing | |
| CN109667576A (en) | A kind of high salinity origin cause of formation low-resistivity reservoir Logging Identification Method | |
| CN113775326A (en) | Method and device for evaluating movable water saturation, electronic device, and medium | |
| Buddo et al. | Joint interpretation of seismic and TEM data from the Kovykta gas-condensate field, East Siberia | |
| MacLennan | Monitoring CO2 storage using well-casing source electromagnetics | |
| CN104749647A (en) | Geological interface structure detection method, device and system | |
| Dutta et al. | Modeling tools for drilling, reservoir navigation, and formation evaluation | |
| CN109598049B (en) | Method for drilling rock fracture development degree and regional rock fracture development rule | |
| Gong et al. | Estimating net sand from borehole images in laminated deepwater reservoirs with a neural network | |
| Smith et al. | Case histories illustrating the characteristics of the HeliGEOTEM system | |
| Gohari et al. | DAS/DTS/DSS/DPS/DxS-do we measure what adds value? | |
| Gong et al. | Improving the estimation of fracture length using resistivity and borehole image logs | |
| Sun et al. | 3D spatial characterization of sand body for uranium reservoirs based on geostatistical resistivity inversion | |
| Lun et al. | A procedure for integrating geologic concepts into history matching |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |