CN101892838B - Method and device for obtaining high-resolution well logging curve - Google Patents

Abstract

The embodiment of the invention provides a method and a device for obtaining a high-resolution well logging curve, and relates to the oil and gas exploration technical field. The method comprises the following steps: collecting a target geological parameter of a target core through a core experiment; acquiring multi-stage calibration curves corresponding to correction formulas according to multi-stage deconvolution correction formulas and the original well logging curve; and judging the multi-stage calibration curve closest to the collected target geological parameter according to the acquired multi-stage calibration curves. The method and the device provided by the embodiment of the invention can correct the well logging curves at different strata, thus improving the resolution of the well logging curve and acquiring the high-resolution well logging curve.

Description

Obtain method and the device of high-resolution well logging curve

Technical field

The embodiment of the invention relates to the oil-gas exploration technical field, particularly a kind of method and device that obtains high-resolution well logging curve.

Background technology

Along with the fast development of human civilization, the demand of the energy is grown with each passing day, wherein the demand to oil, natural gas is more strong.In order effectively to survey oil reservoir, can carry out the geologic parameter collection to formation at target locations when usually logging well on the stratum, these geologic parameters comprise formation resistivity, radioactivity etc., be linked to be the corresponding geological condition that a curve (being log) is described formation at target locations by a large amount of geologic parameters that will get access to same alike result, thereby help the geologic structure that the researcher can understanding more directly perceived stratum, judge.

Yet the structure on stratum is very complicated, and such as for some thinner stratum, the resolution ratio of log is lower, can not describe accurately the stratum.Although can proofread and correct log in the prior art, this correction does not have specific aim, namely proofreaies and correct the log of Different Strata different characteristic with unified formula, so calibration result is not good.

Summary of the invention

The purpose of the embodiment of the invention is to provide a kind of method and device that obtains high-resolution well logging curve, in analyzing earth formation, the method that the embodiment of the invention provides and device can be proofreaied and correct for the log of Different Strata, improve the resolution ratio of log, thereby obtain high-resolution log.

For achieving the above object, the embodiment of the invention provides a kind of method of obtaining high-resolution well logging curve, and described method comprises:

Gather the target geologic parameter of target rock core by core experiment;

Obtain the multistage calibration curve of corresponding each updating formula according to each multistage Deconvolution Correction formula and original log, described multistage Deconvolution Correction formula is 3+(n-1) * 2 grades of Deconvolution Correction formula, wherein n is natural number;

Judge and the immediate multistage calibration curve of the described target geologic parameter that collects according to each the multistage calibration curve that gets access to.

In order to gather more accurately geologic data, when gathering the target geologic parameter of described target rock core, sampling interval is not more than 0.125 meter, such as 0.1 meter.

Wherein each multistage Deconvolution Correction formula comprises three grades updating formula when n is 1 at least:

XR(z _i)＝M[-X(z _i-1)+N·X(z _i)-X(z _i+1)]

XR (Z wherein _i) corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ²

Each multistage Deconvolution Correction formula also comprises the updating formula of the Pyatyi when n is 2:

$\mathrm{XR}\left(z_i\right)=\frac{1}{4}M[(\mathrm{\λ}-1)X\left(z_{i+2}\right)-4\mathrm{\λ}\·X\left(z_{i+1}\right)+(4N+6\mathrm{\λ}+2)X\left(z_i\right)-4\mathrm{\λ}\·X\left(z_{i-1}\right)+(\mathrm{\λ}-1)X\left(Z_{i-2}\right)]$

XR (Z wherein _i) corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ², λ and μ are weight coefficient.

Each multistage Deconvolution Correction formula also comprises seven grades updating formula when n is 3:

$\mathrm{XR}\left(z_i\right)=\frac{1}{36}M[-4\mathrm{\λ}\·X\left(z_{i+3}\right)-9\mathrm{\μ}\·X\left(z_{i+2}\right)-36(1-\mathrm{\λ}-\mathrm{\μ})X\left(z_{i+1}\right)+(36N-64\mathrm{\λ}-64\mathrm{\μ})X\left(z_i\right)$

$-36(1-\mathrm{\λ}-\mathrm{\μ})X\left(z_{i-1}\right)-9\mathrm{\μ}\·X\left(Z_{i-2}\right)-4\mathrm{\λ}\·X\left(Z_{i-3}\right)]$

