CN105137495A - Oil gas detection method and oil gas detection system - Google Patents

Abstract

An embodiment of the invention provides an oil gas detection method and an oil gas detection system. The oil gas detection method comprises the steps of acquiring a resistivity curve and a polarizability curve of an observation point in a target stratum; calculating the resistivity curve and the polarizability curve of a deep point below the observation point; calculating a correlation coefficient between the resistivity curve and the polarizability curve of the deep point; drawing a correlation coefficient isopleth section diaphragm of the target stratum, and identifying an isopleth area of which the correlation coefficient is larger than a first threshold in the isopleth section diaphragm. Based on calculation of resistivity curvature and polarizability curvature and correlation analysis in electromagnetic exploration, a false anomaly is effectively eliminated, and effective information which reflects whether the oil gas is contained is extracted. An effective oil gas containing area is identified according to correlation coefficient analysis. An ability for discriminating an effective abnormity and an oil gas containing area identification precision are improved.

Description

A kind of method and system of oil and gas detection

Technical field

The invention belongs to geophysics petroleum exploration domain, particularly a kind of method and system of oil and gas detection.

Background technology

In petroleum exploration domain, be no matter the exploration of subtle reservoir or structural deposit, gas-oil detecting method always praise highly by oil-gas exploration expert and pay attention to.In seismic prospecting, often utilize the exception that the seismic properties such as the seismic amplitude of hydrocarbon-bearing formation, frequency, speed occur, identify or judge hydrocarbon-bearing pool.Than " bright spot " technology as everyone knows, this technology is applied to identification of hydrocarbon target always, but seismic properties research also can accurately not identify at every turn, oil gas drilling success ratio is not still improved significantly, and is mainly also wave impedance and the filling oil gas of oil gas and water and does not have filling wave impedance difference not obvious.Geochemical method is called unique direct oil prospecting method, and the method directly measures hydrocarbon component content, but owing to just gathering on earth's surface, and with the hydrocarbon-bearing pool of deep under ground apart from too far away, effect is often not as people's will.

Because rock electricity and Hydrocarbon Relationship are close, oil gas and resistivity of water difference can reach 1000 times, filling oil gas and do not have filling resistivity difference clearly, and therefore, electrical method accounts for critical role in oil-gas recognition.As the most basic parameter of identification and evaluation oil-bearing reservoir, resistivity receives much concern in electric logging technology always, is no matter the oily situation that resistivity that well logging or ground electromagnetic are measured can reflect reservoir targets to a certain extent.

Realizing in the application's process, inventor finds that in prior art, at least there are the following problems: be mixed in hydrocarbon zone containing the non-hydrocarbon formations of high resistivity that fiery rock stratum is such in stratum, if directly evaluated the oily situation of reservoir targets by resistivity value, the abnormal resistivity value so obtained probably is produced by the non-hydrocarbon formations of high resistivity, and this will cause the error of evaluation result.

By the method for resistivity value anomaly evaluation hydrocarbon-bearing pool in prior art, there is larger error, for above-mentioned situation, the present invention proposes a kind of method and system of oil and gas detection, and the method and system of described oil and gas detection are specifically achieved in that

A kind of gas-oil detecting method, described method comprises:

Obtain resistivity curve and the polarizability curve of observation station in formation at target locations;

According to described resistivity curve and polarizability curve, calculate resistivity curvature and the polarizability curvature of the depth point in described formation at target locations under observation station;

Calculate the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature;

According to the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature, draw the related coefficient isoline sectional drawing of described formation at target locations; Identify that the isoline region that correlation coefficient value in described isoline sectional drawing is greater than first threshold is hydrocarbon-zone.

Optionally, in an embodiment of the present invention, the computing formula of the resistivity curvature of the depth point in the described formation at target locations of described calculating under observation station comprises:

$k_{\mathrm{\ρ}}=\left(\frac{{\mathrm{\ρ}}^{\′\′}}{{(1+{\mathrm{\ρ}}^{\′2})}^{\frac{3}{2}}}\right)$

