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

Method and device for calculating logging parameters of low-porosity and compact limestone rock structure components Download PDF

Info

Publication number
CN105549114B
CN105549114B CN201510876133.XA CN201510876133A CN105549114B CN 105549114 B CN105549114 B CN 105549114B CN 201510876133 A CN201510876133 A CN 201510876133A CN 105549114 B CN105549114 B CN 105549114B
Authority
CN
China
Prior art keywords
rock
mrow
number parameter
fabric number
arfn
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510876133.XA
Other languages
Chinese (zh)
Other versions
CN105549114A (en
Inventor
李昌
乔占峰
沈安江
李扬
潘立银
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201510876133.XA priority Critical patent/CN105549114B/en
Publication of CN105549114A publication Critical patent/CN105549114A/en
Application granted granted Critical
Publication of CN105549114B publication Critical patent/CN105549114B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

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

Low hole, the computational methods of compact limestone texture constituent of rock log parameter and device
Technical field
The present invention relates to petroleum geology exploration and logging technique field, more particularly to a kind of low hole, compact limestone rock knot The computational methods and device of structure component log parameter.
Background technology
At present, recognize that limestone texture constituent of rock is always a world-famous puzzle using well-log information, and how to be had Effect log parameter is the key for improving accuracy of identification.Currently obtain effective log parameter research method and thinking be broadly divided into Under several steps：1) rock-electricity relationship analysis：First to the core hole of core hole, sample is obtained, sample is ground, obtained Thin slice is taken, thin section identification is carried out using microscope, texture constituent of rock information is obtained, and the playback of thin slice depth is deep to well logging Degree, then carries out texture constituent of rock and well logging correlation analysis, it is determined that the log parameter best with texture constituent of rock correlation (Logging Curves, synthetic parameters curve etc.), sets up the logging characteristic parameters Sample Storehouse of texture constituent of rock.2) technology is known Not：It is main to use the mathematical algorithms, sample of the neural network to foundation such as neutral net, clustering, discriminant analysis, fuzzy mathematics This log parameter storehouse is trained, and obtains artificial intelligence；Clustering methodology sets up well logging 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 different texture constituent of rock The membership function of well logging classification mode.3) technology application：By the research to studying area's core hole, it is determined that optimal well logging ginseng Number, chooses most suitable identification technology, carries out texture constituent of rock identification to non-core hole, and tested with other core holes in work area Demonstrate,prove recognition effect.
Study less at home on texture constituent of rock Logging Identification Method in the prior art.And abroad specifically designed for The research of texture constituent of rock well logging recognition is usually to carry out rock physicses three according to limestone rock particle or crystallite dimension to classify, the Size is defined as grain structure more than 100 μm by one class, and size is defined as micrite grain structure by Equations of The Second Kind at 20 μm -100 μm, It is main structure (including part micrite grain structure, particle micrite structure, micrite that size is defined as micrite less than 20 μm by the 3rd class Structure), rock-fabric number (Rock Fabric are established by core porosity, permeability, water saturation parameter Number) calculation formula, calculates porosity, permeability and water saturation by logging well and replaces core parameters, realize rock knot Structure component well logging recognition.This method is set up under following two geological conditions：(1) stratum is mainly with intergranular, intracrystalline pore or dissolution pore Based on, crack agensis；(2) stratum, in the absence of pure water layer, meets formation rock knot under the conditions of the two on gas-water interface Structure component well logging recognition can obtain fine recognition effect.
This method of prior art effectively solves the technical bottleneck of texture constituent of rock well logging recognition, but what it was applied There is also deficiency for rock-fabric number formula：Its set up basis be formation core porosity ranges between 4%-40%, and for Domestic low hole, compact limestone rock formation, core porosity scope pass through the rock knot of prior art between 0.5%-10% Structure number formula can not accurately determine low hole, the structure number parameter of compact limestone rock.
The content of the invention
The embodiment of the present invention provides a kind of low hole, the computational methods of compact limestone texture constituent of rock log parameter and dress Put, core porosity scope has not been suitable between 0.5%-10% to solve rock-fabric number formula of the prior art Low hole, compact limestone rock formation, the problem of causing accurately determine low hole, the structure number parameter of compact limestone rock.
In order to achieve the above object, the present invention is adopted the following technical scheme that：
A kind of low hole, the computational methods of compact limestone texture constituent of rock log parameter, including：
Obtain the log data of limestone rock formation to be measured；The log data includes：Deep resistivity log curve, flushing Strip resistance rate log, mud filtrate resistivity, gamma ray curve and interval transit time curve；
According to the deep resistivity log curve and flushed zone resistivity log definitely layer fracture porosity；
Deep resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, The size of interval transit time value in the size and interval transit time curve of formation fracture porosity, determines that the first rock-fabric number is joined Number, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter；
By the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and Four rock-fabric number parameters are superimposed, and determine rock-fabric number parameter summation；
Deep resistivity size in the rock-fabric number parameter summation and deep resistivity log curve, determines rock Structure number parameter end value.
Specifically, according to the deep resistivity log curve and flushed zone resistivity log definitely layer fracture pore Degree, including：
Judge in the deep resistivity and flushed zone resistivity log under synchronization in deep resistivity log curve The size of flushed zone resistivity；
If deep resistivity is more than flushed zone resistivity, according to formula：
Calculate formation fracture porosity Φfrc；Wherein RXO is flushed zone resistivity；RT is deep resistivity；Rmf filters for mud Liquid resistivity；
If deep resistivity is less than flushed zone resistivity, according to formula：
Calculate formation fracture porosity Φfrc。
Specifically, the natural gamma in deep resistivity size, gamma ray curve in deep resistivity log curve It is worth the size of size, the size of formation fracture porosity and the interval transit time value in interval transit time curve, determines the first rock Structure number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, including：
According to the natural gamma value size and the size of formation fracture porosity, the first rock-fabric number ginseng is determined Number；
It is described according to the natural gamma value size and the size of formation fracture porosity, determine first rock texture Number parameter, including：
If natural gamma value is more than 50API, and formation fracture porosity is less than 0.08%, according to formula：
Calculate the first rock-fabric number parameter ARFN1；Wherein, RT is deep resistivity；GR is natural gamma value；
If natural gamma value is more than 50API, and formation fracture porosity is more than 0.08%, according to formula：
Calculate the first rock-fabric number parameter ARFN1。
In addition, the natural gamma value in deep resistivity size, gamma ray curve in deep resistivity log curve The size of interval transit time value in size, the size of formation fracture porosity and interval transit time curve, determines the first rock knot Structure number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, including：
According to the size and natural gamma value size of the interval transit time value, the second rock-fabric number parameter is determined；
The size and natural gamma value size according to the interval transit time value, determines the second rock-fabric number ginseng Number, including：
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：
Calculate the second rock-fabric number parameter ARFN2；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is Deep resistivity.
In addition, the natural gamma value in deep resistivity size, gamma ray curve in deep resistivity log curve The size of interval transit time value in size, the size of formation fracture porosity and interval transit time curve, determines the first rock knot Structure number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, including：
According to the size of the interval transit time value, natural gamma value size and deep resistivity size, the 3rd rock is determined Stone structure number parameter；
It is described according to the size of the interval transit time value, natural gamma value size and deep resistivity size, determine described the Three rock-fabric number parameters, including：
If interval transit time value is more than 51us/f, natural gamma value is less than or equal to 50API, and deep resistivity is more than 200ohmm, According to formula：
Calculate the 3rd rock-fabric number parameter ARFN3；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is Deep resistivity.
In addition, the natural gamma value in deep resistivity size, gamma ray curve in deep resistivity log curve The size of interval transit time value in size, the size of formation fracture porosity and interval transit time curve, determines the first rock knot Structure number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, including：
According to the big of the size of the interval transit time value, natural gamma value size, deep resistivity and formation fracture porosity It is small, determine the 4th rock-fabric number parameter；
It is described according to the size of the interval transit time value, natural gamma value size, deep resistivity and formation fracture porosity Size, determine the 4th rock-fabric number parameter, including：
If interval transit time value be more than 51us/f, natural gamma value be less than or equal to 40API, deep resistivity be less than 200ohmm and Formation fracture porosity is less than 0.08%, according to formula：
Calculate the 4th rock-fabric number parameter ARFN4；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is Deep resistivity；GR is natural gamma value；
If interval transit time value is more than 51us/f, natural gamma value is more than 40API and less than or equal to 50API, and deep resistivity is small It is more than 0.08% in 200ohmm and formation fracture porosity, according to formula：
Calculate the 4th rock-fabric number parameter ARFN4。
Specifically, 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 superposition, determine rock-fabric number parameter summation, including：
According to formula：
ARFN=ARFN1+ARFN2+ARFN3+ARFN4
Determine rock-fabric number parameter summation ARFN.
Specifically, the deep resistivity size in the rock-fabric number parameter summation and deep resistivity log curve, Rock-fabric number parameter end value is determined, including：
If ARFN is less than or equal to the first preset value, determine that the rock-fabric number parameter end value ARFN ' is equal to described ARFN；
If ARFN is more than first preset value, and deep resistivity is more than 9000ohmm, determines the rock-fabric number ginseng Number end value ARFN ' is equal to the second preset value；
If ARFN is more than first preset value, and deep resistivity is less than or equal to 9000ohmm, determines the rock texture Number parameter end value ARFN ' is equal to first preset value；
If ARFN is less than or equal to the 3rd preset value, determine that the rock-fabric number parameter end value ARFN ' is equal to described the Three preset values.
A kind of low hole, the computing device of compact limestone texture constituent of rock log parameter, including：
Log data acquiring unit, the log data for obtaining limestone rock formation to be measured；The log data includes： Deep resistivity log curve, flushed zone resistivity log, mud filtrate resistivity, gamma ray curve and interval transit time Curve；
Formation fracture porosity determining unit, for according to the deep resistivity log curve and flushed zone resistivity well logging Curve definitely layer fracture porosity；
Rock-fabric number parameter determination unit, for the deep resistivity size in deep resistivity log curve, nature The interval transit time value in natural gamma value size, the size of formation fracture porosity and interval transit time curve in gamma curve Size, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock Stone structure number parameter；
Rock-fabric number parameter superpositing unit, for the first rock-fabric number parameter, the second rock-fabric number to be joined Number, the 3rd rock-fabric number parameter and the superposition of the 4th rock-fabric number parameter, determine rock-fabric number parameter summation；
Rock-fabric number parameter end value determining unit, for according to the rock-fabric number parameter summation and deep resistivity Deep resistivity size in log, determines rock-fabric number parameter end value.
Specifically, the formation fracture porosity determining unit, including：
Judge module, for judging deep resistivity and flushed zone resistivity under synchronization in deep resistivity log curve The size of flushed zone resistivity in log；
First computing module, for when deep resistivity is more than flushed zone resistivity, according to formula：
Calculate formation fracture porosity Φfrc；Wherein RXO is flushed zone resistivity；RT is deep resistivity；Rmf filters for mud Liquid resistivity；
Second computing module, for when deep resistivity is less than flushed zone resistivity, according to formula：
Calculate formation fracture porosity Φfrc。
Specifically, the rock-fabric number parameter determination unit, including：
First rock-fabric number parameter calculating module, for according to the natural gamma value size and formation fracture porosity Size, determine the first rock-fabric number parameter；
The first rock-fabric number parameter calculating module, specifically for：
When natural gamma value is less than 0.08% more than 50API, and formation fracture porosity, according to formula：
Calculate the first rock-fabric number parameter ARFN1；Wherein, RT is deep resistivity；GR is natural gamma value；
When natural gamma value is more than 0.08% more than 50API, and formation fracture porosity, according to formula：
Calculate the first rock-fabric number parameter ARFN1。
In addition, the rock-fabric number parameter determination unit, including：
Second rock-fabric number parameter calculating module is big for the size and natural gamma value according to the interval transit time value It is small, determine the second rock-fabric number parameter；
The second rock-fabric number parameter calculating module, specifically for：
When interval transit time value is less than or equal to 50API less than or equal to 51us/f, and natural gamma value, according to formula：
Calculate the second rock-fabric number parameter ARFN2；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is Deep resistivity.
In addition, the rock-fabric number parameter determination unit, including：
3rd rock-fabric number parameter calculating module is big for size, the natural gamma value according to the interval transit time value Small and deep resistivity size, determines the 3rd rock-fabric number parameter；
The 3rd rock-fabric number parameter calculating module, specifically for：
It is more than 51us/f in interval transit time value, natural gamma value is less than or equal to 50API, and deep resistivity is more than 200ohmm When, according to formula：
Calculate the 3rd rock-fabric number parameter ARFN3；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is Deep resistivity.
In addition, the rock-fabric number parameter determination unit, including：
4th rock-fabric number parameter calculating module is big for size, the natural gamma value according to the interval transit time value The size of small, deep resistivity and formation fracture porosity, determines the 4th rock-fabric number parameter；
The 4th rock-fabric number parameter calculating module, specifically for：
Interval transit time value be more than 51us/f, natural gamma value be less than or equal to 40API, deep resistivity be less than 200ohmm and When formation fracture porosity is less than 0.08%, according to formula：
Calculate the 4th rock-fabric number parameter ARFN4；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is Deep resistivity；GR is natural gamma value；
It is more than 51us/f in interval transit time value, natural gamma value is more than 40API and less than or equal to 50API, and deep resistivity is small When 200ohmm and formation fracture porosity are more than 0.08%, according to formula：
Calculate the 4th rock-fabric number parameter ARFN4。
In addition, the rock-fabric number parameter superpositing unit, specifically for：
According to formula：
ARFN=ARFN1+ARFN2+ARFN3+ARFN4
Determine rock-fabric number parameter summation ARFN.
In addition, the rock-fabric number parameter end value determining unit, specifically for：
When ARFN is less than or equal to the first preset value, determine that the rock-fabric number parameter end value ARFN ' is equal to described ARFN；
When ARFN is more than 9000ohmm more than first preset value, and deep resistivity, the rock-fabric number is determined Parameter end value ARFN ' is equal to the second preset value；
When ARFN is less than or equal to 9000ohmm more than first preset value, and deep resistivity, the rock knot is determined Structure number parameter end value ARFN ' is equal to first preset value；
When ARFN is less than or equal to three preset values, determine that the rock-fabric number parameter end value ARFN ' is equal to described 3rd preset value.
A kind of low hole provided in an embodiment of the present invention, the computational methods of compact limestone texture constituent of rock log parameter and dress Put, four rock-fabric number parameters have been asked for respectively, then final superposition determines rock-fabric number parameter summation, and then according to described Deep resistivity size in rock-fabric number parameter summation and deep resistivity log curve, determines that rock-fabric number parameter is final Value.The present invention is mainly low hole, densification suitable for stratum, and formation porosity is less than 10%, or with fracture development, sedimentary facies belt The predominantly pure limestone formation of open-platform or half open-platform sedimentary facies.Rock-fabric number formula compared to prior art is built Vertical basis be formation core porosity ranges between 4%-40%, present invention, avoiding for domestic low hole, compact limestone rock Stone reservoir, core porosity scope, can not be accurately true by the rock-fabric number formula of prior art between 0.5%-10% The problem of fixed low hole, structure number parameter of compact limestone rock.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this Some embodiments of invention, for those of ordinary skill in the art, without having to pay creative labor, may be used also To obtain other accompanying drawings according to these accompanying drawings.
Fig. 1 is a kind of one kind provided in an embodiment of the present invention low hole provided in an embodiment of the present invention, compact limestone rock knot The flow chart one of the computational methods of structure component log parameter；
Fig. 2 is a kind of one kind provided in an embodiment of the present invention low hole provided in an embodiment of the present invention, compact limestone rock knot The flowchart 2 of the computational methods of structure component log parameter；
Fig. 3 is core porosity, rock core water containing saturability and the core wafer analysis rock texture in the embodiment of the present invention Component cross plot；
Fig. 4 is interval transit time-depth resistivity, natural gamma-depth resistivity and the core wafer point in the embodiment of the present invention Analyse texture constituent of rock cross plot；
Fig. 5 is interval transit time-depth resistivity of four kinds of different models in the embodiment of the present invention, natural gamma-depth resistivity Texture constituent of rock cross plot is analyzed with core wafer；
Fig. 6 be interval transit time-depth resistivity of the first model in the embodiment of the present invention, natural gamma-depth resistivity and Core wafer analyzes texture constituent of rock cross plot；
Fig. 7 be interval transit time-depth resistivity of second model in the embodiment of the present invention, natural gamma-depth resistivity and Core wafer analyzes texture constituent of rock cross plot；
Fig. 8 be interval transit time-depth resistivity of the third model in the embodiment of the present invention, natural gamma-depth resistivity and Core wafer analyzes texture constituent of rock cross plot；
Fig. 9 be interval transit time-depth resistivity of the 4th kind of model in the embodiment of the present invention, natural gamma-depth resistivity and Core wafer analyzes texture constituent of rock cross plot；
Figure 10 regards rock-fabric number well logging for 201 wells in 161 wells, the middle ancient times 51 and tower in the tower in the embodiment of the present invention and calculated Design sketch；
Figure 11 is a kind of low hole provided in an embodiment of the present invention, the calculating dress of compact limestone texture constituent of rock log parameter The structural representation one put；
Figure 12 is a kind of low hole provided in an embodiment of the present invention, the calculating dress of compact limestone texture constituent of rock log parameter The structural representation two put.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made Embodiment, belongs to the scope of protection of the invention.
As shown in figure 1, a kind of low hole provided in an embodiment of the present invention, the meter of compact limestone texture constituent of rock log parameter Calculation method, including：
Step 101, the log data for obtaining limestone rock formation to be measured.
Wherein, log data includes：Deep resistivity log curve, flushed zone resistivity log, F resistance Rate, gamma ray curve and interval transit time curve.
Step 102, according to deep resistivity log curve and flushed zone resistivity log definitely layer fracture porosity.
Step 103, the deep resistivity size in deep resistivity log curve, the natural gamma in gamma ray curve It is worth the size of size, the size of formation fracture porosity and the interval transit time value in interval transit time curve, determines the first rock Structure 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 with And the 4th rock-fabric number parameter superposition, determine rock-fabric number parameter summation.
Step 105, the deep resistivity size in rock-fabric number parameter summation and deep resistivity log curve, it is determined that Rock-fabric number parameter end value.
A kind of low hole provided in an embodiment of the present invention, the computational methods of compact limestone texture constituent of rock log parameter, point Four rock-fabric number parameters have not been asked for, then final superposition determines rock-fabric number parameter summation, and then according to rock texture Deep resistivity size in number parameter summation and deep resistivity log curve, determines rock-fabric number parameter end value.The present invention It is mainly low hole, densification suitable for stratum, formation porosity is less than 10%, or with fracture development, sedimentary facies belt is mainly out Wealthy tableland or the pure limestone formation of half open-platform sedimentary facies.Rock-fabric number formula compared to prior art sets up basis and is Formation core porosity ranges are between 4%-40%, and present invention, avoiding for domestic low hole, compact limestone rock formation, rock Heart porosity ranges can not accurately determine low hole, cause between 0.5%-10% by the rock-fabric number formula of prior art The problem of structure number parameter of close limestone rock.
A more detailed embodiment is set forth below, to cause those skilled in the art to be better understood by the present invention. As shown in Fig. 2 a kind of low hole provided in an embodiment of the present invention, the computational methods of compact limestone texture constituent of rock log parameter, Including：
Step 201, the log data for obtaining limestone rock formation to be measured.
Wherein, log data includes：Deep resistivity log curve, flushed zone resistivity log, F resistance Rate, gamma ray curve and interval transit time curve.
The embodiment of the present invention can be combined herein to be applied to be said exemplified by the Ordovician strata of In The Central Tarim Area It is bright：TZ area Ordovician strata is the Marine source rock that a set of open-platform is deposited, and lithology is limestone.The set stratum Texture constituent of rock has 5 kinds：Micrite structure, particle micrite structure, micrite grain structure, grain structure and reef structure, this area Stratum is fine and close, crack comparative development, and above geological conditions is adapted to the application of the present invention.Parameter acquiring：Logged well using CLS-5700 Series instrument measurement stratum natural gamma GR, interval transit time DT and resistivity RT, resistivity measurement, mud are carried out to F The Ω .m (ohm meter)/84 DEG C of filtrate resistivity 0.15.
Deep resistivity log curve herein, flushed zone resistivity log can by Formation Resistivity Measurement come Arrive.
Interval transit time curve can include formation fluid interval transit time value and rock matrix interval transit time value.Wherein, stratum Fluid acoustic time difference value DTfIt is general to use the μ s/f of theoretical value 189, rock matrix interval transit time value DTmaIt is general main according to stratum Lithology is set, and such as stratum is limestone formation, DTmaUsing the μ s/f of theoretical value 47, dolomite stratum then uses the μ s/f of theoretical value 42, If mixing lithology then takes both weighted calculation average values according to geological condition, the value is between 42-47 μ s/ft, typically The interval transit time on stratum is in 40-90 μ s/f.
Natural gamma GR in above-mentioned gamma ray curve is mainly stratum natural gamma value (GRma), and formation fluid Natural gamma is almost 0.Wherein GRmaGeneral to be set according to stratum principal lithologic, such as stratum is limestone formation, then GRmaUsing Theoretical value 8API, mudstone stratum then uses theoretical value 230API, if mixing lithology, then take both to add according to geological condition Power calculates average value, and the value is between 0-70API, and the natural gamma on general stratum can be between 10-200API.
Above-mentioned deep resistivity RT can include formation fluid resistivity (RTf) and rock matrix resistivity (RTma).Its In, formation fluid is mainly water, RTfGenerally 0.01~1ohmm, RTmaGeneral dense, non-conductive (in addition to pyrite), Typically can be in more than 10000ohmm, the resistivity on general stratum is between 20-200000ohmm.
The log of the present invention can be obtained using prior art device, such as CLS-3700 logging instruments, CLS- 5700 logging instruments etc..
Step 202, according to deep resistivity log curve and flushed zone resistivity log definitely layer fracture porosity.
Specifically, needing to judge the deep resistivity and flushing strip resistance under synchronization in deep resistivity log curve herein The size of flushed zone resistivity in rate log；
If deep resistivity is more than flushed zone resistivity, according to formula：
Calculate formation fracture porosity Φfrc；Wherein RXO is flushed zone resistivity；RT is deep resistivity；Rmf filters for mud Liquid resistivity；
If deep resistivity is less than flushed zone resistivity, according to formula：
Calculate formation fracture porosity Φfrc。
For example by taking 161 core holes in tower as an example, depth of stratum is that 4384.1m is micrite grainstone, this layer of log measurement Deep resistivity is 547.ohmm, and flushed zone resistivity is 501.74ohmm, mud filtrate resistivity RmfSubstituted into equal to 0.15ohmm Above formula, calculates fracture porosity and is equal to 0.014%.
In following steps 203- steps 207, the calculating of the rock-fabric number parameter of the present invention will be illustrated.
Rock-fabric number computing technique starts from Lucia in 2005 and carries out rock according to carbonate rock particle or crystallite dimension Physics three is classified, and the definition that size is more than 100 μm by the first kind is grain structure, and Equations of The Second Kind is by size between 20 μm -100 μm Definition be micrite grain structure, it is that micrite is main structure (including part micrite that size is less than 20 μm of definition by the 3rd class Kernel structure, particle micrite structure, micrite structure), pass through core porosity, the foundation (as shown in Figure 3) of water saturation cross plot Regard rock-fabric number (Apparent Rock Fabric Number) calculation formula as follows：
Φ represents formation porosity (decimal unit), S in formulawRepresent stratum water saturation (decimal unit).
The acquisition that wherein formation porosity and water saturation can be calculated using well logging, so as to realize that well logging calculates rock Structure number, according to rock-fabric number size, it is possible to achieve texture constituent of rock well logging recognition.
Work as ARFN>4, represent micrite structure；
When ARFN is between 2.5~4, particle micrite structure and powder micrite structure are represented；
When ARFN is between 1.5~2.5, micrite particle or fine powder crystal structure are represented；
When ARFN is between 0.5~1.5, grain structure or middle coarse structure are represented.
Carbonate formation (porosity is generally more than 4%) application effect is preferable abroad for above-mentioned formula, but for the country Carbonate reservoir hole is low, particularly porosity 0.5%-4% reservoir, in this case Lucia formula do not applied to, such as Fruit directly calculates rock-fabric number using the formula, and rock-fabric number scope is 0-1000, cannot be used for recognizing structure group Point, such as when stratum is micrite, micrite porosity is 1% (decimal unit is 0.01), and micrite contains water saturation Spend for 100% (decimal unit is 1), substitute into above formula and calculate ARFN=0.0032, the numerical value exceeds 0.5-4 scopes, no The identification of texture constituent of rock can be carried out.
Due to the low hole in stratum, densification and fracture development, its texture constituent of rock logging character has very big plyability, different Texture constituent of rock logging character is difficult to distinguish (such as shown in Fig. 4), according to stratum geologic feature, and point 4 kinds of typical models are (as schemed Shown in 5), set up regard rock-fabric number formula respectively, then by the formula combinations of 4 models, obtain final rock-fabric number and calculate Formula.Mode is as follows：
Step 203, according to natural gamma value size and the size of formation fracture porosity, determine the first rock-fabric number join Number.
Herein, if natural gamma value is more than 50API (such as shown in Fig. 6), and formation fracture porosity is less than 0.08%, root According to formula：
Calculate the first rock-fabric number parameter ARFN1；Wherein, RT is deep resistivity；GR is natural gamma value；
If natural gamma value is more than 50API, and formation fracture porosity is more than 0.08%, according to formula：
Calculate the first rock-fabric number parameter ARFN1。
Step 204, the size according to 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 Fig. 7 institutes Show), according to formula：
Calculate the second rock-fabric number parameter ARFN2；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep electricity Resistance rate.
Step 205, the size according to interval transit time value, natural gamma value size and deep resistivity size, determine the 3rd rock Stone structure number parameter.
Herein, if interval transit time value is more than 51us/f, natural gamma value is less than or equal to 50API (such as shown in Fig. 8), and deep Resistivity is more than 200ohmm, according to formula：
Calculate the 3rd rock-fabric number parameter ARFN3；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep electricity Resistance rate.
Step 206, the size according to interval transit time value, natural gamma value size, deep resistivity and formation fracture porosity Size, determine the 4th rock-fabric number parameter.
Herein, if interval transit time value is more than 51us/f, natural gamma value is less than or equal to 40API, and deep resistivity is less than 200ohmm (such as shown in Fig. 9) and formation fracture porosity are less than 0.08%, according to formula：
Calculate the 4th rock-fabric number parameter ARFN4；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep electricity Resistance rate；GR is natural gamma value.
If interval transit time value is more than 51us/f, natural gamma value is more than 40API and less than or equal to 50API, and deep resistivity is small It is more than 0.08% in 200ohmm and formation fracture porosity, according to formula：
Calculate the 4th rock-fabric number parameter ARFN4。
Step 207, by the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter with And the 4th rock-fabric number parameter superposition, determine rock-fabric number parameter summation.
Herein specifically can be according to formula：
ARFN=ARFN1+ARFN2+ARFN3+ARFN4
Determine rock-fabric number parameter summation ARFN.
The interval transit time for for example having stratum is 48.83us/f, and deep resistivity is 547.5ohmm, and natural gamma is 8.7API, ARFN=2.31 then is obtained using aforesaid way, and actual rock core rock-fabric number is 2.1, it is basically identical with rock core structure number (as shown in Figure 10).
Step 208, the deep resistivity size in rock-fabric number parameter summation and deep resistivity log curve, it is determined that Rock-fabric number parameter end value.
Herein, if ARFN is less than or equal to one first preset value, determine that the rock-fabric number parameter end value ARFN ' is equal to The ARFN.
Herein, such as the first preset value can be 4, then when ARFN is less than or equal to 4, such as 3, then rock-fabric number parameter End value ARFN ' is equal to 3.
If ARFN is more than the first preset value, and deep resistivity is more than 9000ohmm, determines rock-fabric number parameter end value ARFN ' is equal to the second preset value.
The first preset value herein can be 4；Second preset value can be 4 or 4.5.The purpose of each preset value is herein The ARFN curve number ranges that limitation is calculated ensure between 1-4.5, have exceeded and have fallen with regard to cut-off.
If ARFN is more than the first preset value, and deep resistivity is less than or equal to 9000ohmm, determines rock-fabric number parameter most Final value ARFN ' is equal to the first preset value.
For example, the first preset value is 4, then ARFN is more than 4, and deep resistivity is when being less than or equal to 9000ohmm, determines rock Structure number parameter end value ARFN ' is equal to 4.
In addition, when ARFN is less than or equal to three preset values, determining that the rock-fabric number parameter end value ARFN ' is equal to 3rd preset value.
For example, the 3rd preset value can be 1, then when ARFN is less than or equal to 1, the rock-fabric number parameter end value ARFN ' is equal to 1.
In embodiments of the present invention, ARFN number range typically can be 1 to 4.5 between：
(1) when 0 ＜ ARFN≤1, rock texture is defined as micrite structural constituent phase, and lithology is shale micrite, micrite Limestone breathes out partial particulate micrite.
(2) when 1 ＜ ARFN≤2, rock texture is defined as particle micrite structural constituent phase, and lithology is mainly particle micrite ash Rock, bondstone and micrite support micrite grainstone,；
(3) when 2 ＜ ARFN≤3, rock texture is defined as micrite grain structure component phase, and lithology is predominantly particle-supported Micrite grainstone；
(4) when 3 ＜ ARFN≤4, rock texture is defined as grain structure component phase, and lithology is mainly grainstone, dissolution pore Coral rock；
(5) when 4 ＜ ARFN≤4.5, rock texture is defined as reef structural constituent phase, and lithology is fine and close biohermal limestone, caused Close grainstone.
The present invention can obtain good effect on In The Central Tarim Area Ordovician system well compactness stratum.
A kind of low hole provided in an embodiment of the present invention, the computational methods of compact limestone texture constituent of rock log parameter, point Four rock-fabric number parameters have not been asked for, then final superposition determines rock-fabric number parameter summation, and then according to rock texture Deep resistivity size in number parameter summation and deep resistivity log curve, determines rock-fabric number parameter end value.The present invention It is mainly low hole, densification suitable for stratum, formation porosity is less than 10%, or with fracture development, sedimentary facies belt is mainly out Wealthy tableland or the pure limestone formation of half open-platform sedimentary facies.Rock-fabric number formula compared to prior art sets up basis and is Formation core porosity ranges are between 4%-40%, and present invention, avoiding for domestic low hole, compact limestone rock formation, rock Heart porosity ranges can not accurately determine low hole, cause between 0.5%-10% by the rock-fabric number formula of prior art The problem of structure number parameter of close limestone rock.
Corresponding to the embodiment of the method shown in above-mentioned Fig. 1, Fig. 2, a kind of low hole provided in an embodiment of the present invention, compact limestone The computing device of texture constituent of rock log parameter, as shown in figure 11, including：
Log data acquiring unit 31, can obtain the log data of limestone rock formation to be measured.
Log data includes：Deep resistivity log curve, flushed zone resistivity log, mud filtrate resistivity, from Right gamma curve and interval transit time curve.
Formation fracture porosity determining unit 32, can log well bent according to deep resistivity log curve and flushed zone resistivity Line definitely layer fracture porosity.
Rock-fabric number parameter determination unit 33, deep resistivity size that can be in deep resistivity log curve, from The interval transit time in natural gamma value size, the size of formation fracture porosity and interval transit time curve in right gamma curve The size of value, determines 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 34, can by the first rock-fabric number parameter, the second rock-fabric number parameter, 3rd rock-fabric number parameter and the superposition of the 4th rock-fabric number parameter, determine rock-fabric number parameter summation.
Rock-fabric number parameter end value determining unit 35, can be surveyed according to rock-fabric number parameter summation and deep resistivity Deep resistivity size in well curve, determines rock-fabric number parameter end value.
Specifically, as shown in figure 12, the formation fracture porosity determining unit 32, including：
Judge module 321, it can be determined that the deep resistivity and flushing under synchronization in deep resistivity log curve are powered The size of flushed zone resistivity in resistance rate log.
First computing module 322, can be when deep resistivity be more than flushed zone resistivity, according to formula：
Calculate formation fracture porosity Φfrc；Wherein RXO is flushed zone resistivity；RT is deep resistivity；Rmf filters for mud Liquid resistivity.
Second computing module 323, can be when deep resistivity be less than flushed zone resistivity, according to formula：
Calculate formation fracture porosity Φfrc。
Specifically, as shown in figure 12, the rock-fabric number parameter determination unit 33 can include：
First rock-fabric number parameter calculating module 331, can be according to natural gamma value size and formation fracture porosity Size, determine the first rock-fabric number parameter.
First rock-fabric number parameter calculating module 331, specifically can be with：
When natural gamma value is less than 0.08% more than 50API, and formation fracture porosity, according to formula：
Calculate the first rock-fabric number parameter ARFN1；Wherein, RT is deep resistivity；GR is natural gamma value.
When natural gamma value is more than 0.08% more than 50API, and formation fracture porosity, according to formula：
Calculate the first rock-fabric number parameter ARFN1。
In addition, as shown in figure 12, rock-fabric number parameter determination unit 33, in addition to：
Second rock-fabric number parameter calculating module 332, can be big according to the size and natural gamma value of interval transit time value It is small, determine the second rock-fabric number parameter.
Second rock-fabric number parameter calculating module 332, specifically can be with：
When interval transit time value is less than or equal to 50API less than or equal to 51us/f, and natural gamma value, according to formula：
Calculate the second rock-fabric number parameter ARFN2；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep electricity Resistance rate.
In addition, as shown in figure 12, rock-fabric number parameter determination unit 33, in addition to：
3rd rock-fabric number parameter calculating module 333, size, natural gamma value size that can be according to interval transit time value With deep resistivity size, the 3rd rock-fabric number parameter is determined.
3rd rock-fabric number parameter calculating module 333, specifically can be with：
It is more than 51us/f in interval transit time value, natural gamma value is less than or equal to 50API, and deep resistivity is more than 200ohmm When, according to formula：
Calculate the 3rd rock-fabric number parameter ARFN3；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep electricity Resistance rate.
In addition, as shown in figure 12, the rock-fabric number parameter determination unit 33, in addition to：
4th rock-fabric number parameter calculating module 334, can be big according to size, the natural gamma value of interval transit time value The size of small, deep resistivity and formation fracture porosity, determines the 4th rock-fabric number parameter.
4th rock-fabric number parameter calculating module 334, specifically can be with：
Interval transit time value be more than 51us/f, natural gamma value be less than or equal to 40API, deep resistivity be less than 200ohmm and When formation fracture porosity is less than 0.08%, according to formula：
Calculate the 4th rock-fabric number parameter ARFN4；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep electricity Resistance rate；GR is natural gamma value.
It is more than 51us/f in interval transit time value, natural gamma value is more than 40API and less than or equal to 50API, and deep resistivity is small When 200ohmm and formation fracture porosity are more than 0.08%, according to formula：
Calculate the 4th rock-fabric number parameter ARFN4。
, specifically can be with addition, rock-fabric number parameter superpositing unit 34：
According to formula：
ARFN=ARFN1+ARFN2+ARFN3+ARFN4
Determine rock-fabric number parameter summation ARFN.
, specifically can be with addition, rock-fabric number parameter end value determining unit 35：
When ARFN is less than or equal to the first preset value, determine that rock-fabric number parameter end value ARFN ' is equal to the ARFN；
When ARFN is more than 9000ohmm more than the first preset value, and deep resistivity, determine that rock-fabric number parameter is final Value ARFN ' is equal to the second preset value；
When ARFN is less than or equal to 9000ohmm more than the first preset value, and deep resistivity, rock-fabric number parameter is determined End value ARFN ' is equal to the first preset value；
When ARFN is less than or equal to three preset values, determine that rock-fabric number parameter end value ARFN ' is equal to the described 3rd Preset value.
What deserves to be explained is, a kind of low hole provided in an embodiment of the present invention, compact limestone texture constituent of rock log parameter The specific implementation of computing device may refer to the corresponding embodiments of the method for above-mentioned Fig. 1 and Fig. 2, here is omitted.
A kind of low hole provided in an embodiment of the present invention, the computing device of compact limestone texture constituent of rock log parameter, point Four rock-fabric number parameters have not been asked for, then final superposition determines rock-fabric number parameter summation, and then according to rock texture Deep resistivity size in number parameter summation and deep resistivity log curve, determines rock-fabric number parameter end value.The present invention It is mainly low hole, densification suitable for stratum, formation porosity is less than 10%, or with fracture development, sedimentary facies belt is mainly out Wealthy tableland or the pure limestone formation of half open-platform sedimentary facies.Rock-fabric number formula compared to prior art sets up basis and is Formation core porosity ranges are between 4%-40%, and present invention, avoiding for domestic low hole, compact limestone rock formation, rock Heart porosity ranges can not accurately determine low hole, cause between 0.5%-10% by the rock-fabric number formula of prior art The problem of structure number parameter of close limestone rock.
Apply specific embodiment in the present invention to be set forth the principle and embodiment of the present invention, above example Explanation be only intended to help to understand the method and its core concept of the present invention；Simultaneously for those of ordinary skill in the art, According to the thought of the present invention, it will change in specific embodiments and applications, in summary, in this specification Appearance should not be construed as limiting the invention.