XR (Z wherein _i) corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ², λ and μ are weight coefficient.

For achieving the above object, the embodiment of the invention also provides a kind of device that obtains high-resolution well logging curve, and described device comprises:

Collecting unit is used for the target geologic parameter by core experiment collection target rock core;

Acquiring unit is used for obtaining according to each multistage Deconvolution Correction formula and original log the multistage calibration curve of corresponding each updating formula, and described multistage Deconvolution Correction formula is 3+(n-1) * 2 grades of Deconvolution Correction formula, wherein n is natural number;

Decision package, the multistage calibration curve of each that is used for getting access to according to described acquiring unit is judged the immediate multistage calibration curve of target geologic parameter that collects with described collecting unit.

Described device also comprises:

Memory cell is used for storage 3+(n-1) * 2 grades of Deconvolution Correction formula, wherein n is natural number.

Wherein said 3+(n-1) * 2 grade Deconvolution Correction formula comprises that at least n is 1 o'clock three grades updating formula:

XR(z _i)＝M[-X(z _i-1)+N·X(z _i)-X(z _i+1)]

XR (Z wherein _i) the geologic parameter corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) the geologic parameter original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ²

In order effectively to obtain multistage calibration curve, described acquiring unit comprises:

Computing module is used for the 3+(n-1 according to cell stores) * 2 grades of Deconvolution Correction formula calculating geologic parameter corrected values;

Curve fitting module is used for the geologic parameter corrected value matched curve that calculates according to computing module.

The advantage of the embodiment of the invention is and can effectively proofreaies and correct the log of Different Strata, improves the resolution ratio of log, thereby according to the geological condition on the high-resolution well logging curve accurate analysis stratum that gets access to.

Description of drawings

In order to be illustrated more clearly in the technical scheme of the embodiment of the invention, the accompanying drawing of required use was done to introduce simply during the below will describe embodiment, apparently, accompanying drawing in the following describes only is some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain according to these accompanying drawings other accompanying drawing.

Fig. 1 is a kind of schematic flow diagram that obtains the method for high-resolution well logging curve of the embodiment of the invention.

Fig. 2 is the comparison diagram of log before and after core experiment gathers target coregamma Value Data and proofreaies and correct.

Fig. 3 is the functional block diagram of obtaining the high-resolution well logging curve device of the embodiment of the invention.

The specific embodiment

Below in conjunction with the accompanying drawing in the embodiment of the invention, the technical scheme in the embodiment of the invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, rather than whole embodiment.Illustrative examples of the present invention and explanation thereof are used for explaining the present invention; but not as a limitation of the invention; based on the embodiment among the present invention; those of ordinary skills belong to the scope of protection of the invention not making the every other embodiment that obtains under the creative work prerequisite.

The embodiment of the invention provides a kind of technology of obtaining high-resolution well logging curve, this technology by to the correction of log so that the log after proofreading and correct meets the true geological condition of formation at target locations more, thereby establish good basis for the geology of accurate description formation at target locations.

Based on this, present embodiment provides a kind of method of obtaining high-resolution well logging curve, and as shown in Figure 1, Fig. 1 is a kind of schematic flow diagram that obtains the method for high-resolution well logging curve of present embodiment, and this flow process comprises the steps:

101. gather the target geologic parameter of target rock core by core experiment.

Present embodiment gets the target geologic parameter that gathers first the target rock core, this step gathers by core experiment, wherein the target rock core is the geological specimen from the sampling of target Log Strata, this sample can represent the geological characteristics on this stratum, if such as proofreading and correct the original log of surveying thinner stratum, that just will gather rock core in this thinner stratum.And the target geologic parameter refers to the geologic parameter that wish gathers, and this parameter is the parameter identical with original log geological property, and such as being what to describe than thin strate resistivity when original log, that target geologic parameter is exactly testing resistance rate parameter; When original log is to describe than thin strate radioactive (curve that namely is comprised of the natural gamma value), that target geologic parameter is exactly natural gamma value to be measured.

And in order to gather more accurately geologic data, this step is when gathering the target geologic parameter of target rock core, and sampling interval is not more than 0.125 meter, such as on the target rock core every 0.1 meter once sampling geologic data.