Wherein, k _ρfor the resistivity curvature of described depth point, ρ " for described resistivity curve is in the second derivative at described depth point place, ρ ' is for described resistivity curve is in the first order derivative of described sight depth.

Optionally, in an embodiment of the present invention, the computing formula of the polarizability curvature of the depth point in the described formation at target locations of described calculating under observation station comprises:

$k_{\mathrm{\η}}=\left(\frac{{\mathrm{\η}}^{\′\′}}{{(1+{\mathrm{\η}}^{\′2})}^{\frac{3}{2}}}\right)$

Wherein, k _ηfor the polarizability curvature of described depth point, η " for described polarizability curve is in the second derivative of described sight depth, η ' is for described polarizability curve is in the first order derivative at described depth point place.

Optionally, in an embodiment of the present invention, the resistivity curvature of depth point under the described observation station of described calculating and the computing formula of the related coefficient of polarizability curvature comprise:

$r^i=\frac{n\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i{k}_{\mathrm{\η}}^i-\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\·\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i}{\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\ρ}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\right)}^2}\·\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\η}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i\right)}^2}}$

Wherein, r ⁱfor the resistivity curvature of described observation station at i-th depth point place and the related coefficient of polarizability curvature, n is the number of the depth point of described observation station, for described observation station is in the resistivity curvature at i-th depth point place, for described observation station is in the polarizability curvature at i-th depth point place.

Optionally, in an embodiment of the present invention, the resistivity curve of observation station and polarizability curve in described acquisition formation at target locations, comprising:

Obtain the electromagnetic field amplitude of observation station and the observation data of electromagnetic field phase in formation at target locations;

Electromagnetic inversion method is adopted described observation data inverting to be obtained to resistivity curve and the polarizability curve of observation station in formation at target locations.

Optionally, in an embodiment of the present invention, described electromagnetic inversion method comprises the Aukma method of inversion.

A kind of oil and gas detection system, comprising:

Acquiring unit, for obtaining resistivity curve and the polarizability curve of observation station in formation at target locations;

Curvature estimation unit, for according to described resistivity curve and polarizability curve, calculates resistivity curvature and the polarizability curvature of the depth point in described formation at target locations under observation station;

Calculation of correlation factor unit, for the related coefficient of the resistivity curvature and polarizability curvature that calculate the depth point under described observation station;

Sectional drawing drawing unit, for according to the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature, draw the related coefficient isoline sectional drawing of described formation at target locations, identify that the isoline region that the sectional related coefficient numerical value of described isoline is greater than first threshold is hydrocarbon-zone.

Optionally, in an embodiment of the present invention, described curvature estimation unit comprises: the computing formula of the resistivity curvature of described depth point comprises:

$k_{\mathrm{\ρ}}=\left(\frac{{\mathrm{\ρ}}^{\′\′}}{{(1+{\mathrm{\ρ}}^{\′2})}^{\frac{3}{2}}}\right)$

Wherein, k _ρfor the resistivity curvature of described depth point, ρ " for described resistivity curve is in the second derivative at described depth point place, ρ ' is for described resistivity curve is in the first order derivative of described sight depth.

Optionally, in an embodiment of the present invention, described curvature estimation unit comprises: the computing formula of the polarizability curvature of described depth point comprises:

$k_{\mathrm{\η}}=\left(\frac{{\mathrm{\η}}^{\′\′}}{{(1+{\mathrm{\η}}^{\′2})}^{\frac{3}{2}}}\right)$

Wherein, k _ηfor the polarizability curvature of described depth point, η " for described polarizability curve is in the second derivative of described sight depth, η ' is for described polarizability curve is in the first order derivative at described depth point place.

Optionally, in an embodiment of the present invention, described Calculation of correlation factor unit comprises: the computing formula of the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature comprises:

