CN105628576A - Calcite dolomite porosity measuring method and apparatus - Google Patents

Abstract

The invention provides a calcite dolomite porosity measuring method and an apparatus. The calcite dolomite porosity measuring method comprises following steps: dolomite content of calcite dolomite in an area to be measured is acquired; it is determined that which target content area the acquired dolomite content belongs to; if the acquired dolomite content belongs to a first content area, the porosity of pure dolomite is taken as the porosity of calcite dolomite; if the acquired dolomite content belongs to a second content area, the porosity of pure limestone is taken as the porosity of calcite dolomite; and if the acquired dolomite content belongs to a third content area, the porosity of calcite dolomite is calculated based on the proportion of pure limestone porosity and pure dolomite porosity in the calcite dolomite porosity. According to the calcite dolomite porosity measuring method, the dolomite content of calcite dolomite is used for determining the relationship of calcite dolomite porosity with pure limestone porosity and pure dolomite porosity, the porosity of calcite dolomite is obtained based on the proportion of the pure limestone porosity and the pure dolomite porosity in the calcite dolomite porosity, and accuracy on prediction on calcite dolomite porosity is increased.

Description

Calcite dolomite porosity measurement method and device

Technical field

The present invention relates to geological exploration field, particularly relate to a kind of calcite dolomite porosity measurement method and device.

Background technology

The porosity of ground determines that the important indicator of oil and gas reserves. Thickness for Lower Paleozoic strata carbonate rock series, Cambrian system and Ordovician Dolostones reaches 1560��1692m, and Exploration Potential is huge. Wherein the lithology on Middle?Lower Ordovician-upper Cambrian series stratum is mainly limestone, dolomite, calcite dolomite, impact by the Forming Mechanism of hole and compaction, dolomite more than 4000 meters of deep layer and buried depth and limestone porosity differ greatly, its rock physics rule difference is also relatively larger, the particularly changeover portion between limestone-cloud rock, how calcite dolomite porosity being predicted is the difficult problem in dolomite exploration.

The porosity type of dolostone reservoirs mainly includes cave type, pore type and slit formation three types. For cave type dolostone reservoirs, reservoir space is large-scale corrosion cave, and Seismic reflection character shows as " beading " reflection. And pore type and slit formation dolostone reservoirs, reservoir space is mainly the little pore dimension such as intracrystalline pore, screen work hole, mould pore, crack, the combined factors such as big by buried depth, anisotropism is strong, geologic(al) factor is complicated, seismic resolution is low affects, pore type dolostone reservoirs Seismic reflection character is indefinite, and reservoir prediction difficulty is big.

The technology being currently used for predicting reservoir specifically includes that Seismic attribute analysis technology, three-dimensional coherent body technique, prestack AVO analytical technology, geophysical inversion technology and Rock physical analysis technology. But when carrying out reservoir measurement by above-mentioned technology, all deposit the problem that the porosity prediction result of calcite dolomite reservoir is poor.

Summary of the invention

The present invention provides a kind of calcite dolomite porosity measurement method and device, for solving there is, when dolomite reservoir porosity is predicted by prior art, the problem that porosity prediction result is poor.

To achieve these goals, a kind of calcite dolomite porosity measurement method provided by the invention, including:

Obtain the dolomite content of region to be measured calcite dolomite;

Determine the target content region being subordinate to of described dolomite content; Wherein concentration areas includes: the second concentration areas less than 20% of the dolomite content the first concentration areas more than 80%, dolomite content and the dolomite content the 3rd concentration areas between 20% to 80%;

If it is determined that described target content region is described first concentration areas, using the pure dolomitic porosity in the described region to be measured porosity as described calcite dolomite;

If it is determined that described target content region is described second concentration areas, using the porosity of the pure limestone in the described region to be measured porosity as described calcite dolomite;

If it is determined that described target content region is described 3rd concentration areas, porosity according to described pure limestone and described pure dolomitic porosity form the shared weight during porosity of described calcite dolomite, calculate the porosity obtaining described calcite dolomite.

To achieve these goals, a kind of calcite dolomite porosity measuring device provided by the invention, including:

Acquisition module, for obtaining the dolomite content of region to be measured calcite dolomite;

Determine module, for determining the target content region being subordinate to of described dolomite content; Wherein concentration areas includes: the second concentration areas less than 20% of the dolomite content the first concentration areas more than 80%, dolomite content and the dolomite content the 3rd concentration areas between 20% to 80%;

Computing module, if determining that module determines that described target content region is described first concentration areas for described, using the pure dolomitic porosity in the described region to be measured porosity as described calcite dolomite, if and described determine that module determines that described target content region is described second concentration areas, using the porosity of the pure limestone in the described region to be measured porosity as described calcite dolomite, if and described determine that module determines that described target content region is described 3rd concentration areas, porosity according to described pure limestone and described pure dolomitic porosity form the shared weight during porosity of described calcite dolomite, calculate the porosity obtaining described calcite dolomite.

Calcite dolomite porosity measurement method provided by the invention and device, obtain the dolomite content of region to be measured calcite dolomite, it is determined that the target content region being subordinate to of dolomite content, wherein concentration areas includes: the dolomite content the first concentration areas more than 80%, the dolomite content the second concentration areas less than 20% and the dolomite content the 3rd concentration areas between 20% to 80%, if it is determined that target content region is the first concentration areas, using the pure dolomitic porosity in the region to be measured porosity as calcite dolomite, if it is determined that target content region is the second concentration areas, using the porosity of the pure limestone in the region to be measured porosity as calcite dolomite, if it is determined that target content region is the 3rd concentration areas, porosity according to pure limestone and pure dolomitic porosity form the shared weight during porosity of calcite dolomite, calculate the porosity obtaining calcite dolomite. the present invention determines the relation between the porosity of calcite dolomite and the porosity of pure limestone and pure dolomitic porosity by the dolomite content of calcite dolomite, again based on the shared weight of the porosity of pure limestone and pure dolomitic porosity, obtain the porosity of calcite dolomite, there is the poor problem that predicts the outcome when can solve calcite dolomite porosity is predicted, improve the accuracy of prediction.

Accompanying drawing explanation

The schematic flow sheet of a kind of calcite dolomite porosity measurement method that Fig. 1 provides for the embodiment of the present invention one;

The schematic flow sheet of the another kind of calcite dolomite porosity measurement method that Fig. 2 provides for the embodiment of the present invention two;

The dolomite that Fig. 3 provides for the present embodiment two and limestone porosity and p-wave impedance cross figure;

The dolomite content that Fig. 4 provides for the present embodiment two and Poisson's ratio cross figure;

The structural representation of a kind of calcite dolomite porosity measuring device that Fig. 5 provides for the embodiment of the present invention three;

The structural representation of the another kind of calcite dolomite porosity measuring device that Fig. 6 provides for the embodiment of the present invention four.

Detailed description of the invention

For making the purpose of the embodiment of the present invention, technical scheme and advantage clearly, below in conjunction with the accompanying drawing in the embodiment of the present invention, technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is a part of embodiment of the present invention, rather than whole embodiments.

Embodiment one

The schematic flow sheet of a kind of calcite dolomite porosity measurement method that Fig. 1 provides for the embodiment of the present invention one. As it is shown in figure 1, this calcite dolomite porosity measurement method comprises the following steps:

101, the dolomite content of region to be measured calcite dolomite is obtained.

Before step 101 obtains the dolomite content of region to be measured calcite dolomite, also include:

Obtaining the related data in region to be measured, wherein related data includes geological data, log data and core data.

Wherein, geological data includes: subangle superposition of data body or pre-stack time migration data, seismic horizon. Log data includes basis well data, if coordinate, well depth, reinforcing the heart height above sea level etc., geological layering data, log and velocity of longitudinal wave, shear wave velocity, density, dolomite are containing discharge curve and well log interpretation porosity curve etc., other rock core and imaging logging data.

After getting above-mentioned related data, log data and core data are analyzed, obtain the corresponding relation between the dolomite content of calcite dolomite and Rock Elastic Parameters and Poisson's ratio. Wherein, the corresponding relation between dolomite content and the Poisson's ratio of calcite dolomite is:

$V_{\mathrm{dol}}=\left\{\begin{array}{cc}1& \mathrm{\σ}\≤0.23\\ -11.11*\mathrm{\σ}+3.56& 0.23\mathrm{p\σp}0.32\\ 0& \mathrm{\σ}\≥0.32\end{array}\right.$

Wherein, V_dolRepresenting dolomite content, �� represents Poisson's ratio.

Further, carry out prestack ripple joint inversion in length and breadth according to geological data and log data, obtain the p-wave impedance of calcite dolomite and the Poisson's ratio of calcite dolomite. After getting the Poisson's ratio of above-mentioned corresponding relation and calcite dolomite, based on this Poisson's ratio and described corresponding relation, it is possible to calculate the dolomite content of calcite dolomite.

Need explanation herein, the constant in calculation expression corresponding as 0.23p �� p0.32 in above-mentioned corresponding relation, can be able to change according to the actual geographic position at dolomite place, constant in the present embodiment is applicable to Xinjiang region, and when other areas use such as Sichuan, need geological data and the log data in other areas, the constant in above-mentioned calculation expression is adjusted.

102, the target content region being subordinate to of dolomite content is determined.

Wherein concentration areas includes: the second concentration areas less than 20% of the dolomite content the first concentration areas more than 80%, dolomite content and the dolomite content the 3rd concentration areas between 20% to 80%.

The present embodiment is provided with three concentration areas for the dolomite content of calcite dolomite in advance: the first concentration areas, the second concentration areas and the 3rd concentration areas. Wherein, the first concentration areas be dolomite content more than 80%, when dolomite content is more than 80%, the cloud rock in region to be measured is equivalent to pure dolomite; Second concentration areas be dolomite content less than 20%, when dolomite content is less than 20%, the cloud rock in region to be measured is equivalent to pure limestone; 3rd concentration areas be dolomite content between 20% to 80%, when dolomite content is between 20% to 80%, the cloud rock in region to be measured is as calcite dolomite.

After getting the dolomite content of calcite dolomite in region to be measured, it is capable of determining that, according to the dolomite content of this calcite dolomite, the target content region that this dolomite content is subordinate to. If target content region is the first concentration areas, perform step 103; If target content region is the second concentration areas, perform step 104; If target content region is the 3rd concentration areas, perform step 105.

103 if it is determined that target content region is the first concentration areas, using the pure dolomitic porosity in the region to be measured porosity as calcite dolomite.

When target content region is the first concentration areas, illustrate that the cloud rock in region to be measured is equivalent to pure dolomite, in such cases, the calcite dolomite in region to be measured processes as pure dolomite, using the pure dolomitic porosity in the region to be measured porosity as calcite dolomite.

104 if it is determined that target content region is the second concentration areas, using the porosity of the pure limestone in the region to be measured porosity as calcite dolomite.

When target content region is the second concentration areas, illustrate that the cloud rock in region to be measured is equivalent to pure limestone, in such cases, the calcite dolomite in region to be measured processes as pure limestone, using the porosity of the pure limestone in the region to be measured porosity as calcite dolomite.

105 if it is determined that target content region is the 3rd concentration areas, forms the shared weight during porosity of calcite dolomite according to the porosity of pure limestone and pure dolomitic porosity, calculates the porosity obtaining calcite dolomite.

When target content region is three concentration areas, the cloud rock that region to be measured is described is calcite dolomite, in such cases, pure limestone and pure dolomitic shared weight during this calcite dolomite of composition can be got, namely form the shared weight during porosity of calcite dolomite according to the porosity of pure limestone and pure dolomitic porosity, calculate the porosity obtaining calcite dolomite.

The calcite dolomite porosity measurement method that the present embodiment provides, by obtaining the dolomite content of region to be measured calcite dolomite, determine the target content region being subordinate to of dolomite content, if it is determined that target content region is the first concentration areas, using the pure dolomitic porosity in the region to be measured porosity as calcite dolomite, if it is determined that target content region is the second concentration areas, using the porosity of the pure limestone in the region to be measured porosity as calcite dolomite, if it is determined that target content region is the 3rd concentration areas, porosity according to pure limestone and pure dolomitic porosity form the shared weight during porosity of calcite dolomite, calculate the porosity obtaining calcite dolomite. the present embodiment determines the relation between the porosity of calcite dolomite and the porosity of pure limestone and pure dolomitic porosity by the dolomite content of calcite dolomite, shared weight based on the porosity of pure limestone and pure dolomitic porosity, obtain the porosity of calcite dolomite, there is the poor problem that predicts the outcome when can solve calcite dolomite porosity is predicted, improve the accuracy to calcite dolomite porosity prediction.

Embodiment two

The schematic flow sheet of the another kind of calcite dolomite porosity measurement method that Fig. 2 provides for the embodiment of the present invention two. As in figure 2 it is shown, this calcite dolomite porosity measurement method comprises the following steps:

201, the related data in described region to be measured is obtained.

Wherein, described related data includes geological data, log data and core data.

Wherein, geological data includes: subangle superposition of data body or pre-stack time migration data, seismic horizon. Log data includes basis well data, if coordinate, well depth, reinforcing the heart height above sea level etc., geological layering data, log and velocity of longitudinal wave, shear wave velocity, density, dolomite are containing discharge curve and well log interpretation porosity curve etc., other rock core and imaging logging data.

202, carry out Rock physical analysis according to log data and core data, obtain the first mapping relations between pure limestone porosity and pure limestone p-wave impedance and the second mapping relations between pure white cloud petrosal foramen porosity and pure dolomite p-wave impedance.

Specifically, the Rock physical analysis method log data to obtaining and core data is adopted to be analyzed, obtain the first mapping relations between pure limestone porosity and pure limestone p-wave impedance and the second mapping relations between pure white cloud petrosal foramen porosity and pure dolomite p-wave impedance.

Wherein, the first mapping relations are: ��₁=f (PI)=-2.12e^-6*PI₁+ 37;

Second mapping relations are: ��₂=g (PI)=-1.85e^-6*PI₂+ 36;

��₁For the porosity of pure limestone, PI₁P-wave impedance for pure limestone.

��₂For pure dolomitic porosity, PI₂For pure dolomitic p-wave impedance.

Herein it is to be appreciated that coefficient-2.12e in above-mentioned first mapping relations^-6With coefficient-1.85e in constant 37 and the second mapping relations^-6With constant 36, changing according to the actual geographic position at dolomite place, the coefficient in the present embodiment is applicable to Xinjiang region, and when other areas use such as Sichuan, need log data and the core data in other areas, the coefficient in above-mentioned mapping relations and constant are adjusted.

203, log data and core data are analyzed, obtain the corresponding relation between dolomite content and the Poisson's ratio of calcite dolomite.

After getting above-mentioned related data, log data and core data are analyzed, obtain the corresponding relation between the dolomite content of calcite dolomite and Rock Elastic Parameters and Poisson's ratio. Wherein, the corresponding relation between dolomite content and the Poisson's ratio of calcite dolomite is:

$V_{\mathrm{dol}}=\left\{\begin{array}{cc}1& \mathrm{\σ}\≤0.23\\ -11.11*\mathrm{\σ}+3.56& 0.23\mathrm{p\σp}0.32\\ 0& \mathrm{\σ}\≥0.32\end{array}\right.$

Wherein, V_dolRepresenting dolomite content, �� represents Poisson's ratio.

Need explanation herein, the coefficient-11.11 in calculation expression corresponding as 0.23p �� p0.32 in above-mentioned corresponding relation and constant 3.56, can be able to change according to the actual geographic position at dolomite place, constant in the present embodiment is applicable to Xinjiang region, and when other areas use such as Sichuan, need geological data and the log data in other areas, the coefficient in above-mentioned calculation expression and constant are adjusted.

204, geological data and log data are carried out prestack ripple joint inversion in length and breadth, obtain the p-wave impedance of calcite dolomite and the Poisson's ratio of calcite dolomite.

205, based on the Poisson's ratio of calcite dolomite and described corresponding relation, the dolomite content of calcite dolomite is got.

Further, carry out prestack ripple joint inversion in length and breadth according to geological data and log data, obtain the p-wave impedance of calcite dolomite and the Poisson's ratio of calcite dolomite. After getting the Poisson's ratio of above-mentioned corresponding relation and calcite dolomite, based on this Poisson's ratio and described corresponding relation, it is possible to calculate the dolomite content of calcite dolomite.

206, the target content region being subordinate to of dolomite content is determined.

Wherein concentration areas includes: the second concentration areas less than 20% of the dolomite content the first concentration areas more than 80%, dolomite content and the dolomite content the 3rd concentration areas between 20% to 80%.

Step 206 referring to the record of related content in step 102 in above-described embodiment one, can repeat no more herein.

If target content region is the first concentration areas, perform step 207; If target content region is the second concentration areas, perform step 208; If target content region is the 3rd concentration areas, perform step 209.

207 if it is determined that target content region is the first concentration areas, using the pure dolomitic porosity in the region to be measured porosity as calcite dolomite.

Specifically, using the first mapping relations of obtaining as the relation between the porosity and calcite dolomite p-wave impedance of calcite dolomite, based on these first mapping relations and calcite dolomite p-wave impedance, calculate the porosity obtaining calcite dolomite.

208 if it is determined that target content region is the second concentration areas, using the porosity of the pure limestone in the region to be measured porosity as calcite dolomite.

Specifically, using the second mapping relations of obtaining as the relation between the porosity and calcite dolomite p-wave impedance of calcite dolomite, based on these second mapping relations and calcite dolomite p-wave impedance, calculate the porosity obtaining calcite dolomite.

209 if it is determined that target content region is the 3rd concentration areas, forms the shared weight during porosity of calcite dolomite according to the porosity of pure limestone and pure dolomitic porosity, calculates the porosity obtaining calcite dolomite.

Specifically, after getting the p-wave impedance of described calcite dolomite, the p-wave impedance according to the p-wave impedance of pure limestone and calcite dolomite, calculate and obtain the p-wave impedance of pure limestone and the shared weight of pure dolomitic p-wave impedance. Wherein, this shared weight is the shared weight when porosity of pure limestone forms the porosity of calcite dolomite with pure dolomitic porosity. After getting described shared weight, described shared weight is multiplied with the second mapping relations with the first corresponding mapping relations respectively, obtain relation corresponding between the porosity of calcite dolomite and the p-wave impedance of described calcite dolomite, based on relation corresponding between porosity and the dolomitic p-wave impedance of matter of the calcite dolomite obtained and calcite dolomite p-wave impedance, calculate the porosity obtaining calcite dolomite.

Based on the record of above-mentioned steps 207��step 209, the expression formula of the porosity that can get calcite dolomite in the present embodiment is as follows:

$\mathrm{\φ}=\left\{\begin{array}{cc}f\left(\mathrm{PI}\right)={-2.12e}^{-6}*\mathrm{PI}+37& V_{\mathrm{dol}}p0.2\\ g\left(\mathrm{PI}\right)={-1.85e}^{-6}*\mathrm{PI}+36& V_{\mathrm{dol}}f0.8\\ q*g\left(\mathrm{PI}\right)+(1-q)*f\left(\mathrm{PI}\right)& 0.2\≤V_{\mathrm{dol}}\≤0.8\end{array}\right.$

Wherein, �� represents the porosity of calcite dolomite, and PI represents the p-wave impedance of calcite dolomite, and q represents weight shared by pure dolomitic porosity, and (1-q) represents weight shared by pure dolomitic porosity, V_dolRepresent dolomite content.

Within 2012, certain field drilling explores Lower Ordovician Series-upper Cambrian series dolomite oil-gas possibility risk well flatly, and according to poststack Seismic Reservoir Prediction result, this well dolomite High-quality Reservoir is grown. Actual results of drilling shows, owing to deep layer dolomite grey matter content is big, causes this wellhole gap aplasia, causes drilling well to defeat. Use up-to-date related data subsequently, adopt the method that the present embodiment provides that this well is re-started reservoir prediction, prediction porosity and well log interpretation porosity are in the identical rate more than 95% of purpose, and its porosity error is less than 1%, thus proving the method effectiveness in actual applications. The dolomite that Fig. 3 provides for the present embodiment two and limestone porosity and p-wave impedance cross figure; The dolomite content that Fig. 4 provides for the present embodiment two and Poisson's ratio cross figure.

In the present embodiment, the relation between the porosity of calcite dolomite and the porosity of pure limestone and pure dolomitic porosity is determined by the dolomite content of calcite dolomite, again based on the porosity of pure limestone and pure dolomitic porosity and shared weight, obtain the porosity of calcite dolomite, there is the poor problem that predicts the outcome when can solve calcite dolomite porosity is predicted, improve the accuracy to calcite dolomite porosity prediction.

Embodiment three

The structural representation of a kind of calcite dolomite porosity measuring device that Fig. 5 provides for the embodiment of the present invention three. As it is shown in figure 5, this calcite dolomite porosity measuring device includes: acquisition module 21, determine module 22 and computing module 23.

Wherein, acquisition module 21, for obtaining the dolomite content of region to be measured calcite dolomite.

Specifically, obtain the dolomite content of region to be measured calcite dolomite at acquisition module 21 before, it is necessary to obtaining the related data in region to be measured, wherein related data includes geological data, log data and core data. Wherein, geological data includes: subangle superposition of data body or pre-stack time migration data, seismic horizon. Log data includes basis well data, if coordinate, well depth, reinforcing the heart height above sea level etc., geological layering data, log and velocity of longitudinal wave, shear wave velocity, density, dolomite are containing discharge curve and well log interpretation porosity curve etc., other rock core and imaging logging data.

After getting above-mentioned related data, log data and core data are analyzed, obtain the corresponding relation between the dolomite content of calcite dolomite and Rock Elastic Parameters and Poisson's ratio. Wherein, the corresponding relation between dolomite content and the Poisson's ratio of calcite dolomite is:

$V_{\mathrm{dol}}=\left\{\begin{array}{cc}1& \mathrm{\σ}\≤0.23\\ -11.11*\mathrm{\σ}+3.56& 0.23\mathrm{p\σp}0.32\\ 0& \mathrm{\σ}\≥0.32\end{array}\right.$

Wherein, V_dolRepresenting dolomite content, �� represents Poisson's ratio.

Further, carry out prestack ripple joint inversion in length and breadth according to geological data and log data, obtain the p-wave impedance of calcite dolomite and the Poisson's ratio of calcite dolomite. After getting the Poisson's ratio of above-mentioned corresponding relation and calcite dolomite, based on this Poisson's ratio and described corresponding relation, it is possible to calculate the dolomite content of calcite dolomite.

Need explanation herein, the constant in calculation expression corresponding as 0.23p �� p0.32 in above-mentioned corresponding relation, can be able to change according to the actual geographic position at dolomite place, constant in the present embodiment is applicable to Xinjiang region, and when other areas use such as Sichuan, need geological data and the log data in other areas, the constant in above-mentioned calculation expression is adjusted.

Cover half block 22 really it is connected, for determining the target content region being subordinate to of described dolomite content with acquisition module 21.

Wherein concentration areas includes: the second concentration areas less than 20% of the dolomite content the first concentration areas more than 80%, dolomite content and the dolomite content the 3rd concentration areas between 20% to 80%.

The present embodiment is provided with three concentration areas for the dolomite content of calcite dolomite in advance: the first concentration areas, the second concentration areas and the 3rd concentration areas. Wherein, the first concentration areas be dolomite content more than 80%, when dolomite content is more than 80%, the cloud rock in region to be measured is equivalent to pure dolomite; Second concentration areas be dolomite content less than 20%, when dolomite content is less than 20%, the cloud rock in region to be measured is equivalent to pure limestone; 3rd concentration areas be dolomite content between 20% to 80%, when dolomite content is between 20% to 80%, the cloud rock in region to be measured is as calcite dolomite.

Get the dolomite content of calcite dolomite in region to be measured at acquisition module 21 after, it is determined that module 22 is capable of determining that, according to the dolomite content of this calcite dolomite, the target content region that this dolomite content is subordinate to.

With the computing module 23 determining that module 22 is connected, for if it is determined that module 22 determines that target content region is the first concentration areas, using the pure dolomitic porosity in the region to be measured porosity as calcite dolomite, and if it is determined that module 22 determines that target content region is the second concentration areas, using the porosity of the pure limestone in the region to be measured porosity as calcite dolomite, and if it is determined that module 22 determines that target content region is the 3rd concentration areas, porosity according to pure limestone and pure dolomitic porosity form the shared weight during porosity of calcite dolomite, calculate the porosity obtaining calcite dolomite.

The calcite dolomite porosity measuring device that the present embodiment provides, by obtaining the dolomite content of region to be measured calcite dolomite, determine the target content region being subordinate to of dolomite content, if it is determined that target content region is the first concentration areas, using the pure dolomitic porosity in the region to be measured porosity as calcite dolomite, if it is determined that target content region is the second concentration areas, using the porosity of the pure limestone in the region to be measured porosity as calcite dolomite, if it is determined that target content region is the 3rd concentration areas, porosity according to pure limestone and pure dolomitic porosity form the shared weight during porosity of calcite dolomite, calculate the porosity obtaining calcite dolomite. the present embodiment determines the relation between the porosity of calcite dolomite and the porosity of pure limestone and pure dolomitic porosity by the dolomite content of calcite dolomite, based on weight shared by the porosity of pure limestone and pure dolomitic porosity, obtain the porosity of calcite dolomite, there is the poor problem that predicts the outcome when can solve calcite dolomite porosity is predicted, improve the accuracy to calcite dolomite porosity prediction.

Embodiment four

Fig. 6 is the structural representation of the another kind of calcite dolomite porosity measuring device of example four offer for the present invention. As shown in Figure 6, this device is except including the acquisition module 21 in above-described embodiment three, determining module 22 and except computing module 23, also include: data acquisition module 24 and analysis module 25.

Wherein, data acquisition module 24, for obtaining the related data in described region to be measured, described related data includes geological data, log data and core data. Introduction about related data referring to the record of related content in above-described embodiment, can repeat no more herein.

The analysis module 25 being connected with data acquisition module 24, for carrying out Rock physical analysis according to log data and core data, obtain the first mapping relations between pure limestone porosity and described pure limestone p-wave impedance and the second mapping relations between pure white cloud petrosal foramen porosity and pure dolomite p-wave impedance.

Wherein, the first mapping relations are: ��₁=f (PI)=-2.12e^-6*PI₁+ 37;

Second mapping relations are: ��₂=g (PI)=-1.85e^-6*PI₂+ 36;

��₁For the porosity of pure limestone, PI₁P-wave impedance for pure limestone;

��₂For pure dolomitic porosity, PI₂For pure dolomitic p-wave impedance.

In the present embodiment, a kind of optional frame mode of acquisition module 21 includes: analytic unit 211, inverting unit 212 and acquiring unit 213.

Wherein, analytic unit 211, for log data and core data are analyzed, obtain the corresponding relation between dolomite content and the Poisson's ratio of calcite dolomite.

Wherein, the corresponding relation between dolomite content and the Poisson's ratio of calcite dolomite is:

$V_{\mathrm{dol}}=\left\{\begin{array}{cc}1& \mathrm{\σ}\≤0.23\\ -11.11*\mathrm{\σ}+3.56& 0.23\mathrm{p\σp}0.32\\ 0& \mathrm{\σ}\≥0.32\end{array}\right.$

Wherein, V_dolRepresenting dolomite content, �� represents Poisson's ratio.

Inverting unit 212, for geological data and log data are carried out prestack ripple joint inversion in length and breadth, gets the p-wave impedance of calcite dolomite and the Poisson's ratio of calcite dolomite.

The acquiring unit 213 being connected with analytic unit 211 and inverting unit 212 respectively, is used for the Poisson's ratio based on calcite dolomite and described corresponding relation, gets the dolomite content of calcite dolomite.

Further, carry out prestack ripple joint inversion in length and breadth according to geological data and log data, obtain the p-wave impedance of calcite dolomite and the Poisson's ratio of calcite dolomite. After getting the Poisson's ratio of above-mentioned corresponding relation and calcite dolomite, based on this Poisson's ratio and described corresponding relation, it is possible to calculate the dolomite content of calcite dolomite.

Further, described acquisition module 21 also includes: weight calculation unit 214.

Described weight calculation unit 214, for after inverting unit 213 gets the p-wave impedance of calcite dolomite, the p-wave impedance according to the p-wave impedance of pure limestone, pure dolomitic p-wave impedance and calcite dolomite, calculating obtains described shared weight.

In the present embodiment, a kind of optional frame mode of computing module 23 includes: map unit 231 and computing unit 232.

Wherein, map unit 231, for when determining that target content region is the first concentration areas, using the first mapping relations of obtaining as the relation between the porosity and described calcite dolomite p-wave impedance of calcite dolomite.

The computing unit 232 being connected with map unit 231, for based on the first mapping relations and calcite dolomite p-wave impedance, calculating the porosity obtaining calcite dolomite.

Further, map unit 231, it is additionally operable to when determining that target content region is the second concentration areas, using the second mapping relations of obtaining as the relation between the porosity and calcite dolomite p-wave impedance of calcite dolomite.

Then computing unit 232, are additionally operable to, based on the second mapping relations and calcite dolomite p-wave impedance, calculate the porosity obtaining calcite dolomite.

Further, map unit 231, it is additionally operable to, when determining that target content region is three concentration areas, shared weight be multiplied with the second mapping relations with the first corresponding mapping relations respectively, obtain relation corresponding between the porosity of calcite dolomite and the p-wave impedance of calcite dolomite.

Then computing unit 232, are additionally operable to relation corresponding between the porosity based on the calcite dolomite obtained and the p-wave impedance of calcite dolomite and calcite dolomite p-wave impedance, calculate the porosity obtaining calcite dolomite.

The expression formula of the porosity that can get calcite dolomite in the present embodiment is as follows:

$\mathrm{\φ}=\left\{\begin{array}{cc}f\left(\mathrm{PI}\right)={-2.12e}^{-6}*\mathrm{PI}+37& V_{\mathrm{dol}}p0.2\\ g\left(\mathrm{PI}\right)={-1.85e}^{-6}*\mathrm{PI}+36& V_{\mathrm{dol}}f0.8\\ q*g\left(\mathrm{PI}\right)+(1-q)*f\left(\mathrm{PI}\right)& 0.2\≤V_{\mathrm{dol}}\≤0.8\end{array}\right.$

Wherein, �� represents the porosity of calcite dolomite, and PI represents the p-wave impedance of calcite dolomite, and q represents weight shared by pure dolomitic porosity, and (1-q) represents weight shared by pure dolomitic porosity, V_dolRepresent dolomite content.

In the present embodiment, the relation between the porosity of calcite dolomite and the porosity of pure limestone and pure dolomitic porosity is determined by the dolomite content of calcite dolomite, shared weight based on the porosity of pure limestone and pure dolomitic porosity, obtain the porosity of calcite dolomite, solve and when calcite dolomite porosity is predicted, there is the poor problem that predicts the outcome, improve the accuracy to calcite dolomite porosity prediction.

Last it is noted that various embodiments above is only in order to illustrate technical scheme, it is not intended to limit; Although the present invention being described in detail with reference to foregoing embodiments, it will be understood by those within the art that: the technical scheme described in foregoing embodiments still can be modified by it, or wherein some or all of technical characteristic is carried out equivalent replacement; And these amendments or replacement, do not make the essence of appropriate technical solution depart from the scope of various embodiments of the present invention technical scheme.

Claims (10)

1. a calcite dolomite porosity measurement method, it is characterised in that including:

Obtain the dolomite content of region to be measured calcite dolomite;

Determine the target content region being subordinate to of described dolomite content; Wherein concentration areas includes: the second concentration areas less than 20% of the dolomite content the first concentration areas more than 80%, dolomite content and the dolomite content the 3rd concentration areas between 20% to 80%;

If it is determined that described target content region is described first concentration areas, using the pure dolomitic porosity in the described region to be measured porosity as described calcite dolomite;

If it is determined that described target content region is described second concentration areas, using the porosity of the pure limestone in the described region to be measured porosity as described calcite dolomite;

If it is determined that described target content region is described 3rd concentration areas, porosity according to described pure limestone and described pure dolomitic porosity form the shared weight during porosity of described calcite dolomite, calculate the porosity obtaining described calcite dolomite.

2. calcite dolomite porosity measurement method according to claim 1, it is characterised in that before the dolomite content in described acquisition region to be measured, including:

Obtaining the related data in described region to be measured, described related data includes geological data, log data and core data;

Carry out Rock physical analysis according to described log data and described core data, obtain the first mapping relations between described pure limestone porosity and described pure limestone p-wave impedance and the second mapping relations between described pure white cloud petrosal foramen porosity and described pure dolomite p-wave impedance;

Wherein, the first mapping relations are: ��₁=f (PI)=-2.12e^-6*PI₁+ 37;

Second mapping relations are: ��₂=g (PI)=-1.85e^-6*PI₂+ 36;

��₁For the porosity of described pure limestone, PI₁P-wave impedance for described pure limestone;

��₂For described pure dolomitic porosity, PI₂For described pure dolomitic p-wave impedance.

3. calcite dolomite porosity measurement method according to claim 2, it is characterised in that the dolomite content of described acquisition region to be measured calcite dolomite, including:

Described log data and described core data are analyzed, obtain the corresponding relation between described dolomite content and described Poisson's ratio;

Described geological data and described log data are carried out prestack ripple joint inversion in length and breadth, obtains the p-wave impedance of described calcite dolomite and described Poisson's ratio;

Based on described Poisson's ratio and described corresponding relation, get described dolomite content.

4. calcite dolomite porosity measurement method according to claim 3, it is characterized in that, also include: after getting the p-wave impedance of described calcite dolomite, the p-wave impedance of the p-wave impedance according to described pure limestone, described pure dolomitic p-wave impedance and described calcite dolomite, calculates and obtains described shared weight.

5. calcite dolomite porosity measurement method according to claim 4, it is characterized in that, described if it is determined that described target content region is described first concentration areas, using the pure dolomitic porosity in the described region to be measured dolomitic porosity as described region to be measured, including:

Using described first mapping relations that obtain as the relation between the porosity and described calcite dolomite p-wave impedance of described calcite dolomite;

Based on described first mapping relations and described calcite dolomite p-wave impedance, calculate the porosity obtaining described calcite dolomite;

If it is determined that described target content region is described second concentration areas, using the porosity of the pure limestone in the described region to be measured porosity as described calcite dolomite, including:

Using described second mapping relations that obtain as the relation between the porosity and described calcite dolomite p-wave impedance of described calcite dolomite;

Based on described second mapping relations and described calcite dolomite p-wave impedance, calculate the porosity obtaining described calcite dolomite;

If it is determined that described target content region is described 3rd concentration areas, porosity according to described pure limestone and described pure dolomitic porosity form the shared weight during porosity of described calcite dolomite, calculate the porosity obtaining described calcite dolomite, including:

Described shared weight is multiplied with described second mapping relations with corresponding described first mapping relations respectively, obtains relation corresponding between the porosity of described calcite dolomite and the p-wave impedance of described calcite dolomite;

Based on relation corresponding between porosity and the p-wave impedance of described calcite dolomite of the described calcite dolomite obtained and described calcite dolomite p-wave impedance, calculate the porosity obtaining described calcite dolomite.

6. a calcite dolomite porosity measuring device, it is characterised in that including:

Acquisition module, for obtaining the dolomite content of region to be measured calcite dolomite;

Determine module, for determining the target content region being subordinate to of described dolomite content; Wherein concentration areas includes: the second concentration areas less than 20% of the dolomite content the first concentration areas more than 80%, dolomite content and the dolomite content the 3rd concentration areas between 20% to 80%;

Computing module, if determining that module determines that described target content region is described first concentration areas for described, using the pure dolomitic porosity in the described region to be measured porosity as described calcite dolomite, if and described determine that module determines that described target content region is described second concentration areas, using the porosity of the pure limestone in the described region to be measured porosity as described calcite dolomite, if and described determine that module determines that described target content region is described 3rd concentration areas, porosity according to described pure limestone and described pure dolomitic porosity form the shared weight during porosity of described calcite dolomite, calculate the porosity obtaining described calcite dolomite.

7. calcite dolomite porosity measuring device according to claim 6, it is characterised in that also bag:

Data acquisition module, for obtaining the related data in described region to be measured, described related data includes geological data, log data and core data;

Analysis module, for carrying out Rock physical analysis according to described log data and described core data, obtain the first mapping relations between described pure limestone porosity and described pure limestone p-wave impedance and the second mapping relations between described pure white cloud petrosal foramen porosity and described pure dolomite p-wave impedance;

Wherein, the first mapping relations are: ��₁=f (PI)=-2.12e^-6*PI₁+ 37;

Second mapping relations are: ��₂=g (PI)=-1.85e^-6*PI₂+ 36;

��₁For the porosity of described pure limestone, PI₁P-wave impedance for described pure limestone;

��₂For described pure dolomitic porosity, PI₂For described pure dolomitic p-wave impedance.

8. calcite dolomite porosity measuring device according to claim 7, it is characterised in that described acquisition module, specifically includes:

Analytic unit, for described log data and described core data are analyzed, obtains the corresponding relation between described dolomite content and the Poisson's ratio of described calcite dolomite;

Inverting unit, for described geological data and described log data are carried out prestack ripple joint inversion in length and breadth, gets the p-wave impedance of described calcite dolomite and described Poisson's ratio;

Acquiring unit, for based on described Poisson's ratio and described corresponding relation, getting described dolomite content.

9. calcite dolomite porosity measuring device according to claim 8, it is characterised in that described acquisition module, also includes: weight calculation unit;

Described weight calculation unit, for get described calcite dolomite at described inverting unit p-wave impedance after, the p-wave impedance of the p-wave impedance according to described pure limestone, described pure dolomitic p-wave impedance and described calcite dolomite, calculates and obtains described shared weight.

10. calcite dolomite porosity measuring device according to claim 9, it is characterised in that described computing module, including: map unit and computing unit;

Described map unit, for when determining that described target content region is described first concentration areas, using described first mapping relations that obtain as the relation between the porosity and described calcite dolomite p-wave impedance of described calcite dolomite;

Described computing unit, for based on described first mapping relations and described calcite dolomite p-wave impedance, calculating the porosity obtaining described calcite dolomite;

Described map unit, is additionally operable to when determining that described target content region is described second concentration areas, using described second mapping relations that obtain as the relation between the porosity and described calcite dolomite p-wave impedance of described calcite dolomite;

Described computing unit, is additionally operable to, based on described second mapping relations and described calcite dolomite p-wave impedance, calculate the porosity obtaining described calcite dolomite;

Described map unit, it is additionally operable to when determining that described target content region is described three concentration areas, described shared weight is multiplied with described second mapping relations with corresponding described first mapping relations respectively, obtains relation corresponding between the porosity of described calcite dolomite and the p-wave impedance of described calcite dolomite;

Described computing unit, is additionally operable to relation corresponding between the porosity based on the described calcite dolomite obtained with the p-wave impedance of described calcite dolomite and described calcite dolomite p-wave impedance, calculates the porosity obtaining described calcite dolomite.