CN111911140A - Method for simulating continental facies sedimentary reservoir - Google Patents
Method for simulating continental facies sedimentary reservoir Download PDFInfo
- Publication number
- CN111911140A CN111911140A CN202010685610.5A CN202010685610A CN111911140A CN 111911140 A CN111911140 A CN 111911140A CN 202010685610 A CN202010685610 A CN 202010685610A CN 111911140 A CN111911140 A CN 111911140A
- Authority
- CN
- China
- Prior art keywords
- layer
- sand
- sand body
- interlayer
- body group
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 208000035126 Facies Diseases 0.000 title claims abstract description 24
- 239000010410 layer Substances 0.000 claims abstract description 286
- 239000004576 sand Substances 0.000 claims abstract description 226
- 239000011229 interlayer Substances 0.000 claims abstract description 84
- 238000004088 simulation Methods 0.000 claims abstract description 64
- 238000004364 calculation method Methods 0.000 claims abstract description 44
- 238000012360 testing method Methods 0.000 claims abstract description 42
- 239000011435 rock Substances 0.000 claims abstract description 10
- 235000006679 Mentha X verticillata Nutrition 0.000 claims description 11
- 235000002899 Mentha suaveolens Nutrition 0.000 claims description 11
- 235000001636 Mentha x rotundifolia Nutrition 0.000 claims description 11
- 125000006850 spacer group Chemical group 0.000 claims description 9
- 238000005192 partition Methods 0.000 claims description 2
- 239000003921 oil Substances 0.000 abstract description 19
- 239000010779 crude oil Substances 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007788 roughening Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The disclosure provides a simulation method of a continental facies sedimentary reservoir, belonging to the field of oil and gas exploitation. The simulation method comprises the following steps: acquiring information of a test well in an area to be researched; determining a sand body layer in the area to be researched according to the information of the test well, wherein the sand body layer is a sand rock layer containing oil in a reservoir; dividing the reservoir into sand body group layers and an interlayer according to the distribution of each sand body layer, wherein each sand body group layer comprises at least one sand body layer, and the interlayer is a part between two adjacent sand body group layers; and carrying out coarsening calculation on the sand body group layer, and carrying out coarsening calculation on the interlayer so as to obtain a simulation model of the reservoir. Not only ensures the model precision, but also avoids the existence of a large number of invalid grids. The method can simulate the structural characteristics of the continental facies sedimentary reservoir with high precision, and further guide the exploitation of crude oil.
Description
Technical Field
The disclosure belongs to the field of oil and gas exploitation, and particularly relates to a simulation method of a continental facies sedimentary reservoir.
Background
Since the continental facies sedimentary reservoir is relatively complex, in order to enable more efficient exploitation of the continental facies sedimentary reservoir, it is generally necessary to perform digital simulation of the continental facies sedimentary reservoir. That is, the geological structure change of the reservoir is truly reflected by establishing a simulation model correspondingly matched with the reservoir, so that the extraction of crude oil in the reservoir is guided through the simulation result of the geological model. The establishment of the simulation model generally requires that based on the relevant data of a test well in a research area, the reservoir of the research area is subjected to grid division through computer numerical simulation software, calculation is performed according to the divided grids to establish a grid model, then simulation calculation is performed according to the grid model, and finally the simulation model is obtained.
In the related art, when a simulation model is established, the optimal mode is that no mesh is coarsened, however, if direct calculation is performed without coarsening, the problem that the calculation speed is slow or the calculation cannot be performed due to too large calculation amount of the mesh occurs. If the coarsening is carried out, the prior technical means is mainly to carry out the coarsening in equal proportion or the coarsening in equal thickness through modeling software, the problems that the geological model is distorted and the oil reservoir cannot be accurately described due to the heterogeneity of reservoir development can occur.
Disclosure of Invention
The embodiment of the disclosure provides a method for simulating a continental facies sedimentary reservoir, which can simulate the structural characteristics of the continental facies sedimentary reservoir with high precision so as to guide the exploitation of crude oil. The technical scheme is as follows:
the embodiment of the disclosure provides a simulation method of a continental facies sedimentary reservoir, which comprises the following steps:
acquiring information of a test well in an area to be researched;
determining a sand body layer in the area to be researched according to the information of the test well, wherein the sand body layer is a sand rock layer containing oil in a reservoir;
according to the distribution of each sand body layer, carrying out horizon division on a reservoir in the area to be researched so as to divide the reservoir into sand body group layers and an interlayer, wherein each sand body group layer comprises at least one sand body layer, and the interlayer is a part between two adjacent sand body group layers;
and carrying out coarsening calculation on the sand body group layer, and carrying out coarsening calculation on the interlayer so as to obtain the simulation model of the reservoir stratum.
In an implementation manner of the present disclosure, before the calculating while maintaining the actual model precision for the sand layer group and performing the coarsening calculation for the interlayer, so as to obtain the simulation model of the reservoir, the simulation method further includes:
calibrating the top boundary and the bottom boundary of all the sand body group layers;
determining a non-roughened thickness boundary in each sand bank layer and a roughened thickness boundary in each interlayer based on the top and bottom boundaries of each sand bank layer.
In another implementation of the present disclosure, the determining a non-coarsened thickness limit in each sand bank layer according to a top boundary and a bottom boundary of each sand bank layer includes:
respectively determining the highest top boundary position in each sand body group layer;
respectively determining the lowest bottom boundary position in each sand body group layer;
and calculating the thickness limit which is not coarsened in each sand body group layer according to the highest top boundary position and the lowest bottom boundary position of each sand body group layer.
In yet another implementation of the present disclosure, the determining the highest top boundary position in each sand group layer respectively includes:
and determining the minimum depth of the top boundary of each sand body group layer, and determining the position corresponding to the minimum depth of the top boundary of each sand body group layer as the highest top boundary position in the sand body group layers.
In yet another implementation of the present disclosure, the determining the lowest bottom-boundary position in each sand group layer respectively includes:
and determining the maximum depth of the bottom boundary of each sand group layer, and determining the position corresponding to the maximum depth of the bottom boundary of each sand group layer as the position of the lowest bottom boundary in the sand group layers.
In yet another implementation of the present disclosure, the calculating the thickness boundary in each sand bank layer without coarsening according to the highest top boundary position and the lowest bottom boundary position of each sand bank layer includes:
the thickness limit of the sand group layer without coarsening meets the following formula:
hn＝MaxBn-MinTn；
wherein h isnA non-roughened thickness limit in the nth sand bank layer; MaxBnIs the maximum depth of the bottom boundary, MinT, in the nth sand bank layernIs the minimum depth of the top boundary in the nth sand bank layer.
In yet another implementation of the present disclosure, the determining the roughened thickness limit in each of the spacer layers includes:
determining the highest top boundary position and the lowest bottom boundary position of each interlayer according to the highest top boundary position and the lowest bottom boundary position of the sand body group layer;
and calculating a coarsening thickness limit in the interlayer by using the highest top boundary position and the lowest bottom boundary position of the interlayer.
In yet another implementation of the present disclosure, the determining the highest top boundary position and the lowest bottom boundary position of each of the interlayer layers includes:
determining an upper sand body group layer and a lower sand body group layer which are adjacent to the interlayer, wherein the upper sand body group layer is the sand body group layer which is adjacent to the interlayer and is positioned above the interlayer, and the lower sand body group layer is the sand body group layer which is adjacent to the interlayer and is positioned below the interlayer;
determining the lowest bottom boundary position of the upper sand body group layer and the highest top boundary position of the lower sand body group layer;
and determining the highest top boundary position and the lowest bottom boundary position of each interlayer according to the lowest bottom boundary position of the upper sand body group layer and the highest top boundary position of the lower sand body group layer.
In another implementation of the present disclosure, the calculating the coarsened thickness limit in the interlayer using the highest top boundary position and the lowest bottom boundary position of the interlayer includes:
the roughened thickness margin in the spacer layer satisfies the following formula:
hJn＝MinTn+1-MaxBn；
wherein h isJnA roughened thickness limit in the nth spacer layer; MaxBnThe maximum depth, MinT, of the bottom boundary of the upper sand body group layer corresponding to the nth partition interlayern+1The minimum depth of the top boundary in the lower sand group layer corresponding to the nth interlayer is obtained.
In another implementation manner of the present disclosure, the performing a non-coarsening calculation on the sand layer group and performing a coarsening calculation on the interlayer layer to obtain a grid model of the reservoir includes:
according to the thickness boundary which is not coarsened in each sand body group layer, node division is carried out on the sand body group layers, and sand body grids are established;
coarsening the interlayer according to the coarsened thickness limit in each interlayer, thereby obtaining coarsened grids;
and calculating to obtain a simulation model of the reservoir according to the sand body grid and the coarsening grid.
The technical scheme provided by the embodiment of the disclosure has the following beneficial effects:
when the simulation method of the continental facies sedimentary reservoir provided by the embodiment of the disclosure is used for three-dimensional simulation of the continental facies reservoir, information of a test well in an area to be researched is firstly obtained, so that data support can be provided for subsequent determination of a sand body layer. And then, determining a sand body layer in the area to be researched according to the information of the test well, so that a corresponding sand rock layer in the area to be researched can be obtained, and preparing for subsequently establishing a simulation model.
Then, according to the distribution of each sand body layer, the reservoir stratum of the area to be researched is subjected to layer position division to obtain a plurality of sand body group layers and a plurality of interlayer layers, so that the layer position distribution of the area to be researched can be obtained preliminarily, and the approximate distribution of the sand body group layers and the non-sand body group layers (namely the interlayer layers) can be determined approximately. Then, according to the division result, the sand body group layer is subjected to non-coarsening calculation (namely calculation is carried out while maintaining the actual model precision), and the interlayer is subjected to coarsening calculation. Because the corresponding sand body layers in the sand body group layers belong to oil-containing sandstone layers. In order to ensure that a geological model of a reservoir obtained by subsequent simulation can accurately reflect the actual oil reservoir of the reservoir, a sand body group layer is required to be as fine as possible and cannot be subjected to rough calculation. And the interlayer belongs to a non-oil-bearing rock stratum and is generally a soil layer with relatively low porosity. Therefore, in order to reduce a large amount of operations in the simulation process, the interlayer can be used as an invalid part for carrying out coarsening calculation, and a high-precision reservoir simulation model can be obtained in a targeted manner through the division determination, and meanwhile, the problem of too large calculation amount caused by too many grid nodes can be avoided.
According to the simulation method of the continental facies sedimentary reservoir, when the simulation model is established, the calculation is carried out by keeping the actual model precision on the oil-containing sand body group layer in a targeted manner, and the coarsening calculation is carried out on the non-oil-containing interlayer, so that the high precision of the simulation model can be ensured when the coarsening calculation is carried out. That is to say, the number of invalid grid nodes of the simulation model can be reduced, the operation speed is increased, the calculation of the number of valid grid nodes can be ensured, and the precision of the simulation model is ensured.
The method for simulating the continental facies sedimentary reservoir provided by the embodiment has the advantages of simple steps and convenience in implementation, and can be widely used.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present disclosure, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present disclosure, and it is obvious for those skilled in the art to obtain other drawings based on the drawings without creative efforts.
Fig. 1 is a flow chart of a method for simulating a continental facies sedimentary reservoir according to an embodiment of the present disclosure;
FIG. 2 is a flow chart of another method for simulating a continental facies sedimentary reservoir provided by an embodiment of the present disclosure;
fig. 3 is a schematic diagram of distribution of a test well and a sand layer set provided in the embodiment of the present disclosure.
Detailed Description
To make the objects, technical solutions and advantages of the present disclosure more apparent, embodiments of the present disclosure will be described in detail with reference to the accompanying drawings.
The disclosed embodiment provides a simulation method of a continental facies sedimentary reservoir, as shown in fig. 1, the simulation method includes:
s101: information is obtained for a test well in an area to be investigated.
S102: and determining a sand body layer in the area to be researched according to the information of the test well, wherein the sand body layer is a sand rock layer containing oil in the reservoir.
S103: according to the distribution of each sand body layer, the reservoir in the area to be researched is subjected to layer position division so as to divide the reservoir into sand body group layers and an interlayer, wherein each sand body group layer comprises at least one sand body layer, and the interlayer is a position between two adjacent sand body group layers.
S104: and carrying out coarsening-free calculation on the sand body group layer and carrying out coarsening calculation on the interlayer so as to obtain a simulation model of the reservoir.
When the simulation method of the continental facies sedimentary reservoir provided by the embodiment of the disclosure is used for three-dimensional simulation of the continental facies reservoir, information of a test well in an area to be researched is firstly obtained, so that data support can be provided for subsequent determination of a sand body layer. And then, determining a sand body layer in the area to be researched according to the information of the test well, so that a corresponding sand rock layer in the area to be researched can be obtained, and preparing for subsequently establishing a simulation model.
Then, according to the distribution of each sand body layer, the reservoir stratum of the area to be researched is subjected to layer position division to obtain a plurality of sand body group layers and a plurality of interlayer layers, so that the layer position distribution of the area to be researched can be obtained preliminarily, and the approximate distribution of the sand body group layers and the non-sand body group layers (namely the interlayer layers) can be determined approximately. Then, according to the division result, the sand body group layer is subjected to non-coarsening calculation (namely calculation is carried out while maintaining the actual model precision), and the interlayer is subjected to coarsening calculation. Because the corresponding sand body layers in the sand body group layers belong to oil-containing sandstone layers. In order to ensure that a geological model of a reservoir obtained by subsequent simulation can accurately reflect the actual oil reservoir of the reservoir, a sand body group layer is required to be as fine as possible and cannot be subjected to rough calculation. And the interlayer belongs to a non-oil-bearing rock stratum and is generally a soil layer with relatively low porosity. Therefore, in order to reduce a large amount of operations in the simulation process, the interlayer can be used as an invalid part for carrying out coarsening calculation, and a high-precision reservoir simulation model can be obtained in a targeted manner through the division determination, and meanwhile, the problem of too large calculation amount caused by too many grid nodes can be avoided.
According to the simulation method of the continental facies sedimentary reservoir, when the simulation model is established, the calculation is carried out by keeping the actual model precision on the oil-containing sand body group layer in a targeted manner, and the coarsening calculation is carried out on the non-oil-containing interlayer, so that the high precision of the simulation model can be ensured when the coarsening calculation is carried out. That is to say, the number of invalid grid nodes of the simulation model can be reduced, the operation speed is increased, the calculation of the number of valid grid nodes can be ensured, and the precision of the simulation model is ensured.
The method for simulating the continental facies sedimentary reservoir provided by the embodiment has the advantages of simple steps and convenience in implementation, and can be widely used.
Fig. 2 is a flow chart of another method for simulating a terrestrial sedimentary reservoir according to an embodiment of the present disclosure, and in conjunction with fig. 2, the method for simulating a terrestrial sedimentary reservoir includes:
s201: information is obtained for a test well in an area to be investigated.
In the above implementation, the test wells in the area to be studied are selected test wells within the test area. The information for the test well typically includes formation electrical depth test data, porosity curves, depth resistivity (RT/RXO, RLLD/RLLS) curves, and the like.
It is easy to understand that by acquiring information of the test well, the information can be used as basic data of the geological structure simulation of the region to be researched, and further the geological structure simulation can be analyzed and calculated.
S202: and determining a sand body layer in the area to be researched according to the information of the test well, wherein the sand body layer is a sand rock layer containing oil in the reservoir.
In the above implementation, by determining the sand layer of the region to be researched, the approximate distribution of the oil reservoirs in the reservoir in the region to be researched can be known, so as to prepare for the subsequent simulation of the reservoir in the region to be researched.
S203: and according to the distribution of each sand body layer, performing layer division on the reservoir in the region to be researched so as to divide the reservoir into a sand body group layer and an interlayer. The sand body group layer comprises at least one sand body layer, and the interlayer is a part between two adjacent sand body group layers.
In the implementation manner, the reservoirs in the area to be researched can be divided into two categories by dividing the reservoir in the area to be researched into the horizon, namely, the reservoir in the area to be researched is a sand group layer, and the reservoir in the area to be researched is a non-sand group layer, namely, an interlayer. The compartments are non-oil bearing rock formations and are typically layers of relatively low porosity earth. By dividing the reservoir in this way, different grid divisions can be performed on the reservoir in a targeted manner in the subsequent steps, so that a high-precision simulation model can be obtained through calculation.
FIG. 3 is a schematic diagram of the distribution of the test wells and sand groups provided by the embodiment of the present disclosure, and in conjunction with FIG. 3, the reservoir shown in the diagram can be divided into n sand groups, i.e., S1、S2、S3……SnAnd each sand body group layer is provided with a plurality of sand body layers. The blank areas between the sand group layers are corresponding interlayer layers.
S204: calibrating the top boundary and the bottom boundary of all the sand body group layers; and determining the thickness limit of the sand body group layers without coarsening and the thickness limit of the coarsening in the interlayer according to the top boundary and the bottom boundary of each sand body group layer.
In the implementation mode, the depth of each sand body group layer relative to the ground can be roughly determined by calibrating the top boundary and the bottom boundary of all the sand body group layers in the test well. The thickness limits of the sand groups, which are not coarsened, and the coarsened thickness limits of the separation layers can be correspondingly determined based on the depth of each sand group layer.
The term "sand bank" as used herein refers to a portion of the sand bank dispersed in the test well. The top boundary of the sand group layer refers to the interface of the sand group layer closest to the ground when each sand group layer is positioned in all the test wells. The bottom boundary of the sand group layer refers to the interface of the sand layer farthest from the ground when each sand group layer is positioned in all the test wells.
Step S204 is implemented by:
4.1: and respectively determining the highest top boundary position in each sand body group layer.
Exemplarily, step 4.1 is achieved by;
and determining the minimum depth of the top boundary of each sand group layer, and determining the position corresponding to the minimum depth of the top boundary of each sand group layer as the highest top boundary position in the sand group layers.
Illustratively, the minimum depth of the top boundary in each sand bank layer satisfies the following formula:
MinTn＝min(W1Tn，W2Tn，W3Tn……WmTn)； (1)
wherein, MinTnIs the minimum depth of the top boundary of the nth sand bank layer; wmTnThe depth corresponding to the top boundary in the nth sand group layer in the mth test well.
In the implementation mode, the depths corresponding to the top boundaries of each sand body group layer in each test well can be compared through the formula (1), and then the minimum depth of the top boundaries of each sand body group layer is rapidly calculated.
Continuing with FIG. 3, illustratively, there are m test wells (3 shown) in FIG. 3, labeled W1, W2, W3 … … W one by onemWherein, according to the related information of the test well, the corresponding sand body layer (the shaded part in the figure is the corresponding sand body layer) can be determined, and according to the related information of the test well, the sand body layer can be determinedThe sand body layer divides the layer of the reservoir in the area to be researched, then the top boundary and the bottom boundary of the sand body group layer in each test well are respectively calibrated, namely the depths corresponding to the top boundary and the bottom boundary of the sand body group layer in all the test wells in each sand body group layer are calibrated, namely the calibration positions can comprise the depths corresponding to the top boundary and the bottom boundary of all the sand body group layers in each test well, namely W1T1、W1B1，W2T1、W2B1……WmT1、WmB1；W1T2、W1B2，W2T2、W2B2……WmT2、WmB2；……；W1Tn、W1Bn，W2Tn、W2Bn……WmTn、WmBn(WmTnFor the depth, W, corresponding to the top boundary in the nth sand bank layer located in the mth test wellmBnThe depth corresponding to the bottom boundary in the nth sand bank layer located in the mth test well).
Then, by using the above formula (1), the minimum value of the depths corresponding to the top boundaries of all the sand bank layers in each sand bank layer is obtained, so that the minimum depth of the top boundaries in each sand bank layer can be determined.
4.2: and respectively determining the position of the lowest bottom boundary in each sand body group layer.
Exemplarily, step 4.2 is implemented by:
determining the maximum depth of the bottom boundary of each sand group layer, and determining the position corresponding to the maximum depth of the bottom boundary of each sand group layer as the lowest bottom boundary position in the sand group layers;
the maximum depth of the bottom boundary in each sand group layer satisfies the following formula:
MaxBn＝max(W1Bn，W2Bn，W3Bn……WmBn)； (2)
wherein, MaxBnThe maximum depth of the bottom boundary of the nth sand group layer; wmBnThe depth corresponding to the bottom boundary in the nth sand group layer in the mth test well.
In the implementation manner, the depths corresponding to the bottom boundaries of all the sand body group layers in each sand body group layer can be compared through the formula (2), and then the maximum depth of the bottom boundaries of each sand body group layer is rapidly calculated.
With continued reference to fig. 3, the maximum depth of the bottom boundary of each sand bank layer can be determined by calibrating the depths corresponding to the top boundary and the bottom boundary in all the test wells in each sand bank layer, and then calculating the maximum value of the depths corresponding to the bottom boundary in each sand bank layer by using the above formula (2).
4.3: and calculating the thickness limit which is not coarsened in each sand body group layer according to the highest top boundary position and the lowest bottom boundary position of each sand body group layer.
Illustratively, step 4.3 is implemented by:
the thickness boundary in the sand bank that is not coarsened satisfies the following formula:
hn＝MaxBn-MinTn； (3)
wherein h isnA thickness limit of the nth sand group layer which is not coarsened; MaxBnDepth of bottom boundary, MinT, in the nth sand groupnIs the depth of the top boundary in the nth sand bank layer.
In the implementation manner, the thickness limit which is not coarsened in each sand body group layer can be quickly and simply calculated through the formula (3).
Continuing with FIG. 3, the thickness margin of each sand section that is not coarsened, i.e., h, is obtained by subtracting the minimum depth of the top boundary from the maximum depth of the bottom boundary in each sand sectionn。
4.4: and determining the highest top boundary position and the lowest bottom boundary position of each interlayer according to the highest top boundary position and the lowest bottom boundary position of the sand body group layer.
Step 4.4 is realized by the following method:
(1) determining an upper sand body group layer and a lower sand body group layer which are adjacent to the interlayer, wherein the upper sand body group layer is a sand body group layer which is adjacent to the interlayer and is positioned above the interlayer, and the lower sand body group layer is a sand body group layer which is adjacent to the interlayer and is positioned below the interlayer;
(2) determining the lowest bottom boundary position of the upper sand body group layer and the highest top boundary position of the lower sand body group layer;
(3) and determining the highest top boundary position and the lowest bottom boundary position of each interlayer through the lowest bottom boundary position of the upper sand body group layer and the highest top boundary position of the lower sand body group layer.
4.5: and calculating the coarsening thickness limit in the interlayer by using the highest top boundary position and the lowest bottom boundary position of the interlayer.
Exemplarily, step 4.5 is achieved by;
the thickness margin of the roughening in the spacer layer satisfies the following formula:
hJn＝MinTn+1-MaxBn；
wherein h isJnA roughened thickness limit in the nth spacer layer; MaxBnIs the maximum depth, MinT, of the bottom boundary of the upper sand group layer corresponding to the nth interlayern+1The minimum depth of the top boundary in the lower sand group layer corresponding to the nth interlayer.
Referring to FIG. 3, it is shown from FIG. 3 that the thickness of each spacer layer is limited to a distance h from two adjacent sand bank layersJn。
S205: and carrying out coarsening-free calculation on the sand body group layer and carrying out coarsening calculation on the interlayer so as to obtain a simulation model of the reservoir.
Step S205 is implemented by:
according to the thickness boundary which is not coarsened in each sand body group layer, node division is carried out on the sand body group layers, and a sand body grid is established;
coarsening the interlayer according to the coarsened thickness limit in each interlayer, thereby obtaining coarsened grids;
and calculating to obtain a simulation model of the reservoir according to the sand body grid and the coarsening grid.
For example, according to fig. 3, for sand group layer No. n:
according to the calibration of the sand body group layers in each test well, the minimum depth of the top boundary of each sand body group layer in all the test wells is firstly solved: MinT (minimum time to live)n＝min(W1Tn，W2Tn，W3Tn……WmTn)；
And (3) solving the maximum depth of the bottom boundary of each sand body group layer in all the test wells: MaxBn＝max(W1Bn，W2Bn，W3Bn……WmBn)；
Subtracting the corresponding maximum depth of the bottom boundary and the minimum depth of the top boundary to obtain the non-coarsening thickness boundary of each sand body group layer: h isn＝MaxBn-MinTnAnd a sand body grid model (namely a sand body grid) is established by keeping certain high precision in the thickness range, so that high-precision description of a reservoir is realized, namely, the high-precision requirement of the simulation model is ensured by calculating the original actual model precision of a sand body layer.
Then, the coarsened thickness limit of each interlayer is obtained, and a coarsened grid model (namely a coarsened grid) of the interlayer is established. Aiming at the interlayer between the n-th sand body group layer and the n + 1-th sand body group layer: using MinT in n +1 sand group layern+1Subtracting MaxB in the n-th sand body group layernAnd obtaining the coarsening thickness limit of the coarsening grid model corresponding to the interlayer between the n-th sand body and the n + 1-th sand body: h isJn＝MinTn+1-MaxBnIn the thickness range, a plurality of longitudinal grid nodes of the geological model can be combined into one grid, the number of invalid grid nodes of the geological model is reduced, and the operation speed is increased.
The specific operation of the simulation method for a continental facies sedimentary reservoir provided by the embodiment of the present disclosure is illustrated below by specific examples:
example (c):
reservoir fault block in certain areaThe total number of the inner 3 wells is W1、W2、W3And if the reservoir simulation is carried out in the fault block, establishing a simulation model to carry out a selective coarsening process, wherein the fault block is divided into 3 series of layers (namely comprises 3 sand body group layers):
(1) in the first layer system, W1Depth W corresponding to top boundary of the sand body group layer in well1T11612 m, depth W corresponding to the bottom boundary1B1Is 1621 m, W2Depth W corresponding to top boundary in the sand body group layer in well2T11608 meters, depth W corresponding to the bottom boundary2B1Is 1616 m, W3Depth W corresponding to top boundary of the sand body group layer in well3T11617 meters, depth W corresponding to the bottom boundary2B11629 meters, the minimum depth of the top boundary of the first set of layer sand group layers is calculated by the formula (1): MinT (minimum time to live)1＝min(W1T1，W2T1，W3T1) The maximum depth of the bottom boundary of the sand group layer is 1608 meters, and is calculated by the formula (2): MaxB1＝max(W1Bn，W2Bn，W3Bn……WmBn) 1629 m.
(2) In the second layer system, W1Depth W corresponding to top boundary of well sand body group layer1T11642 m, depth W corresponding to the bottom boundary1B1Is 1649 m, W2Depth W corresponding to top boundary of the sand body group layer in well2T11644 m, depth W corresponding to bottom boundary2B11657 m, W3Depth W corresponding to top boundary of the sand body group layer in well3T11643 m, depth W corresponding to the bottom boundary2B1At 1652 meters, the minimum depth of the top boundary in the second casing sand group layer is calculated by the formula (1): MinT (minimum time to live)2＝min(W1T1，W2T1，W3T1) 1642 m, the maximum depth of the bottom boundary in the sand group layer is calculated by formula (2): MaxB2＝max(W1Bn，W2Bn，W3Bn……WmBn) 1657 m.
(3) The third layer system is inner, W1Depth W corresponding to top boundary of well sand body group layer1T11689 m, depth W corresponding to bottom boundary1B1Is 1709 m, W2Depth W corresponding to top boundary of the sand body group layer in well2T11694 m, depth W corresponding to bottom boundary2B1Is 1706 m, W3Depth W corresponding to top boundary of the sand body group layer in well3T11690 m, depth W corresponding to the bottom boundary2B1The minimum depth of the top boundary of the sand group layer of the third casing system is 1711 m, and is calculated by the formula (1): MinT (minimum time to live)3＝min(W1T1，W2T1，W3T1) The maximum depth of the bottom boundary of the sand group layer is 1689 meters, and is calculated by the formula (2): MaxB3＝max(W1Bn，W2Bn，W3Bn……WmBn) 1711 m;
(4) solving the non-coarsening thickness limit of the longitudinal sand body grid model of each sand body group layer: h is1＝MaxB1-MinT11629-2＝MaxB2-MinT21657-15 m, h3＝MaxB3-MinT31711-;
(5) obtaining a coarsening thickness limit of the interlayer longitudinal coarsening grid model: h isJ1＝MinT2-MaxB11642 and 1629 13 m, hJ2＝MinT3-MaxB21689-.
The above description is meant to be illustrative of the principles of the present disclosure and not to be taken in a limiting sense, and any modifications, equivalents, improvements and the like that are within the spirit and scope of the present disclosure are intended to be included therein.
Claims (10)
1. A method of simulating a continental facies depositional reservoir, the method comprising:
acquiring information of a test well in an area to be researched;
determining a sand body layer in the area to be researched according to the information of the test well, wherein the sand body layer is a sand rock layer containing oil in a reservoir;
according to the distribution of each sand body layer, carrying out horizon division on a reservoir in the area to be researched so as to divide the reservoir into sand body group layers and an interlayer, wherein each sand body group layer comprises at least one sand body layer, and the interlayer is a part between two adjacent sand body group layers;
and carrying out coarsening calculation on the sand body group layer, and carrying out coarsening calculation on the interlayer so as to obtain the simulation model of the reservoir stratum.
2. The simulation method of claim 1, wherein prior to performing the non-upscaling calculation on the sand bank layer and the upscaling calculation on the interbedded layer to obtain the simulation model of the reservoir, the simulation method further comprises:
calibrating the top boundary and the bottom boundary of all the sand body group layers;
determining a non-roughened thickness boundary in each sand bank layer and a roughened thickness boundary in each interlayer based on the top and bottom boundaries of each sand bank layer.
3. The simulation method of claim 2, wherein determining the non-coarsened thickness limits in each sand bank layer according to the top and bottom boundaries of each sand bank layer comprises:
respectively determining the highest top boundary position in each sand body group layer;
respectively determining the lowest bottom boundary position in each sand body group layer;
and calculating the thickness limit which is not coarsened in each sand body group layer according to the highest top boundary position and the lowest bottom boundary position of each sand body group layer.
4. The simulation method of claim 3, wherein the separately determining a highest top boundary position in each of the sand bank layers comprises:
and determining the minimum depth of the top boundary of each sand body group layer, and determining the position corresponding to the minimum depth of the top boundary of each sand body group layer as the highest top boundary position in the sand body group layers.
5. The simulation method of claim 4, wherein the separately determining a lowest bottom boundary position in each of the sand bank layers comprises:
and determining the maximum depth of the bottom boundary of each sand group layer, and determining the position corresponding to the maximum depth of the bottom boundary of each sand group layer as the position of the lowest bottom boundary in the sand group layers.
6. The simulation method of claim 5, wherein calculating the non-coarsened thickness limits in each sand bank layer according to the highest top boundary position and the lowest bottom boundary position of each sand bank layer comprises:
the thickness limit of the sand group layer without coarsening meets the following formula:
hn＝MaxBn-MinTn；
wherein h isnA non-roughened thickness limit in the nth sand bank layer; MaxBnIs the maximum depth of the bottom boundary, MinT, in the nth sand bank layernIs the minimum depth of the top boundary in the nth sand bank layer.
7. The simulation method of claim 3, wherein said determining a roughened thickness limit in each of said spacers comprises:
determining the highest top boundary position and the lowest bottom boundary position of each interlayer according to the highest top boundary position and the lowest bottom boundary position of the sand body group layer;
and calculating a coarsening thickness limit in the interlayer by using the highest top boundary position and the lowest bottom boundary position of the interlayer.
8. The simulation method of claim 7, wherein the determining a highest top boundary position and a lowest bottom boundary position for each of the compartments comprises:
determining an upper sand body group layer and a lower sand body group layer which are adjacent to the interlayer, wherein the upper sand body group layer is the sand body group layer which is adjacent to the interlayer and is positioned above the interlayer, and the lower sand body group layer is the sand body group layer which is adjacent to the interlayer and is positioned below the interlayer;
determining the lowest bottom boundary position of the upper sand body group layer and the highest top boundary position of the lower sand body group layer;
and determining the highest top boundary position and the lowest bottom boundary position of each interlayer according to the lowest bottom boundary position of the upper sand body group layer and the highest top boundary position of the lower sand body group layer.
9. The simulation method of claim 8, wherein calculating the coarsened thickness limit in the interlayer using the highest top boundary position and the lowest bottom boundary position of the interlayer comprises:
the roughened thickness margin in the spacer layer satisfies the following formula:
hJn＝MinTn+1-MaxBn；
wherein h isJnA roughened thickness limit in the nth spacer layer; MaxBnThe maximum depth, MinT, of the bottom boundary of the upper sand body group layer corresponding to the nth partition interlayern+1The minimum depth of the top boundary in the lower sand group layer corresponding to the nth interlayer is obtained.
10. The simulation method of claim 1, wherein the performing a non-coarsening calculation on the sand bank layer and a coarsening calculation on the interlayer layer to obtain the grid model of the reservoir comprises:
according to the thickness boundary which is not coarsened in each sand body group layer, node division is carried out on the sand body group layers, and sand body grids are established;
coarsening the interlayer according to the coarsened thickness limit in each interlayer, thereby obtaining coarsened grids;
and calculating to obtain a simulation model of the reservoir according to the sand body grid and the coarsening grid.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010685610.5A CN111911140B (en) | 2020-07-16 | 2020-07-16 | Simulation method of land phase sedimentary reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010685610.5A CN111911140B (en) | 2020-07-16 | 2020-07-16 | Simulation method of land phase sedimentary reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111911140A true CN111911140A (en) | 2020-11-10 |
| CN111911140B CN111911140B (en) | 2023-10-27 |
Family
ID=73280999
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010685610.5A Active CN111911140B (en) | 2020-07-16 | 2020-07-16 | Simulation method of land phase sedimentary reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111911140B (en) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100017181A1 (en) * | 2008-06-27 | 2010-01-21 | Thibaud Mouton | Method for Constructing a Hybrid Grid From a CPG Type Grid |
| US20100076738A1 (en) * | 2008-09-19 | 2010-03-25 | Chevron U.S.A. Inc. | Computer-implemented systems and methods for use in modeling a geomechanical reservoir system |
| US20100185428A1 (en) * | 2008-12-31 | 2010-07-22 | Vink Jeroen Cornelis | Method and system for simulating fluid flow in an underground formation with uncertain properties |
| CN103679540A (en) * | 2013-12-31 | 2014-03-26 | 中国石油大学(北京) | Reservoir model coarsening method and system based on reservoir architecture element interfaces |
| CN104331925A (en) * | 2014-09-28 | 2015-02-04 | 长江大学 | Permeability coarsening method considering interlayer influence |
| CN105298484A (en) * | 2015-10-28 | 2016-02-03 | 中国石油化工股份有限公司 | Remaining oil description method based on reservoir architecture |
| US20160131800A1 (en) * | 2014-11-07 | 2016-05-12 | Schlumberger Technology Corporation | Modeling fluid-conducting fractures in reservoir simulation grids |
| CN107366534A (en) * | 2017-08-10 | 2017-11-21 | 中国石油天然气股份有限公司 | It is roughened the determination method and apparatus of permeability |
| CN108131134A (en) * | 2017-11-22 | 2018-06-08 | 中国石油天然气股份有限公司 | The determining method and apparatus of remaining oil |
| CN109426687A (en) * | 2017-08-29 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of reservoir modeling grid conversion method and system |
| GB201907855D0 (en) * | 2019-06-03 | 2019-07-17 | Cognitive Geology Ltd | Apparatus method and computer-program product for processing geological data |
| CN110838175A (en) * | 2019-11-07 | 2020-02-25 | 成都北方石油勘探开发技术有限公司 | Geological model building method for gas injection oil reservoir development |
-
2020
- 2020-07-16 CN CN202010685610.5A patent/CN111911140B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100017181A1 (en) * | 2008-06-27 | 2010-01-21 | Thibaud Mouton | Method for Constructing a Hybrid Grid From a CPG Type Grid |
| US20100076738A1 (en) * | 2008-09-19 | 2010-03-25 | Chevron U.S.A. Inc. | Computer-implemented systems and methods for use in modeling a geomechanical reservoir system |
| US20100185428A1 (en) * | 2008-12-31 | 2010-07-22 | Vink Jeroen Cornelis | Method and system for simulating fluid flow in an underground formation with uncertain properties |
| CN103679540A (en) * | 2013-12-31 | 2014-03-26 | 中国石油大学(北京) | Reservoir model coarsening method and system based on reservoir architecture element interfaces |
| CN104331925A (en) * | 2014-09-28 | 2015-02-04 | 长江大学 | Permeability coarsening method considering interlayer influence |
| US20160131800A1 (en) * | 2014-11-07 | 2016-05-12 | Schlumberger Technology Corporation | Modeling fluid-conducting fractures in reservoir simulation grids |
| CN105298484A (en) * | 2015-10-28 | 2016-02-03 | 中国石油化工股份有限公司 | Remaining oil description method based on reservoir architecture |
| CN107366534A (en) * | 2017-08-10 | 2017-11-21 | 中国石油天然气股份有限公司 | It is roughened the determination method and apparatus of permeability |
| CN109426687A (en) * | 2017-08-29 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of reservoir modeling grid conversion method and system |
| CN108131134A (en) * | 2017-11-22 | 2018-06-08 | 中国石油天然气股份有限公司 | The determining method and apparatus of remaining oil |
| GB201907855D0 (en) * | 2019-06-03 | 2019-07-17 | Cognitive Geology Ltd | Apparatus method and computer-program product for processing geological data |
| CN110838175A (en) * | 2019-11-07 | 2020-02-25 | 成都北方石油勘探开发技术有限公司 | Geological model building method for gas injection oil reservoir development |
Non-Patent Citations (1)
| Title |
|---|
| 叶俊: "松辽盆地州葡 32-X50 井区储层地质建模及数值模拟", pages 1 - 45 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111911140B (en) | 2023-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8095349B2 (en) | Dynamic updating of simulation models | |
| CN104765065B (en) | Method and device for determining sedimentary stratum ancient landform | |
| CN105651676B (en) | A kind of reservoir heterogeneity characterizing method under horizontal well rule well pattern | |
| CN107808413A (en) | A kind of three-dimensional geological modeling method based on GOCAD | |
| US7340385B2 (en) | Method and program storage device for generating grids representing the architecture of fluvial reservoirs | |
| CN110056346B (en) | Oil reservoir three-dimensional original water saturation simulation method based on trend change function | |
| CN109102564B (en) | Coupling modeling method for numerical model of complex geologic body | |
| CN110244021B (en) | Stratum layering method based on anisotropic interpolation | |
| CN111383336A (en) | Three-dimensional geological model construction method | |
| van Dijk et al. | Reduced-complexity probabilistic reconstruction of alluvial aquifer stratigraphy, and application to sedimentary fans in northwestern India | |
| CN114549774A (en) | Three-dimensional stratum modeling method based on drilling data | |
| CN108121009A (en) | Development late stage complex fault block velocity-variant depth mapping method | |
| CN110390154A (en) | A method of improving Complex reservoir reservoir numerical simulation efficiency | |
| CN110863818A (en) | Description method and device for residual oil/gas distribution | |
| Zhang et al. | Method and application of urban 3D rapid modeling of geology based on CAD Borehole Logs | |
| Popovs et al. | Bedrock surface topography of Latvia | |
| CN111911140A (en) | Method for simulating continental facies sedimentary reservoir | |
| CN116449423A (en) | Method, device, equipment and storage medium for analyzing single sand body configuration of braided river | |
| CN105089659A (en) | Conglomerate oil reservoir seepage unit recognizing method | |
| CN115035258A (en) | Efficient urban three-dimensional geological modeling method based on CAD (computer-aided design) drilling histogram | |
| CN115019004A (en) | Oil-gas reserve calculation method based on reservoir geological modeling | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| CN112649893B (en) | Multi-data and multi-parameter fusion modeling method and system for thin reservoir | |
| CN113514879A (en) | 'ear' layer identification method | |
| CN113361111A (en) | Partitioning method of low-permeability reservoir characteristic model |
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 |