Claims (6)

1. a kind of low hole, the computational methods of compact limestone texture constituent of rock log parameter, it is characterised in that including：
Obtain the log data of limestone rock formation to be measured；The log data includes：Deep resistivity log curve, flushing are powered Resistance rate log, mud filtrate resistivity, gamma ray curve and interval transit time curve；
According to the deep resistivity log curve and flushed zone resistivity log definitely layer fracture porosity；
Deep resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, stratum The size of interval transit time value in the size and interval transit time curve of fracture porosity, determine the first rock-fabric number parameter, Second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter；
By the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock Stone structure number parameter is superimposed, and determines rock-fabric number parameter summation；
Deep resistivity size in the rock-fabric number parameter summation and deep resistivity log curve, determines rock texture Number parameter end value；
According to the deep resistivity log curve and flushed zone resistivity log definitely layer fracture porosity, including：
Judge the flushing in the deep resistivity and flushed zone resistivity log under synchronization in deep resistivity log curve The size of strip resistance rate；
If deep resistivity is more than flushed zone resistivity, according to formula：
<mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>f</mi> <mi>r</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <mn>8.52253</mn> <mrow> <mi>R</mi> <mi>X</mi> <mi>O</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mn>8.242778</mn> <mrow> <mi>R</mi> <mi>T</mi> </mrow> </mfrac> <mo>+</mo> <mn>0.00071236</mn> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mi>R</mi> <mi>m</mi> <mi>f</mi> </mrow>
Calculate formation fracture porosity Φfrc；Wherein RXO is flushed zone resistivity；RT is deep resistivity；Rmf is F electricity Resistance rate；
If deep resistivity is less than flushed zone resistivity, according to formula：
<mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>f</mi> <mi>r</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <mn>1.99247</mn> <mrow> <mi>R</mi> <mi>T</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mn>0.992719</mn> <mrow> <mi>R</mi> <mi>X</mi> <mi>O</mi> </mrow> </mfrac> <mo>+</mo> <mn>0.00031829</mn> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mi>R</mi> <mi>m</mi> <mi>f</mi> </mrow>
Calculate formation fracture porosity Φfrc；
Deep resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, stratum The size of interval transit time value in the size and interval transit time curve of fracture porosity, determine the first rock-fabric number parameter, Second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, including：
According to the natural gamma value size and the size of formation fracture porosity, the first rock-fabric number parameter is determined；
It is described according to the natural gamma value size and the size of formation fracture porosity, determine first rock-fabric number ginseng Number, including：
If natural gamma value is more than 50API, and formation fracture porosity is less than 0.08%, according to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> </mrow> <mn>100</mn> </mfrac> <mo>+</mo> <mn>0.0225</mn> <mo>&amp;times;</mo> <mi>G</mi> <mi>R</mi> <mo>-</mo> <mn>1.125</mn> </mrow>
Calculate the first rock-fabric number parameter ARFN1；Wherein, RT is deep resistivity；GR is natural gamma value；
If natural gamma value is more than 50API, and formation fracture porosity is more than 0.08%, according to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> </mrow> <mn>100</mn> </mfrac> <mo>+</mo> <mn>0.0225</mn> <mo>&amp;times;</mo> <mi>G</mi> <mi>R</mi> <mo>-</mo> <mn>0.625</mn> </mrow>
Calculate the first rock-fabric number parameter ARFN1；
Deep resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, stratum The size of interval transit time value in the size and interval transit time curve of fracture porosity, determine the first rock-fabric number parameter, Second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, including：
According to the size and natural gamma value size of the interval transit time value, the second rock-fabric number parameter is determined；
The size and natural gamma value size according to the interval transit time value, determines the second rock-fabric number parameter, Including：
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：
<mrow> <msub> <mi>ARFN</mi> <mn>2</mn> </msub> <mo>=</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mo>(</mo> <mn>2.9</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>+</mo> <mn>3.2</mn> <mo>&amp;times;</mo> <mi>lg</mi> <mo>(</mo> <mi>R</mi> <mi>T</mi> <mo>&amp;times;</mo> <msup> <mi>e</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> <mo>-</mo> <mn>200</mn> </mrow> <mn>1800</mn> </mfrac> <mo>)</mo> </mrow> </msup> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>6.9</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>)</mo> </mrow> </mfrac> </msup> <mo>-</mo> <mn>3</mn> </mrow>
Calculate the second rock-fabric number parameter ARFN2；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep resistance Rate；
Deep resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, stratum The size of interval transit time value in the size and interval transit time curve of fracture porosity, determine the first rock-fabric number parameter, Second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, including：
According to the size of the interval transit time value, natural gamma value size and deep resistivity size, the 3rd rock knot is determined Structure number parameter；
It is described according to the size of the interval transit time value, natural gamma value size and deep resistivity size, determine the 3rd rock Stone structure number parameter, including：
If interval transit time value is more than 51us/f, natural gamma value is less than or equal to 50API, and deep resistivity is more than 200ohmm, according to Formula：
<mrow> <msub> <mi>ARFN</mi> <mn>3</mn> </msub> <mo>=</mo> <mn>9.5</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mo>(</mo> <mn>3.2</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>+</mo> <mn>2</mn> <mo>*</mo> <mi>lg</mi> <mo>(</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> </mrow> <mrow> <mi>A</mi> <mi>C</mi> </mrow> </mfrac> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>3.1</mn> <mo>+</mo> <mi>l</mi> <mi>g</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>)</mo> </mrow> </mfrac> </msup> </mrow>
Calculate the 3rd rock-fabric number parameter ARFN3；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep resistance Rate；
Deep resistivity size in deep resistivity log curve, the natural gamma value size in gamma ray curve, stratum The size of interval transit time value in the size and interval transit time curve of fracture porosity, determine the first rock-fabric number parameter, Second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock-fabric number parameter, including：
According to the size of the size of the interval transit time value, natural gamma value size, deep resistivity and formation fracture porosity, really The fixed 4th rock-fabric number parameter；
It is described according to the big of the size of the interval transit time value, natural gamma value size, deep resistivity and formation fracture porosity It is small, the 4th rock-fabric number parameter is determined, including：
If interval transit time value is more than 51us/f, natural gamma value is less than or equal to 40API, and deep resistivity is less than 200ohmm and stratum Fracture porosity is less than 0.08%, according to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>4</mn> </msub> <mo>=</mo> <mn>7.2</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mo>(</mo> <mn>3.2</mn> <mo>+</mo> <mn>1.1</mn> <mo>*</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> <mo>*</mo> <mi>R</mi> <mi>T</mi> </mrow> <mrow> <mi>G</mi> <mi>R</mi> <mo>*</mo> <mi>G</mi> <mi>R</mi> </mrow> </mfrac> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>3.1</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>G</mi> <mi>R</mi> <mo>)</mo> <mo>)</mo> </mrow> </mfrac> </msup> </mrow>
Calculate the 4th rock-fabric number parameter ARFN4；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep resistance Rate；GR is natural gamma value；
If interval transit time value is more than 51us/f, natural gamma value is more than 40API and less than or equal to 50API, and deep resistivity is less than 200ohmm and formation fracture porosity are more than 0.08%, according to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>4</mn> </msub> <mo>=</mo> <mn>7.2</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mo>(</mo> <mn>3.2</mn> <mo>+</mo> <mn>1.1</mn> <mo>*</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> <mo>*</mo> <mi>R</mi> <mi>T</mi> </mrow> <mrow> <mi>G</mi> <mi>R</mi> <mo>*</mo> <mi>G</mi> <mi>R</mi> </mrow> </mfrac> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>3.1</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>G</mi> <mi>R</mi> <mo>)</mo> <mo>)</mo> </mrow> </mfrac> </msup> <mo>+</mo> <mn>0.5</mn> </mrow>
Calculate the 4th rock-fabric number parameter ARFN4。
2. low hole according to claim 1, the computational methods of compact limestone texture constituent of rock log parameter, its feature exist In by the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock Structure number parameter is superimposed, and determines rock-fabric number parameter summation, including：
According to formula：
ARFN=ARFN1+ARFN2+ARFN3+ARFN4
Determine rock-fabric number parameter summation ARFN.
3. low hole according to claim 2, the computational methods of compact limestone texture constituent of rock log parameter, its feature exist In the deep resistivity size in the rock-fabric number parameter summation and deep resistivity log curve determines rock texture Number parameter end value, including：
If ARFN is less than or equal to the first preset value, determine that the rock-fabric number parameter end value ARFN ' is equal to the ARFN；
If ARFN is more than first preset value, and deep resistivity is more than 9000ohmm, determines the rock-fabric number parameter most Final value ARFN ' is equal to the second preset value；
If ARFN is more than first preset value, and deep resistivity is less than or equal to 9000ohmm, determines the rock-fabric number ginseng Number end value ARFN ' is equal to first preset value；
If ARFN is less than or equal to the 3rd preset value, determine that the rock-fabric number parameter end value ARFN ' is equal to the described 3rd pre- If value.
4. a kind of low hole, the computing device of compact limestone texture constituent of rock log parameter, it is characterised in that including：
Log data acquiring unit, the log data for obtaining limestone rock formation to be measured；The log data includes：Deep electricity Resistance rate log, flushed zone resistivity log, mud filtrate resistivity, gamma ray curve and interval transit time are bent Line；
Formation fracture porosity determining unit, for according to the deep resistivity log curve and flushed zone resistivity log Definitely layer fracture porosity；
Rock-fabric number parameter determination unit, for the deep resistivity size in deep resistivity log curve, natural gamma Natural gamma value size, the size of formation fracture porosity in curve and the interval transit time value in interval transit time curve it is big It is small, determine the first rock-fabric number parameter, the second rock-fabric number parameter, the 3rd rock-fabric number parameter and the 4th rock knot Structure number parameter；
Rock-fabric number parameter superpositing unit, for by the first rock-fabric number parameter, the second rock-fabric number parameter, Three rock-fabric number parameters and the superposition of the 4th rock-fabric number parameter, determine rock-fabric number parameter summation；
Rock-fabric number parameter end value determining unit, for according to the rock-fabric number parameter summation and deep resistivity log Deep resistivity size in curve, determines rock-fabric number parameter end value；
Formation fracture porosity determining unit, including：
Judge module, for judging the deep resistivity under synchronization in deep resistivity log curve and flushed zone resistivity well logging The size of flushed zone resistivity in curve；
First computing module, for when deep resistivity is more than flushed zone resistivity, according to formula：
<mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>f</mi> <mi>r</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <mn>8.52253</mn> <mrow> <mi>R</mi> <mi>X</mi> <mi>O</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mn>8.242778</mn> <mrow> <mi>R</mi> <mi>T</mi> </mrow> </mfrac> <mo>+</mo> <mn>0.00071236</mn> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mi>R</mi> <mi>m</mi> <mi>f</mi> </mrow>
Calculate formation fracture porosity Φfrc；Wherein RXO is flushed zone resistivity；RT is deep resistivity；Rmf is F electricity Resistance rate；
Second computing module, for when deep resistivity is less than flushed zone resistivity, according to formula：
<mrow> <msub> <mi>&amp;Phi;</mi> <mrow> <mi>f</mi> <mi>r</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <mn>1.99247</mn> <mrow> <mi>R</mi> <mi>T</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mn>0.992719</mn> <mrow> <mi>R</mi> <mi>X</mi> <mi>O</mi> </mrow> </mfrac> <mo>+</mo> <mn>0.00031829</mn> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mi>R</mi> <mi>m</mi> <mi>f</mi> </mrow>
Calculate formation fracture porosity Φfrc；
The rock-fabric number parameter determination unit, including：
First rock-fabric number parameter calculating module, for according to the big of the natural gamma value size and formation fracture porosity It is small, determine the first rock-fabric number parameter；
The first rock-fabric number parameter calculating module, specifically for：
When natural gamma value is less than 0.08% more than 50API, and formation fracture porosity, according to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> </mrow> <mn>100</mn> </mfrac> <mo>+</mo> <mn>0.0225</mn> <mo>&amp;times;</mo> <mi>G</mi> <mi>R</mi> <mo>-</mo> <mn>1.125</mn> </mrow>
Calculate the first rock-fabric number parameter ARFN1；Wherein, RT is deep resistivity；GR is natural gamma value；
When natural gamma value is more than 0.08% more than 50API, and formation fracture porosity, according to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> </mrow> <mn>100</mn> </mfrac> <mo>+</mo> <mn>0.0225</mn> <mo>&amp;times;</mo> <mi>G</mi> <mi>R</mi> <mo>-</mo> <mn>0.625</mn> </mrow>
Calculate the first rock-fabric number parameter ARFN1；
The rock-fabric number parameter determination unit, including：
Second rock-fabric number parameter calculating module, for the size and natural gamma value size according to the interval transit time value, Determine the second rock-fabric number parameter；
The second rock-fabric number parameter calculating module, specifically for：
When interval transit time value is less than or equal to 50API less than or equal to 51us/f, and natural gamma value, according to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>2</mn> </msub> <mo>=</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mo>(</mo> <mn>2.9</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>+</mo> <mn>3.2</mn> <mo>&amp;times;</mo> <mi>lg</mi> <mo>(</mo> <mi>R</mi> <mi>T</mi> <mo>&amp;times;</mo> <msup> <mi>e</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> <mo>-</mo> <mn>200</mn> </mrow> <mn>1800</mn> </mfrac> <mo>)</mo> </mrow> </msup> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>6.9</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>)</mo> </mrow> </mfrac> </msup> <mo>-</mo> <mn>3</mn> </mrow>
Calculate the second rock-fabric number parameter ARFN2；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep resistance Rate；
The rock-fabric number parameter determination unit, including：
3rd rock-fabric number parameter calculating module, for the size according to the interval transit time value, natural gamma value size and Deep resistivity size, determines the 3rd rock-fabric number parameter；
The 3rd rock-fabric number parameter calculating module, specifically for：
It is more than 51us/f in interval transit time value, when natural gamma value is more than 200ohmm less than or equal to 50API, and deep resistivity, root According to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>3</mn> </msub> <mo>=</mo> <mn>9.5</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mo>(</mo> <mn>3.2</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>+</mo> <mn>2</mn> <mo>*</mo> <mi>lg</mi> <mo>(</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> </mrow> <mrow> <mi>A</mi> <mi>C</mi> </mrow> </mfrac> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>3.1</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>)</mo> </mrow> </mfrac> </msup> </mrow>
Calculate the 3rd rock-fabric number parameter ARFN3；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep resistance Rate；
The rock-fabric number parameter determination unit, including：
4th rock-fabric number parameter calculating module, for the size according to the interval transit time value, natural gamma value size, depth The size of resistivity and formation fracture porosity, determines the 4th rock-fabric number parameter；
The 4th rock-fabric number parameter calculating module, specifically for：
It is more than 51us/f in interval transit time value, natural gamma value is less than or equal to 40API, and deep resistivity is less than 200ohmm and stratum When fracture porosity is less than 0.08%, according to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>4</mn> </msub> <mo>=</mo> <mn>7.2</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mo>(</mo> <mn>3.2</mn> <mo>+</mo> <mn>1.1</mn> <mo>*</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> <mo>*</mo> <mi>R</mi> <mi>T</mi> </mrow> <mrow> <mi>G</mi> <mi>R</mi> <mo>*</mo> <mi>G</mi> <mi>R</mi> </mrow> </mfrac> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>3.1</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>G</mi> <mi>R</mi> <mo>)</mo> <mo>)</mo> </mrow> </mfrac> </msup> </mrow>
Calculate the 4th rock-fabric number parameter ARFN4；Wherein, e is the nature truth of a matter；AC is interval transit time value；RT is deep resistance Rate；GR is natural gamma value；
It is more than 51us/f in interval transit time value, natural gamma value is more than 40API and less than or equal to 50API, and deep resistivity is less than When 200ohmm and formation fracture porosity are more than 0.08%, according to formula：
<mrow> <msub> <mi>ARFN</mi> <mn>4</mn> </msub> <mo>=</mo> <mn>7.2</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mo>(</mo> <mn>3.2</mn> <mo>+</mo> <mn>1.1</mn> <mo>*</mo> <mi>lg</mi> <mo>(</mo> <mi>A</mi> <mi>C</mi> <mo>)</mo> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mfrac> <mrow> <mi>R</mi> <mi>T</mi> <mo>*</mo> <mi>R</mi> <mi>T</mi> </mrow> <mrow> <mi>G</mi> <mi>R</mi> <mo>*</mo> <mi>G</mi> <mi>R</mi> </mrow> </mfrac> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>3.1</mn> <mo>+</mo> <mi>lg</mi> <mo>(</mo> <mi>G</mi> <mi>R</mi> <mo>)</mo> <mo>)</mo> </mrow> </mfrac> </msup> <mo>+</mo> <mn>0.5</mn> </mrow>
Calculate the 4th rock-fabric number parameter ARFN4。
5. low hole according to claim 4, the computing device of compact limestone texture constituent of rock log parameter, its feature exist In, the rock-fabric number parameter superpositing unit, specifically for：
According to formula：
ARFN=ARFN1+ARFN2+ARFN3+ARFN4
Determine rock-fabric number parameter summation ARFN.
6. low hole according to claim 5, the computing device of compact limestone texture constituent of rock log parameter, its feature exist In, the rock-fabric number parameter end value determining unit, specifically for：
When ARFN is less than or equal to the first preset value, determine that the rock-fabric number parameter end value ARFN ' is equal to the ARFN；
When ARFN is more than 9000ohmm more than first preset value, and deep resistivity, the rock-fabric number parameter is determined End value ARFN ' is equal to the second preset value；
When ARFN is less than or equal to 9000ohmm more than first preset value, and deep resistivity, the rock-fabric number is determined Parameter end value ARFN ' is equal to first preset value；
When ARFN is less than or equal to three preset values, determine that the rock-fabric number parameter end value ARFN ' is equal to the described 3rd Preset value.
CN201510876133.XA 2015-12-03 2015-12-03 Method and device for calculating logging parameters of low-porosity and compact limestone rock structure components Active CN105549114B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510876133.XA CN105549114B (en) 2015-12-03 2015-12-03 Method and device for calculating logging parameters of low-porosity and compact limestone rock structure components

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510876133.XA CN105549114B (en) 2015-12-03 2015-12-03 Method and device for calculating logging parameters of low-porosity and compact limestone rock structure components