102. obtain the multistage calibration curve of corresponding each updating formula according to each multistage Deconvolution Correction formula and original log.

Present embodiment obtains calibration curve by multistage Deconvolution Correction formula, and this multistage deconvolution formula can be 3+(n-1) * 2 grades of Deconvolution Correction formula, wherein n is natural number.This multistage Deconvolution Correction formula is at least 3 grades of deconvolution like this, so this each multistage Deconvolution Correction formula comprises 3 grades updating formula when n is 1 at least in the present embodiment:

XR(z _i)＝M[-X(z _i-1)+N·X(z _i)-X(z _i+1)]

XR (Z wherein _i) the geologic parameter corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) the geologic parameter original value of corresponding target collection point on the original log of expression, this parameter can obtain Z by original log _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ²

This formula can be derived by following method, when having determined collection point Z _iAfter, can determine apart from this collection point Z _iThe neighbouring sample point Z of a sampling interval Δ z _I+1And Z _I-1, because this geologic parameter corrected value XR (Z _i) can be expressed as:

$\mathrm{XR}\left(z_i\right)=X\left(z_i\right)-\frac{1}{{\mathrm{\α}}^2}\frac{d^{2}X}{{\mathrm{dz}}^2}{|}_{z=z_i}$

So in the discretization second dervative, but by second order center difference coefficient derived expression:

$\frac{d^{2}X}{{\mathrm{dz}}^2}{|}_{z=z_i}\≈\frac{X(z_i+\mathrm{\Δz})-2X\left(z_i\right)+X(z_i-\mathrm{\Δz})}{{\left(\mathrm{\Δz}\right)}^2}$

With derived expression substitution geologic parameter corrected value XR (Z _i) behind the formula, can obtain 3 grades of Deconvolution Correction formula:

XR(z _i)＝M[-X(z _i-1)+N·X(z _i)-X(z _i+1)]

In order to proofread and correct more accurately original log, determine and the immediate calibration curve of target geologic parameter that present embodiment can obtain calibration curve by more multistage Deconvolution Correction formula.Such as 5 grades updating formula when n is 2:

$\mathrm{XR}\left(z_i\right)=\frac{1}{4}M[(\mathrm{\λ}-1)X\left(z_{i+2}\right)-4\mathrm{\λ}\·X\left(z_{i+1}\right)+(4N+6\mathrm{\λ}+2)X\left(z_i\right)-4\mathrm{\λ}\·X\left(z_{i-1}\right)+(\mathrm{\λ}-1)X\left(Z_{i-2}\right)]$

When obtaining formula according to weighted average:

$\frac{d^{2}X}{d_{z^2}}{|}_{z=z_i}\≈\mathrm{\λ}\frac{X\left(z_{i+1}\right)-2X\left(z_i\right)+X\left(z_{i-1}\right)}{{\left(\mathrm{\Δz}\right)}^2}+(1-\mathrm{\λ})\frac{X\left(z_{i+2}\right)-2X\left(z_i\right)+X\left(z_{i-2}\right)}{{\left(2\mathrm{\Δz}\right)}^2}$

After, with its substitution geologic parameter corrected value XR (Z _i) formula can obtain 5 grades updating formula.Wherein the implication of each parameter is identical with 3 grades updating formula, thus repeat no more, and Z _I+2And Z _I-2The sampled point that has two sampling intervals for the target collection point.

In like manner, multistage Deconvolution Correction formula also comprises seven grades updating formula when n is 3:

$\mathrm{XR}\left(z_i\right)=\frac{1}{36}M[-4\mathrm{\λ}\·X\left(z_{i+3}\right)-9\mathrm{\μ}\·X\left(z_{i+2}\right)-36(1-\mathrm{\λ}-\mathrm{\μ})X\left(z_{i+1}\right)+(36N-64\mathrm{\λ}-64\mathrm{\μ})X\left(z_i\right)$

$-36(1-\mathrm{\λ}-\mathrm{\μ})X\left(z_{i-1}\right)-9\mathrm{\μ}\·X\left(Z_{i-2}\right)-4\mathrm{\λ}\·X\left(Z_{i-3}\right)]$

