CN102798896B - A kind of logging signal synthesis disposal route of array induction logging instrument and system thereof - Google Patents

Abstract

The invention discloses a kind of logging signal synthetic method of array induction logging instrument, by obtaining logging instrumentation signal data pair array measurement of inductance signal carries out borehole environment correction, and by the Signal transmissions after this correction to the true resolution synthesis unit in signal syntheses processing module and effective background conductance rate computing unit; Described system comprises array induction instrument, hole correction module and signal syntheses processing module, wherein signal syntheses processing module comprises true resolution synthesis unit and resolution match unit, and be respectively equipped with true resolution composite filter storehouse and resolution match filter bank, true resolution synthesis unit and resolution match unit are provided with effective background conductance rate computing unit jointly, and effective background conductance rate computing unit comprises effective background conductance rate calculating filter storehouse; By building each filter bank, overcome the defect of prior art median filter difficult design, the easy measurement by magnification noise of result, and without the need to carrying out SEC in advance.

Description

A kind of logging signal synthesis disposal route of array induction logging instrument and system thereof

Technical field

The invention belongs to Geophysical Logging, relating to a kind of logging signal synthesis disposal route and system thereof of the array induction logging instrument for detecting uncased hole electrical resistivity of earth layer around.

Background technology

Array induction logging instrument is a kind of surveying instrument measuring uncased hole electrical resistivity of earth layer around, by down-hole, it arranges that multiple coil system carrys out the circumocular resistivity information of production wells, then ground is passed in measuring-signal digitizing, the formation information needed for user is extracted by software process measuring-signal in ground, and signal syntheses process is one of core of array induction logging instrument.The object of signal syntheses process eliminates the unnecessary environmental impact such as well, intrusion, country rock and skin effect impact, extracts the formation information needed for user.Have three kinds of array induction logging instruments both at home and abroad at present, be the HRAI of HDIL and the Halliburton company of AIT, BakerAtlas company of Schlumberger company respectively, different instrument adopts different signal syntheses disposal routes.

The coil system of array induction logging instrument AIT is launched 8 group of subarrays of arranging with both sides by one and is become, each subarray comprises a transmitting coil, a potted coil and a main receiving coil, transmitting coil launches 3 kinds of frequency of operation, 8 subarrays receive 1 or 2 frequency signals selectively by bandpass filtering, be specially subarray 1 and 2 and receive only high-frequency signal, subarray 3 and 4 receives high frequency and intermediate-freuqncy signal, subarray 5,6,7 and 8 receives intermediate frequency and low frequency signal, and real part and imaginary part have 28 barss altogether.Measuring-signal, through hole correction, SEC and signal syntheses process, provides 3 groups of resolution (1ft, 2ft, 4ft), 5 kinds of investigation depths (10inch, 20inch, 30inch, 60inch, 90inch) totally 15 resistivity curves.

The coil system of high resolution induction logging tool device HDIL is become with 7 group of subarrays of single-sided arrangement by a transmitting, and each subarray structure is consistent with AIT, is all a transmitting coil, a potted coil and a main receiving coil.Transmitting coil launches the composite signal of 8 kinds of frequencies, and 7 subarrays all receive 8 frequency signals, and real part and imaginary part be totally 112 barss.Real part measuring-signal is first by multi-frequency SEC, then borehole environment correction is carried out, the synthesis of true resolution and resolution match process afterwards, result provides 3 groups of resolution (1ft, 2ft, 4ft), 6 kinds of investigation depths (10inch, 20inch, 30inch, 60inch, 90inch, 120inch) totally 18 curves, comprises middle 6 true resolution resultant curves.

The coil system of high definition induction log logging instrumentation HRAI is launched 10 group of subarrays of arranging with both sides by one and is become.Two short subarrays are asymmetric, all the other 4 subarrays full symmetrics.Transmitting coil launches two kinds of frequency 8kHz and 32kHz, and 10 subarrays receive 2 frequency signals simultaneously, real part and imaginary part totally 40.Measuring-signal is first through bifrequency SEC, then be the hole correction based on geometric factor, be finally the synthesis of true resolution and resolution match, produce 3 groups of resolution (1ft, 2ft, 4ft), 6 kinds of investigation depths (10inch, 20inch, 30inch, 60inch, 90inch, 120inch) totally 18 curves.

Realizing in process of the present invention, inventor finds that prior art at least exists following shortcoming:

For the AIT instrument of Schlumberger company, SEC wherein adopts the pointwise correction method of use in traditional double induction or utilizes imaginary to correct.For last method, can not discriminating measurement signal be due to country rock and intrusion effect lower than true value, or skin effect impact.In practical logging, imaginary signals influence factor is complicated, measures unstable, because after this method can not use during middle high resistivity.In addition, original 28 measuring-signals (or 14 real parts) are directly treated to the curve of 3 groups of resolution, 5 kinds of investigation depths by signal syntheses process, filter design difficult, the easy measurement by magnification noise of result.

For the HDIL instrument of BakerAtlas company, the signal processing method due to this instrument utilizes the approximately linear characteristic of measuring-signal and frequency to efficiently solve skin effect and the inseparable problem of environmental impact, and result is stablized.Limitation is that all subarrays of requirement receive all frequency signals.

For the HRAI instrument of Halliburton company, signal processing method due to this instrument make use of the thought of HDIL, requires that all subarrays receive 2 frequency signals simultaneously, when measuring-signal has error, only use 2 frequency signals, SEC is incorrect.

Summary of the invention

In order to realize the two dimension impact eliminating skin effect and environment, overcome the defect of prior art median filter difficult design, the easy measurement by magnification noise of result, embodiments provide a kind of logging signal synthesis disposal route of array induction logging instrument, described technical scheme is as follows: the acquisition being carried out well-log information by array induction logging instrument, obtain apparatus measures signal, then borehole environment correction process is carried out to surveyed array induction measuring-signal, and Signal transmissions to the signal synthesis unit after being corrected by borehole environment carries out the synthesis process of signal

Described signal syntheses disposal route comprises the steps:

Step one: build effective background conductance rate calculating filter storehouse, true resolution composite filter storehouse and resolution match filter bank respectively based on some background conductance rates;

Step 2: the calculating according to constructed effective background conductance rate calculating filter storehouse, all measuring-signals after hole correction being carried out to background conductance rate, complete the calculating of all signal background conductivity of given depth, and all background conductance rates of given depth are weighted on average, obtain effective background conductance rate σ _aevb;

Step 3: according to the effective background conductance rate σ obtained in step 2 _aevbcalculate with true resolution composite filter storehouse constructed in step one and generate new composite filter, and the measuring-signal convolution after new composite filter and borehole environment are corrected in certain limit, obtain the true resolution resultant curve of all investigation depths;

Step 4: according to the effective background conductance rate σ obtained in step 2 _aevbcalculate with resolution match filter bank constructed in step one and generate new resolution match wave filter, new resolution match wave filter is carried out convolution from the resultant curve of the different investigation depths obtained in step 3, calculate the resolution difference information between adjacent detector depth curve, according to determined coupling aim curve resolution, depth curve is surveyed to institute's probing and carry out information compensation process, thus complete the resolution synthesis process that depth curve is surveyed in institute's probing.

The calculating of resolution match wave filter new in the calculating generation of composite filter new in described step 3 and step 4 generates and all adopts method of interpolation to calculate.

The structure in the effective background conductance rate calculating filter storehouse in described step one is realized by following formula $t_{\mathrm{VGF}}\left(z\right)=\frac{1.65}{w_z\sqrt{2\mathrm{\π}}}{e}^{-\frac{1}{2}{\left[\frac{1.65\×z}{w_z}\right]}^2}$

In formula: parameter w _zfor the width that target Gaussian function vertical response 90% main lobe signal is corresponding;

Wherein w _zdetermining step is as follows:

1. longitudinal infinitesimal geometry factor of selected subarray in segmentation place of background conductance rate is calculated, and normalized, contribute corresponding width as the initial parameter of Gaussian function with longitudinal infinitesimal geometry factor main lobe 90%

2. by the measuring-signal convolution on selected Gaussian function and stratum, background conductance rate must be looked;