$r^i=\frac{n\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i{k}_{\mathrm{\η}}^i-\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\·\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i}{\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\ρ}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\right)}^2}\·\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\η}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i\right)}^2}}$

Wherein, r ⁱfor the resistivity curvature of described observation station at i-th depth point place and the related coefficient of polarizability curvature, n is the number of the depth point of described observation station, for described observation station is in the resistivity curvature at i-th depth point place, for described observation station is in the polarizability curvature at i-th depth point place.

Optionally, in an embodiment of the present invention, described acquiring unit comprises:

Obtain subelement, for the observation data of the electromagnetic field amplitude and electromagnetic field phase that obtain observation station in formation at target locations;

Inverting unit, obtains resistivity curve and the polarizability curve of observation station in formation at target locations for adopting electromagnetic inversion method to described observation data inverting.

Optionally, in an embodiment of the present invention, described electromagnetic inversion method comprises the Aukma method of inversion.

The technique scheme of embodiment of the present invention introduction has following beneficial effect: based on calculating and the correlation analysis of resistivity curvature and polarizability curvature in electromagnetic survey, effectively reject spurious anomaly, extracts the effective information reflecting whether oily.According to the effective hydrocarbon-zone of correlation analysis identification, enhance the ability differentiating effective anomaly, improve the accuracy of identification of hydrocarbon-zone.

Summary of the invention

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

The process flow diagram of the gas-oil detecting method that Fig. 1 provides for one embodiment of the invention;

Fig. 2 obtains the resistivity curve of observation station and the method flow diagram of polarizability curve in one embodiment of the invention;

Fig. 3 is observation station and depth point structural representation in the earth formation in one embodiment of the invention;

Fig. 4 is the related coefficient isoline sectional drawing in certain work area in one embodiment of the invention;

The structural representation of the oil and gas detection system that Fig. 5 provides for one embodiment of the invention;

The structural representation of the acquiring unit that Fig. 6 provides for one embodiment of the invention.

Embodiment

For making the object of the embodiment of the present invention, technical scheme and advantage clearly, below in conjunction with the accompanying drawing in the embodiment of the present invention, technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.

As two basic parameters of identification and evaluation oil-bearing reservoir, resistivity and polarizability can reflect the oily situation of target reservoir all to a certain extent.Resistivity method is based on the conductivity difference of rock and ore, analyzes stratum oily situation by the distribution observing and study the artificial underground current field set up.Polarizability is different from resistivity, and it is not the intrinsic physical quantity of material.Polarizability method feeds excitation current field to rock ore bed, and a series of charge movement process will occur in rock ore deposit, therefore produces polarization electromotive force, observe and study this polarization electromotive force and can analyze stratum oily situation.Under normal circumstances, in oil-bearing reservoir, resistivity and polarizability all can produce exception to a certain degree, and in non-oil and gas reservoir, resistivity and polarizability can not occur all to occur abnormal situation, but these of oil-bearing reservoir are not obvious especially extremely, need just can be extracted by some method.In view of this, the gas-oil detecting method that the present invention proposes combines the curvature correlation properties of resistivity and polarizability, by the correlativity of both analyses, judges the oily situation of target reservoir.

As shown in Figure 1, concrete grammar step is as follows for the method flow diagram of the gas-oil detecting method that the present invention proposes:

Step S1: the resistivity curve and the polarizability curve that obtain observation station in formation at target locations.

The resistivity curve of formation at target locations or polarizability curve generally can be obtained by electromagnetic inversion method.Fig. 2 obtains the resistivity curve of formation at target locations and the method flow diagram of polarizability curve, specifically can comprise following two steps:

Step S101: the observation data obtaining electromagnetic field amplitude and electromagnetic field phase in formation at target locations.

