CN105549114A - Method and device for calculating logging parameters of low-porosity and compact limestone rock structure components - Google Patents

Abstract

The invention provides a method and a device for calculating logging parameters of low-porosity and compact limestone rock structure components, which relate to the technical field of petroleum geological exploration and logging, and the method comprises the following steps: acquiring logging data of a limestone rock reservoir to be detected; determining the porosity of the formation fracture according to the deep resistivity logging curve and the flushing zone resistivity logging curve; determining a first rock structure number parameter, a second rock structure number parameter, a third rock structure number parameter and a fourth rock structure number parameter according to the logging data, and then superposing to determine the sum of the rock structure number parameters; and determining the final value of the rock structure number parameter according to the sum of the rock structure number parameter and the depth resistivity in the depth resistivity logging curve. The method can solve the problem that the rock structure number formula in the prior art is not suitable for a low-porosity and compact limestone rock reservoir with the core porosity range of 0.5-10%, so that the structure number parameters of the low-porosity and compact limestone rock cannot be accurately determined.

Description

The computing method of low hole, vaughanite texture constituent of rock log parameter and device

Technical field

The present invention relates to petroleum geology exploration and logging technology field, particularly relate to a kind of low hole, the computing method of vaughanite texture constituent of rock log parameter and device.

Background technology

At present, utilize well-log information identification limestone texture constituent of rock to be a world-famous puzzle always, and how to obtain the key that effective log parameter is raising accuracy of identification.Research method and the thinking of the effective log parameter of current acquisition are mainly divided into following step: 1) rock-electric relationship analysis: first to the core hole of core hole, obtain sample, abrasive disc is carried out to sample, obtain thin slice, microscope is utilized to carry out thin section identification, obtain texture constituent of rock information, and depth-logger that the thin slice degree of depth is playbacked to, then texture constituent of rock and well logging correlation analysis is carried out, determine the log parameter (Logging Curves best with texture constituent of rock correlativity, synthetic parameters curve etc.), set up the logging characteristic parameters Sample Storehouse of texture constituent of rock.2) technology identification: mainly adopt the mathematical algorithms such as neural network, cluster analysis, discriminatory analysis, fuzzy mathematics, neural network is trained the sample log parameter storehouse set up, and obtains artificial intelligence; Clustering methodology sets up electrofacies and texture constituent of rock corresponding relation; Discriminant analysis method sets up the well logging discrimination formula of different texture constituent of rock.Mathematics method sets up the membership function of different texture constituent of rock well logging classification mode.3) technology application: by the research to study area core hole, determines best log parameter, chooses most suitable recognition technology, carry out texture constituent of rock identification to non-core hole, and verifies recognition effect with other core holes in work area.

Study less at home about texture constituent of rock Logging Identification Method in prior art.And special be abroad carry out rock physics three classification according to limestone rock particle or crystallite dimension for texture constituent of rock well logging recognition study general, size is greater than 100 μm and is defined as grain pattern by the first kind, size is defined as micrite grain pattern at 20 μm-100 μm by Equations of The Second Kind, size is less than 20 μm by the 3rd class, and to be defined as micrite be that main structure (comprises part micrite grain pattern, particle micrite structure, micrite structure), pass through core porosity, permeability, water saturation parameter establishes rock-fabric number (RockFabricNumber) computing formula, factor of porosity is calculated by well logging, permeability and water saturation replace core parameters, realize texture constituent of rock well logging recognition.Under the method is based upon following two geologic conditions: (1) stratum mainly based on intergranular, intracrystalline pore or dissolution pore, crack agensis; (2) stratum is on gas-water interface, there is not pure water layer, meets these two condition sub-surface texture constituent of rock well logging recognition and can obtain fine recognition effect.

The method of prior art effectively solves the technical bottleneck of texture constituent of rock well logging recognition, but its rock-fabric number formula also Shortcomings applied: it sets up basis is that formation core porosity ranges is between 4%-40%, and for domestic low hole, vaughanite rock formation, core porosity scope, between 0.5%-10%, accurately cannot determine the structure number parameter of low hole, vaughanite rock by the rock-fabric number formula of prior art.

Summary of the invention

The embodiment of the present invention provides a kind of low hole, the computing method of vaughanite texture constituent of rock log parameter and device, be not suitable for the low hole of core porosity scope between 0.5%-10%, vaughanite rock formation to solve rock-fabric number formula of the prior art, cause the problem of structure number parameter accurately cannot determining low hole, vaughanite rock.

In order to achieve the above object, the present invention adopts following technical scheme:

Computing method for low hole, vaughanite texture constituent of rock log parameter, comprising:

Obtain the log data of limestone rock formation to be measured; Described log data comprises: deep resistivity log curve, flushed zone resistivity logging trace, mud filtrate resistivity, gamma ray curve and interval transit time curve;

According to described deep resistivity log curve and flushed zone resistivity logging trace layer fracture porosity definitely;

According to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter;

Described first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter are superposed, determines rock-fabric number parameter summation;

According to the dark resistivity size in described rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value.

Concrete, according to described deep resistivity log curve and flushed zone resistivity logging trace layer fracture porosity definitely, comprising:

Judge the size of the dark resistivity under synchronization in deep resistivity log curve and the flushed zone resistivity in flushed zone resistivity logging trace;

If dark resistivity is greater than flushed zone resistivity, then according to formula:

${\mathrm{\Φ}}_{frc}=(\frac{8.52253}{RXO}-\frac{8.242778}{RT}+0.00071236)\×Rmf$

Calculate formation fracture factor of porosity Φ _frc; Wherein RXO is flushed zone resistivity; RT is dark resistivity; Rmf is mud filtrate resistivity;

If dark resistivity is less than flushed zone resistivity, then according to formula:

${\mathrm{\Φ}}_{frc}=(\frac{1.99247}{RT}-\frac{0.992719}{RXO}+0.00031829)\×Rmf$

Calculate formation fracture factor of porosity Φ _frc.

Concrete, according to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, comprising:

According to the size of described natural gamma value size and formation fracture factor of porosity, determine described first rock-fabric number parameter;

The described size according to described natural gamma value size and formation fracture factor of porosity, determine described first rock-fabric number parameter, comprising:

If natural gamma value is greater than 50API, and formation fracture factor of porosity is less than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\×GR-1.125$

Calculate described first rock-fabric number parameter ARFN ₁; Wherein, RT is dark resistivity; GR is natural gamma value;

If natural gamma value is greater than 50API, and formation fracture factor of porosity is greater than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\×GR-0.625$

Calculate described first rock-fabric number parameter ARFN ₁.

In addition, according to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, comprising:

According to size and the natural gamma value size of described interval transit time value, determine described second rock-fabric number parameter;

The described size according to described interval transit time value and natural gamma value size, determine described second rock-fabric number parameter, comprising:

If interval transit time value is less than or equal to 51us/f, and natural gamma value is less than or equal to 50API, according to formula:

${\mathrm{ARFN}}_2=e^{\frac{\left(2.9+\lg \left(AC\right)+3.2\×\lg \left(RT\×e^{\left(\frac{RT-200}{1800}\right)}\right)\right)}{(6.9+\lg (AC\left)\right)}}-3$

Calculate described second rock-fabric number parameter ARFN ₂; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

In addition, according to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, comprising:

According to the size of described interval transit time value, natural gamma value size and dark resistivity size, determine described 3rd rock-fabric number parameter;

The described size according to described interval transit time value, natural gamma value size and dark resistivity size, determine described 3rd rock-fabric number parameter, comprising:

If interval transit time value is greater than 51us/f, natural gamma value is less than or equal to 50API, and dark resistivity is greater than 200ohmm, according to formula:

${\mathrm{ARFN}}_3=9.5-e^{\frac{(3.2+\lg (AC)+2*\lg (\frac{RT}{AC}\left)\right)}{(3.1+lg(AC\left)\right)}}$

Calculate described 3rd rock-fabric number parameter ARFN ₃; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

In addition, according to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, comprising:

According to the size of described interval transit time value, the size of natural gamma value size, deeply resistivity and formation fracture factor of porosity, determine described 4th rock-fabric number parameter;

The size of the described size according to described interval transit time value, natural gamma value size, deeply resistivity and formation fracture factor of porosity, determine described 4th rock-fabric number parameter, comprising:

If interval transit time value is greater than 51us/f, natural gamma value is less than or equal to 40API, and dark resistivity is less than 200ohmm and formation fracture factor of porosity is less than 0.08%, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{(3.2+1.1*\lg (AC)+\lg (\frac{RT*RT}{GR*GR}\left)\right)}{(3.1+\lg (GR\left)\right)}}$

Calculate described 4th rock-fabric number parameter ARFN ₄; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity; GR is natural gamma value;

If interval transit time value is greater than 51us/f, natural gamma value is greater than 40API and is less than or equal to 50API, and dark resistivity is less than 200ohmm and formation fracture factor of porosity is greater than 0.08%, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{(3.2+1.1*\lg (AC)+\lg (\frac{RT*RT}{GR*GR}\left)\right)}{(3.1+\lg (GR\left)\right)}}+0.5$

Calculate described 4th rock-fabric number parameter ARFN ₄.

Concrete, described first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter are superposed, determine rock-fabric number parameter summation, comprising:

According to formula:

ARFN＝ARFN ₁+ARFN ₂+ARFN ₃+ARFN ₄

Determine rock-fabric number parameter summation ARFN.

Concrete, according to the dark resistivity size in described rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value, comprising:

If ARFN is less than or equal to the first preset value, determine that described rock-fabric number parameter end value ARFN ' equals described ARFN;

If ARFN is greater than described first preset value, and dark resistivity is greater than 9000ohmm, determines that described rock-fabric number parameter end value ARFN ' equals the second preset value;

If ARFN is greater than described first preset value, and dark resistivity is less than or equal to 9000ohmm, determines that described rock-fabric number parameter end value ARFN ' equals described first preset value;

If ARFN is less than or equal to the 3rd preset value, determine that described rock-fabric number parameter end value ARFN ' equals described 3rd preset value.

A calculation element for low hole, vaughanite texture constituent of rock log parameter, comprising:

Log data acquiring unit, for obtaining the log data of limestone rock formation to be measured; Described log data comprises: deep resistivity log curve, flushed zone resistivity logging trace, mud filtrate resistivity, gamma ray curve and interval transit time curve;

Formation fracture factor of porosity determining unit, for according to described deep resistivity log curve and flushed zone resistivity logging trace layer fracture porosity definitely;

Rock-fabric number parameter determination unit, for the size according to the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter;

Rock-fabric number parameter superpositing unit, for described first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter being superposed, determines rock-fabric number parameter summation;

Rock-fabric number parameter end value determining unit, for according to the dark resistivity size in described rock-fabric number parameter summation and deep resistivity log curve, determines rock-fabric number parameter end value.

Concrete, this formation fracture factor of porosity determining unit, comprising:

Judge module, for judging the size of the dark resistivity under synchronization in deep resistivity log curve and the flushed zone resistivity in flushed zone resistivity logging trace;

First computing module, for when dark resistivity is greater than flushed zone resistivity, according to formula:

${\mathrm{\Φ}}_{frc}=(\frac{8.52253}{RXO}-\frac{8.242778}{RT}+0.00071236)\×Rmf$

Calculate formation fracture factor of porosity Φ _frc; Wherein RXO is flushed zone resistivity; RT is dark resistivity; Rmf is mud filtrate resistivity;

Second computing module, for when dark resistivity is less than flushed zone resistivity, according to formula:

${\mathrm{\Φ}}_{frc}=(\frac{1.99247}{RT}-\frac{0.992719}{RXO}+0.00031829)\×Rmf$

Calculate formation fracture factor of porosity Φ _frc.

Concrete, described rock-fabric number parameter determination unit, comprising:

First rock-fabric number parameter calculating module, for the size according to described natural gamma value size and formation fracture factor of porosity, determines described first rock-fabric number parameter;

Described first rock-fabric number parameter calculating module, specifically for:

50API is greater than in natural gamma value, and when formation fracture factor of porosity is less than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\×GR-1.125$

Calculate described first rock-fabric number parameter ARFN ₁; Wherein, RT is dark resistivity; GR is natural gamma value;

50API is greater than in natural gamma value, and when formation fracture factor of porosity is greater than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\×GR-0.625$

Calculate described first rock-fabric number parameter ARFN ₁.

In addition, described rock-fabric number parameter determination unit, comprising:

Second rock-fabric number parameter calculating module, for according to the size of described interval transit time value and natural gamma value size, determines described second rock-fabric number parameter;

Described second rock-fabric number parameter calculating module, specifically for:

51us/f is less than or equal in interval transit time value, and when natural gamma value is less than or equal to 50API, according to formula:

${\mathrm{ARFN}}_2=e^{\frac{\left(2.9+\lg \left(AC\right)+3.2\×\lg \left(RT\×e^{\left(\frac{RT-200}{1800}\right)}\right)\right)}{(6.9+\lg (AC\left)\right)}}-3$

Calculate described second rock-fabric number parameter ARFN ₂; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

In addition, described rock-fabric number parameter determination unit, comprising:

3rd rock-fabric number parameter calculating module, for the size according to described interval transit time value, natural gamma value size and dark resistivity size, determines described 3rd rock-fabric number parameter;

Described 3rd rock-fabric number parameter calculating module, specifically for::

Be greater than 51us/f in interval transit time value, natural gamma value is less than or equal to 50API, and when dark resistivity is greater than 200ohmm, according to formula:

${\mathrm{ARFN}}_3=9.5-e^{\frac{(3.2+\lg (AC)+2*\lg (\frac{RT}{AC}\left)\right)}{(3.1+lg(AC\left)\right)}}$

Calculate described 3rd rock-fabric number parameter ARFN ₃; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

In addition, described rock-fabric number parameter determination unit, comprising:

4th rock-fabric number parameter calculating module, for the size of the size according to described interval transit time value, natural gamma value size, dark resistivity and formation fracture factor of porosity, determines described 4th rock-fabric number parameter;

Described 4th rock-fabric number parameter calculating module, specifically for:

Be greater than 51us/f in interval transit time value, natural gamma value is less than or equal to 40API, dark resistivity be less than 200ohmm and formation fracture factor of porosity is less than 0.08% time, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{(3.2+1.1*\lg (AC)+\lg (\frac{RT*RT}{GR*GR}\left)\right)}{(3.1+\lg (GR\left)\right)}}$

Calculate described 4th rock-fabric number parameter ARFN ₄; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity; GR is natural gamma value;

Be greater than 51us/f in interval transit time value, natural gamma value is greater than 40API and is less than or equal to 50API, dark resistivity be less than 200ohmm and formation fracture factor of porosity is greater than 0.08% time, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{(3.2+1.1*\lg (AC)+\lg (\frac{RT*RT}{GR*GR}\left)\right)}{(3.1+\lg (GR\left)\right)}}+0.5$

Calculate described 4th rock-fabric number parameter ARFN ₄.

In addition, described rock-fabric number parameter superpositing unit, specifically for:

According to formula:

ARFN＝ARFN ₁+ARFN ₂+ARFN ₃+ARFN ₄

Determine rock-fabric number parameter summation ARFN.

In addition, described rock-fabric number parameter end value determining unit, specifically for:

When ARFN is less than or equal to the first preset value, determine that described rock-fabric number parameter end value ARFN ' equals described ARFN;

Be greater than described first preset value at ARFN, and when dark resistivity is greater than 9000ohmm, determine that described rock-fabric number parameter end value ARFN ' equals the second preset value;

Be greater than described first preset value at ARFN, and when dark resistivity is less than or equal to 9000ohmm, determine that described rock-fabric number parameter end value ARFN ' equals described first preset value;

When ARFN is less than or equal to the 3rd preset value, determine that described rock-fabric number parameter end value ARFN ' equals described 3rd preset value.

The computing method of a kind of low hole that the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter and device, four rock-fabric number parameters are asked for respectively, rock-fabric number parameter summation is determined in final superposition again, and then according to the dark resistivity size in described rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value.The present invention is applicable to stratum and is mainly low hole, densification, and formation porosity is less than 10%, or with fracture development, sedimentary facies belt is mainly the pure limestone formation of open-platform or half open-platform sedimentary facies.Setting up basis compared to the rock-fabric number formula of prior art is that formation core porosity ranges is between 4%-40%, present invention, avoiding for domestic low hole, vaughanite rock formation, core porosity scope, between 0.5%-10%, accurately cannot determine the problem of structure number parameter of low hole, vaughanite rock by the rock-fabric number formula of prior art.

Accompanying drawing explanation

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

The process flow diagram one of the computing method of a kind of low hole that a kind of embodiment of the present invention that Fig. 1 provides for the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter;

The flowchart 2 of the computing method of a kind of low hole that a kind of embodiment of the present invention that Fig. 2 provides for the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter;

Fig. 3 is that core porosity, rock core water containing saturability and the core wafer in the embodiment of the present invention analyzes texture constituent of rock X plot;

Fig. 4 is that interval transit time-dark resistivity, natural gamma-dark resistivity and the core wafer in the embodiment of the present invention analyzes texture constituent of rock X plot;

Fig. 5 is the interval transit time-dark resistivity of four kinds of different models in the embodiment of the present invention, natural gamma-dark resistivity and core wafer analyze texture constituent of rock X plot;

Fig. 6 is the interval transit time-dark resistivity of the first model in the embodiment of the present invention, natural gamma-dark resistivity and core wafer analyze texture constituent of rock X plot;

Fig. 7 is the interval transit time-dark resistivity of second model in the embodiment of the present invention, natural gamma-dark resistivity and core wafer analyze texture constituent of rock X plot;

Fig. 8 is the interval transit time-dark resistivity of the third model in the embodiment of the present invention, natural gamma-dark resistivity and core wafer analyze texture constituent of rock X plot;

Fig. 9 is the interval transit time-dark resistivity of the 4th kind of model in the embodiment of the present invention, natural gamma-dark resistivity and core wafer analyze texture constituent of rock X plot;

Figure 10 is that in 161 wells in the tower in the embodiment of the present invention, the middle ancient times 51 and tower, rock-fabric number well logging calculating design sketch looked by 201 wells;

The structural representation one of the calculation element of a kind of low hole that Figure 11 provides for the embodiment of the present invention, vaughanite texture constituent of rock log parameter;

The structural representation two of the calculation element of a kind of low hole that Figure 12 provides for the embodiment of the present invention, vaughanite texture constituent of rock log parameter.

Embodiment

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

As shown in Figure 1, the computing method of a kind of low hole that the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter, comprising:

Step 101, obtain the log data of limestone rock formation to be measured.

Wherein, log data comprises: deep resistivity log curve, flushed zone resistivity logging trace, mud filtrate resistivity, gamma ray curve and interval transit time curve.

Step 102, according to deep resistivity log curve and flushed zone resistivity logging trace layer fracture porosity definitely.

Step 103, size according to the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter.

Step 104, by the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter superpose, determine rock-fabric number parameter summation.

Step 105, according to the dark resistivity size in rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value.

The computing method of a kind of low hole that the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter, four rock-fabric number parameters are asked for respectively, rock-fabric number parameter summation is determined in final superposition again, and then according to the dark resistivity size in rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value.The present invention is applicable to stratum and is mainly low hole, densification, and formation porosity is less than 10%, or with fracture development, sedimentary facies belt is mainly the pure limestone formation of open-platform or half open-platform sedimentary facies.Setting up basis compared to the rock-fabric number formula of prior art is that formation core porosity ranges is between 4%-40%, present invention, avoiding for domestic low hole, vaughanite rock formation, core porosity scope, between 0.5%-10%, accurately cannot determine the problem of structure number parameter of low hole, vaughanite rock by the rock-fabric number formula of prior art.

Enumerate an embodiment specifically below, better understand the present invention to make those skilled in the art.As shown in Figure 2, the computing method of a kind of low hole that the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter, comprising:

Step 201, obtain the log data of limestone rock formation to be measured.

Wherein, log data comprises: deep resistivity log curve, flushed zone resistivity logging trace, mud filtrate resistivity, gamma ray curve and interval transit time curve.

Can be applied to In The Central Tarim Area Ordovician strata in conjunction with the embodiment of the present invention is herein that example is described: TZ area Ordovician strata is the Marine source rock of a set of open-platform deposition, and lithology is limestone.This cover formation rock structural constituent has 5 kinds: micrite structure, particle micrite structure, micrite grain pattern, grain pattern and reef structure, and stratum, this area is fine and close, and crack comparative development, above geologic condition is applicable to application of the present invention.Parameter acquiring: utilize CLS-5700 logging suite apparatus measures stratum natural gamma GR, interval transit time DT and resistivity RT, resistivity measurement is carried out to mud filtrate, mud filtrate resistivity 0.15 Ω .m (ohm meter)/84 DEG C.

Deep resistivity log curve herein, flushed zone resistivity logging trace can be obtained by Formation Resistivity Measurement.

Interval transit time curve can comprise resident fluid interval transit time value and rock skeleton interval transit time value.Wherein, resident fluid interval transit time value DT _fgeneral employing theoretical value 189 μ s/f, rock skeleton interval transit time value DT _magenerally according to the setting of stratum principal lithologic, such as stratum is limestone formation, DT _maadopt theoretical value 47 μ s/f, pierite stratum then adopts theoretical value 42 μ s/f, if mixing lithology then takes both weighted calculation mean value according to geological condition, this value is between 42-47 μ s/ft, and the interval transit time on general stratum is at 40-90 μ s/f.

Natural gamma GR in above-mentioned gamma ray curve is mainly stratum natural gamma value (GR _ma), and resident fluid natural gamma is almost 0.Wherein GR _magenerally according to the setting of stratum principal lithologic, such as stratum is limestone formation, then GR _maadopt theoretical value 8API, mudstone stratum then adopts theoretical value 230API, if mixing lithology, then take both weighted calculation mean value according to geological condition, this value is between 0-70API, and the natural gamma on general stratum can between 10-200API.

Above-mentioned dark resistivity RT can comprise resident fluid resistivity (RT _f) and rock skeleton resistivity (RT _ma).Wherein, resident fluid is mainly water, RT _fbe generally 0.01 ~ 1ohmm, RT _magenerally dense, non-conductive (except pyrite), generally can at more than 10000ohmm, and the resistivity on general stratum is between 20-200000ohmm.

Logging trace of the present invention can adopt prior art equipment to obtain, such as CLS-3700 well logger, CLS-5700 well logger etc.

Step 202, according to deep resistivity log curve and flushed zone resistivity logging trace layer fracture porosity definitely.

Concrete, need to judge the size of the dark resistivity under synchronization in deep resistivity log curve and the flushed zone resistivity in flushed zone resistivity logging trace herein;

If dark resistivity is greater than flushed zone resistivity, then according to formula:

${\mathrm{\Φ}}_{frc}=(\frac{8.52253}{RXO}-\frac{8.242778}{RT}+0.00071236)\×Rmf$

Calculate formation fracture factor of porosity Φ _frc; Wherein RXO is flushed zone resistivity; RT is dark resistivity; Rmf is mud filtrate resistivity;

If dark resistivity is less than flushed zone resistivity, then according to formula:

${\mathrm{\Φ}}_{frc}=(\frac{1.99247}{RT}-\frac{0.992719}{RXO}+0.00031829)\×Rmf$

Calculate formation fracture factor of porosity Φ _frc.

Such as 161 core holes in tower, depth of stratum is 4384.1m is micrite grainstone, and this layer of dark resistivity of well log measurement is 547.ohmm, and flushed zone resistivity is 501.74ohmm, mud filtrate resistivity R _mfequal 0.15ohmm and substitute into above formula, calculate fracture porosity and equal 0.014%.

In following steps 203-step 207, the calculating of rock-fabric number parameter of the present invention will be described.

Rock-fabric number computing technique starts from Lucia in 2005 and carries out rock physics three classification according to carbonatite particle or crystallite dimension, size is greater than 100 μm and is defined as grain pattern by the first kind, size is defined as micrite grain pattern by Equations of The Second Kind between 20 μm-100 μm, size is less than 20 μm by the 3rd class, and to be defined as micrite be that main structure (comprises part micrite grain pattern, particle micrite structure, micrite structure), pass through core porosity, water saturation X plot (as shown in Figure 3) establishes that to look rock-fabric number (ApparentRockFabricNumber) computing formula as follows:

$ARFN={10}^{\left(\frac{3.11+1.88\log \left(\mathrm{\Φ}\right)+\log \left(S_w\right)}{3.06+1.4\log \left(\mathrm{\Φ}\right)}\right)}$

In formula, Φ represents formation porosity (decimal unit), S _wrepresent stratum water saturation (decimal unit).

Wherein formation porosity and water saturation can utilize the acquisition that well logging calculates, thus realize well logging calculating rock-fabric number, according to rock-fabric number size, can realize texture constituent of rock well logging recognition.

Work as ARFN>4, represent micrite structure;

When ARFN is between 2.5 ~ 4, represent particle micrite structure and powder micrite structure;

When ARFN is between 1.5 ~ 2.5, represent micrite particle or fine powder crystal structure;

When ARFN is between 0.5 ~ 1.5, represent grain pattern or middle coarse structure.

Above-mentioned formula abroad carbonate formation (factor of porosity is generally greater than 4%) effect is better, but it is low for domestic carbonate reservoir hole, particularly factor of porosity is at the reservoir of 0.5%-4%, in this situation, Lucia formula is inapplicable, if directly utilize this formulae discovery rock-fabric number, rock-fabric number scope is 0-1000, can not be used for identifying structural constituent, such as when stratum is cryptite, cryptite factor of porosity is 1% (decimal unit is 0.01), cryptite water saturation is 100% (decimal unit is 1), substitute into above formula and calculate ARFN=0.0032, this numerical value has exceeded 0.5-4 scope, the identification of texture constituent of rock can not be carried out.

Due to low hole, stratum, densification and fracture development, its texture constituent of rock logging character has very large plyability, different texture constituent of rock logging character is difficult to distinguish (such as shown in Fig. 4), according to stratum geologic feature, divide 4 kinds of typical models (as shown in Figure 5), set up respectively and look rock-fabric number formula, then by the formula combinations of 4 models, obtain final rock-fabric number computing formula.Mode is as follows:

Step 203, size according to natural gamma value size and formation fracture factor of porosity, determine the first rock-fabric number parameter.

Herein, if natural gamma value is greater than 50API (such as shown in Fig. 6), and formation fracture factor of porosity is less than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\×GR-1.125$

Calculate the first rock-fabric number parameter ARFN ₁; Wherein, RT is dark resistivity; GR is natural gamma value;

If natural gamma value is greater than 50API, and formation fracture factor of porosity is greater than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\×GR-0.625$

Calculate the first rock-fabric number parameter ARFN ₁.

Step 204, according to the size of interval transit time value and natural gamma value size, determine the second rock-fabric number parameter.

Herein, if interval transit time value is less than or equal to 51us/f, and natural gamma value is less than or equal to 50API (such as shown in Fig. 7), according to formula:

${\mathrm{ARFN}}_2=e^{\frac{\left(2.9+\lg \left(AC\right)+3.2\×\lg \left(RT\×e^{\left(\frac{RT-200}{1800}\right)}\right)\right)}{(6.9+\lg (AC\left)\right)}}-3$

Calculate the second rock-fabric number parameter ARFN ₂; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

Step 205, size according to interval transit time value, natural gamma value size and dark resistivity size, determine the 3rd rock-fabric number parameter.

Herein, if interval transit time value is greater than 51us/f, natural gamma value is less than or equal to 50API (such as shown in Fig. 8), and dark resistivity is greater than 200ohmm, according to formula:

${\mathrm{ARFN}}_3=9.5-e^{\frac{(3.2+\lg (AC)+2*\lg (\frac{RT}{AC}\left)\right)}{(3.1+lg(AC\left)\right)}}$

Calculate the 3rd rock-fabric number parameter ARFN ₃; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

The size of step 206, size according to interval transit time value, natural gamma value size, dark resistivity and formation fracture factor of porosity, determines the 4th rock-fabric number parameter.

Herein, if interval transit time value is greater than 51us/f, natural gamma value is less than or equal to 40API, and dark resistivity is less than 200ohmm (such as shown in Fig. 9) and formation fracture factor of porosity is less than 0.08%, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{(3.2+1.1*\lg (AC)+\lg (\frac{RT*RT}{GR*GR}\left)\right)}{(3.1+\lg (GR\left)\right)}}$

Calculate the 4th rock-fabric number parameter ARFN ₄; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity; GR is natural gamma value.

If interval transit time value is greater than 51us/f, natural gamma value is greater than 40API and is less than or equal to 50API, and dark resistivity is less than 200ohmm and formation fracture factor of porosity is greater than 0.08%, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{(3.2+1.1*\lg (AC)+\lg (\frac{RT*RT}{GR*GR}\left)\right)}{(3.1+\lg (GR\left)\right)}}+0.5$

Calculate the 4th rock-fabric number parameter ARFN ₄.

Step 207, by the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter superpose, determine rock-fabric number parameter summation.

Herein specifically can according to formula:

ARFN＝ARFN ₁+ARFN ₂+ARFN ₃+ARFN ₄

Determine rock-fabric number parameter summation ARFN.

The interval transit time on stratum is such as had to be 48.83us/f, dark resistivity is 547.5ohmm, and natural gamma is 8.7API, then utilize aforesaid way to obtain ARFN=2.31, and actual rock core rock-fabric number is 2.1, basically identical with rock core structure number (as shown in Figure 10).

Step 208, according to the dark resistivity size in rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value.

Herein, if ARFN is less than or equal to one first preset value, determine that described rock-fabric number parameter end value ARFN ' equals described ARFN.

Herein, such as the first preset value can be 4, then when ARFN is less than or equal to 4, and such as 3, then rock-fabric number parameter end value ARFN ' equals 3.

If ARFN is greater than the first preset value, and dark resistivity is greater than 9000ohmm, determines that rock-fabric number parameter end value ARFN ' equals the second preset value.

The first preset value herein can be 4; This second preset value can be 4 or 4.5.The object of each preset value is that the ARFN curve numerical range of limit calculation ensures between 1-4.5 herein, has exceeded just cut-off and has fallen.

If ARFN is greater than the first preset value, and dark resistivity is less than or equal to 9000ohmm, determines that rock-fabric number parameter end value ARFN ' equals the first preset value.

Such as, the first preset value is 4, then ARFN is greater than 4, and when dark resistivity is less than or equal to 9000ohmm, determines that rock-fabric number parameter end value ARFN ' equals 4.

In addition, when ARFN is less than or equal to the 3rd preset value, determine that described rock-fabric number parameter end value ARFN ' equals described 3rd preset value.

Such as, the 3rd preset value can be 1, then, when ARFN is less than or equal to 1, described rock-fabric number parameter end value ARFN ' equals 1.

In embodiments of the present invention, the numerical range of ARFN can be generally between 1 to 4.5:

(1) when 0 < ARFN≤1, rock texture is defined as micrite structural constituent phase, and lithology is shale cryptite, cryptite breathes out partial particulate cryptite.

(2) when 1 < ARFN≤2, rock texture is defined as particle micrite structural constituent phase, and lithology is mainly particle cryptite, and bondstone and micrite support micrite grainstone;

(3) when 2 < ARFN≤3, rock texture is defined as micrite grain pattern component phase, and lithology is mainly particle-supported micrite grainstone;

(4) when 3 < ARFN≤4, rock texture is defined as grain pattern component phase, and lithology is mainly grainstone, dissolution pore reef limestone;

(5) when 4 < ARFN≤4.5, rock texture is defined as reef structural constituent phase, and lithology is fine and close biohermal limestone, dense granule limestone.

The present invention can obtain good effect on Ordovician system well compactness stratum, In The Central Tarim Area.

The computing method of a kind of low hole that the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter, four rock-fabric number parameters are asked for respectively, rock-fabric number parameter summation is determined in final superposition again, and then according to the dark resistivity size in rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value.The present invention is applicable to stratum and is mainly low hole, densification, and formation porosity is less than 10%, or with fracture development, sedimentary facies belt is mainly the pure limestone formation of open-platform or half open-platform sedimentary facies.Setting up basis compared to the rock-fabric number formula of prior art is that formation core porosity ranges is between 4%-40%, present invention, avoiding for domestic low hole, vaughanite rock formation, core porosity scope, between 0.5%-10%, accurately cannot determine the problem of structure number parameter of low hole, vaughanite rock by the rock-fabric number formula of prior art.

Corresponding to the embodiment of the method shown in above-mentioned Fig. 1, Fig. 2, the calculation element of a kind of low hole that the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter, as shown in figure 11, comprising:

Log data acquiring unit 31, can obtain the log data of limestone rock formation to be measured.

Log data comprises: deep resistivity log curve, flushed zone resistivity logging trace, mud filtrate resistivity, gamma ray curve and interval transit time curve.

Formation fracture factor of porosity determining unit 32, can according to deep resistivity log curve and flushed zone resistivity logging trace layer fracture porosity definitely.

Rock-fabric number parameter determination unit 33, according to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter can be determined.

Rock-fabric number parameter superpositing unit 34, can superpose the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, determine rock-fabric number parameter summation.

Rock-fabric number parameter end value determining unit 35, according to the dark resistivity size in rock-fabric number parameter summation and deep resistivity log curve, can determine rock-fabric number parameter end value.

Concrete, as shown in figure 12, this formation fracture factor of porosity determining unit 32, comprising:

Judge module 321, can judge the size of the dark resistivity under synchronization in deep resistivity log curve and the flushed zone resistivity in flushed zone resistivity logging trace.

First computing module 322, can when dark resistivity be greater than flushed zone resistivity, according to formula:

${\mathrm{\&Phi;}}_{frc}=(\frac{8.52253}{RXO}-\frac{8.242778}{RT}+0.00071236)\&times;Rmf$

Calculate formation fracture factor of porosity Φ _frc; Wherein RXO is flushed zone resistivity; RT is dark resistivity; Rmf is mud filtrate resistivity.

Second computing module 323, can when dark resistivity be less than flushed zone resistivity, according to formula:

${\mathrm{\&Phi;}}_{frc}=(\frac{1.99247}{RT}-\frac{0.992719}{RXO}+0.00031829)\&times;Rmf$

Calculate formation fracture factor of porosity Φ _frc.

Concrete, as shown in figure 12, this rock-fabric number parameter determination unit 33, can comprise:

First rock-fabric number parameter calculating module 331, according to the size of natural gamma value size and formation fracture factor of porosity, can determine the first rock-fabric number parameter.

First rock-fabric number parameter calculating module 331, specifically can:

50API is greater than in natural gamma value, and when formation fracture factor of porosity is less than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\&times;GR-1.125$

Calculate the first rock-fabric number parameter ARFN ₁; Wherein, RT is dark resistivity; GR is natural gamma value.

50API is greater than in natural gamma value, and when formation fracture factor of porosity is greater than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\&times;GR-0.625$

Calculate the first rock-fabric number parameter ARFN ₁.

In addition, as shown in figure 12, rock-fabric number parameter determination unit 33, also comprises:

Second rock-fabric number parameter calculating module 332, according to the size of interval transit time value and natural gamma value size, can determine the second rock-fabric number parameter.

Second rock-fabric number parameter calculating module 332, specifically can:

51us/f is less than or equal in interval transit time value, and when natural gamma value is less than or equal to 50API, according to formula:

${\mathrm{ARFN}}_2=e^{\frac{\left(2.9+\lg \left(AC\right)+3.2\&times;\lg \left(RT\&times;e^{\left(\frac{RT-200}{1800}\right)}\right)\right)}{(6.9+\lg (AC\left)\right)}}-3$

Calculate the second rock-fabric number parameter ARFN ₂; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

In addition, as shown in figure 12, rock-fabric number parameter determination unit 33, also comprises:

3rd rock-fabric number parameter calculating module 333, according to the size of interval transit time value, natural gamma value size and dark resistivity size, can determine the 3rd rock-fabric number parameter.

3rd rock-fabric number parameter calculating module 333, specifically can:

Be greater than 51us/f in interval transit time value, natural gamma value is less than or equal to 50API, and when dark resistivity is greater than 200ohmm, according to formula:

${\mathrm{ARFN}}_3=9.5-e^{\frac{(3.2+\lg (AC)+2*\lg (\frac{RT}{AC}\left)\right)}{(3.1+lg(AC\left)\right)}}$

Calculate the 3rd rock-fabric number parameter ARFN ₃; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

In addition, as shown in figure 12, this rock-fabric number parameter determination unit 33, also comprises:

4th rock-fabric number parameter calculating module 334, according to the size of the size of interval transit time value, natural gamma value size, dark resistivity and formation fracture factor of porosity, can determine the 4th rock-fabric number parameter.

4th rock-fabric number parameter calculating module 334, specifically can:

Be greater than 51us/f in interval transit time value, natural gamma value is less than or equal to 40API, dark resistivity be less than 200ohmm and formation fracture factor of porosity is less than 0.08% time, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{(3.2+1.1*\lg (AC)+\lg (\frac{RT*RT}{GR*GR}\left)\right)}{(3.1+\lg (GR\left)\right)}}$

Calculate the 4th rock-fabric number parameter ARFN ₄; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity; GR is natural gamma value.

Be greater than 51us/f in interval transit time value, natural gamma value is greater than 40API and is less than or equal to 50API, dark resistivity be less than 200ohmm and formation fracture factor of porosity is greater than 0.08% time, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{(3.2+1.1*\lg (AC)+\lg (\frac{RT*RT}{GR*GR}\left)\right)}{(3.1+\lg (GR\left)\right)}}+0.5$

Calculate the 4th rock-fabric number parameter ARFN ₄.

In addition, rock-fabric number parameter superpositing unit 34, specifically can:

According to formula:

ARFN＝ARFN ₁+ARFN ₂+ARFN ₃+ARFN ₄

Determine rock-fabric number parameter summation ARFN.

In addition, rock-fabric number parameter end value determining unit 35, specifically can:

When ARFN is less than or equal to the first preset value, determine that rock-fabric number parameter end value ARFN ' equals described ARFN;

Be greater than the first preset value at ARFN, and when dark resistivity is greater than 9000ohmm, determine that rock-fabric number parameter end value ARFN ' equals the second preset value;

Be greater than the first preset value at ARFN, and when dark resistivity is less than or equal to 9000ohmm, determine that rock-fabric number parameter end value ARFN ' equals the first preset value;

When ARFN is less than or equal to the 3rd preset value, determine that rock-fabric number parameter end value ARFN ' equals described 3rd preset value.

What deserves to be explained is, the specific implementation of the calculation element of a kind of low hole that the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter see embodiment of the method corresponding to above-mentioned Fig. 1 and Fig. 2, can repeat no more herein.

The calculation element of a kind of low hole that the embodiment of the present invention provides, vaughanite texture constituent of rock log parameter, four rock-fabric number parameters are asked for respectively, rock-fabric number parameter summation is determined in final superposition again, and then according to the dark resistivity size in rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value.The present invention is applicable to stratum and is mainly low hole, densification, and formation porosity is less than 10%, or with fracture development, sedimentary facies belt is mainly the pure limestone formation of open-platform or half open-platform sedimentary facies.Setting up basis compared to the rock-fabric number formula of prior art is that formation core porosity ranges is between 4%-40%, present invention, avoiding for domestic low hole, vaughanite rock formation, core porosity scope, between 0.5%-10%, accurately cannot determine the problem of structure number parameter of low hole, vaughanite rock by the rock-fabric number formula of prior art.

Apply specific embodiment in the present invention to set forth principle of the present invention and embodiment, the explanation of above embodiment just understands method of the present invention and core concept thereof for helping; Meanwhile, for one of ordinary skill in the art, according to thought of the present invention, all will change in specific embodiments and applications, in sum, this description should not be construed as limitation of the present invention.

Claims (16)

1. computing method for low hole, vaughanite texture constituent of rock log parameter, is characterized in that, comprising:

Obtain the log data of limestone rock formation to be measured; Described log data comprises: deep resistivity log curve, flushed zone resistivity logging trace, mud filtrate resistivity, gamma ray curve and interval transit time curve;

According to described deep resistivity log curve and flushed zone resistivity logging trace layer fracture porosity definitely;

According to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter;

Described first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter are superposed, determines rock-fabric number parameter summation;

According to the dark resistivity size in described rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value.

2. the computing method of low hole according to claim 1, vaughanite texture constituent of rock log parameter, is characterized in that, according to described deep resistivity log curve and flushed zone resistivity logging trace layer fracture porosity definitely, comprising:

Judge the size of the dark resistivity under synchronization in deep resistivity log curve and the flushed zone resistivity in flushed zone resistivity logging trace;

If dark resistivity is greater than flushed zone resistivity, then according to formula:

${\mathrm{\&Phi;}}_{frc}=\left(\frac{8.52253}{RXO}-\frac{8.242778}{RT}+0.00071236\right)\&times;Rmf$

Calculate formation fracture factor of porosity Φ _frc; Wherein RXO is flushed zone resistivity; RT is dark resistivity; Rmf is mud filtrate resistivity;

If dark resistivity is less than flushed zone resistivity, then according to formula:

${\mathrm{\&Phi;}}_{frc}=(\frac{1.99247}{RT}-\frac{0.992719}{RXO}+0.00031829)\&times;Rmf$

Calculate formation fracture factor of porosity Φ _frc.

3. the computing method of low hole according to claim 2, vaughanite texture constituent of rock log parameter, it is characterized in that, according to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, comprising:

According to the size of described natural gamma value size and formation fracture factor of porosity, determine described first rock-fabric number parameter;

The described size according to described natural gamma value size and formation fracture factor of porosity, determine described first rock-fabric number parameter, comprising:

If natural gamma value is greater than 50API, and formation fracture factor of porosity is less than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\&times;GR-1.125$

Calculate described first rock-fabric number parameter ARFN ₁; Wherein, RT is dark resistivity; GR is natural gamma value;

If natural gamma value is greater than 50API, and formation fracture factor of porosity is greater than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\&times;GR-0.625$

Calculate described first rock-fabric number parameter ARFN ₁.

4. the computing method of low hole according to claim 3, vaughanite texture constituent of rock log parameter, it is characterized in that, according to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, comprising:

According to size and the natural gamma value size of described interval transit time value, determine described second rock-fabric number parameter;

The described size according to described interval transit time value and natural gamma value size, determine described second rock-fabric number parameter, comprising:

If interval transit time value is less than or equal to 51us/f, and natural gamma value is less than or equal to 50API, according to formula:

${\mathrm{ARFN}}_2=e^{\frac{(2.9+\lg (AC)+3.2\&times;\lg (RT\&times;e^{\left(\frac{RT-200}{1800}\right)}\left)\right)}{(6.9+\lg (AC\left)\right)}}-3$

Calculate described second rock-fabric number parameter ARFN ₂; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

5. the computing method of low hole according to claim 4, vaughanite texture constituent of rock log parameter, it is characterized in that, according to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, comprising:

According to the size of described interval transit time value, natural gamma value size and dark resistivity size, determine described 3rd rock-fabric number parameter;

The described size according to described interval transit time value, natural gamma value size and dark resistivity size, determine described 3rd rock-fabric number parameter, comprising:

If interval transit time value is greater than 51us/f, natural gamma value is less than or equal to 50API, and dark resistivity is greater than 200ohmm, according to formula:

${\mathrm{ARFN}}_3=9.5-e^{\frac{(3.2+\lg (AC)+2*\lg (\frac{RT}{AC}\left)\right)}{(3.1+\lg (AC\left)\right)}}$

Calculate described 3rd rock-fabric number parameter ARFN ₃; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

6. the computing method of low hole according to claim 5, vaughanite texture constituent of rock log parameter, it is characterized in that, according to the size of the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, comprising:

According to the size of described interval transit time value, the size of natural gamma value size, deeply resistivity and formation fracture factor of porosity, determine described 4th rock-fabric number parameter;

The size of the described size according to described interval transit time value, natural gamma value size, deeply resistivity and formation fracture factor of porosity, determine described 4th rock-fabric number parameter, comprising:

If interval transit time value is greater than 51us/f, natural gamma value is less than or equal to 40API, and dark resistivity is less than 200ohmm and formation fracture factor of porosity is less than 0.08%, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{\left(3.2+1.1*\lg \left(AC\right)+\lg \left(\frac{RT*RT}{GR*GR}\right)\right)}{(3.1+\lg (GR\left)\right)}}$

Calculate described 4th rock-fabric number parameter ARFN ₄; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity; GR is natural gamma value;

If interval transit time value is greater than 51us/f, natural gamma value is greater than 40API and is less than or equal to 50API, and dark resistivity is less than 200ohmm and formation fracture factor of porosity is greater than 0.08%, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{\left(3.2+1.1*\lg \left(AC\right)+\lg \left(\frac{RT*RT}{GR*GR}\right)\right)}{(3.1+\lg (GR\left)\right)}}+0.5$

Calculate described 4th rock-fabric number parameter ARFN ₄.

7. the computing method of low hole according to claim 6, vaughanite texture constituent of rock log parameter, it is characterized in that, described first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter are superposed, determine rock-fabric number parameter summation, comprising:

According to formula:

ARFN＝ARFN ₁+ARFN ₂+ARFN ₃+ARFN ₄

Determine rock-fabric number parameter summation ARFN.

8. the computing method of low hole according to claim 7, vaughanite texture constituent of rock log parameter, it is characterized in that, according to the dark resistivity size in described rock-fabric number parameter summation and deep resistivity log curve, determine rock-fabric number parameter end value, comprising:

If ARFN is less than or equal to the first preset value, determine that described rock-fabric number parameter end value ARFN ' equals described ARFN;

If ARFN is greater than described first preset value, and dark resistivity is greater than 9000ohmm, determines that described rock-fabric number parameter end value ARFN ' equals the second preset value;

If ARFN is greater than described first preset value, and dark resistivity is less than or equal to 9000ohmm, determines that described rock-fabric number parameter end value ARFN ' equals described first preset value;

If ARFN is less than or equal to the 3rd preset value, determine that described rock-fabric number parameter end value ARFN ' equals described 3rd preset value.

9. a calculation element for low hole, vaughanite texture constituent of rock log parameter, is characterized in that, comprising:

Log data acquiring unit, for obtaining the log data of limestone rock formation to be measured; Described log data comprises: deep resistivity log curve, flushed zone resistivity logging trace, mud filtrate resistivity, gamma ray curve and interval transit time curve;

Formation fracture factor of porosity determining unit, for according to described deep resistivity log curve and flushed zone resistivity logging trace layer fracture porosity definitely;

Rock-fabric number parameter determination unit, for the size according to the interval transit time value in the size of the dark resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, formation fracture factor of porosity and interval transit time curve, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter;

Rock-fabric number parameter superpositing unit, for described first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter being superposed, determines rock-fabric number parameter summation;

Rock-fabric number parameter end value determining unit, for according to the dark resistivity size in described rock-fabric number parameter summation and deep resistivity log curve, determines rock-fabric number parameter end value.

10. the calculation element of low hole according to claim 9, vaughanite texture constituent of rock log parameter, is characterized in that, formation fracture factor of porosity determining unit, comprising:

Judge module, for judging the size of the dark resistivity under synchronization in deep resistivity log curve and the flushed zone resistivity in flushed zone resistivity logging trace;

First computing module, for when dark resistivity is greater than flushed zone resistivity, according to formula:

${\mathrm{\&Phi;}}_{frc}=(\frac{8.52253}{RXO}-\frac{8.242778}{RT}+0.00071236)\&times;Rmf$

Calculate formation fracture factor of porosity Φ _frc; Wherein RXO is flushed zone resistivity; RT is dark resistivity; Rmf is mud filtrate resistivity;

Second computing module, for when dark resistivity is less than flushed zone resistivity, according to formula:

${\mathrm{\&Phi;}}_{frc}=(\frac{1.99247}{RT}-\frac{0.992719}{RXO}+0.00031829)\&times;Rmf$

Calculate formation fracture factor of porosity Φ _frc.

The calculation element of 11. low holes according to claim 10, vaughanite texture constituent of rock log parameter, is characterized in that, described rock-fabric number parameter determination unit, comprising:

First rock-fabric number parameter calculating module, for the size according to described natural gamma value size and formation fracture factor of porosity, determines described first rock-fabric number parameter;

Described first rock-fabric number parameter calculating module, specifically for:

50API is greater than in natural gamma value, and when formation fracture factor of porosity is less than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\&times;GR-1.125$

Calculate described first rock-fabric number parameter ARFN ₁; Wherein, RT is dark resistivity; GR is natural gamma value;

50API is greater than in natural gamma value, and when formation fracture factor of porosity is greater than 0.08%, according to formula:

${\mathrm{ARFN}}_1=\frac{RT}{100}+0.0225\&times;GR-0.625$

Calculate described first rock-fabric number parameter ARFN ₁.

The calculation element of 12. low holes according to claim 11, vaughanite texture constituent of rock log parameter, is characterized in that, described rock-fabric number parameter determination unit, comprising:

Second rock-fabric number parameter calculating module, for according to the size of described interval transit time value and natural gamma value size, determines described second rock-fabric number parameter;

Described second rock-fabric number parameter calculating module, specifically for:

51us/f is less than or equal in interval transit time value, and when natural gamma value is less than or equal to 50API, according to formula:

${\mathrm{ARFN}}_2=e^{\frac{\left(2.9+\lg \left(AC\right)+3.2\&times;\lg \left(RT\&times;e^{\left(\frac{RT-200}{1800}\right)}\right)\right)}{(6.9+\lg (AC\left)\right)}}-3$

Calculate described second rock-fabric number parameter ARFN ₂; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

The calculation element of 13. low holes according to claim 12, vaughanite texture constituent of rock log parameter, is characterized in that, described rock-fabric number parameter determination unit, comprising:

3rd rock-fabric number parameter calculating module, for the size according to described interval transit time value, natural gamma value size and dark resistivity size, determines described 3rd rock-fabric number parameter;

Described 3rd rock-fabric number parameter calculating module, specifically for::

Be greater than 51us/f in interval transit time value, natural gamma value is less than or equal to 50API, and when dark resistivity is greater than 200ohmm, according to formula:

${\mathrm{ARFN}}_3=9.5-e^{\frac{\left(3.2+\lg \left(AC\right)+2*\lg \left(\frac{RT}{AC}\right)\right)}{(3.1+\lg (AC\left)\right)}}$

Calculate described 3rd rock-fabric number parameter ARFN ₃; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity.

The calculation element of 14. low holes according to claim 13, vaughanite texture constituent of rock log parameter, is characterized in that, described rock-fabric number parameter determination unit, comprising:

4th rock-fabric number parameter calculating module, for the size of the size according to described interval transit time value, natural gamma value size, dark resistivity and formation fracture factor of porosity, determines described 4th rock-fabric number parameter;

Described 4th rock-fabric number parameter calculating module, specifically for:

Be greater than 51us/f in interval transit time value, natural gamma value is less than or equal to 40API, dark resistivity be less than 200ohmm and formation fracture factor of porosity is less than 0.08% time, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{\left(3.2+1.1*\lg \left(AC\right)+\lg \left(\frac{RT*RT}{GR*GR}\right)\right)}{(3.1+\lg (GR\left)\right)}}$

Calculate described 4th rock-fabric number parameter ARFN ₄; Wherein, e is the nature truth of a matter; AC is interval transit time value; RT is dark resistivity; GR is natural gamma value;

Be greater than 51us/f in interval transit time value, natural gamma value is greater than 40API and is less than or equal to 50API, dark resistivity be less than 200ohmm and formation fracture factor of porosity is greater than 0.08% time, according to formula:

${\mathrm{ARFN}}_4=7.2-e^{\frac{\left(3.2+1.1*\lg \left(AC\right)+\lg \left(\frac{RT*RT}{GR*GR}\right)\right)}{(3.1+\lg (GR\left)\right)}}+0.5$

Calculate described 4th rock-fabric number parameter ARFN ₄.

The calculation element of 15. low holes according to claim 14, vaughanite texture constituent of rock log parameter, is characterized in that, described rock-fabric number parameter superpositing unit, specifically for:

According to formula:

ARFN＝ARFN ₁+ARFN ₂+ARFN ₃+ARFN ₄

Determine rock-fabric number parameter summation ARFN.

The calculation element of 16. low holes according to claim 15, vaughanite texture constituent of rock log parameter, is characterized in that, described rock-fabric number parameter end value determining unit, specifically for:

When ARFN is less than or equal to the first preset value, determine that described rock-fabric number parameter end value ARFN ' equals described ARFN;

Be greater than described first preset value at ARFN, and when dark resistivity is greater than 9000ohmm, determine that described rock-fabric number parameter end value ARFN ' equals the second preset value;

Be greater than described first preset value at ARFN, and when dark resistivity is less than or equal to 9000ohmm, determine that described rock-fabric number parameter end value ARFN ' equals described first preset value;

When ARFN is less than or equal to the 3rd preset value, determine that described rock-fabric number parameter end value ARFN ' equals described 3rd preset value.