3. according to the log response storehouse calculated in uniform dielectric, by linear interpolation, counter pushing away looks true background conductance rate corresponding to background conductance rate;

4. the anti-true background conductance rate pushed away compared with given background conductance rate, if both are consistent, then selected Gaussian function is the wave filter of selected background conductance rate, otherwise adjustment gaussian function is recalculated, until obtain the wave filter of selected background conductance rate from the 1. step.

5. repeat above-mentioned steps, complete the design of effective background conductance rate filter bank of all background conductance rates and all subarrays.

First determine true resolution synthesis aim curve investigation depth and the segmentation of background conductance rate when true resolution composite filter storehouse builds in described step one, then realized the structure of true resolution composite filter by following formula,

$t^k(\mathrm{\ρ},z)=\underset{j=1}{\overset{J}{\mathrm{\Σ}}}\underset{z^{\′}=z_{\min }}{\overset{z_{\max }}{\mathrm{\Σ}}}{w}_j^k(z^{\′},{\mathrm{\σ}}_b)g_j(\mathrm{\ρ},z-z^{\′},{\mathrm{\σ}}_b)$

In formula: t ^k(ρ, z) synthesizes the two dimension target function of kth kind investigation depth for focusing on, for the true resolution composite filter of jth subarray during synthesis kth kind investigation depth, g _j(ρ, z-z ', σ _b) be the two-dimentional Born geometric factor of a jth subarray; J=1 ~ J is synthesis subarray number used, z '=z _min~ z _maxfor focusing on the window ranges of synthesis process;

By above formula discretize, least square optimization method is utilized to realize the design in true resolution composite filter storehouse, the t after discrete ^k(ρ, z) is represented by following expression formula:

$t^k(\mathrm{\ρ},z)=t_{\mathrm{VGF}}^k\left(z\right)t_{\mathrm{RGF}}^k\left(\mathrm{\ρ}\right)$

$t_{\mathrm{VGF}}^k\left(z\right)=\frac{2\sqrt{\ln 2}}{w_z\sqrt{\mathrm{\π}}}{e}^{-4\ln 2{\left(\frac{z-c_z}{w_z}\right)}^2}$

$t_{\mathrm{RGF}}^k\left(\mathrm{\ρ}\right)=\frac{\ln w_{\mathrm{\ρ}}}{\sqrt{\mathrm{\π}}}{e}^{-{\left(\frac{\ln \mathrm{\ρ}-\ln c_{\mathrm{\ρ}}}{\ln w_{\mathrm{\ρ}}}\right)}^2}$

In formula, w _zand c _zfull width at half maximum width and the center of longitudinal objective function respectively; w _ρand c _ρfor the shape parameters of radial object function, the requirement according to different investigation depths and radial shape is determined;

Concrete steps are as follows:

1) two-dimentional object function parameter is determined, by adjustment w _zmake the resolution of composite function equal the resolution of object function, then determine radial object function parameter w according to the requirement of institute's investigation depth resolution _ρand c _ρ;

2) two-dimentional object function and the cross correlation function measuring sub-array responds function is calculated for the discrete one dimension of above-mentioned generation, two-dimentional response function, and according to allowing the length of degree of correlation determination wave filter;

3) alignment error transmission controe parameter and optimal control parameter;

4) adopt least square method optimization to obtain one group and focus on composite filter coefficient, if not required optimum focusing composite filter, repeat the 3rd step, until find out optimum focusing composite filter;

5) the true resolution composite filter storehouse that step 1-4 obtains the segmentation of all background conductance rates is repeated.

Wherein, the adjustment of described error propagation controling parameters is zero for making the coefficient attenuation of wave filter both sides, and it is realized by following expression formula, and this is added in two-dimentional object function,

$s\left(z\right)=\frac{1}{\mathrm{\α}+y\left(z\right)}-\frac{1}{\mathrm{\α}+y\left(0\right)}$

$y\left(z\right)=e^{-\frac{1}{2}{\left(\frac{z-z_0}{\mathrm{\β}}\right)}^2}$

Wherein function y (z) is error distribution function; S (z) is error complementary function; β is the controling parameters of soft-constraint width; α is the precision controlling of wave filter constraint both sides; Force the both sides value of wave filter close to zero by regulating parameter α and β.

First determine the segmentation of background conductance rate when described step one intermediate-resolution matched filter storehouse builds, then realized the structure of resolution match filter bank (9) by following formula,

g _pvgf(i+1)(z)＝h _pvgfi(z)*g _pvgfi(z)

In formula: g _{pvgf (i+1)}z () is for focusing on longitudinal Born geometric factor of synthesis low resolution curve, h _pvgfifor matched filter, g _pvgfiz () synthesizes longitudinal Born geometric factor of high resolving power curve for focusing on, i=1, and 2 ..., n-1, n are true resolution resultant curve number;

By above formula discretize, least square optimization method is then utilized to realize the design of resolution match filter bank.

Information compensation process in described step 4 is by new matched filter and high resolving power curve convolution, form adjacent low resolution curve, then deduct the low resolution curve of generation with high resolving power curve, obtain the difference information of adjacent resolution curve, its mapping mode comprises:

Low resolution curvilinear transformation is high resolving power curve, is added to the curve that low resolution curve obtains by difference information;

High resolving power curvilinear transformation is low resolution curve, is added to the curve that high resolving power curve obtains by difference information.

On the other hand, provide a kind of logging signal synthesis disposal system of array induction logging instrument, described system comprises:

Apparatus measures device, comprises array induction instrument and instrument truck, for providing apparatus measures signal;

Hole correction module, the measuring-signal for pair array induction instrument aratus carries out borehole environment correction, eliminates the impact of borehole environment;

Signal syntheses processing module, the measuring-signal after described hole correction resume module transfers to the synthesis process that described signal syntheses processing module carries out signal, it is characterized in that:

Described signal syntheses processing module comprises true resolution synthesis unit and resolution match unit; Described true resolution synthesis unit is provided with true resolution composite filter storehouse; Described resolution match unit is provided with resolution match filter bank; True resolution synthesis unit and resolution match unit are provided with effective background conductance rate computing unit jointly, and wherein, described effective background conductance rate computing unit comprises effective background conductance rate calculating filter storehouse;

The measuring-signal of correction is transferred to true resolution synthesis unit and effective background conductance rate computing unit by described hole correction module respectively, the structure value in the background conductance rate score after the measuring-signal that hole correction module corrects by described true resolution synthesis unit and described effective background conductance rate computing unit calculate and true resolution composite filter storehouse carries out true resolution synthesis, Signal transmissions after synthesis is to described resolution match unit, background conductance rate score after effective background conductance rate computing unit calculates by resolution match unit and the structure value of resolution match filter bank carry out resolution match, thus complete the logging signal synthesis process of array induction logging instrument.

Further, described effective background conductance rate computing unit and described true resolution composite filter storehouse build new composite filter, and new composite filter and measuring-signal after described hole correction module corrects carry out true resolution and synthesize and process.

Further, described effective background conductance rate computing unit and described resolution match filter bank build new matched filter, and new matched filter carries out the coupling of resolution to the composite signal after the process of described true resolution synthesis unit.

The true resolution synthesis of described array induction apparatus measures signal realizes according to the following equation:

${\mathrm{\σ}}_p({\mathrm{\ρ}}_k,z)=\underset{j=1}{\overset{J}{\mathrm{\Σ}}}\underset{z^{\′}=z_{\min }}{\overset{z_{\max }}{\mathrm{\Σ}}}{w}_j^k(z^{\′},\mathrm{\σ}){\mathrm{\σ}}_a^j(z-z^{\′})$

Wherein, σ _p(ρ _k, z) be radial depth of investigetion ρ after synthesis _k, the conductivity at longitudinal degree of depth z place; focus on to be synthesized to investigation depth ρ _ktime, the true resolution composite filter of a jth subarray, for a jth subarray is at the measured value of z point through hole correction; K=1,2 ..., K, K are synthesis investigation depth numbers, j=1,2 ..., J, J are synthesis subarray number used; Z=z _min~ z _maxfor focusing on the window ranges of synthesis process.