Before carrying out electromagnetic inversion, need to obtain the electromagnetic field amplitude of observation station and the observation data of electromagnetic field phase described in formation at target locations.To obtain resistivity observation data instance, two supply terminals specifically can be adopted to hit to the power supply of formation at target locations the earth, thus form underground artificial electric field.Afterwards the survey line of Resistivity testing instrument is passed into underground, measure the potential difference (PD) that two electrodes pass into each depth point, underground respectively, calculate and record electromagnetic field amplitude and the electromagnetic field phase value of each Frequency point on survey line.Choose that the multiple observation station of formation at target locations carries out measuring, record, define electromagnetic field amplitude and the electromagnetic field phase observation data of formation at target locations.Described observation station can along a line of observation, and the distance between the number of observation station and observation station is chosen according to the actual requirements.

Step S102: adopt electromagnetic inversion method described observation data inverting to be obtained to resistivity curve and the polarizability curve of observation station in formation at target locations.

Electromagnetic inversion method is adopted described observation data inverting to be obtained to resistivity curve and the polarizability curve of observation station in formation at target locations.Described electromagnetic inversion method can comprise one dimension, two dimension and the 3 D electromagnetic method of inversion.Described one dimension electromagnetic inversion method can comprise the Bostick method of inversion, the gradient method of inversion, the Gauss-Newton method of inversion, least square method and simulated annealing.The described two-dimensional electromagnetic method of inversion can comprise Aukma (hereinafter referred to as the OCCAM) method of inversion, the fast relaxation method of inversion (RRI method) and Nonlinear Conjugate Gradient Methods (NLCG method).Described one dimension electromagnetic inversion method can comprise OCCAM method, the fast relaxation method of inversion (RRI method), Nonlinear Conjugate Gradient Methods (NLCG method) and the artificial neural network method of inversion (ANN method).Resistivity-polarizability related coefficient isoline the sectional drawing of the degree of depth-survey line length two dimensions can be obtained because the inventive method is follow-up, therefore can adopt the two-dimensional electromagnetic method of inversion in embodiments of the present invention.

The OCCAM method of inversion can all observation data of matching, and computing is stablized, and step is simple, fast convergence rate, and therefore the embodiment of the present invention adopts the OCCAM method of inversion to carry out inverting to described resistivity observation data and polarizability observation data.The described OCCAM method of inversion adopts Lagrange multiplier to carry out the smooth and objective function that is data fitting degree of balance model:

$U=\left|\right|Rm|{|}^2+{\mathrm{\μ}}^{-1}\left(\right||Wd-WF\[m\]|{|}^2{X}_{*}^2)---\left(1\right)$

Wherein, μ ^-1for Lagrange multiplier, d be resistivity or polarizability observation data vector, F for just to calculate son, || Rm|| is the roughness of model, and R is roughness matrix, || Wd-WF [m] || ²for the second order norm of standard, represent that observation data d and forward modeling respond the matching difference X of F [m] ², for X ²expectation value, W is the data normalization diagonal matrix of n × n.

Inverting target to obtain a smooth enough resistivity curve or polarizability curve model, specifically can be realized by interative computation repeatedly, and the iteration expression formula of kth time is:

m _k+1＝[μR ^TR+(WJ _k) ^TWJ _k] ^-1(WJ _k) ^TWd _k(2)

In an interative computation, first to calculate the partial derivative matrix J of this iteration _kthen given a series of μ value, iterative model m (μ) is obtained by Cholesky breakdown (2), adopt the forward response F [m] of Finite Element computation model, the relative matching of employing formula (3) computation model is poor, determines μ value and the model of current iteration optimum according to matching difference.

$R_{ms}=\sqrt{\frac{\underset{j=1}{\overset{n}{\Σ}}{(d_j-F_j\[m\])}^2/d_j^2}{n}}---\left(3\right)$

Step S2: according to resistivity curve and the polarizability curve of described observation station, calculates resistivity curvature and the polarizability curvature of the depth point in described formation at target locations under observation station.

Step S102 inverting obtains resistivity curve and the polarizability curve of each observation station, according to resistivity curve and the polarizability curve of described observation station, calculates resistivity curvature and the polarizability curvature of the depth point in described formation at target locations under observation station.Fig. 3 is observation station and depth point structural representation in the earth formation, and as shown in the figure, described observation station A, B are along line of observation direction, and the line of observation can be positioned at ground.Described depth point is the Inversion Calculation point along depth direction under described observation station, and as shown in Figure 3, some A1, A2, A3 and A4 are the Inversion Calculation depth points of observation station A, and some B1, B2, B3 and B4 are the Inversion Calculation depth points with observation station B.The computing formula of the resistivity curvature of described depth point is:

$k_{\mathrm{\ρ}}=\left(\frac{{\mathrm{\ρ}}^{\′\′}}{{(1+{\mathrm{\ρ}}^{\′2})}^{\frac{3}{2}}}\right)---\left(4\right)$

Wherein, k _ρfor the resistivity curvature of described depth point, ρ " for described resistivity curve is in the second derivative at described depth point place, ρ ' is for described resistivity curve is in the first order derivative of described sight depth.

The computing formula of the polarizability curvature of described depth point is:

$k_{\mathrm{\η}}=\left(\frac{{\mathrm{\η}}^{\′\′}}{{(1+{\mathrm{\η}}^{\′2})}^{\frac{3}{2}}}\right)---\left(5\right)$

Wherein, k _ηfor the polarizability curvature of described depth point, η " for described polarizability curve is in the second derivative of described sight depth, η ' is for described polarizability curve is in the first order derivative at described depth point place.

In step s 2, by the calculating of resistivity curvature and polarizability curvature in electromagnetic survey, extract the effective information reflecting whether oily

Step S3: calculate the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature.

There are resistivity curvature corresponding with it and polarizability curvature in each depth point under formation at target locations observation station.Calculate the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature.The computing formula of the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature comprises:

$r^i=\frac{n\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i{k}_{\mathrm{\η}}^i-\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\·\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i}{\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\ρ}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\right)}^2}\·\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\η}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i\right)}^2}}---\left(6\right)$

Wherein, r ⁱfor the resistivity curvature of described observation station at i-th depth point place and the related coefficient of polarizability curvature, n is the number of the depth point of described observation station, for described observation station is in the resistivity curvature at i-th depth point place, for described observation station is in the polarizability curvature at i-th depth point place.

The electromagnetism rate curvature that step S3 calculates based on step S2 and polarizability curvature, the correlativity of both analyses, effectively rejects spurious anomaly, enhances the ability differentiating effective anomaly, improves the accuracy of identification of hydrocarbon-zone.

Step S4: according to the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature, draw the related coefficient isoline sectional drawing of described objective area, when described related coefficient isoline numerical value is greater than first threshold, described isoline region is hydrocarbon-zone.

Described related coefficient isoline sectional drawing to draw along the result of the related coefficient of depth direction multiple depth point resistivity curvature and polarizability curvature according to each observation station to form.Described related coefficient isoline sectional drawing can reflect the change of the related coefficient along line direction (direction of observation station when straight line) and depth direction resistivity curvature and polarizability curvature.The sectional horizontal ordinate of described related coefficient isoline can be survey line length, and ordinate can be the degree of depth.Generally, hydrocarbon-zone has the advantages that resistivity curvature anomalies is high and polarizability curvature anomalies is high, and time extremely high while of therefore, both related coefficients are also relatively high.And exception can not appear in the resistivity curvature of hydrocarbon-zone and polarizability curvature simultaneously, therefore both related coefficients are less.In sum, the related coefficient of resistivity curvature and polarizability curvature is directly proportional to the probability of formation at target locations oily, and related coefficient is larger, and the probability of formation at target locations oily is larger.When described related coefficient isoline numerical value is greater than first threshold, described isoline region is hydrocarbon-zone.Described first threshold is arranged according to the data of known prospect pit.

Fig. 4 is the related coefficient isoline sectional drawing in certain work area, and being 6 kms in survey line length, is that in the stratum of 5 kms, the distribution of related coefficient isoline as shown in the figure deeply.Can judge that the degree of depth is 3 kms, survey line length is 3.5 km places is oil and gas reservoir center.According to the data of known prospect pit, arranging first threshold is 0.2, and so related coefficient is greater than the isoline region of 0.2 is hydrocarbon-zone.Hydrocarbon-zone in Fig. 4 is the region that isoline 0.2 surrounds.According to the determination of above-mentioned hydrocarbon-zone, can judge that W1 is dry-well, W2 is producing oil well.

Introduce the oil and gas detection system corresponding with above-described embodiment method below, Fig. 5 is the structural representation of described oil and gas detection system, as shown in the figure, described system 50 comprises: acquiring unit 51, curvature estimation unit 52, Calculation of correlation factor unit 53, sectional drawing drawing unit 54, wherein

Acquiring unit 51, for obtaining resistivity curve and the polarizability curve of observation station in formation at target locations.

Curvature estimation unit 52, for according to described resistivity curve and polarizability curve, calculates resistivity curvature and the polarizability curvature of the depth point in described formation at target locations under observation station.

Calculation of correlation factor unit 53, for the related coefficient of the resistivity curvature and polarizability curvature that calculate the depth point under described observation station.

Sectional drawing drawing unit 54, for according to the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature, draw the related coefficient isoline sectional drawing of described objective area, when described related coefficient isoline numerical value is greater than first threshold, described isoline region is hydrocarbon-zone.

Fig. 6 is the structural representation of described acquiring unit 51, and as shown in the figure, described acquiring unit 51 also comprises: obtain subelement 61, inverting unit 62, wherein,

Obtain subelement 61, for the observation data of the electromagnetic field amplitude and electromagnetic field phase that obtain observation station in formation at target locations.

Inverting unit 62, obtains resistivity curve and the polarizability curve of observation station in formation at target locations for adopting electromagnetic inversion method to described observation data inverting.

The technique scheme of embodiment of the present invention introduction has following beneficial effect: based on calculating and the correlation analysis of resistivity curvature and polarizability curvature in electromagnetic survey, effectively reject spurious anomaly, extracts the effective information reflecting whether oily.According to the effective hydrocarbon-zone of correlation analysis identification, enhance the ability differentiating effective anomaly, improve the accuracy of identification of hydrocarbon-zone.

Those skilled in the art can also recognize the various illustrative components, blocks (illustrativelogicalblock) that the embodiment of the present invention is listed, unit, and step can pass through electronic hardware, computer software, or both combinations realize.For the replaceability (interchangeability) of clear displaying hardware and software, above-mentioned various illustrative components (illustrativecomponents), unit and step have universally described their function.Such function is the designing requirement realizing depending on specific application and whole system by hardware or software.Those skilled in the art for often kind of specifically application, can use the function described in the realization of various method, but this realization can should not be understood to the scope exceeding embodiment of the present invention protection.

Various illustrative logical block described in the embodiment of the present invention, or unit can pass through general processor, digital signal processor, special IC (ASIC), field programmable gate array (FPGA) or other programmable logic device, discrete gate or transistor logic, discrete hardware components, or the design of above-mentioned any combination realizes or operates described function.General processor can be microprocessor, and alternatively, this general processor also can be any traditional processor, controller, microcontroller or state machine.Processor also can be realized by the combination of calculation element, such as digital signal processor and microprocessor, multi-microprocessor, and a Digital Signal Processor Core combined by one or more microprocessor, or other similar configuration any realizes.

The software module that method described in the embodiment of the present invention or the step of algorithm directly can embed hardware, processor performs or the combination of both.Software module can be stored in the storage medium of other arbitrary form in RAM storer, flash memory, ROM storer, eprom memory, eeprom memory, register, hard disk, moveable magnetic disc, CD-ROM or this area.Exemplarily, storage medium can be connected with processor, with make processor can from storage medium reading information, and write information can be deposited to storage medium.Alternatively, storage medium can also be integrated in processor.Processor and storage medium can be arranged in ASIC, and ASIC can be arranged in user terminal.Alternatively, processor and storage medium also can be arranged in the different parts in user terminal.

In one or more exemplary design, the above-mentioned functions described by the embodiment of the present invention can realize in the combination in any of hardware, software, firmware or this three.If realized in software, these functions can store on the medium with computer-readable, or are transmitted on the medium of computer-readable with one or more instruction or code form.Computer readable medium comprises computer storage medium and is convenient to make to allow computer program transfer to the telecommunication media in other place from a place.Storage medium can be that any general or special computer can the useable medium of access.Such as, such computer readable media can include but not limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage device, or other anyly may be used for carrying or store the medium that can be read the program code of form with instruction or data structure and other by general or special computer or general or special processor.In addition, any connection can be properly termed computer readable medium, such as, if software is by a concentric cable, fiber optic cables, twisted-pair feeder, Digital Subscriber Line (DSL) or being also comprised in defined computer readable medium with wireless way for transmittings such as such as infrared, wireless and microwaves from a web-site, server or other remote resource.Described video disc (disk) and disk (disc) comprise Zip disk, radium-shine dish, CD, DVD, floppy disk and Blu-ray Disc, and disk is usually with magnetic duplication data, and video disc carries out optical reproduction data with laser usually.Above-mentioned combination also can be included in computer readable medium.

The above is only the preferred embodiment of the present invention; it should be pointed out that for those skilled in the art, under the premise without departing from the principles of the invention; can also make some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.

Claims (12)

1. a gas-oil detecting method, is characterized in that, described method comprises:

Obtain resistivity curve and the polarizability curve of observation station in formation at target locations;

According to described resistivity curve and polarizability curve, calculate resistivity curvature and the polarizability curvature of the depth point in described formation at target locations under observation station;

Calculate the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature;

According to the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature, draw the related coefficient isoline sectional drawing of described formation at target locations; Identify that the isoline region that correlation coefficient value in described isoline sectional drawing is greater than first threshold is hydrocarbon-zone.

2. gas-oil detecting method according to claim 1, is characterized in that, the computing formula of the resistivity curvature of the depth point in the described formation at target locations of described calculating under observation station comprises:

$k_{\mathrm{\ρ}}=\left(\frac{{\mathrm{\ρ}}^{\′\′}}{{(1+{\mathrm{\ρ}}^{\′2})}^{\frac{3}{2}}}\right)$

Wherein, k _ρfor the resistivity curvature of described depth point, ρ " for described resistivity curve is in the second derivative at described depth point place, ρ ' is for described resistivity curve is in the first order derivative of described sight depth.

3. gas-oil detecting method according to claim 1, is characterized in that, the computing formula of the polarizability curvature of the depth point in the described formation at target locations of described calculating under observation station comprises:

$k_{\mathrm{\η}}=\left(\frac{{\mathrm{\η}}^{\′\′}}{{(1+{\mathrm{\η}}^{\′2})}^{\frac{3}{2}}}\right)$

Wherein, k _ηfor the polarizability curvature of described depth point, η " for described polarizability curve is in the second derivative of described sight depth, η ' is for described polarizability curve is in the first order derivative at described depth point place.

4. gas-oil detecting method according to claim 1, is characterized in that, the computing formula of the resistivity curvature of the depth point under the described observation station of described calculating and the related coefficient of polarizability curvature comprises:

$r^i=\frac{n\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i{k}_{\mathrm{\η}}^i-\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\·\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i}{\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\ρ}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\right)}^2}\·\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\η}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i\right)}^2}}$

Wherein, r ⁱfor the resistivity curvature of described observation station at i-th depth point place and the related coefficient of polarizability curvature, n is the number of the depth point of described observation station, for described observation station is in the resistivity curvature at i-th depth point place, for described observation station is in the polarizability curvature at i-th depth point place.