Wherein the implication of each parameter is identical with 3 grades updating formula, thus repeat no more, and Z _I+3And Z _I-3The sampled point that has three sampling intervals for the target collection point.

So just can obtain the corresponding discrete geology parameter correction value of each Deconvolution Correction formula, each connection of discrete geology parameter correction value that these are adjacent can obtain the corresponding multistage calibration curve of each Deconvolution Correction formula.

203. judge and the immediate multistage calibration curve of the described target geologic parameter that collects according to each the multistage calibration curve that gets access to.

Because the geologic parameter by the core experiment collection is more accurately, so the closer to the geologic parameter that gathers by core experiment, then the calibration result of calibration curve is better.This step is judged and the immediate multistage calibration curve of the described target geologic parameter that collects according to each multistage calibration curve that step 202 gets access to.Such as coincideing (namely the geologic parameter by the core experiment collection nearly all drops on the curve) when 3 grades of nearly all corresponding with 3 grades of updating formulas calibration curves of the geologic parameter that gathers by core experiment, these 3 grades of calibration curves are exactly near the curve of target geologic parameter so.For oil well corresponding to target rock core, the geological condition of this oil well can be described by 3 grades of calibration curves.And for other oil wells on this stratum, then can come log is proofreaied and correct by 3 grades of Deconvolution Correction formula.Can effectively proofread and correct the log of Different Strata like this, improve the resolution ratio of log, thereby according to the geological condition on the high-resolution well logging curve accurate analysis stratum that gets access to.

For those skilled in the art can clearerly understand the described method of the embodiment of the invention; the below quotes an example and illustrates; it should be noted that numeral and geologic parameter in this example only are exemplary; can not be with its explanation as restriction protection domain of the present invention, those skilled in the art should know that understanding the described method of the embodiment of the invention can also express by other concrete examples.

As shown in Figure 2, Fig. 2 is the comparison diagram of log before and after core experiment gathers target coregamma Value Data and proofreaies and correct, can know from Fig. 2 and to find out that the target rock core is 2313～2320 meters rock core sample, coregamma value in this section degree of depth obtains and has been reflected among Fig. 2 by the core experiment collection, and wherein sampling interval Δ z is 0.125 meter.Full curve among Fig. 2 is original log, and this log is the radioactive log in this stratum of description of field survey.

For this original log is proofreaied and correct, this example calculations three multistage Deconvolution Correction curves: 3 grades of Deconvolution Correction curves, 5 grades of Deconvolution Correction curves and 7 grades of Deconvolution Correction curves, wherein get α=2 in 3 grades of Deconvolution Correction formula, Δ z=0.125, so M=16, N=2.0625 is with these parameter substitution formula, namely

XR(z _i)＝16[-X(z _i-1)+2.0625·X(z _i)-X(z _i+1)]

According to the geologic parameter corrected value that can calculate each collection point after the geologic parameter original value substitution of original log with the collection point, will be 3 grades of Deconvolution Correction curves after the adjacent corrected value connection.In like manner, when getting λ=0.6, μ=0.4 o'clock, 5 grades of Deconvolution Correction formula are:

$\mathrm{XR}\left(z_i\right)=\frac{1}{4}\×16[(0.6-1)X\left(z_{i+2}\right)-4\×0.6X\left(z_{i+1}\right)+(4\×2.0625+6\×0.6+2)X\left(z_i\right)-4\×0.6X\left(z_{i-1}\right)+(0.6-1)X\left(Z_{i-2}\right)]$

Thereby also can calculate the geologic parameter corrected value of each collection point, be 5 grades of Deconvolution Correction curves after adjacent corrected value is connected.When getting λ=0.2, μ=0.3 o'clock, 7 grades of Deconvolution Correction formula are:

$\mathrm{XR}\left(z_i\right)=\frac{1}{36}\×16[-4\×0.2\mathrm{GR}\left(z_{i+3}\right)-9\×0.3\mathrm{GR}\left(z_{i+2}\right)-36\×0.5\mathrm{GR}\left(z_{i+1}\right)+(36\×2.0625-64\×0.5)\mathrm{GR}\left(z_i\right)$

$-36\×0.5\mathrm{GR}\left(z_{i-1}\right)-9\×0.3\mathrm{GR}\left(Z_{i-2}\right)-4\×0.2\mathrm{GR}\left(Z_{i-3}\right)]$