Wherein, coil in described array induction instrument is made up of multiple subarray, each described subarray comprises a transmitting coil, a main receiving coil and 1-2 potted coil, subarray wherein adopts single-sided arrangement and/or both sides arrangement, and the quantity of frequency signal that each subarray receives is less than the quantity of the frequency signal that described transmitting coil is launched.

The beneficial effect of the technical scheme that the embodiment of the present invention provides is:

The present invention constructs effective background conductance rate calculating filter storehouse, and proposes a kind of method for designing of true resolution composite filter storehouse, this system overcomes the defect of prior art median filter difficult design, the easy measurement by magnification noise of result; And the subarray of the method is measured when not seeking common ground and is received all frequency signals, does not require to carry out SEC in advance.

Array induction logging signal is after process of the present invention, and eliminate skin effect and environment (comprising intrusion and country rock) two dimension impact, result provides 5 kinds or 6 kinds of investigation depth curves of 3 groups of resolution (0.3m, 0.6m and 1.2m); Provide middle true resolution resultant curve simultaneously.

Accompanying drawing explanation

In order to be illustrated more clearly in the technical scheme in the embodiment of the present invention, below the accompanying drawing used required in describing embodiment is briefly described, apparently, accompanying drawing in the following describes is only 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.

Fig. 1 is array induction signal syntheses disposal system schematic diagram;

Fig. 2 is array induction signal syntheses processing flow chart;

Fig. 3 is a kind of two-dimentional response function of array induction logging;

Fig. 4 is one group of true resolution composite filter waveform schematic diagram;

Fig. 5 is that one group of true resolution synthesizes two-dimentional response function;

Fig. 6 is one group of resolution matched filter waveform schematic diagram;

Fig. 7 is the two-dimentional response function after one group of resolution match;

Fig. 8 a is Oklahoma stratum measurement data;

Fig. 8 b is the true resolution synthesis result of Oklahoma stratum measurement data;

Fig. 8 c is the 0.6m resolution match result of Oklahoma stratum measurement data.

Embodiment

For making the object, technical solutions and advantages of the present invention clearly, below in conjunction with accompanying drawing, embodiment of the present invention is described further in detail.

Measuring-signal below by way of array induction logging instrument AIT structure of coil system is that example is described.

As shown in Figure 1, array induction signal syntheses disposal system as shown in the figure, present system forms primarily of three links: apparatus measures device 1, hole correction module 2 and signal syntheses processing module 3, the array induction measuring-signal obtained of being logged well by apparatus measures device 1 is delivered to the borehole environment correction that hole correction module 2 carries out measured value, and the Signal transmissions after this correction is carried out the synthesis process of signal to signal syntheses processing module 3, thus complete the synthesis process of array induction apparatus measures signal;

Wherein, apparatus measures device 1 is made up of, for providing apparatus measures signal array induction instrument 11 and surface seismic records instrument car 12;

Wherein, hole correction module 2 carries out borehole environment correction for the measuring-signal of array induction instrument, eliminates the impact of borehole environment;

Wherein, signal syntheses processing module 3 comprises true resolution synthesis unit 4 and resolution match unit 5; And true resolution synthesis unit 4 is provided with and true resolution composite filter storehouse 8; Resolution match unit 5 is provided with resolution match filter bank 9; True resolution synthesis unit 4 and resolution match unit 5 are provided with effective background conductance rate computing unit 6 jointly, wherein, effective background conductance rate computing unit 6 comprises effective background conductance rate calculating filter storehouse 7, for carrying out the calculating of effective background conductance rate to measuring-signal; The measuring-signal of correction is transferred to true resolution synthesis unit 4 and effective background conductance rate computing unit 6 by hole correction module 2 respectively; The structure value in the background conductance rate score after the measuring-signal that hole correction module 2 corrects by true resolution synthesis unit 4 and effective background conductance rate computing unit 6 calculate and true resolution composite filter storehouse 8 carries out true resolution synthesis; After this synthesis, Signal transmissions is to resolution match unit 5, background conductance rate score after effective background conductance rate computing unit 6 calculates by resolution match unit 5 and the structure value of resolution match filter bank 9 carry out resolution match, thus complete the logging signal synthesis process of array induction logging instrument.

In the present invention further, effective background conductance rate computing unit 6 and true resolution composite filter storehouse 8 build new composite filter, and carry out true resolution and synthesize with the measuring-signal after correcting through borehole environment correction module 2 and process.

Further, effective background conductance rate computing unit 6 builds new matched filter with resolution match filter bank 9, the composite signal after the process of synthesis processing unit 4 is carried out to the coupling of resolution.

Further, effective background conductance rate calculating filter storehouse 7 is the filter set set up based on some background conductance rates; Defining method wherein for specific background conductivity and its wave filter of measuring-signal realizes according to the following equation:

$t_{\mathrm{VGF}}\left(z\right)=\frac{1.65}{w_z\sqrt{2\mathrm{\π}}}{e}^{-\frac{1}{2}{\left[\frac{1.65\×z}{w_z}\right]}^2}$ ①

In formula: parameter w _zfor the width that target Gaussian function vertical response 90% main lobe signal is corresponding, select initial value according to different background conductivity and measuring-signal.

Further, true resolution composite filter storehouse 8 is the filter set set up based on some background conductance rates; Defining method wherein for specific background conductivity wave filter realizes according to the following equation:

$t^k(\mathrm{\ρ},z)=\underset{j=1}{\overset{J}{\mathrm{\Σ}}}\underset{z^{\′}=z_{\min }}{\overset{z_{\max }}{\mathrm{\Σ}}}{w}_j^k(z^{\′},{\mathrm{\σ}}_b)g_j(\mathrm{\ρ},z-z^{\′},{\mathrm{\σ}}_b)$ ………….②