5. gas-oil detecting method according to claim 1, is characterized in that, the resistivity curve of observation station and polarizability curve in described acquisition formation at target locations, comprising:

Obtain the electromagnetic field amplitude of observation station and the observation data of electromagnetic field phase in formation at target locations;

Electromagnetic inversion method is adopted described observation data inverting to be obtained to resistivity curve and the polarizability curve of observation station in formation at target locations.

6. gas-oil detecting method according to claim 5, is characterized in that, described electromagnetic inversion method comprises the Aukma method of inversion.

7. an oil and gas detection system, is characterized in that, comprising:

Acquiring unit, for obtaining resistivity curve and the polarizability curve of observation station in formation at target locations;

Curvature estimation unit, for according to described resistivity curve and polarizability curve, calculates resistivity curvature and the polarizability curvature of the depth point in described formation at target locations under observation station;

Calculation of correlation factor unit, for the related coefficient of the resistivity curvature and polarizability curvature that calculate the depth point under described observation station;

Sectional drawing drawing unit, for according to the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature, draw the related coefficient isoline sectional drawing of described formation at target locations, identify that the isoline region that the sectional related coefficient numerical value of described isoline is greater than first threshold is hydrocarbon-zone.

8. oil and gas detection system according to claim 7, is characterized in that, described curvature estimation unit comprises: the computing formula of the resistivity curvature of described depth point comprises:

$k_{\mathrm{\ρ}}=\left(\frac{{\mathrm{\ρ}}^{\′\′}}{{(1+{\mathrm{\ρ}}^{\′2})}^{\frac{3}{2}}}\right)$

Wherein, k _ρfor the resistivity curvature of described depth point, ρ " for described resistivity curve is in the second derivative at described depth point place, ρ ' is for described resistivity curve is in the first order derivative of described sight depth.

9. oil and gas detection system according to claim 7, is characterized in that, described curvature estimation unit comprises: the computing formula of the polarizability curvature of described depth point comprises:

$k_{\mathrm{\η}}=\left(\frac{{\mathrm{\η}}^{\′\′}}{{(1+{\mathrm{\η}}^{\′2})}^{\frac{3}{2}}}\right)$

Wherein, k _ηfor the polarizability curvature of described depth point, η " for described polarizability curve is in the second derivative of described sight depth, η ' is for described polarizability curve is in the first order derivative at described depth point place.

10. oil and gas detection system according to claim 7, is characterized in that, described Calculation of correlation factor unit comprises: the computing formula of the resistivity curvature of the depth point under described observation station and the related coefficient of polarizability curvature comprises:

$r^i=\frac{n\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i{k}_{\mathrm{\η}}^i-\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\·\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i}{\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\ρ}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\ρ}}^i\right)}^2}\·\sqrt{n\underset{i}{\overset{n}{\Σ}}{\left(k_{\mathrm{\η}}^i\right)}^2-{\left(\underset{i}{\overset{n}{\Σ}}{k}_{\mathrm{\η}}^i\right)}^2}}$

Wherein, r ⁱfor the resistivity curvature of described observation station at i-th depth point place and the related coefficient of polarizability curvature, n is the number of the depth point of described observation station, for described observation station is in the resistivity curvature at i-th depth point place, for described observation station is in the polarizability curvature at i-th depth point place.

11. oil and gas detection systems according to claim 7, it is characterized in that, described acquiring unit comprises:

Obtain subelement, for the observation data of the electromagnetic field amplitude and electromagnetic field phase that obtain observation station in formation at target locations;

Inverting unit, obtains resistivity curve and the polarizability curve of observation station in formation at target locations for adopting electromagnetic inversion method to described observation data inverting.

12. oil and gas detection systems according to claim 11, is characterized in that, described electromagnetic inversion method comprises the Aukma method of inversion.