CN111561312A - Saturation calculation method and device based on conventional logging data - Google Patents

Saturation calculation method and device based on conventional logging data Download PDF

Info

Publication number
CN111561312A
CN111561312A CN202010315304.2A CN202010315304A CN111561312A CN 111561312 A CN111561312 A CN 111561312A CN 202010315304 A CN202010315304 A CN 202010315304A CN 111561312 A CN111561312 A CN 111561312A
Authority
CN
China
Prior art keywords
saturation
height
function
relational expression
logging
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.)
Granted
Application number
CN202010315304.2A
Other languages
Chinese (zh)
Other versions
CN111561312B (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 CN202010315304.2A priority Critical patent/CN111561312B/en
Publication of CN111561312A publication Critical patent/CN111561312A/en
Application granted granted Critical
Publication of CN111561312B publication Critical patent/CN111561312B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Environmental & Geological Engineering (AREA)
  • Computational Mathematics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a saturation calculation method and a saturation calculation device based on conventional logging data, wherein the saturation calculation method based on the conventional logging data comprises the following steps: determining a relational expression between the saturation and the height of the oil column according to the logging data of the target block; generating a saturation height function of the target block according to a relational expression between the saturation and the height of the oil column; and calculating the saturation of the target block according to the saturation height function. The method can quickly and effectively determine the saturation height function based on the conventional logging data so as to accurately determine the saturation, and lays a solid foundation for fine evaluation of oil and gas field exploration and development reservoirs, reserve calculation and establishment of a saturation field attribute model in a three-dimensional geological model.

Description

Saturation calculation method and device based on conventional logging data
Technical Field
The invention relates to the technical field of petroleum and natural gas development, in particular to saturation calculation based on conventional logging data.
Background
In geological modeling, the establishment of a saturation field is an important work in attribute modeling, because the well logging interpretation saturation changes greatly due to the reservoir per se in the longitudinal direction, an ideal result cannot be obtained by utilizing the well point well logging interpretation saturation to perform transverse interpolation, and the establishment of the saturation field through a saturation height function is often the most effective means; reserves calculated using the saturation field established in the geological model are closer to the actual condition of the reservoir. The saturation height function, namely Sw ═ f (H, POR, K) (Sw is the water saturation, H is the height of the oil column above the free water interface, POR is the porosity, K is the permeability), is a very important function for reservoir evaluation, reserve accurate calculation and geological modeling in oil and gas exploration and development. In reservoir evaluation, the water saturation of the logging interpretation can be checked whether the function is reasonable or not so as to complete the function in time;
at present, the saturation height function is generally obtained by a capillary pressure curve obtained by measuring a rock core in a laboratory by using a high-pressure mercury intrusion method, a diaphragm method or a centrifugal method. The capillary pressure curves obtained in the laboratories can obtain the saturation height function through various conversions and multiple regression, but some uncertain factors exist in the process of determining the saturation height function under the oil reservoir condition by using laboratory data, which are mainly shown in the following steps:
(1) the interfacial tension and contact angle used in the conversion of indoor capillary pressure to capillary pressure at reservoir conditions are fixed empirical constants, and these parameters are varied by the reservoir fluid properties.
(2) The difference between room conditions and reservoir conditions (mainly reflected in temperature and pressure) results in some differences in porosity and permeability, which are not taken into account in the capillary pressure transition process.
(3) Permeability is used in the saturation height function and is difficult to determine accurately, especially in carbonate reservoirs (e.g., the middle east).
Therefore, the saturation height function established by the indoor rock core experimental data is limited in the practical application process due to the factors. On the other hand, the capillary pressure measurement cost in the laboratory is high, the experiment period is long, and particularly, the capillary pressure experiment by a partition plate method causes that the number of cores for capillary pressure measurement is small; furthermore, some old oil field exploration and development do not have the condition to carry out capillary pressure rock core experimental data and the like in the initial stage. Therefore, the saturation height function established based on the core data cannot be realized at all.
Disclosure of Invention
Aiming at the problems in the prior art, the method can quickly and effectively determine the height function of the saturation degree based on the conventional logging data, so as to accurately determine the saturation degree, and lay a solid foundation for the fine evaluation of the exploration and development reservoir of the oil and gas field, the calculation of the reserve volume and the establishment of a saturation field attribute model in a three-dimensional geological model.
In order to solve the technical problems, the invention provides the following technical scheme:
in a first aspect, the present invention provides a saturation calculation method based on conventional well logging data, including:
determining a relational expression between the saturation and the height of the oil column according to the logging data of the target block;
generating a saturation height function of the target block according to a relational expression between the saturation and the height of the oil column;
and calculating the saturation of the target block according to the saturation height function.
In one embodiment, the determining the relation between the saturation and the height of the oil column according to the logging data of the target block includes:
determining a relational expression between the saturation and the logging depth according to the logging data;
converting the log depth to the oil column height to generate a relationship between the saturation and oil column height.
In an embodiment, the generating a saturation height function of the target block according to the relation between the saturation and the oil column height includes:
converting the relation between the saturation and the height of the oil column into a relation between the saturation and capillary pressure;
generating a relation between a J function and the saturation according to the relation between the saturation and capillary pressure;
and generating a saturation height function of the target block according to a relation between the J function and the saturation.
In an embodiment, the generating a saturation height function of the target block according to the relation between the saturation and the oil column height includes:
and generating the saturation height function according to the height of the oil column, the water saturation and the porosity by utilizing a multiple regression method.
In one embodiment, the well data in the transition zone of the target block with the altitude within a preset range is selected as the well logging data.
In a second aspect, the present invention provides a saturation calculation apparatus based on conventional well log data, the apparatus comprising:
the first relational expression determining unit is used for determining a relational expression between the saturation and the height of the oil column according to the logging data of the target block;
the saturation height function generating unit is used for generating a saturation height function of the target block according to a relational expression between the saturation and the height of the oil column;
and the saturation calculation unit is used for calculating the saturation of the target block according to the saturation height function.
In one embodiment, the first relational expression determining unit includes:
the second relational expression determining module is used for determining a relational expression between the saturation and the logging depth according to the logging data;
and the logging depth conversion module is used for converting the logging depth into the oil column height so as to generate a relational expression between the saturation and the oil column height.
In one embodiment, the saturation height function generating unit includes:
the relational expression conversion module is used for converting the relational expression between the saturation and the height of the oil column into the relational expression between the saturation and the capillary pressure;
the relational expression generating module is used for generating a relational expression between a J function and the saturation according to the relational expression between the saturation and the capillary pressure;
and the saturation height function generation first module is used for generating a saturation height function of the target block according to a relation between the J function and the saturation.
In one embodiment, the saturation height function generating unit includes:
and the saturation height function generation second module is used for generating the saturation height function according to the height of the oil column, the water saturation and the porosity by utilizing a multiple regression method.
In one embodiment, the saturation calculation device based on conventional well log data further comprises:
and the well logging data selection unit is used for selecting the well data of which the altitude in the target block transition zone is within a preset range as the well logging data.
In a third aspect, the present invention provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the method for saturation calculation based on conventional well log data when executing the program.
In a fourth aspect, the present invention provides a computer readable storage medium having stored thereon a computer program which, when being executed by a processor, carries out the steps of the method for saturation calculation based on conventional logging data.
As can be seen from the above description, the saturation calculation method and apparatus based on conventional logging data according to the embodiments of the present invention can quickly determine the saturation height function by directly using the conventional logging interpretation result, porosity, permeability, and saturation, logging depth, and free water interface. Specifically, firstly, the basic data is sorted, the quality is checked and optimized, and the method comprises the following steps: the conventional logging information explains the quality check of porosity, permeability and water saturation of the result, the quality check of a free water interface and the optimization of the logging information so as to obtain more ideal basic logging information; then, the conversion of the logging data interpretation result parameters comprises the following steps: converting the logging depth into an oil column height, a saturation-logging depth relation graph, a saturation-oil column height (H) relation graph, a saturation-capillary pressure relation graph, a J function and a saturation relation graph; obtaining coefficients a and b through multivariate regression, and obtaining a saturation height function through proper conversion. Compared with the laboratory core data, the function establishing process hardly involves intermediate parameters (such as interfacial tension, contact angle and the like), and the conventional logging data are directly the real reflection of the reservoir characteristics under the oil reservoir condition, so that the saturation height function established by the logging data can reflect the actual characteristics of the reservoir more, and the function is more reliable and more effective when used in reservoir evaluation and reserve accurate calculation. And the actual characteristics of the reservoir can be reflected better by using the saturation height function established by the logging information. In addition, the method only utilizes the existing logging information, can obtain good effect, greatly reduces the cost and has great practical value, thereby creatively solving a great practical problem in the field. The invention can quickly determine the saturation height function. The method aims to lay a foundation for improving reservoir quantitative evaluation, reserve accurate calculation and three-dimensional geological model saturation field attribute modeling.
In conclusion, the method is easy to operate and the obtained result is reliable after practical application and inspection. Reliable technical support and guarantee are provided for improving the saturation field attribute modeling and reserve calculation in the reservoir fine evaluation and three-dimensional geological modeling in the exploration and development period of the oil and gas field.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a first flowchart illustrating a saturation calculation method based on conventional well log data according to an embodiment of the present invention;
FIG. 2 is a flow chart illustrating step 100 according to an embodiment of the present invention;
FIG. 3 is a first flowchart illustrating a step 200 according to an embodiment of the present invention;
FIG. 4 is a second flowchart illustrating a step 200 according to an embodiment of the present invention;
FIG. 5 is a second flowchart illustrating a saturation calculation method based on conventional well log data according to an embodiment of the present invention;
FIG. 6 is a flow chart illustrating a method for saturation calculation based on conventional well log data in an embodiment of the present invention;
FIG. 7 is a mental diagram of a saturation calculation method based on conventional well logging data in an embodiment of the present invention;
FIG. 8 is a graph of well logging interpretation saturation vs. TVD in an embodiment of the present invention;
FIG. 9 is a graph of an example of an embodiment of the present invention illustrating the relationship between saturation and TVD for an oilfield A logging;
FIG. 10 is a graph of column height versus saturation for different porosities and permeabilities in an example embodiment of the present invention;
FIG. 11 is a graph of oil column height versus saturation for an oil field A in accordance with an embodiment of the present invention;
FIG. 12 is a graph of capillary pressure versus saturation for an exemplary embodiment of the present invention;
FIG. 13 is a graph of oilfield tubular pressure A versus saturation for an example embodiment of the present invention;
FIG. 14 is a graph of J function versus saturation for an exemplary embodiment of the present invention;
FIG. 15 is a graph of the relationship between the J function and the saturation of the A field in an embodiment of the present invention;
FIG. 16 is a first schematic structural diagram of a saturation calculation apparatus based on conventional well log data according to an embodiment of the present invention;
fig. 17 is a block diagram showing the configuration of a first relational expression determining unit in the embodiment of the present invention;
FIG. 18 is a block diagram of a first structure of a saturation height function generating unit according to an embodiment of the present invention;
FIG. 19 is a block diagram of a saturation height function generating unit according to an embodiment of the present invention;
FIG. 20 is a second schematic structural diagram of a saturation calculation apparatus based on conventional well log data according to an embodiment of the present invention;
fig. 21 is a schematic structural diagram of an electronic device in an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In view of the prior art that some old oil fields lack laboratory core experimental data (especially capillary pressure experimental data) for economic and technical reasons in the early exploration period, the embodiment of the invention provides a specific implementation mode of a saturation calculation method based on conventional logging data, and the method specifically comprises the following contents, referring to fig. 1:
step 100: and determining a relational expression between the saturation and the height of the oil column according to the logging data of the target block.
It is understood that in the oil and gas industry, the saturation of a certain fluid refers to: the volume percentage of a certain fluid in the pores of the reservoir rock. It indicates the degree of pore space occupied by a certain fluid. If the rock is filled with several phases of fluid, the sum of the saturation levels of these phases of fluid is 1 (100%). The height of the oil column refers to the vertical distance from the top of the oil-gas reservoir to the oil (gas) water interface.
Step 200: and generating a saturation height function of the target block according to the relation between the saturation and the height of the oil column.
It can be understood that the saturation height function in step 200 refers to a relation between saturation and oil column height, the function establishment process hardly involves intermediate parameters (such as interfacial tension, contact angle and the like) compared with laboratory core data, and conventional logging data utilized by the function is directly a real reflection of reservoir characteristics under reservoir conditions, so that the saturation height function established by the logging data can reflect the actual characteristics of the reservoir more.
Step 300: and calculating the saturation of the target block according to the saturation height function.
The saturation calculated in step 300 can be very effective and cost reduced by using only existing well log data.
As can be seen from the above description, the saturation calculation method based on conventional logging data according to the embodiments of the present invention can quickly determine the saturation height function by directly using the conventional logging interpretation result, porosity, permeability, and saturation, logging depth, and free water interface. Specifically, firstly, the basic data is sorted, the quality is checked and optimized, and the method comprises the following steps: the conventional logging information explains the quality check of porosity, permeability and water saturation of the result, the quality check of a free water interface and the optimization of the logging information so as to obtain more ideal basic logging information; then, the conversion of the logging data interpretation result parameters comprises the following steps: converting the logging depth into an oil column height, a saturation-logging depth relation graph, a saturation-oil column height (H) relation graph, a saturation-capillary pressure relation graph, a J function and a saturation relation graph; obtaining coefficients a and b through multivariate regression, and obtaining a saturation height function through proper conversion. Compared with the laboratory core data, the function establishing process hardly involves intermediate parameters (such as interfacial tension, contact angle and the like), and the conventional logging data are directly the real reflection of the reservoir characteristics under the oil reservoir condition, so that the saturation height function established by the logging data can reflect the actual characteristics of the reservoir more, and the function is more reliable and more effective when used in reservoir evaluation and reserve accurate calculation. And the actual characteristics of the reservoir can be reflected better by using the saturation height function established by the logging information. In addition, the method only utilizes the existing logging information, can obtain good effect, greatly reduces the cost and has great practical value, thereby creatively solving a great practical problem in the field. The invention can quickly determine the saturation height function. The method aims to lay a foundation for improving reservoir quantitative evaluation, reserve accurate calculation and three-dimensional geological model saturation field attribute modeling.
In one embodiment, referring to fig. 2, step 100 specifically includes:
step 101: and determining a relational expression between the saturation and the logging depth according to the logging data.
It should be noted that the logging depth in step 101 refers to the vertical depth of the logging data point, i.e. tvd (true vertical depth).
Step 102: converting the log depth to the oil column height to generate a relationship between the saturation and oil column height.
Specifically, the logging depth is converted into the oil column height with reference to equation (1).
H=TVDFWL-TVD (1)
In the above formula: h is the height of the oil column above the free water interface, m; TVDFWLIs the vertical depth of the free water interface, m; TVD is the vertical depth of logging, m.
In one embodiment, referring to fig. 3, step 200 specifically includes:
step 201: and converting the relation between the saturation and the height of the oil column into the relation between the saturation and the capillary pressure.
Specifically, the relation (diagram) of the saturation and the oil column height (H) is converted into a relation (diagram) of the saturation and the capillary pressure with reference to the formula (2).
Figure BDA0002459352910000071
In the formula: pcRCapillary pressure under reservoir conditions, MPa; h is the height of the oil column above the free water interface, m; RHOwAnd RHOoDensity of formation water and crude oil in pore system, g/cm3
Step 202: and generating a relation between the J function and the saturation according to the relation between the saturation and the capillary pressure.
When the step 202 is implemented, it specifically includes: and (4) converting the saturation-capillary pressure relation graph into a J function-saturation relation graph according to the formulas (3) and (4), wherein the regression coefficients a and b in the formulas can be obtained through multiple regression.
Figure BDA0002459352910000081
Figure BDA0002459352910000082
In the above formula, J (Sw) -Leverett J function; pc is capillary pressure, Psi; sigma is surface tension, dynes/cm; k is the permeability, md; phi is porosity, v/v; a. b is a regression coefficient.
Step 203: and generating a saturation height function of the target block according to a relation between the J function and the saturation.
Specifically, the saturation height function can be generated using equations (5) and (6).
Figure BDA0002459352910000083
Figure BDA0002459352910000084
In one embodiment, referring to fig. 4, step 200 may also be implemented in another way:
step 201 a: and generating the saturation height function according to the height of the oil column, the water saturation and the porosity by utilizing a multiple regression method.
And (3) directly performing multiple regression on the parameters H, Sw and POR to obtain a saturation height function, namely:
Figure BDA0002459352910000085
it will be appreciated that if the permeability interpreted using the well log data is of high accuracy, then the saturation height function is calculated using steps 201 to 203, otherwise the saturation height function is calculated using step 201 a.
In an embodiment, referring to fig. 5, the method for calculating saturation based on conventional logging data further includes:
step 400: and selecting the well data of which the altitude is within a preset range in the transition zone of the target block as the logging data.
In selecting a well, typical well data with significant transitions and similar elevations are selected (as much as possible, all well data are not used in their entirety, so as not to cause undesirable results). In addition, the basic data needs to be sorted and quality checked, specifically: the porosity, permeability and water saturation of the interpretation result of the conventional logging data are checked by quality to obtain more ideal basic logging data; and (4) checking the quality of the free water interface, and comprehensively checking the reasonability of the free water interface by using data such as well logging, geology, oil deposit and the like, and perfecting.
As can be seen from the above description, the saturation calculation method based on conventional logging data according to the embodiments of the present invention can quickly determine the saturation height function by directly using the conventional logging interpretation result, porosity, permeability, and saturation, logging depth, and free water interface. Specifically, firstly, the basic data is sorted, the quality is checked and optimized, and the method comprises the following steps: the conventional logging information explains the quality check of porosity, permeability and water saturation of the result, the quality check of a free water interface and the optimization of the logging information so as to obtain more ideal basic logging information; then, the conversion of the logging data interpretation result parameters comprises the following steps: converting the logging depth into an oil column height, a saturation-logging depth relation graph, a saturation-oil column height (H) relation graph, a saturation-capillary pressure relation graph, a J function and a saturation relation graph; obtaining coefficients a and b through multivariate regression, and obtaining a saturation height function through proper conversion. Compared with the laboratory core data, the function establishing process hardly involves intermediate parameters (such as interfacial tension, contact angle and the like), and the conventional logging data are directly the real reflection of the reservoir characteristics under the oil reservoir condition, so that the saturation height function established by the logging data can reflect the actual characteristics of the reservoir more, and the function is more reliable and more effective when used in reservoir evaluation and reserve accurate calculation. And the actual characteristics of the reservoir can be reflected better by using the saturation height function established by the logging information. In addition, the method only utilizes the existing logging information, can obtain good effect, greatly reduces the cost and has great practical value, thereby creatively solving a great practical problem in the field. The invention can quickly determine the saturation height function. The method aims to lay a foundation for improving reservoir quantitative evaluation, reserve accurate calculation and three-dimensional geological model saturation field attribute modeling.
In conclusion, the method is easy to operate and the obtained result is reliable after practical application and inspection. Reliable technical support and guarantee are provided for improving the saturation field attribute modeling and reserve calculation in the reservoir fine evaluation and three-dimensional geological modeling in the exploration and development period of the oil and gas field.
To further illustrate the present solution, the present invention provides a specific application example of the saturation calculation method based on conventional logging data, taking the a oil field as an example, and the specific application example specifically includes the following contents, see fig. 6 and fig. 7.
S1: collecting and sorting basic data.
The basic data includes: conventional well log data interprets the results, such as porosity, permeability, and water saturation, and performs a quality check on the well log interpretation results. And comprehensively utilizing relevant data such as well logging, geology, oil deposit and the like to check the reasonability of the free water interface and perfect the free water interface.
S2: well log data is preferred.
In selecting a well, a typical well profile with a distinct transition zone and similar elevation in the a field is selected (as much as possible, all the well profiles are not used in their entirety so as not to cause undesirable results).
S3: and determining a relational expression between the saturation and the height of the oil column according to the logging data.
Specifically, a saturation-logging depth (TVD) relationship map is created from the preferred logging data in step S2 (see fig. 8 and 9). Then, the logging depth is converted into the oil column height according to the formula (1). And converting the established saturation-logging depth relation graph into a saturation-oil column height (H) relation graph (FIG. 10 and FIG. 11), and counting the porosity and permeability values of the data points.
S4: and generating a saturation height function according to the relation between the saturation and the height of the oil column.
In practice, step S4 includes two parallel schemes: if the precision of the permeability of the well logging interpretation is high, the saturation height function is established by using the steps S41 to S43, and if the permeability interpretation result is unreliable, the saturation height function is established by directly using the step S41 a.
S41: referring to fig. 12 and 13, the relationship between the saturation and the oil column height is converted into a relationship between the saturation and the capillary pressure.
S42: referring to fig. 14 and 15, the relationship between the J function and the saturation is generated from the relationship between the saturation and the capillary pressure.
S43: and generating a saturation height function of the target block according to a relation between the J function and the saturation.
S41 a: and generating the saturation height function according to the height of the oil column, the water saturation and the porosity by utilizing a multiple regression method.
Note that, in the above steps, attention needs to be paid to the unit of each parameter to select a reasonable conversion coefficient.
S5: and calculating the saturation of the oil field A according to the saturation height function.
As can be seen from the above description, the saturation calculation method based on conventional logging data according to the embodiments of the present invention can quickly determine the saturation height function by directly using the conventional logging interpretation result, porosity, permeability, and saturation, logging depth, and free water interface. Specifically, firstly, the basic data is sorted, the quality is checked and optimized, and the method comprises the following steps: the conventional logging information explains the quality check of porosity, permeability and water saturation of the result, the quality check of a free water interface and the optimization of the logging information so as to obtain more ideal basic logging information; then, the conversion of the logging data interpretation result parameters comprises the following steps: converting the logging depth into an oil column height, a saturation-logging depth relation graph, a saturation-oil column height (H) relation graph, a saturation-capillary pressure relation graph, a J function and a saturation relation graph; obtaining coefficients a and b through multivariate regression, and obtaining a saturation height function through proper conversion. Compared with the laboratory core data, the function establishing process hardly involves intermediate parameters (such as interfacial tension, contact angle and the like), and the conventional logging data are directly the real reflection of the reservoir characteristics under the oil reservoir condition, so that the saturation height function established by the logging data can reflect the actual characteristics of the reservoir more, and the function is more reliable and more effective when used in reservoir evaluation and reserve accurate calculation. And the actual characteristics of the reservoir can be reflected better by using the saturation height function established by the logging information. In addition, the method only utilizes the existing logging information, can obtain good effect, greatly reduces the cost and has great practical value, thereby creatively solving a great practical problem in the field. The invention can quickly determine the saturation height function. The method aims to lay a foundation for improving reservoir quantitative evaluation, reserve accurate calculation and three-dimensional geological model saturation field attribute modeling.
In conclusion, the method is easy to operate and the obtained result is reliable after practical application and inspection. Reliable technical support and guarantee are provided for improving the saturation field attribute modeling and reserve calculation in the reservoir fine evaluation and three-dimensional geological modeling in the exploration and development period of the oil and gas field.
Based on the same inventive concept, the embodiments of the present application further provide a saturation calculation apparatus based on conventional well logging data, which can be used to implement the methods described in the above embodiments, as described in the following embodiments. Because the principle of solving the problem of the conventional logging data-based saturation calculation device is similar to that of the conventional logging data-based saturation calculation method, the conventional logging data-based saturation calculation device can be implemented by the conventional logging data-based saturation calculation method, and repeated parts are not described again. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. While the system described in the embodiments below is preferably implemented in software, implementations in hardware, or a combination of software and hardware are also possible and contemplated.
The embodiment of the present invention provides a specific implementation of a conventional logging data-based saturation calculation apparatus capable of implementing a conventional logging data-based saturation calculation method, and referring to fig. 16, the conventional logging data-based saturation calculation apparatus specifically includes the following contents:
the first relational expression determining unit 10 is used for determining a relational expression between the saturation and the height of the oil column according to the logging data of the target block;
a saturation height function generating unit 20, configured to generate a saturation height function of the target block according to a relation between the saturation and the height of the oil column;
a saturation calculating unit 30, configured to calculate the saturation of the target block according to the saturation height function.
In one embodiment, referring to fig. 17, the first relation determination unit 10 includes:
a second relation determining module 101, configured to determine a relation between the saturation and the logging depth according to the logging data;
a logging depth conversion module 102, configured to convert the logging depth into the oil column height to generate a relation between the saturation and the oil column height.
In one embodiment, referring to fig. 18, the saturation height function generating unit 20 includes:
a relational expression conversion module 201, configured to convert a relational expression between the saturation and the oil column height into a relational expression between the saturation and the capillary pressure;
a relation generation module 202, configured to generate a relation between a J function and the saturation according to the relation between the saturation and the capillary pressure;
a saturation height function generating first module 203, configured to generate a saturation height function of the target block according to a relation between the J function and the saturation.
In one embodiment, referring to fig. 19, the saturation height function generating unit 20 includes:
a saturation height function generating second module 201a for generating the saturation height function according to the oil column height, the water saturation and the porosity by using a multiple regression method.
In one embodiment, referring to fig. 20, the saturation calculation apparatus based on conventional logging data further includes:
and the well logging data selecting unit 40 is used for selecting the well data of which the altitude in the target block transition zone is within a preset range as the well logging data.
As can be seen from the above description, the saturation calculation apparatus based on conventional logging data according to the embodiments of the present invention can quickly determine the saturation height function by directly using the conventional logging interpretation result, porosity, permeability, and saturation, logging depth, and free water interface. Specifically, firstly, the basic data is sorted, the quality is checked and optimized, and the method comprises the following steps: the conventional logging information explains the quality check of porosity, permeability and water saturation of the result, the quality check of a free water interface and the optimization of the logging information so as to obtain more ideal basic logging information; then, the conversion of the logging data interpretation result parameters comprises the following steps: converting the logging depth into an oil column height, a saturation-logging depth relation graph, a saturation-oil column height (H) relation graph, a saturation-capillary pressure relation graph, a J function and a saturation relation graph; obtaining coefficients a and b through multivariate regression, and obtaining a saturation height function through proper conversion. Compared with the laboratory core data, the function establishing process hardly involves intermediate parameters (such as interfacial tension, contact angle and the like), and the conventional logging data are directly the real reflection of the reservoir characteristics under the oil reservoir condition, so that the saturation height function established by the logging data can reflect the actual characteristics of the reservoir more, and the function is more reliable and more effective when used in reservoir evaluation and reserve accurate calculation. And the actual characteristics of the reservoir can be reflected better by using the saturation height function established by the logging information. In addition, the method only utilizes the existing logging information, can obtain good effect, greatly reduces the cost and has great practical value, thereby creatively solving a great practical problem in the field. The invention can quickly determine the saturation height function. The method aims to lay a foundation for improving reservoir quantitative evaluation, reserve accurate calculation and three-dimensional geological model saturation field attribute modeling.
In conclusion, the method is easy to operate and the obtained result is reliable after practical application and inspection. Reliable technical support and guarantee are provided for improving the saturation field attribute modeling and reserve calculation in the reservoir fine evaluation and three-dimensional geological modeling in the exploration and development period of the oil and gas field.
The embodiment of the present application further provides a specific implementation manner of an electronic device capable of implementing all steps in the saturation calculation method based on conventional well logging data in the foregoing embodiment, and referring to fig. 21, the electronic device specifically includes the following contents:
a processor (processor)1201, a memory (memory)1202, a communication interface 1203, and a bus 1204;
the processor 1201, the memory 1202 and the communication interface 1203 complete communication with each other through the bus 1204; the communication interface 1203 is configured to implement information transmission between related devices, such as a server-side device, a logging instrument, and a client device.
The processor 1201 is used to call the computer program in the memory 1202, and the processor executes the computer program to implement all the steps in the conventional logging data-based saturation calculation method in the above-described embodiments, for example, the processor executes the computer program to implement the following steps:
step 100: and determining a relational expression between the saturation and the height of the oil column according to the logging data of the target block.
Step 200: and generating a saturation height function of the target block according to the relation between the saturation and the height of the oil column.
Step 300: and calculating the saturation of the target block according to the saturation height function.
From the above description, it can be seen that the electronic device in the embodiment of the present application can quickly determine the saturation height function by directly using the conventional well logging interpretation result, the porosity, the permeability, and the saturation, the well logging depth, and the free water interface. Specifically, firstly, the basic data is sorted, the quality is checked and optimized, and the method comprises the following steps: the conventional logging information explains the quality check of porosity, permeability and water saturation of the result, the quality check of a free water interface and the optimization of the logging information so as to obtain more ideal basic logging information; then, the conversion of the logging data interpretation result parameters comprises the following steps: converting the logging depth into an oil column height, a saturation-logging depth relation graph, a saturation-oil column height (H) relation graph, a saturation-capillary pressure relation graph, a J function and a saturation relation graph; obtaining coefficients a and b through multivariate regression, and obtaining a saturation height function through proper conversion. Compared with the laboratory core data, the function establishing process hardly involves intermediate parameters (such as interfacial tension, contact angle and the like), and the conventional logging data are directly the real reflection of the reservoir characteristics under the oil reservoir condition, so that the saturation height function established by the logging data can reflect the actual characteristics of the reservoir more, and the function is more reliable and more effective when used in reservoir evaluation and reserve accurate calculation. And the actual characteristics of the reservoir can be reflected better by using the saturation height function established by the logging information. In addition, the method only utilizes the existing logging information, can obtain good effect, greatly reduces the cost and has great practical value, thereby creatively solving a great practical problem in the field. The invention can quickly determine the saturation height function. The method aims to lay a foundation for improving reservoir quantitative evaluation, reserve accurate calculation and three-dimensional geological model saturation field attribute modeling.
In conclusion, the method is easy to operate and the obtained result is reliable after practical application and inspection. Reliable technical support and guarantee are provided for improving the saturation field attribute modeling and reserve calculation in the reservoir fine evaluation and three-dimensional geological modeling in the exploration and development period of the oil and gas field.
Embodiments of the present application also provide a computer-readable storage medium capable of implementing all steps in the conventional well logging data-based saturation calculation method in the above embodiments, where the computer-readable storage medium stores thereon a computer program, and the computer program implements all steps of the conventional well logging data-based saturation calculation method in the above embodiments when being executed by a processor, for example, the processor implements the following steps when executing the computer program:
step 100: and determining a relational expression between the saturation and the height of the oil column according to the logging data of the target block.
Step 200: and generating a saturation height function of the target block according to the relation between the saturation and the height of the oil column.
Step 300: and calculating the saturation of the target block according to the saturation height function.
From the above description, it can be seen that the computer readable storage medium in the embodiments of the present application can quickly determine the saturation height function by directly using the conventional well logging interpretation, porosity, permeability, and saturation, and well logging depth, free water interface. Specifically, firstly, the basic data is sorted, the quality is checked and optimized, and the method comprises the following steps: the conventional logging information explains the quality check of porosity, permeability and water saturation of the result, the quality check of a free water interface and the optimization of the logging information so as to obtain more ideal basic logging information; then, the conversion of the logging data interpretation result parameters comprises the following steps: converting the logging depth into an oil column height, a saturation-logging depth relation graph, a saturation-oil column height (H) relation graph, a saturation-capillary pressure relation graph, a J function and a saturation relation graph; obtaining coefficients a and b through multivariate regression, and obtaining a saturation height function through proper conversion. Compared with the laboratory core data, the function establishing process hardly involves intermediate parameters (such as interfacial tension, contact angle and the like), and the conventional logging data are directly the real reflection of the reservoir characteristics under the oil reservoir condition, so that the saturation height function established by the logging data can reflect the actual characteristics of the reservoir more, and the function is more reliable and more effective when used in reservoir evaluation and reserve accurate calculation. And the actual characteristics of the reservoir can be reflected better by using the saturation height function established by the logging information. In addition, the method only utilizes the existing logging information, can obtain good effect, greatly reduces the cost and has great practical value, thereby creatively solving a great practical problem in the field. The invention can quickly determine the saturation height function. The method aims to lay a foundation for improving reservoir quantitative evaluation, reserve accurate calculation and three-dimensional geological model saturation field attribute modeling.
In conclusion, the method is easy to operate and the obtained result is reliable after practical application and inspection. Reliable technical support and guarantee are provided for improving the saturation field attribute modeling and reserve calculation in the reservoir fine evaluation and three-dimensional geological modeling in the exploration and development period of the oil and gas field.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the hardware + program class embodiment, since it is substantially similar to the method embodiment, the description is simple, and the relevant points can be referred to the partial description of the method embodiment.
The foregoing description has been directed to specific embodiments of this disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.
Although the present application provides method steps as described in an embodiment or flowchart, additional or fewer steps may be included based on conventional or non-inventive efforts. The order of steps recited in the embodiments is merely one manner of performing the steps in a multitude of orders and does not represent the only order of execution. When an actual apparatus or client product executes, it may execute sequentially or in parallel (e.g., in the context of parallel processors or multi-threaded processing) according to the embodiments or methods shown in the figures.
Although embodiments of the present description provide method steps as described in embodiments or flowcharts, more or fewer steps may be included based on conventional or non-inventive means. The order of steps recited in the embodiments is merely one manner of performing the steps in a multitude of orders and does not represent the only order of execution. When an actual apparatus or end product executes, it may execute sequentially or in parallel (e.g., parallel processors or multi-threaded environments, or even distributed data processing environments) according to the method shown in the embodiment or the figures. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, the presence of additional identical or equivalent elements in a process, method, article, or apparatus that comprises the recited elements is not excluded.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.

Claims (12)

1. A method of saturation calculation based on conventional well log data, comprising:
determining a relational expression between the saturation and the height of the oil column according to the logging data of the target block;
generating a saturation height function of the target block according to a relational expression between the saturation and the height of the oil column;
and calculating the saturation of the target block according to the saturation height function.
2. The method of claim 1, wherein determining the relationship between the saturation and the height of the oil column from the well log data for the target zone comprises:
determining a relational expression between the saturation and the logging depth according to the logging data;
converting the log depth to the oil column height to generate a relationship between the saturation and oil column height.
3. The method of claim 1, wherein generating the saturation height function for the target block according to the relationship between the saturation and the oil column height comprises:
converting the relation between the saturation and the height of the oil column into a relation between the saturation and capillary pressure;
generating a relation between a J function and the saturation according to the relation between the saturation and capillary pressure;
and generating a saturation height function of the target block according to a relation between the J function and the saturation.
4. The method of claim 1, wherein generating the saturation height function for the target block according to the relationship between the saturation and the oil column height comprises:
and generating the saturation height function according to the height of the oil column, the water saturation and the porosity by utilizing a multiple regression method.
5. The computing method of claim 1, further comprising: and selecting the well data of which the altitude is within a preset range in the transition zone of the target block as the logging data.
6. A saturation calculation apparatus based on conventional well log data, comprising:
the first relational expression determining unit is used for determining a relational expression between the saturation and the height of the oil column according to the logging data of the target block;
the saturation height function generating unit is used for generating a saturation height function of the target block according to a relational expression between the saturation and the height of the oil column;
and the saturation calculation unit is used for calculating the saturation of the target block according to the saturation height function.
7. The computing apparatus of claim 6, wherein the first relational expression determination unit comprises:
the second relational expression determining module is used for determining a relational expression between the saturation and the logging depth according to the logging data;
and the logging depth conversion module is used for converting the logging depth into the oil column height so as to generate a relational expression between the saturation and the oil column height.
8. The computing apparatus of claim 6, wherein the saturation height function generating unit comprises:
the relational expression conversion module is used for converting the relational expression between the saturation and the height of the oil column into the relational expression between the saturation and the capillary pressure;
the relational expression generating module is used for generating a relational expression between a J function and the saturation according to the relational expression between the saturation and the capillary pressure;
and the saturation height function generation first module is used for generating a saturation height function of the target block according to a relation between the J function and the saturation.
9. The computing apparatus of claim 6, wherein the saturation height function generating unit comprises:
and the saturation height function generation second module is used for generating the saturation height function according to the height of the oil column, the water saturation and the porosity by utilizing a multiple regression method.
10. The computing apparatus of claim 6, further comprising:
and the well logging data selection unit is used for selecting the well data of which the altitude in the target block transition zone is within a preset range as the well logging data.
11. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the steps of the method for saturation calculation based on conventional well log data according to any one of claims 1 to 5.
12. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method for saturation calculation based on conventional well log data of any one of claims 1 to 5.
CN202010315304.2A 2020-04-21 2020-04-21 Saturation calculation method and device based on conventional logging data Active CN111561312B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010315304.2A CN111561312B (en) 2020-04-21 2020-04-21 Saturation calculation method and device based on conventional logging data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010315304.2A CN111561312B (en) 2020-04-21 2020-04-21 Saturation calculation method and device based on conventional logging data

Publications (2)

Publication Number Publication Date
CN111561312A true CN111561312A (en) 2020-08-21
CN111561312B CN111561312B (en) 2023-09-26

Family

ID=72070570

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010315304.2A Active CN111561312B (en) 2020-04-21 2020-04-21 Saturation calculation method and device based on conventional logging data

Country Status (1)

Country Link
CN (1) CN111561312B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103210180A (en) * 2010-11-24 2013-07-17 雪佛龙美国公司 System and method for estimating fluid distribution in a subterranean reservoir
CN104179498A (en) * 2013-05-22 2014-12-03 中国石油化工股份有限公司 Zone-dividing saturation simulation method
US20160187532A1 (en) * 2013-08-06 2016-06-30 Schlumberger Technology Corporation Methods for determining a saturation-height function in oil and gas reservoirs
CN105781539A (en) * 2016-03-15 2016-07-20 中国石油大学(华东) Saturability well logging calculation method of tight oil and gas reservoir
WO2016164507A1 (en) * 2015-04-09 2016-10-13 Schlumberger Technology Corporation Oilfield reservior saturation and permeability modeling
CN107701180A (en) * 2017-09-19 2018-02-16 中国石油天然气集团公司 A kind of the original reservoir water saturation computational methods based on sealing core drilling
CN108241785A (en) * 2018-01-11 2018-07-03 中国海洋石油集团有限公司 A kind of fine characterizing method of Heterogeneous reservoir saturation field
CN108252710A (en) * 2018-02-26 2018-07-06 中国石油天然气股份有限公司 Method and device for determining water saturation of oil reservoir
CN110056346A (en) * 2019-04-17 2019-07-26 中国石油大港油田勘探开发研究院 A kind of three-dimensional original water saturation analogy method of the oil reservoir based on Long-term change trend function

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103210180A (en) * 2010-11-24 2013-07-17 雪佛龙美国公司 System and method for estimating fluid distribution in a subterranean reservoir
CN104179498A (en) * 2013-05-22 2014-12-03 中国石油化工股份有限公司 Zone-dividing saturation simulation method
US20160187532A1 (en) * 2013-08-06 2016-06-30 Schlumberger Technology Corporation Methods for determining a saturation-height function in oil and gas reservoirs
WO2016164507A1 (en) * 2015-04-09 2016-10-13 Schlumberger Technology Corporation Oilfield reservior saturation and permeability modeling
CN105781539A (en) * 2016-03-15 2016-07-20 中国石油大学(华东) Saturability well logging calculation method of tight oil and gas reservoir
CN107701180A (en) * 2017-09-19 2018-02-16 中国石油天然气集团公司 A kind of the original reservoir water saturation computational methods based on sealing core drilling
CN108241785A (en) * 2018-01-11 2018-07-03 中国海洋石油集团有限公司 A kind of fine characterizing method of Heterogeneous reservoir saturation field
CN108252710A (en) * 2018-02-26 2018-07-06 中国石油天然气股份有限公司 Method and device for determining water saturation of oil reservoir
CN110056346A (en) * 2019-04-17 2019-07-26 中国石油大港油田勘探开发研究院 A kind of three-dimensional original water saturation analogy method of the oil reservoir based on Long-term change trend function

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
曾大乾等: "确定砂岩油藏原始含油饱和度的毛管-测井法", 石油天然气学报 *
杨永亮: "一种储层三维地质建模技术方法——以哈萨克斯坦K 油田为例" *
杨永亮: "一种储层三维地质建模技术方法——以哈萨克斯坦K油田为例", 承德石油高等专科学校学报 *

Also Published As

Publication number Publication date
CN111561312B (en) 2023-09-26

Similar Documents

Publication Publication Date Title
US8359184B2 (en) Method, program and computer system for scaling hydrocarbon reservoir model data
CN103485770B (en) A kind of method and system asking for oil saturation based on artificial neural network
CN112946783B (en) Hydrate saturation determination method, device and equipment
CN108873065B (en) Sandstone high-quality reservoir prediction method and device
CN108629459B (en) Method and device for detecting hydrocarbon-containing pore of reservoir
Ebeltoft et al. Parameter based scal-analysing relative permeability for full field application
WO2015038162A1 (en) Pseudo phase production simulation: a signal processing approach to assess quasi-multiphase flow production via successive analogous step-function relative permeability controlled models in reservoir flow simulation in order to rank multiple petro-physical realizations
WO2018102199A1 (en) Method for prediction of live oil interfacial tension at reservoir conditions from dead oil measurements
CN110344822B (en) Determination method of water saturation of calcareous compact reservoir
CN113552036B (en) Method and device for determining oil-water relative permeability of medium-low pore permeability reservoir
CN111561312A (en) Saturation calculation method and device based on conventional logging data
US20220342105A1 (en) Method for validating rock formations compaction parameters using geomechanical modeling
CN113051697B (en) Method and device for generating evolution history curve of pore pressure of ancient gas reservoir stratum
CN104675391A (en) Method for calculating oil saturation of stratum
Maschio et al. A methodology to reduce uncertainty constrained to observed data
Amanipoor Providing a subsurface reservoir quality maps in oil fields by geostatistical methods
CN111648766A (en) Method and device for determining free water interface
CN113914853B (en) Method for determining water saturation of sandstone reservoir of deep water deposition system sand shale thin interbed
CN113138412B (en) Deep shale gas porosity earthquake prediction method and device
US11905807B2 (en) Workflow to model and condition water distribution in a hydrocarbon reservoir
CN112487599B (en) Method and device for determining single-well liquid production capacity of oil field in different water-containing periods
CN112746839B (en) Method and device for determining volume coefficient of crude oil
Johnson Upscaling of Saturation Height Functions
WO2023187452A1 (en) System and self-learning method for the interpretation of petrophysical parameters
CN115440312A (en) Method and device for quickly calculating crude oil density of high-pressure oil reservoir

Legal Events

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