In formula: t ^k(ρ, z) synthesizes the two dimension target function of kth kind investigation depth for focusing on, for the true resolution composite filter of jth subarray during synthesis kth kind investigation depth, g _j(ρ, z-z ', σ _b) be the two-dimentional Born geometric factor of a jth subarray; J=1 ~ J is synthesis subarray number used, z '=z _min~ z _maxfor focusing on the window ranges of synthesis process;

By above formula discretize, then realize the design in true resolution composite filter storehouse by least square optimization method;

T ^k(ρ, z) is represented by following expression formula:

$t^k(\mathrm{\ρ},z)=t_{\mathrm{VGF}}^k\left(z\right)t_{\mathrm{RGF}}^k\left(\mathrm{\ρ}\right)...\left(2.1\right)$

$t_{\mathrm{VGF}}^k\left(z\right)=\frac{2\sqrt{\ln 2}}{w_z\sqrt{\mathrm{\π}}}{e}^{-4\ln 2{\left(\frac{z-c_z}{w_z}\right)}^2}...\left(2.2\right)$

$t_{\mathrm{RGF}}^k\left(\mathrm{\ρ}\right)=\frac{\ln w_{\mathrm{\ρ}}}{\sqrt{\mathrm{\π}}}{e}^{-{\left(\frac{\ln \mathrm{\ρ}-\ln c_{\mathrm{\ρ}}}{\ln w_{\mathrm{\ρ}}}\right)}^2}...\left(2.3\right)$

In formula, w _zand c _zfull width at half maximum width and the center of longitudinal objective function respectively; w _ρand c _ρfor the shape parameters of radial objective function, the requirement according to different investigation depths and radial shape is determined.

Further, resolution match filter bank 9 is the filter set set up based on some background conductance rates; Matched filter wherein for specific background conductivity realizes according to the following equation:

g _pvgf(i+1)(z)＝h _pvgfi(z)*g _pvgfi(z)③

In formula: g _{pvgf (i+1)}z () is for focusing on longitudinal Born geometric factor of synthesis low resolution curve, h _pvgfifor matched filter, g _pvgfiz () synthesizes longitudinal Born geometric factor of high resolving power curve for focusing on, i=1, and 2 ..., n-1, n are true resolution resultant curve number;

By above formula discretize, then realize the design of resolution match filter bank by least square optimization method.

In the present invention, the true resolution synthesis process of array induction apparatus measures signal realizes according to the following equation:

${\mathrm{\σ}}_p({\mathrm{\ρ}}_k,z)=\underset{j=1}{\overset{J}{\mathrm{\Σ}}}\underset{z^{\′}=z_{\min }}{\overset{z_{\max }}{\mathrm{\Σ}}}{w}_j^k(z^{\′},\mathrm{\σ}){\mathrm{\σ}}_a^j(z-z^{\′})$ ④

Wherein, σ _p(ρ _k, z) be radial depth of investigetion ρ after synthesis _k, the conductivity at longitudinal degree of depth z place; focus on to be synthesized to investigation depth ρ _ktime, the true resolution composite filter of a jth subarray, for a jth subarray is at the measured value of z point through hole correction; K=1,2 ..., K, K are synthesis investigation depth numbers, and j=1 ~ J is synthesis subarray number used; Z '=z _min~ z _maxfor focusing on the window ranges of synthesis process.

In the present invention, the resolution match of array induction apparatus measures signal realizes in the following manner:

Information compensation principle is adopted to carry out the resolution match process of true resolution resultant curve; First, effective background conductance rate computing unit 6 builds new matched filter with resolution match filter bank 9, by new matched filter and high resolving power curve convolution, form adjacent low resolution curve, then deduct the low resolution curve of generation with high resolving power curve, obtain the difference information of adjacent resolution curve; If be high resolving power curve by low resolution curvilinear transformation, then difference information is added on low resolution curve; Otherwise, if wanting high resolving power curvilinear transformation is low resolution curve, only difference information need be deducted from high resolving power curve.

As shown in Figure 2, the array induction provided below by the present invention synthesizes disposal route to describe embodiments of systems of the invention in detail.

This array induction instrument signal synthesis disposal route comprises the steps:

One, effective background conductance rate calculating filter storehouse 7, true resolution composite filter storehouse 8 and resolution match filter bank 9 is first built:

1, build effective background conductance rate calculating filter storehouse 7, realized by following method:

Described effective background conductance rate calculating filter storehouse 7 is the filter set set up based on some background conductance rates; Defining method wherein for specific background conductivity and its wave filter of measuring-signal realizes according to the following equation: $t_{\mathrm{VGF}}\left(z\right)=\frac{1.65}{w_z\sqrt{2\mathrm{\π}}}{e}^{-\frac{1}{2}{\left[\frac{1.65\×z}{w_z}\right]}^2}$ ①

In formula: parameter w _zfor the width that target Gaussian function vertical response 90% main lobe signal is corresponding, according to different conductivity and He Ne laser initial value.

(1) segmentation of background conductance rate is determined.According to measurement range and the skin effect non-linear effects situation considering each subarray of formation conductivity, conductivity is divided into some sections.In the present invention, formation conductivity scope is 0-10S/m, and conductivity is divided into 13 sections: 0.001,0.01,0.02,0.05,0.1,0.2,0.5,1.0,2.0,3.5,5.0,7.0 and 10.0S/m.

(2) target Gaussian function is determined.

(3) for selected background conductance rate and measurement subarray, parameter w is determined _znumerical value.W _zdetermining step is as follows:

1. longitudinal infinitesimal geometry factor of selected subarray in segmentation place of background conductance rate is calculated, and normalized, contribute corresponding width as the initial parameter of Gaussian function with longitudinal infinitesimal geometry factor main lobe 90%

2. by the measuring-signal convolution on Gaussian function and stratum, its result is looks background conductance rate;

3. according to the log response storehouse calculated in uniform dielectric, by linear interpolation, counter pushing away looks true background conductance rate corresponding to background conductance rate;

4. the anti-true background conductance rate pushed away compared with given background conductance rate, if both are consistent, then corresponding Gaussian function is the wave filter of selected background conductance rate, otherwise adjustment gaussian function is recalculated from (1) step.

(4) repeat step (3), complete the design of all background conductance rates and all subarrays effective background conductance rate filter bank.

Here uniform dielectric refers to electrically, the completely uniform Single Medium of physical property, refers to without well without country rock, be a kind of lithology in infinitely great scope in well logging.

2, build true resolution composite filter storehouse 8, realized by following method:

True resolution composite filter storehouse of the present invention is the filter set set up based on some background conductance rates.Defining method wherein for specific background conductivity wave filter realizes according to the following equation:

$t^k(\mathrm{\ρ},z)=\underset{j=1}{\overset{J}{\mathrm{\Σ}}}\underset{z^{\′}=z_{\min }}{\overset{z_{\max }}{\mathrm{\Σ}}}{w}_j^k(z^{\′},{\mathrm{\σ}}_b)g_j(\mathrm{\ρ},z-z^{\′},{\mathrm{\σ}}_b)$ ……….②

In formula: t ^k(ρ, z) synthesizes the two dimension target function of kth kind investigation depth for focusing on, for the true resolution composite filter of jth subarray during synthesis kth kind investigation depth, g _j(ρ, z-z ', σ _b) be the two-dimentional Born geometric factor of a jth subarray; J=1 ~ J is synthesis subarray number used, z '=z _min~ z _maxfor focusing on the window ranges of synthesis process;

By above formula discretize, then realize the design in true resolution composite filter storehouse by least square optimization method;

T ^k(ρ, z) is represented by following expression formula:

$t^k(\mathrm{\ρ},z)=t_{\mathrm{VGF}}^k\left(z\right)t_{\mathrm{RGF}}^k\left(\mathrm{\ρ}\right)...\left(2.1\right)$

$t_{\mathrm{VGF}}^k\left(z\right)=\frac{2\sqrt{\ln 2}}{w_z\sqrt{\mathrm{\π}}}{e}^{-4\ln 2{\left(\frac{z-c_z}{w_z}\right)}^2}...\left(2.2\right)$

$t_{\mathrm{RGF}}^k\left(\mathrm{\ρ}\right)=\frac{\ln w_{\mathrm{\ρ}}}{\sqrt{\mathrm{\π}}}{e}^{-{\left(\frac{\ln \mathrm{\ρ}-\ln c_{\mathrm{\ρ}}}{\ln w_{\mathrm{\ρ}}}\right)}^2}...\left(2.3\right)$

In formula, w _zand c _zfull width at half maximum width and the center of longitudinal objective function respectively.W _ρand c _ρfor the shape parameters of radial objective function, the requirement according to different investigation depths and radial shape is determined.

Concrete implementation step is:

(1) true resolution synthesis aim curve investigation depth is determined.

(2) segmentation of background conductance rate is determined.According to measurement range and the skin effect non-linear effects situation considering each subarray of formation conductivity, conductivity is divided into some sections.In the present invention, formation conductivity scope is 0-10S/m, and background conductance rate is divided into 16 sections: 0.001,0.01,0.02,0.05,0.1,0.2,0.4,0.5,0.6,0.8,1.0,2.0,3.5,5.0,7.0 and 10.0S/m.

(3) according to the target detection degree of depth and segmentation background conductance rate, true resolution composite filter storehouse is designed:

1. two-dimentional object function parameter is determined:

First the controling parameters w of longitudinal object function is determined _zand c _z.By adjustment w _zthe resolution of composite function is made to equal the resolution of object function; c _zget 0.In example of the present invention, determine that corresponding resolution is 0.3,0.6,1.2,1.5 and 1.8m.

According to the investigation depth requirement of a certain resolution, determine radial object function parameter w _ρand c _ρ.Wherein w _ρand c _ρrequirement according to the different depth of focus and radial shape is determined.Various investigation depth 0.25,0.50,0.75,1.50 and parameter (c corresponding to 2.25m is determined in the present invention _ρ, w _ρ) be respectively (0.261,0.404), (0.515,0.97), (0.76,1.4) and (2.1,4.2);

2. discrete one dimension, two-dimentional response function is 2. produced for formula;

3. two-dimentional object function and the cross correlation function measuring sub-array responds function is calculated, and according to allowing the length of degree of correlation determination wave filter;

4. alignment error transmission controe parameter and optimal control parameter;

5. adopt least square method optimization to obtain one group and focus on composite filter coefficient, if not required optimum focusing composite filter, repeat the 4th step, until find out optimum focusing composite filter.

(4) the true resolution composite filter storehouse that step (3) obtains the segmentation of all background conductance rates is repeated.

In optimized design, because each subarray two-dimensional space response function exists a large amount of redundant informations near wellbore; The external signal of launching and receiving coil is decayed all slowly.The redundancy of information causes the existence having multiple optimum solution under given object function and constraint condition, shows as filter coefficient both sides and effectively can not decay to zero.Introduce a kind of soft-constraint control function, soft-constraint is added in two-dimentional object function as new regulating error item for this reason.Soft-constraint function is

$s\left(z\right)=\frac{1}{\mathrm{\α}+y\left(z\right)}-\frac{1}{\mathrm{\α}+y\left(0\right)}$

$y\left(z\right)=e^{-\frac{1}{2}{\left(\frac{z-z_0}{\mathrm{\β}}\right)}^2}$

Wherein function y is the function being similar to the distribution of Gaussian error; S also claims error complementary function; β is the controling parameters of soft-constraint width; α is the precision controlling of wave filter constraint both sides.Soft-constraint is equivalent to the penalty of a change in optimization, in main region, punishes little, and in other region, punishment is large, thus forces the both sides value of wave filter close to zero.When optimized design focusing filter, undesirable optimum solution can be eliminated by regulating parameter α and β.

3, build resolution match filter bank 9, realized by following method:

Resolution match filter bank of the present invention is the filter set set up based on some background conductance rates.Matched filter wherein for specific background conductivity realizes according to the following equation:

g _pvgf(i+1)(z)＝h _pvgfi(z)*g _pvgfi(z)…………………③

In formula: g _{pvgf (i+1)}z () is for focusing on longitudinal Born geometric factor of synthesis low resolution curve, h _pvgfifor matched filter, g _pvgfiz () synthesizes longitudinal Born geometric factor of high resolving power curve for focusing on, i=1, and 2 ..., n-1, n are true resolution resultant curve number.

By above formula discretize, then realize the design of resolution match filter bank by least square optimization method;

Concrete implementation step is:

(1) segmentation of background conductance rate is determined.According to measurement range and the skin effect non-linear effects situation considering each subarray of formation conductivity, conductivity is divided into some sections.In the present invention, formation conductivity scope is 0-10S/m, and background conductance rate is divided into 16 sections: 0.001,0.01,0.02,0.05,0.1,0.2,0.4,0.5,0.6,0.8,1.0,2.0,3.5,5.0,7.0 and 10.0S/m.

(2) selected adjacent two true resolution resultant curves, according to selected background conductance rate, calculating resolution matched filter.According to longitudinal Born geometric factor of selected curve, 5. formula is separated into matrix form, and determines the length of wave filter.Then least square method optimization is adopted to obtain one group of resolution matched filter.Determination filter length of the present invention is 30 meters.

(3) repeat step (2), obtain the lower 4 kinds of resolution match wave filters of all background conductance rate segmentation situations.Two, the process of array induction apparatus measures signal syntheses is carried out:

1) carried out the acquisition of well-log information by array induction logging instrument and surface laboratory, obtain apparatus measures signal.

2) hole correction module 2 realizes the correction process of pair array measurement of inductance signal borehole environment.