Publications (2)

Publication Number Publication Date
CN105549114A CN105549114A (en) 2016-05-04
CN105549114B true CN105549114B (en) 2017-10-17

Family

ID=55828413

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510876133.XA Active CN105549114B (en) 2015-12-03 2015-12-03 Method and device for calculating logging parameters of low-porosity and compact limestone rock structure components

Country Status (1)

Country Link
CN (1) CN105549114B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109116440B (en) * 2017-06-23 2020-05-22 中国石油化工股份有限公司 Dense limestone reservoir fracture identification method
WO2024077538A1 (en) * 2022-10-13 2024-04-18 Saudi Arabian Oil Company Methods and systems for predicting lithology and formation boundary ahead of the bit

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6462549B1 (en) * 1999-04-21 2002-10-08 Schlumberger Technology Corporation Method and system for electroseismic monitoring of microseismicity
CN103968811A (en) * 2013-04-18 2014-08-06 常州华达科捷光电仪器有限公司 Adjusting mechanism and laser collimator with same
CN104865614A (en) * 2014-02-20 2015-08-26 中国石油化工股份有限公司 Complicated reservoir fluid identification method based on variable skeleton parameter

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2700666C (en) * 2007-11-27 2016-07-12 Exxonmobil Upstream Research Company Method for determining the properties of hydrocarbon reservoirs from geophysical data

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6462549B1 (en) * 1999-04-21 2002-10-08 Schlumberger Technology Corporation Method and system for electroseismic monitoring of microseismicity
CN103968811A (en) * 2013-04-18 2014-08-06 常州华达科捷光电仪器有限公司 Adjusting mechanism and laser collimator with same
CN104865614A (en) * 2014-02-20 2015-08-26 中国石油化工股份有限公司 Complicated reservoir fluid identification method based on variable skeleton parameter

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Carbonate reservoir characterization;F.Jerry Lucia 等;《Technology Today Series》;20030630;第70-72页 *