Thereby also can calculate the geologic parameter corrected value of each collection point, be 7 grades of Deconvolution Correction curves after adjacent corrected value is connected.The target coregamma value and each the multistage calibration curve that collect are compared, and can judge with the immediate multistage calibration curve of the target coregamma value that collects is 7 grades of Deconvolution Correction curves.From Fig. 2, also can know and find out that GRR7 curve (i.e. 7 grades of Deconvolution Correction curves) and coregamma can coincide preferably.

The device of realizing said method embodiment can be an autonomous device, and the below describes the device of realizing above-mentioned functions.As shown in Figure 3, Fig. 3 is the functional block diagram of obtaining the high-resolution well logging curve device of the embodiment of the invention, and as can be known from Fig. 3, this device can comprise:

Collecting unit 301 is used for the target geologic parameter by core experiment collection target rock core;

Acquiring unit 302 is used for obtaining according to each multistage Deconvolution Correction formula and original log the multistage calibration curve of corresponding each updating formula;

Decision package 303, the multistage calibration curve of each that is used for getting access to according to acquiring unit 302 is judged the immediate multistage calibration curve of target geologic parameter that collects with collecting unit 301.

Memory cell 304 is used for storage 3+(n-1) * 2 grades of Deconvolution Correction formula, wherein n is natural number.Wherein said 3+(n-1) * 2 grade Deconvolution Correction formula comprises that at least n is 1 o'clock three grades updating formula:

XR(z _i)＝M[-X(z _i-1)+N·X(z _i)-X(z _i+1)]

XR (Z wherein _i) the geologic parameter corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) the geologic parameter original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ²Each multistage Deconvolution Correction formula also comprises the updating formula of the Pyatyi when n is 2:

$\mathrm{XR}\left(z_i\right)=\frac{1}{4}M[(\mathrm{\λ}-1)X\left(z_{i+2}\right)-4\mathrm{\λ}\·X\left(z_{i+1}\right)+(4N+6\mathrm{\λ}+2)X\left(z_i\right)-4\mathrm{\λ}\·X\left(z_{i-1}\right)+(\mathrm{\λ}-1)X\left(Z_{i-2}\right)]$

XR (Z wherein _i) corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ², λ and μ are weight coefficient.

Each multistage Deconvolution Correction formula also comprises seven grades updating formula when n is 3:

$\mathrm{XR}\left(z_i\right)=\frac{1}{36}M[-4\mathrm{\λ}\·X\left(z_{i+3}\right)-9\mathrm{\μ}\·X\left(z_{i+2}\right)-36(1-\mathrm{\λ}-\mathrm{\μ})X\left(z_{i+1}\right)+(36N-64\mathrm{\λ}-64\mathrm{\μ})X\left(z_i\right)$

$-36(1-\mathrm{\λ}-\mathrm{\μ})X\left(z_{i-1}\right)-9\mathrm{\μ}\·X\left(Z_{i-2}\right)-4\mathrm{\λ}\·X\left(Z_{i-3}\right)]$

XR (Z wherein _i) corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ², λ and μ are weight coefficient.

In order effectively to obtain multistage calibration curve, acquiring unit 302 comprises:

Computing module 321 is used for the 3+(n-1 according to memory cell 304 storages) * 2 grades of Deconvolution Correction formula calculating geologic parameter corrected values;

Curve fitting module 322 is used for the geologic parameter corrected value matched curve that calculates according to computing module 321.

The described device of present embodiment can be the equipment that computer etc. has information processing capability, effectively proofread and correct by the log to Different Strata, improve the resolution ratio of log, thereby according to the geological condition on the high-resolution well logging curve accurate analysis stratum that gets access to.

Certainly; above-described specific embodiment; purpose of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the above is specific embodiments of the invention only, and the protection domain that is not intended to limit the present invention is all within logic of the present invention and principle; any modification of making, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (9)

1. a method of obtaining high-resolution well logging curve is characterized in that, described method comprises:

Gather the target geologic parameter of target rock core by core experiment;

Obtain the multistage calibration curve of corresponding each updating formula according to each multistage Deconvolution Correction formula and original log, described multistage Deconvolution Correction formula is 3+(n-1) * 2 grades of Deconvolution Correction formula, wherein n is natural number;

Judge and the immediate multistage calibration curve of the described target geologic parameter that collects according to each the multistage calibration curve that gets access to.

2. method according to claim 1 is characterized in that, when gathering the target geologic parameter of described target rock core, sampling interval is not more than 0.125 meter.

3. method according to claim 2 is characterized in that, described each multistage Deconvolution Correction formula comprises three grades updating formula when n is 1 at least:

XR(z _i)＝M[-X(z _i-1)+N·X(z _i)-X(z _i+1)]

XR (Z wherein _i) the geologic parameter corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) the geologic parameter original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ²

4. method according to claim 3 is characterized in that, described each multistage Deconvolution Correction formula also comprises the updating formula of the Pyatyi when n is 2:

$\mathrm{XR}\left(z_i\right)=\frac{1}{4}M[(\mathrm{\λ}-1)X\left(z_{i+2}\right)-4\mathrm{\λ}\·X\left(z_{i+1}\right)+(4N+6\mathrm{\λ}+2)X\left(z_i\right)-4\mathrm{\λ}\·X\left(z_{i-1}\right)+(\mathrm{\λ}-1)X\left(Z_{i-2}\right)]$

XR (Z wherein _i) corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+2And Z _I-2For the target collection point has the sampled point of two sampling intervals, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ², λ and μ are weight coefficient.

5. method according to claim 4 is characterized in that, described each multistage Deconvolution Correction formula also comprises seven grades updating formula when n is 3:

$\mathrm{XR}\left(z_i\right)=\frac{1}{36}M[-4\mathrm{\λ}\·X\left(z_{i+3}\right)-9\mathrm{\μ}\·X\left(z_{i+2}\right)-36(1-\mathrm{\λ}-\mathrm{\μ})X\left(z_{i+1}\right)+(36N-64\mathrm{\λ}-64\mathrm{\μ})X\left(z_i\right)$

$-36(1-\mathrm{\λ}-\mathrm{\μ})X\left(z_{i-1}\right)-9\mathrm{\μ}\·X\left(Z_{i-2}\right)-4\mathrm{\λ}\·X\left(Z_{i-3}\right)]$

XR (Z wherein _i) corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+3And Z _I-3For the target collection point has the sampled point of three sampling intervals, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ², λ and μ are weight coefficient.

6. a device that obtains high-resolution well logging curve is characterized in that, described device comprises:

Collecting unit is used for the target geologic parameter by core experiment collection target rock core;

Acquiring unit is used for obtaining according to each multistage Deconvolution Correction formula and original log the multistage calibration curve of corresponding each updating formula, and described multistage Deconvolution Correction formula is 3+(n-1) * 2 grades of Deconvolution Correction formula, wherein n is natural number;

Decision package, the multistage calibration curve of each that is used for getting access to according to described acquiring unit is judged the immediate multistage calibration curve of target geologic parameter that collects with described collecting unit.

7. device according to claim 6 is characterized in that, described device also comprises:

Memory cell is used for storage 3+(n-1) * 2 grades of Deconvolution Correction formula, wherein n is natural number.

8. device according to claim 7 is characterized in that, described 3+(n-1) * 2 grades of Deconvolution Correction formula comprise that at least n is 1 o'clock three grades updating formula:

XR(z _i)＝M[-X(z _i-1)+N·X(z _i)-X(z _i+1)]

XR (Z wherein _i) the geologic parameter corrected value of target collection point on the corresponding target rock core of expression, X (Z _i) the geologic parameter original value of corresponding target collection point on the original log of expression, Z _iBe target collection point, Z _I+1And Z _I-1For the target collection point has the neighbouring sample point of a sampling interval, M=1/(α Δ z) ², α is geologic characteristic parameter α, Δ z is sampling interval, N=2+(α Δ z) ²

9. device according to claim 8 is characterized in that, described acquiring unit comprises:

Computing module is used for the 3+(n-1 according to cell stores) * 2 grades of Deconvolution Correction formula calculating geologic parameter corrected values;

Curve fitting module is used for the geologic parameter corrected value matched curve that calculates according to computing module.

Images