3) after effective background conductance rate computing unit 6 pairs of hole corrections, measuring-signal realizes the calculating of effective background conductance rate:

1. smothing filtering is done to the measuring-signal after hole correction, as looking background conductance rate estimated value

2. background conductance rate estimated value will be looked compare with uniform dielectric the calculated results, the estimated value of background conductance rate will be converted into depending on background conductance rate estimated value

3. for the estimated value of background conductance rate in conjunction with the set of background conductance rate and effective background conductance rate calculating filter storehouse 7, calculate and generate new background conductance rate calculating filter;

For given depth, if meet then correspond to wave filter can be tried to achieve by following interpolation formula

$w_{\mathrm{ba}}^j=a_2{w}_{\mathrm{bai}}^j+a_1{w}_{\mathrm{ba}(i+1)}^j...\left(1.1\right)$

Wherein a ₁and a ₂by background conductance rate with the weights determined, namely

$a_1=\frac{\ln {\mathrm{\σ}}_{b0}^j-\ln {\mathrm{\σ}}_{\mathrm{bi}}^j}{\ln {\mathrm{\σ}}_{b(i+1)}^j-\ln {\mathrm{\σ}}_{\mathrm{bi}}^j}...\left(1.2\right)$

$a_2=\frac{\ln {\mathrm{\σ}}_{b(i+1)}^j-\ln {\mathrm{\σ}}_{b0}^j}{\ln {\mathrm{\σ}}_{b(i+1)}^j-\ln {\mathrm{\σ}}_{\mathrm{bi}}^j}...(1.3)$

4. the measurement signal value convolution after hole correction in the wave filter 3. calculated and certain limit is obtained the background conductance rate of given depth;

5. to all measuring-signals after hole correction repeat steps 1. ~ 4., complete the calculating of all signal background conductivity;

6. all background conductance rates that 5. upper step obtains are weighted on average, obtain unique effective background conductance rate σ _aevb.

8 sub-array signal in example of the present invention, when background conductance rate is less than 1S/m, with the effective background conductance rate σ of subarray 6,7 and 8 weighted average calculation _aevb, namely

${\mathrm{\σ}}_{\mathrm{aveb}}\left(z\right)=d_9{\mathrm{\σ}}_b^9+d_{10}{\mathrm{\σ}}_b^{10}+d_{11}{\mathrm{\σ}}_b^{11}+d_{12}{\mathrm{\σ}}_b^{12}+d_{13}{\mathrm{\σ}}_b^{13}+d_{14}{\mathrm{\σ}}_b^{14}$

In formula, j=9 ..., the 14 low and medium frequency signals being respectively subarray 6,7 and 8; d _j: j=9 ..., 14 is weighting factor, meets normalizing condition

$\underset{j=9}{\overset{14}{\mathrm{\Σ}}}{d}_j=1$

4) after true resolution synthesis unit 4 pairs of hole corrections, measuring-signal realizes the synthesis process of true resolution:

Self-adaptation nonlinear of the present invention true resolution synthesis process, can realize normalization SEC and true resolution and synthesize and synchronously complete.

1. true resolution synthesis aim curve investigation depth is determined;

In example of the present invention, true resolution synthesis aim curve investigation depth determines 5 kinds, is respectively 0.25m, 0.5m, 0.75m, 1.50m and 2.25m.

2. for step 3) effective background conductance rate σ of calculating _aevb, according to fixed true resolution composite filter storehouse 8, calculate and produce new composite filter.

Judge that effective background conductance rate belongs to the piecewise interval of background conductance rate, determine with look-up interpolations the true resolution composite filter being with skin effect; Work as σ _avebmeet inequality σ _bl≤ σ _aveb< σ _{b (l+1)}time, then σ _avebcorresponding new wave filter is:

$w_j^k\left(i\right)=c_2{w}_{\mathrm{jbl}}^k\left(i\right)+c_1{w}_{\mathrm{jb}(l+1)}^k\left(i\right)...\left(2.1\right)$

In formula, the true resolution composite filter coefficient of a jth subarray at degree of depth i place when being synthesis kth investigation depth, by the adjacent filter in filter bank with weighted calculation obtains; c ₁and c ₂for according to σ _aveb, σ _bland σ _{b (l+1)}the weighting coefficient calculated, namely

$c_1=\frac{\ln {\mathrm{\&sigma;}}_{\mathrm{aveb}}-\ln {\mathrm{\&sigma;}}_{\mathrm{bl}}}{\ln {\mathrm{\&sigma;}}_{b(l+1)}-\ln {\mathrm{\&sigma;}}_{\mathrm{bl}}}...\left(2.2\right)$

c ₂＝1-c ₁...........................(2.3)

3. the measuring-signal convolution after the true resolution composite filter 2. calculated and borehole environment being corrected in certain limit, just obtains true resolution resultant curve;

4. repeat 1. ~ 3., obtain the true resolution resultant curve of all investigation depths.

The present invention gives accompanying drawing for further illustrating its signal syntheses processing procedure.

Fig. 3 to Fig. 5 is respectively array induction subarray two dimension response function, true resolution composite filter and true resolution resultant curve two dimension response function.Wherein:

Fig. 3 is the two-dimentional response function of array induction the 5th subarray.

Fig. 4 is the true resolution composite filter being synthesized to 0.75m investigation depth; Fig. 4 illustrates and fluctuates the true resolution composite filter that the present invention designs smooth, do not have high frequency noise, and main interval fluctuation is large, fluctuates at a distance little, efficiently solves redundant information problem; Original two dimensional response is after true resolution synthesis process, and the smooth symmetry of figure, wellbore effect are little, and information concentrates on certain area near investigation depth.

Fig. 5 is the 0.75m investigation depth function that 5 kinds of true resolution synthesize in two-dimentional response function.

5) resolution match unit 5 realizes resolution match process to true resolution resultant curve:

Adaptive resolution matching treatment of the present invention, its implementation procedure is as follows:

1. for step 3) effective background conductance rate σ of calculating _aevb, according to fixed resolution match filter bank 9, calculate new resolution match wave filter;

Judge effective background conductance rate σ _avebbelong to the piecewise interval of background conductance rate.With method of interpolation determination resolution match wave filter.Work as σ _avebmeet inequality σ _bl≤ σ _aveb< σ _{b (l+1)}time, then σ _avebcorresponding resolution match wave filter is

h _k(i)＝c ₂h _kbl(i)+c ₁h _kb(l+1)(i)................(3.1)

In formula, h _ki () (k=1,2,3,4) are the resolution match wave filter of a kth investigation depth to kth+1, by the adjacent filter h in filter bank _kbl(i) and h _{kb (l+1)}i () weighted calculation obtains; c ₁and c ₂for according to σ _aveb, σ _bland σ _{b (l+1)}the weighting coefficient calculated, namely

$c_1=\frac{\ln {\mathrm{\&sigma;}}_{\mathrm{aveb}}-\ln {\mathrm{\&sigma;}}_{\mathrm{bl}}}{\ln {\mathrm{\&sigma;}}_{b(l+1)}-\ln {\mathrm{\&sigma;}}_{\mathrm{bl}}}...\left(3.2\right)$

c ₂＝1-c ₁...............................(3.3)

2. by new matched filter from through step 4) the different investigation depth resultant curve convolutions that obtain;

3. the resolution difference information between adjacent detector depth curve is calculated;

4. coupling aim curve resolution is determined;

5. require to carry out information compensation process according to resolution match, obtain the curve of 3 kinds of different resolutions.

In example of the present invention, the 0.6m resolution match formula of five kinds of investigation depths is:

${\mathrm{\&sigma;}}_{p1}^{0.6m}\left(n\right)={\mathrm{\&sigma;}}_p^1\left(n\right)-\mathrm{\&Delta;}{\mathrm{\&sigma;}}_p^{\mathrm{1,2}}$

${\mathrm{\&sigma;}}_{p2}^{0.6m}\left(n\right)={\mathrm{\&sigma;}}_p^2\left(n\right)$

${\mathrm{\&sigma;}}_{p3}^{0.6m}\left(n\right)={\mathrm{\&sigma;}}_p^3\left(n\right)+\mathrm{\&Delta;}{\mathrm{\&sigma;}}_p^{\mathrm{2,3}}$

${\mathrm{\&sigma;}}_{p4}^{0.6m}\left(n\right)={\mathrm{\&sigma;}}_p^4\left(n\right)+\mathrm{\&Delta;}{\mathrm{\&sigma;}}_p^{\mathrm{2,3}}+\mathrm{\&Delta;}{\mathrm{\&sigma;}}_p^{\mathrm{3,4}}$

${\mathrm{\&sigma;}}_{p5}^{0.6m}\left(n\right)={\mathrm{\&sigma;}}_p^5\left(n\right)+\mathrm{\&Delta;}{\mathrm{\&sigma;}}_p^{\mathrm{2,3}}+\mathrm{\&Delta;}{\mathrm{\&sigma;}}_p^{\mathrm{3,4}}+\mathrm{\&Delta;}{\mathrm{\&sigma;}}_p^{\mathrm{4,5}}$

In formula: be respectively 5 kinds of investigation depth curves that resolution match is 0.6 meter, be respectively the true resolution resultant curve of 5 different investigation depths, i=1,2,3, the 4 resolution differences being respectively adjacent two kinds of investigation depth curves.

The present invention produces 3 groups of resolution curve of 5 kinds of investigation depths, totally 15.For the array induction logging instrument of 6 kinds of investigation depths, 3 groups of resolution, produce 18 curves.Wherein Fig. 6 is one group of resolution matched filter; Fig. 7 is the two-dimentional response function after 1.50m investigation depth Curve Matching to 0.6m.Fig. 7 display resolution matching treatment improves the resolution of deep investigation curve.

Fig. 8 is to the result of Oklahoma stratum measurement data with numerical simulation the present invention.Fig. 8 a is raw measurement results, and Fig. 8 b is true resolution synthesis result, and Fig. 8 c is 0.6m resolution match result.Fig. 8 shows raw measured signal, and (Fig. 8 a) affects by skin effect and shoulder effect simultaneously, and different subarray measured value is different; After true resolution synthesis (Fig. 8 b), skin effect is corrected, but because different investigation depth resolution is different, zone of interest is still separated; After resolution match (Fig. 8 c), different investigation depth curve co-insides, there is identical resolution, can differentiate 0.6m thin layer (see 150m place), the simulation of this example describes the validity that invention instrument processes various resistivity contrast and thickness in without invaded formation analogue measurement.

The invention described above embodiment sequence number, just to describing, does not represent the quality of embodiment.

The foregoing is only preferred embodiment of the present invention, not in order to limit the present invention, within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (8)

1. the logging signal synthesis disposal route of an array induction logging instrument, the acquisition of well-log information is carried out by array induction logging instrument, obtain apparatus measures signal, then borehole environment correction process is carried out to surveyed array induction measuring-signal, and Signal transmissions to the signal synthesis unit after being corrected by borehole environment carries out the synthesis process of signal

It is characterized in that: described signal syntheses disposal route comprises the steps:

Step one: build effective background conductance rate calculating filter storehouse (7), true resolution composite filter storehouse (8) and resolution match filter bank (9) respectively based on some background conductance rates;

Step 2: the calculating according to constructed effective background conductance rate calculating filter storehouse (7), all measuring-signals after hole correction being carried out to background conductance rate, complete the calculating of all signal background conductivity of given depth, and all background conductance rates of given depth are weighted on average, obtain effective background conductance rate σ _aevb;

Step 3: according to the effective background conductance rate σ obtained in step 2 _aevbnew composite filter is generated with true resolution composite filter storehouse (8) calculating constructed in step one, and the measuring-signal convolution after new composite filter and borehole environment are corrected within the scope of certain depth, obtain the true resolution resultant curve of all investigation depths;

Step 4: according to the effective background conductance rate σ obtained in step 2 _aevbcalculate with resolution match filter bank (9) constructed in step one and generate new resolution match wave filter, new resolution match wave filter is carried out convolution from the resultant curve of the different investigation depths obtained in step 3, calculate the resolution difference information between adjacent detector depth curve, according to determined coupling aim curve resolution, depth curve is surveyed to institute's probing and carry out information compensation process, complete the resolution synthesis process that depth curve is surveyed in institute's probing;

The structure of the effective background conductance rate calculating filter storehouse (7) in described step one is realized by following formula $t_{VGF}\left(z\right)=\frac{1.65}{w_z\sqrt{2\mathrm{\&pi;}}}{e}^{-\frac{1}{2}{\&lsqb;\frac{1.65\&times;z}{w_z}\&rsqb;}^2}$

In formula: t _vGFz () is area normalization Gaussian function corresponding to effective background conductance rate calculating filter storehouse; Z is independent variable; Parameter w _zfor the width that target Gaussian function vertical response 90% main lobe signal is corresponding;

Wherein w _zdetermining step is as follows:

1. longitudinal infinitesimal geometry factor of selected subarray in segmentation place of background conductance rate is calculated, and normalized, contribute corresponding width as the initial parameter of Gaussian function with longitudinal infinitesimal geometry factor main lobe 90%

2. by the measuring-signal convolution on selected Gaussian function and stratum, background conductance rate must be looked;

3. according to the log response storehouse calculated in uniform dielectric, by linear interpolation, counter pushing away looks true background conductance rate corresponding to background conductance rate;

4. the anti-true background conductance rate pushed away compared with given background conductance rate, if both are consistent, then selected Gaussian function is the wave filter of given background conductance rate, otherwise adjustment gaussian function is recalculated, until obtain the wave filter of given background conductance rate from the 1. step;

5. repeat above-mentioned steps, complete the design of effective background conductance rate filter bank of all background conductance rates and all subarrays;

First determine true resolution synthesis aim curve investigation depth and the segmentation of background conductance rate when true resolution composite filter storehouse (8) builds in described step one, then realized the structure of true resolution composite filter by following formula,

$t^k(\mathrm{\&rho;},z)=\underset{j=1}{\overset{J}{\&Sigma;}}\underset{z^{\&prime;}=z_{\min }}{\overset{z_{\max }}{\&Sigma;}}{w}_j^k(z^{\&prime;},{\mathrm{\&sigma;}}_b)g_j(\mathrm{\&rho;},z-z^{\&prime;},{\mathrm{\&sigma;}}_b)$ ①

In formula: ρ is radial depth of investigetion coordinate, z is longitudinal depth coordinate, σ _bfor background conductance rate, t ^k(ρ, z) synthesizes the two dimension target function of kth kind investigation depth for focusing on, for the true resolution composite filter of jth subarray during synthesis kth kind investigation depth, g _j(ρ, z-z ', σ _b) be the two-dimentional Born geometric factor of a jth subarray; J=1 ~ J is synthesis subarray number used, z '=z _min~ z _maxfor focusing on the window ranges of synthesis process;

By formula 1. discretize, least square optimization method is utilized to realize the design in true resolution composite filter storehouse, the t after discrete ^k(ρ, z) is represented by following expression formula:

$t^k(\mathrm{\&rho;},z)=t_{VGF}^k\left(z\right)t_{RGF}^k\left(\mathrm{\&rho;}\right)$

$t_{VGF}^k\left(z\right)=\frac{2\sqrt{ln2}}{w_z\sqrt{\mathrm{\&pi;}}}{e}^{-4\ln 2{\left(\frac{z-c_z}{w_z}\right)}^2}$

$t_{RGF}^k\left(\mathrm{\&rho;}\right)=\frac{{\mathrm{lnw}}_{\mathrm{\&rho;}}}{\sqrt{\mathrm{\&pi;}}}{e}^{-{\left(\frac{ln\mathrm{\&rho;}-{\mathrm{lnc}}_{\mathrm{\&rho;}}}{{\mathrm{lnw}}_{\mathrm{\&rho;}}}\right)}^2}$

In formula, w _zand c _zfull width at half maximum width and the center of longitudinal objective function respectively; w _ρand c _ρfor the shape parameters of radial objective function, the requirement according to different investigation depths and radial shape is determined;

Concrete steps are as follows:

1) two dimension target function parameter is determined, by adjustment w _zmake the resolution of composite function equal the resolution of objective function, then determine radial objective function parameters w according to the requirement of institute's investigation depth resolution _ρand c _ρ;

2) discrete one dimension, two-dimentional response function is 1. produced for formula;

3) two dimension target function and the cross correlation function measuring sub-array responds function is calculated, and according to allowing the length of degree of correlation determination wave filter;

4) alignment error transmission controe parameter and optimal control parameter;

5) adopt least square method optimization to obtain one group and focus on composite filter coefficient, if not required optimum focusing composite filter, repeat the 4th) step, until find out optimum focusing composite filter;

6) step 1 is repeated)-5) obtain the true resolution composite filter storehouse of all background conductance rates segmentation;

The adjustment of described error propagation controling parameters is zero for making the coefficient attenuation of wave filter both sides, and it is realized by following expression formula, and this is added in two dimension target function,

$s\left(z\right)=\frac{1}{\mathrm{\&alpha;}+y\left(z\right)}-\frac{1}{\mathrm{\&alpha;}+y\left(0\right)}$

$y\left(z\right)=e^{-\frac{1}{2}{\left(\frac{z-z_0}{\mathrm{\&beta;}}\right)}^2}$

Wherein function y (z) is error distribution function; Z is independent variable; z ₀for the value mid point of z; Y (0) is error distribution function value when independent variable z value is 0; S (z) is error complementary function; β is the controling parameters of soft-constraint width; α is the precision controlling of wave filter constraint both sides; Force the both sides value of wave filter close to zero by regulating parameter α and β;

First determine the segmentation of background conductance rate when described step one intermediate-resolution matched filter storehouse (9) builds, then realized the structure of resolution match filter bank (9) by following formula,

g _pvgf(i+1)(z)＝h _pvgfi(z)*g _pvgfi(z)

In formula: g _{pvgf (i+1)}z () is for focusing on longitudinal Born geometric factor of synthesis low resolution curve, h _pvgfifor matched filter, g _pvgfiz () synthesizes longitudinal Born geometric factor of high resolving power curve for focusing on, i=1, and 2 ..., n-1, n are true resolution resultant curve number;

By above formula discretize, least square optimization method is then utilized to realize the design of resolution match filter bank.

2. the logging signal synthesis disposal route of array induction logging instrument according to claim 1, is characterized in that:

The calculating of resolution match wave filter new in the calculating generation of composite filter new in described step 3 and step 4 generates and all adopts method of interpolation to calculate.

3. the logging signal synthesis disposal route of array induction logging instrument according to claim 1, is characterized in that:

Information compensation process in described step 4 is by new matched filter and high resolving power curve convolution, form adjacent low resolution curve, then deduct the low resolution curve of generation with high resolving power curve, obtain the difference information of adjacent resolution curve, its mapping mode comprises:

Low resolution curvilinear transformation is high resolving power curve, is added to the curve that low resolution curve obtains by difference information;

High resolving power curvilinear transformation is low resolution curve, is added to the curve that high resolving power curve obtains by difference information.

4. a logging signal synthesis disposal system for array induction logging instrument, comprises

Apparatus measures device (1), comprises array induction instrument (11) and instrument truck (12), for providing apparatus measures signal;

Hole correction module (2), the measuring-signal for pair array induction instrument aratus carries out borehole environment correction, eliminates the impact of borehole environment;

Signal syntheses processing module (3), the measuring-signal after described hole correction resume module transfers to the synthesis process that described signal syntheses processing module (3) carries out signal, it is characterized in that:

Described signal syntheses processing module (3) comprises true resolution synthesis unit (4) and resolution match unit (5); Described true resolution synthesis unit (4) is provided with true resolution composite filter storehouse (8); Described resolution match unit (5) is provided with resolution match filter bank (9); True resolution synthesis unit (4) and resolution match unit (5) are provided with effective background conductance rate computing unit (6) jointly, wherein, described effective background conductance rate computing unit (6) comprises effective background conductance rate calculating filter storehouse (7);

The measuring-signal of correction is transferred to true resolution synthesis unit (4) and effective background conductance rate computing unit (6) by described hole correction module (2) respectively, the structure value in the background conductance rate score after the measuring-signal that hole correction module (2) corrects by described true resolution synthesis unit (4) and described effective background conductance rate computing unit (6) calculate and true resolution composite filter storehouse (8) carries out true resolution synthesis, Signal transmissions after synthesis is to described resolution match unit (5), background conductance rate score after effective background conductance rate computing unit (6) calculates by resolution match unit (5) and the structure value of resolution match filter bank (9) carry out resolution match, thus complete the logging signal synthesis process of array induction logging instrument,

The structure of described effective background conductance rate calculating filter storehouse (7) is realized by following formula

$t_{VGF}\left(z\right)=\frac{1.65}{w_z\sqrt{2\mathrm{\&pi;}}}{e}^{-\frac{1}{2}{\&lsqb;\frac{1.65\&times;z}{w_z}\&rsqb;}^2}$

In formula: t _vGFz () is area normalization Gaussian function corresponding to effective background conductance rate calculating filter storehouse; Z is independent variable; Parameter w _zfor the width that target Gaussian function vertical response 90% main lobe signal is corresponding;

Wherein w _zdetermining step is as follows:

1. longitudinal infinitesimal geometry factor of selected subarray in segmentation place of background conductance rate is calculated, and normalized, contribute corresponding width as the initial parameter of Gaussian function with longitudinal infinitesimal geometry factor main lobe 90%

2. by the measuring-signal convolution on selected Gaussian function and stratum, background conductance rate must be looked;

3. according to the log response storehouse calculated in uniform dielectric, by linear interpolation, counter pushing away looks true background conductance rate corresponding to background conductance rate;

4. the anti-true background conductance rate pushed away compared with given background conductance rate, if both are consistent, then selected Gaussian function is the wave filter of given background conductance rate, otherwise adjustment gaussian function is recalculated, until obtain the wave filter of given background conductance rate from the 1. step;

5. repeat above-mentioned steps, complete the design of effective background conductance rate filter bank of all background conductance rates and all subarrays;

First determine true resolution synthesis aim curve investigation depth and the segmentation of background conductance rate when described true resolution composite filter storehouse (8) builds, then realized the structure of true resolution composite filter by following formula,

$t^k(\mathrm{\&rho;},z)=\underset{j=1}{\overset{J}{\&Sigma;}}\underset{z^{\&prime;}=z_{\min }}{\overset{z_{\max }}{\&Sigma;}}{w}_j^k(z^{\&prime;},{\mathrm{\&sigma;}}_b)g_j(\mathrm{\&rho;},z-z^{\&prime;},{\mathrm{\&sigma;}}_b)$ ①

In formula: ρ is radial depth of investigetion coordinate, z is longitudinal depth coordinate, σ _bfor background conductance rate, t ^k(ρ, z) synthesizes the two dimension target function of kth kind investigation depth for focusing on, for the true resolution composite filter of jth subarray during synthesis kth kind investigation depth, g _j(ρ, z-z ', σ _b) be the two-dimentional Born geometric factor of a jth subarray; J=1 ~ J is synthesis subarray number used, z '=z _min~ z _maxfor focusing on the window ranges of synthesis process;

By formula 1. discretize, least square optimization method is utilized to realize the design in true resolution composite filter storehouse, the t after discrete ^k(ρ, z) is represented by following expression formula:

$t^k(\mathrm{\&rho;},z)=t_{VGF}^k\left(z\right)t_{RGF}^k\left(\mathrm{\&rho;}\right)$

$t_{VGF}^k\left(z\right)=\frac{2\sqrt{ln2}}{w_z\sqrt{\mathrm{\&pi;}}}{e}^{-4\ln 2{\left(\frac{z-c_z}{w_z}\right)}^2}$

$t_{RGF}^k\left(\mathrm{\&rho;}\right)=\frac{{\mathrm{lnw}}_{\mathrm{\&rho;}}}{\sqrt{\mathrm{\&pi;}}}{e}^{-{\left(\frac{ln\mathrm{\&rho;}-{\mathrm{lnc}}_{\mathrm{\&rho;}}}{{\mathrm{lnw}}_{\mathrm{\&rho;}}}\right)}^2}$

In formula, w _zand c _zfull width at half maximum width and the center of longitudinal objective function respectively; w _ρand c _ρfor the shape parameters of radial objective function, the requirement according to different investigation depths and radial shape is determined;

Concrete steps are as follows:

1) two dimension target function parameter is determined, by adjustment w _zmake the resolution of composite function equal the resolution of objective function, then determine radial objective function parameters w according to the requirement of institute's investigation depth resolution _ρand c _ρ;

2) discrete one dimension, two-dimentional response function is 1. produced for formula;

3) two dimension target function and the cross correlation function measuring sub-array responds function is calculated, and according to allowing the length of degree of correlation determination wave filter;

4) alignment error transmission controe parameter and optimal control parameter;

5) adopt least square method optimization to obtain one group and focus on composite filter coefficient, if not required optimum focusing composite filter, repeat the 4th) step, until find out optimum focusing composite filter;

6) step 1 is repeated)-5) obtain the true resolution composite filter storehouse of all background conductance rates segmentation;

The adjustment of described error propagation controling parameters is zero for making the coefficient attenuation of wave filter both sides, and it is realized by following expression formula, and this is added in two dimension target function,

$s\left(z\right)=\frac{1}{\mathrm{\&alpha;}+y\left(z\right)}-\frac{1}{\mathrm{\&alpha;}+y\left(0\right)}$

$y\left(z\right)=e^{-\frac{1}{2}{\left(\frac{z-z_0}{\mathrm{\&beta;}}\right)}^2}$

Wherein function y (z) is error distribution function; Z is independent variable; z ₀for the value mid point of z; Y (0) is error distribution function value when independent variable z value is 0; S (z) is error complementary function; β is the controling parameters of soft-constraint width; α is the precision controlling of wave filter constraint both sides; Force the both sides value of wave filter close to zero by regulating parameter α and β;

First determine the segmentation of background conductance rate when described resolution match filter bank (9) builds, then realized the structure of resolution match filter bank (9) by following formula,

g _pvgf(i+1)(z)＝h _pvgfi(z)*g _pvgfi(z)

In formula: g _{pvgf (i+1)}z () is for focusing on longitudinal Born geometric factor of synthesis low resolution curve, h _pvgfifor matched filter, g _pvgfiz () synthesizes longitudinal Born geometric factor of high resolving power curve for focusing on, i=1, and 2 ..., n-1, n are true resolution resultant curve number;

By above formula discretize, least square optimization method is then utilized to realize the design of resolution match filter bank.

5. the logging signal synthesis disposal system of array induction logging instrument according to claim 4, is characterized in that:

Described effective background conductance rate computing unit (6) and described true resolution composite filter storehouse (8) build new composite filter, and new composite filter and the measuring-signal after described hole correction module (2) corrects carry out true resolution and synthesize and process.

6. the logging signal synthesis disposal system of array induction logging instrument according to claim 4, is characterized in that:

Described effective background conductance rate computing unit (6) and described resolution match filter bank (9) build new matched filter, and new matched filter carries out the coupling of resolution to the composite signal after described true resolution synthesis unit (4) process.

7. the logging signal synthesis disposal system of array induction logging instrument according to claim 4, is characterized in that:

The true resolution synthesis of described array induction instrument (11) measuring-signal realizes according to the following equation:

${\mathrm{\&sigma;}}_p({\mathrm{\&rho;}}_k,z)=\underset{j=1}{\overset{J}{\&Sigma;}}\underset{z^{\&prime;}=z_{\min }}{\overset{z_{\max }}{\&Sigma;}}{w}_j^k(z^{\&prime;},\mathrm{\&sigma;}){\mathrm{\&sigma;}}_a^j(z-z^{\&prime;})$

Wherein, σ _p(ρ _k, z) be radial depth of investigetion ρ after synthesis _k, the conductivity at longitudinal degree of depth z place; focus on to be synthesized to investigation depth ρ _ktime, the true resolution composite filter of a jth subarray, for a jth subarray is at the measured value of z point through hole correction; K=1,2 ..., K, K are synthesis investigation depth numbers, j=1,2 ..., J, J are synthesis subarray number used; Z '=z _min~ z _maxfor focusing on the window ranges of synthesis process.

8., according to the logging signal synthesis disposal system of the arbitrary described array induction logging instrument of claim 4-7, it is characterized in that:

Coil in described array induction instrument (11) is made up of multiple subarray, each described subarray comprises a transmitting coil, a main receiving coil and 1-2 potted coil, subarray wherein adopts single-sided arrangement and/or both sides arrangement, and the quantity of frequency signal that each subarray receives is less than the quantity of the frequency signal that described transmitting coil is launched.