Also Published As

Publication number Publication date
CN105549114A (en) 2016-05-04

Similar Documents

Publication Publication Date Title
Grammer et al. Integration of outcrop and modern analogs in reservoir modeling: Overview with examples from the Bahamas
CN102011583B (en) Method for identifying reservoir by combining electric imaging with reef geologic model
CN109613612A (en) A kind of carbonate rock particle beach meticulous depiction and prediction technique
Worthington Net pay—what is it? What does it do? How do we quantify it? How do we use it?
CN103698811B (en) Carbonate rock structure component well logging quantitative identification method and application thereof
CN106951660A (en) Sea facies clastic rock horizontal well reservoir logging interpretation method and device
Grant et al. Outcrop analog for cyclic-shelf reservoirs, San Andres Formation of Permian Basin: Stratigraphic framework, permeability distribution, geostatistics, and fluid-flow modeling
CN106556863B (en) Porosity prediction method based on Depth Domain prestack angle gathers
Falivene et al. Synthetic seismic models from outcrop-derived reservoir-scale three-dimensional facies models: The Eocene Ainsa turbidite system (southern Pyrenees)
CN104632202B (en) Method and device for determining dry clay three-porosity logging parameter values
CN107589469B (en) Method and device for determining oil-water interface
CN109425900A (en) A kind of Seismic Reservoir Prediction method
Hart Whither seismic stratigraphy?
CN107515957A (en) Mud shale SEQUENCE STRATIGRAPHIC DIVISION method
CN112363242B (en) Reservoir fluid identification method and device based on logging fusion
CN107795320A (en) A kind of computational methods of horizontal well carbonate reservoir parameter
CN107829731A (en) A kind of volcanic rock porosity correction method of clay alteration
Askari et al. A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir
CN105549114B (en) Method and device for calculating logging parameters of low-porosity and compact limestone rock structure components
Doherty et al. Outcrop-based reservoir characterization: A composite phylloid-algal mound, western Orogrande basin (New Mexico)
CN104977612B (en) A kind of method of utilization seismic energy difference attribute forecast reservoir
Fitch Heterogeneity in the petrophysical properties of carbonate reservoirs
Rush et al. Geostatistical facies modeling trends for oolitic tidal sand shoals
CN113219531B (en) Dense sandstone gas-water distribution identification method and device
CN115857047B (en) Comprehensive prediction method for earthquake reservoir

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant