CN111810137B - Drawing method of granite buried hill reservoir development rule chart - Google Patents
Drawing method of granite buried hill reservoir development rule chart Download PDFInfo
- Publication number
- CN111810137B CN111810137B CN202010684890.8A CN202010684890A CN111810137B CN 111810137 B CN111810137 B CN 111810137B CN 202010684890 A CN202010684890 A CN 202010684890A CN 111810137 B CN111810137 B CN 111810137B
- Authority
- CN
- China
- Prior art keywords
- data
- oil
- fracture
- reservoir
- production
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011161 development Methods 0.000 title claims abstract description 81
- 239000010438 granite Substances 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000004519 manufacturing process Methods 0.000 claims abstract description 123
- 238000004088 simulation Methods 0.000 claims abstract description 68
- 239000003921 oil Substances 0.000 claims description 183
- 238000010586 diagram Methods 0.000 claims description 21
- 238000009826 distribution Methods 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000035699 permeability Effects 0.000 claims description 10
- 239000010779 crude oil Substances 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 239000013589 supplement Substances 0.000 claims description 5
- 230000009469 supplementation Effects 0.000 claims description 4
- 239000012530 fluid Substances 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 239000011435 rock Substances 0.000 claims description 3
- 238000011160 research Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000012407 engineering method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 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
- E21B25/00—Apparatus for obtaining or removing undisturbed cores, e.g. core barrels or core extractors
-
- 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/10—Geometric CAD
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Geometry (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geophysics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a drawing method of a granite buried hill oil reservoir development rule plate, which comprises the following steps of firstly obtaining oil reservoir data of M granite oil reservoirs, establishing a discrete fracture numerical simulation model according to the oil reservoir data of each oil reservoir, carrying out numerical simulation under different production pressure differences by using the discrete fracture numerical simulation model, obtaining data of accumulated oil production of different production pressure differences under different oil reservoir data, drawing the development rule plate according to the data of the accumulated oil production, and predicting the productivity of the current oil reservoir according to the actual oil reservoir data by referring to the development rule plate before the specific oil reservoir is developed, thereby formulating a proper development scheme; in the invention, the oil deposit data of a plurality of oil deposits are adopted to generate different discrete fracture models for numerical simulation, so that the development rule prediction chart has stronger applicability, can be suitable for the productivity prediction of most granite buried hill oil deposits, and has convenient use and high accuracy.
Description
Technical Field
The invention relates to the field of a prediction research method for a granite buried hill oil reservoir development rule, in particular to a drawing method for a granite buried hill oil reservoir development rule chart.
Background
With the increase of oil demand and the improvement of industrial exploitation level in recent years, the exploitation of the conventional sandstone oil reservoir does not meet the requirement of the market, and the exploration and development of the buried hill oil reservoir become an important oil production source. The granite buried hill oil reservoir is a buried hill oil reservoir which is widely distributed. The development rule of the granite buried hill reservoir is different from that of the conventional sandstone reservoir, the permeability and porosity of the buried hill matrix are far lower than those of the conventional sandstone reservoir, cracks are main channels for oil-gas migration of the granite buried hill reservoir, the stratum is influenced by weathering, denudation, structure and other actions, pores, holes and cracks in the stratum are randomly distributed, and the reservoir heterogeneity and discontinuity of the stratum are extremely strong. The conventional oil reservoir development rule prediction means is not suitable for the development rule prediction of the granite buried hill oil reservoir. Particularly, in the early production development stage of the oil reservoir with little data or even no available data, the distribution of natural cracks has great influence on the production rule of the granite buried hill oil reservoir, and the existing production rule prediction method is difficult to be applied to the granite buried hill oil reservoir. Therefore, the reservoir production and exploitation can be accurately guided by further research on the aspect of granite buried hill reservoir development rule prediction.
At present, the granite buried hill reservoir development law prediction research is mainly carried out in the following ways. (1) A reservoir engineering method. Such as material balance method, yield decreasing formula method, etc., the practical application has a great limitation because of the high requirement for the precision or data volume of the dynamic and static data of the oil reservoir. The method is characterized in that the theoretical basis is solid and reliable, the physical significance of the expression is definite, but the defects of more required parameters, complicated calculation method and difficult mastering exist. (2) A reservoir numerical simulation method. By establishing an actual geological model and correcting the numerical reservoir model by using historical production data, the production dynamics of the oil field is researched under the production condition, but the numerical simulation method has higher requirement on the oil reservoir data and large workload and is easily influenced by grid division. At present, a dual medium model is mostly adopted for researching the numerical simulation of the granite buried hill oil reservoir, the model highly simplifies the actual crack distribution characteristics, and the large error exists in the production rule of the granite buried hill oil reservoir. (3) A statistical rule method for historical data. The method needs a large amount of accurate actual production data of the granite buried hill reservoir and summarizes the historical production data statistics and development rules, but the method is suitable for the later stage of reservoir production and development and needs a large amount of support of accurate production data information. Therefore, a method for predicting the granite buried hill reservoir development rule chart with low cost and less historical production data is needed to be researched.
Disclosure of Invention
Aiming at the problems in the prior art, the technical problems to be solved by the invention are as follows: the method for drawing the granite buried hill reservoir development rule chart has the advantages of wide application range, accurate productivity data and suitability for field use.
In order to solve the technical problems, the invention adopts the following technical scheme:
a drawing method of a granite buried hill reservoir development rule chart comprises the following steps:
s100: and detecting, coring a rock core, performing conventional imaging logging and seismic data interpretation processing on the M granite buried hill oil reservoirs to correspondingly obtain M groups of oil reservoir data, wherein each group of oil reservoir data comprises oil reservoir fracture data and oil reservoir basic data.
S200: and performing effective characteristic parameter analysis and statistics on fracture distribution characteristics of the oil reservoir according to oil reservoir fracture data in the ith group of oil reservoir data to obtain ith group of fracture data, wherein the ith group of fracture data comprises fracture parameters and fracture distribution characteristics, and i =1,2 and … M.
S300: and (5) generating an ith two-dimensional discrete fracture geometric model by using the matlab program by taking the ith group of fracture data as constraints.
S400: and (3) introducing the ith two-dimensional discrete fracture geometric model into a comsol program, inputting the oil reservoir basic data of the ith group of oil reservoir data into the comsol program for parameter supplement, and inputting a control equation of the reservoir matrix and the fracture area into the comsol program, so as to obtain the ith discrete fracture numerical simulation model.
The input parameters of the ith discrete fracture numerical simulation model comprise production pressure difference, time step and time step length, the output parameters of the discrete fracture numerical simulation model comprise outlet flow rate, and meanwhile, the convergence of the ith discrete fracture numerical simulation model is obtained through a comsol program.
S500: the method comprises the steps of setting the production pressure difference of an ith discrete fracture numerical simulation model according to actual production conditions, setting the time step and the time step length of the ith discrete fracture numerical simulation model according to the convergence of the ith discrete fracture numerical simulation model, and carrying out numerical simulation on an actual production process by using the ith discrete fracture numerical simulation model to obtain the outlet flow rate.
S600: and integrating the outlet flow rate to obtain the daily oil yield, integrating the daily oil yield with a time condition to obtain the accumulated oil yield, and recording the production pressure difference and the accumulated oil yield.
S700: and repeating the steps S500-S600 for multiple times, and setting different production pressure differences at each time, thereby obtaining data of accumulated oil production under multiple different production pressure differences.
S800: repeating the steps S200-S700 for multiple times, sequentially using the oil deposit data of other groups, thereby obtaining the data of the accumulated oil production of each group of oil deposit data under different production pressure differences, and counting the data of the accumulated oil production of each group of oil deposit data under different production pressure differences together.
S900: and establishing a granite buried hill reservoir development rule chart according to the data of the accumulated oil production of each group of reservoir data under different production pressure differences.
The granite buried hill reservoir development rule chart is designed according to the reservoir data of a plurality of granite buried hill reservoirs, so that the obtained chart has strong universality, and when a development scheme is formulated for a specific granite buried hill reservoir, the productivity can be predicted according to the reservoir data of the specific reservoir, so that a reasonable development scheme is formulated according to actual production conditions.
Preferably, the fracture distribution characteristics in S200 include length distribution characteristics of the fractures, and a power law length model is used for statistics of the length distribution characteristics of the fractures, and the formula is as shown in (2-1).
n(l,L)＝αL D l -a ,l∈[l min ,l max ] (2-1)
In the formula (2-1), n (L, L) is the number of cracks with the size of L, the interval is [ L, L + dl ] and dl is far smaller than L, L is a modeling area, D is a fractal dimension, a is a characteristic power law index of length distribution, alpha is a constant related to fracture density, and the only intrinsic characteristic length scales of the power law length model are the minimum fracture length and the maximum fracture length, namely lmin and lmax.
Preferably, the fracture parameters in S200 include fracture strength, and the formula used for analyzing the fracture strength is as shown in (2-2).
Gamma in the formula (2-2) is the fracture Strength, A L Is the modeled zone area, L' is the fracture length in the modeled zone, the L fracture length dimension, the dl fracture length variation dimension.
Through the combination of the two formulas, the calculation result of the fracture strength gamma is more accurate, the finally obtained accumulated oil production data is closer to the real condition, the finally obtained granite buried hill oil reservoir development rule chart is higher in accuracy, and the practicability of the invention is greatly improved.
Preferably, each set of the basic data of the reservoir in the step S100 includes permeability, porosity, crude oil density and dynamic viscosity of the reservoir.
The step S300 of using the ith group of fracture data as the constraint means using the parameter a and the parameter γ obtained from the ith group of fracture data as the constraint.
The step S400 of inputting the reservoir basic data of the ith group of oil reservoir data into the comsol program for parameter supplement refers to performing parameter supplement by using the reservoir permeability, porosity, crude oil density and kinematic viscosity in the reservoir basic data of the ith group of oil reservoir data input into the comsol program.
Preferably, the reservoir matrix control equation in S400 is Darcy' S seepage flow equation, such as (2-3) and (2-4).
In equations (2-3) and (2-4)
Is the matrix porosity, ρ is the fluid density,
is the velocity vector, Q is the source term, k is the matrix permeability, μ dynamic viscosity, p is the pressure, and t is the time.
Preferably, the control equation of the fracture area in S400 adopts fracture flow equations, such as (2-5) and (2-6).
D in equations (2-5) and (2-6) f The f index is the crack opening and the f index is the crack parameter.
The numerical simulation of the discrete fracture numerical simulation model is more accurate by using the Darcy seepage flow equation and the fracture flow equation, so that the finally obtained outlet flow velocity is more in line with the real situation.
Preferably, in the step S800, the step of counting the data of the accumulated oil production of each group of oil deposit data at different production pressure differentials together means to count the a value, the γ value, the values of a plurality of production pressure differentials and the accumulated oil production under the corresponding production pressure differentials in each group of oil deposit data, and obtain a data statistical table of numerical simulation results. By obtaining the numerical simulation result data statistical table, the drawing of the subsequent granite buried hill reservoir development rule chart is more convenient.
Preferably, the step S900 is detailed.
S910: and according to the numerical simulation result data statistical table obtained in the step S800, drawing a scatter diagram with a trend line by using a plurality of fracture strengths γ corresponding to one value a as abscissa and using the accumulated oil production corresponding to each production pressure difference as ordinate, wherein the scatter diagram has a plurality of trend lines with different production pressure differences.
S920: and repeating the step S910 for multiple times, using different a values each time to obtain a plurality of scatter diagrams, and sequentially arranging the scatter diagrams corresponding to the a values according to the a values to obtain a granite buried hill reservoir development rule chart. The granite buried hill reservoir development rule chart is drawn according to the a value, the gamma value, the production pressure difference and the accumulated oil production, so that the data expression of the granite buried hill reservoir development rule chart is simple and clear, and the granite buried hill reservoir development rule chart is convenient to use.
Preferably, the use method of the granite buried hill reservoir development rule chart comprises the following steps.
S1000: and determining that the type of the oil reservoir needing to be developed at present is a granite buried hill oil reservoir, and exploring the oil reservoir needing to be developed at present to obtain the data of the oil reservoir at present.
S1100: and according to the corresponding a value and the gamma value in the current oil reservoir data, the data of the accumulated oil production amount of the current oil reservoir under different production pressure differences can be obtained by referring to a granite buried hill oil reservoir development rule chart.
S1200: and formulating a reasonable oil reservoir development scheme according to the data of the accumulated oil production obtained in the step S1100. By collecting the oil reservoir data of the oil reservoir which needs to be developed currently, the oil production data of the current oil reservoir can be easily obtained by referring to the granite buried hill oil reservoir development rule chart, so that a reasonable exploitation scheme is conveniently formulated, and the practicability of the granite buried hill oil reservoir development rule chart is greatly improved.
Preferably, in the step S500, before the numerical simulation is performed on the actual production process by using the discrete fracture numerical simulation model, the discrete fracture numerical simulation model is subdivided and optimized by using an unstructured network, where the unstructured network is the prior art.
The unstructured network is used for subdividing and optimizing the discrete fracture numerical simulation model, so that outlet flow velocity data obtained by numerical simulation of the discrete fracture numerical simulation model are more accurate, the accuracy of granite buried hill reservoir development rule chart data is improved, and the accuracy and the practicability of the method are greatly improved.
Compared with the prior art, the invention has at least the following advantages:
1. the method has the advantages that different discrete fracture models are generated by adopting oil deposit data of a plurality of oil deposits, the oil deposit generation process is subjected to numerical simulation, and a development rule prediction plate is established, so that the development rule prediction plate has stronger applicability and can be suitable for capacity prediction of most granite buried hill oil deposits, the accuracy of the predicted value of the capacity is high, the oil deposit capacity can be quickly estimated by comparing the field oil deposit data with the development rule prediction plate, and the field use is more convenient and efficient.
2. And predicting the oil reservoir development capacity by combining the established development rule prediction plate with the on-site exploration oil reservoir fracture data, and guiding the development scheme to make, so that the development scheme is more reasonable.
3. The unstructured network is used for subdividing and optimizing the discrete fracture numerical simulation model, so that outlet flow velocity data obtained by numerical simulation of the discrete fracture numerical simulation model are more accurate, the accuracy of granite buried hill reservoir development rule chart data is improved, and the accuracy and the practicability of the method are greatly improved.
Drawings
FIG. 1 is a block flow diagram of the present invention.
FIG. 2 is a geometric model diagram of different fracture parameters in the example.
FIG. 3 is a diagram of an opening distribution of a numerical simulation model of a discrete fracture in an embodiment.
FIG. 4 is a schematic diagram of fracture flow of a numerical simulation model of discrete fractures in an embodiment.
FIG. 5 is a mesh subdivision diagram of a numerical simulation model of discrete fractures in an embodiment.
FIG. 6 is a graph of the effect of production differential pressure and fracture parameters in the examples.
The rule prediction graph is developed in the embodiment of fig. 7.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings.
Referring to fig. 1-7, one embodiment of the present invention is provided:
example 1: a drawing method of a granite buried hill reservoir development rule chart comprises the following steps:
s100: and detecting, coring a rock core, performing conventional imaging logging and seismic data interpretation processing on the M granite buried hill oil reservoirs to correspondingly obtain M groups of oil reservoir data, wherein each group of oil reservoir data comprises oil reservoir fracture data and oil reservoir basic data.
S200: and performing effective characteristic parameter analysis and statistics on the fractures of the oil reservoir according to the oil reservoir fracture data in the ith group of oil reservoir data to obtain the ith group of fracture data, wherein the ith group of fracture data comprises fracture parameters and fracture distribution characteristics, i =1,2 and … M.
S300: and generating an ith two-dimensional discrete fracture geometric model by using the matlab program by taking the ith group of fracture data as constraints.
S400: and (3) introducing the ith two-dimensional discrete fracture geometric model into a comsol program, inputting the oil reservoir basic data of the ith group of oil reservoir data into the comsol program for parameter supplement, and inputting a control equation of the reservoir matrix and the fracture area into the comsol program, so as to obtain the ith discrete fracture numerical simulation model.
The input parameters of the ith discrete fracture numerical simulation model comprise production pressure difference, time step and time step length, the output parameters of the discrete fracture numerical simulation model comprise outlet flow rate, and meanwhile, the convergence of the ith discrete fracture numerical simulation model is obtained through a comsol program. In specific implementation, the fracture parameters further include fracture length, fracture density, and fracture opening.
S500: the method comprises the steps of setting the production pressure difference of an ith discrete fracture numerical simulation model according to actual production conditions, setting the time step and the time step length of the ith discrete fracture numerical simulation model according to the convergence of the ith discrete fracture numerical simulation model, and carrying out numerical simulation on an actual production process by using the ith discrete fracture numerical simulation model to obtain the outlet flow rate.
S600: and integrating the outlet flow rate to obtain the daily oil yield, integrating the daily oil yield with a time condition to obtain the accumulated oil yield, and recording the production pressure difference and the accumulated oil yield.
S700: and repeating the steps S500-S600 for multiple times, and setting different production pressure differences at each time, thereby obtaining data of accumulated oil production under multiple different production pressure differences.
S800: repeating the steps S200-S700 for multiple times, sequentially using the oil deposit data of other groups, thereby obtaining the data of the accumulated oil production of each group of oil deposit data under different production pressure differences, and counting the data of the accumulated oil production of each group of oil deposit data under different production pressure differences together.
S900: and establishing a granite buried hill reservoir development rule chart according to the data of the accumulated oil production of each group of reservoir data under different production pressure differences.
Further, the crack distribution characteristics in the S200 include length distribution characteristics of cracks, and a power law length model is used when statistics is performed on the length distribution characteristics of the cracks, and the formula is shown as (3-1).
n(l,L)＝αL D l -a ,l∈[l min ,l max ] (3-1)
In the formula (3-1), n (L, L) is the number of cracks with the size of L, the interval is [ L, L + dl ] and dl is far smaller than L, L is the modeling area, D is the fractal dimension, a is the characteristic power law index of the length distribution, alpha is a constant related to the fracture density, and the only intrinsic characteristic length scales of the power law length model are the minimum fracture length and the maximum fracture length, lmin and lmax respectively. In particular implementations, D ranges from 1.5 to 2.0 for a natural fracture network.
Further, the fracture parameters in S200 include fracture strength, and the formula used for analyzing the fracture strength is as shown in (3-2).
Gamma in the formula (3-2) is the fracture Strength, A L Is the modeled zone area, L' is the fracture length in the modeled zone, the L fracture length dimension, the dl fracture length variation dimension. In specific implementation, the value of the parameter a is as follows: [1.3,3.5]The value of the parameter γ may be: 1.25, 2.5, 3.75 or 5, the values of parameter a and parameter γ being selected in combination according to the specific reservoir data.
Further, each set of reservoir basic data in the step S100 includes reservoir permeability, porosity, crude oil density, and kinematic viscosity.
The step S300 of using the ith group of fracture data as constraints means using the parameter a and the parameter γ obtained from the ith group of fracture data as constraints.
In the step S400, inputting the oil reservoir basic data of the ith group of oil reservoir data into the comsol program for parameter supplementation refers to inputting the oil reservoir permeability, porosity, crude oil density and dynamic viscosity in the oil reservoir basic data of the ith group of oil reservoir data into the comsol program for parameter supplementation.
Further, the control equation of the reservoir matrix in S400 adopts darcy seepage flow equations, such as (3-3) and (3-4).
In equations (3-3) and (3-4), φ is the porosity of the matrix, ρ is the fluid density, u is → Is the velocity vector, Q is the source term, k is the matrix permeability, μ dynamic viscosity, p is the pressure, and t is the time.
Further, the control equation of the fracture area in S400 adopts a fracture flow equation, and the equations are (3-5) and (3-6).
D in equations (3-5) and (3-6) f The f index is the crack opening and the f index is the crack parameter.
Further, in the step S800, the specific operation of counting the data of the accumulated oil production of each group of oil deposit data at different production pressure differentials together is to count the a value, the γ value, the values of a plurality of production pressure differentials and the accumulated oil production under the corresponding production pressure differentials in each group of oil deposit data to obtain a numerical simulation result data statistical table.
Further, the step S900 is detailed as follows.
S910: and according to the numerical simulation result data statistical table obtained in the step S800, drawing a scatter diagram with a trend line by using a plurality of fracture strengths γ corresponding to one value a as abscissa and using the accumulated oil production corresponding to each production pressure difference as ordinate, wherein the scatter diagram has a plurality of trend lines with different production pressure differences.
S920: and repeating the step S910 for multiple times, using different a values each time to obtain a plurality of scatter diagrams, and sequentially arranging the scatter diagrams corresponding to the a values according to the a values to obtain a granite buried hill reservoir development rule chart.
Further, the use method of the granite buried hill reservoir development rule chart comprises the following steps.
S1000: and determining that the type of the oil reservoir needing to be developed at present is a granite buried hill oil reservoir, and exploring the oil reservoir needing to be developed at present to obtain the data of the oil reservoir at present.
S1100: according to the corresponding a value and gamma value in the current oil deposit data, the data of the accumulated oil production quantity of the current oil deposit under different production pressure differences can be obtained by referring to a granite buried hill oil deposit development rule chart.
S1200: and formulating a reasonable oil reservoir development scheme according to the data of the accumulated oil production obtained in the step S1100.
Further, in the step S500, before the discrete fracture numerical simulation model is used to perform numerical simulation on the actual production process, the unstructured network is used to perform subdivision optimization on the discrete fracture numerical simulation model.
And (3) test analysis:
experimental analysis was performed according to example 1.
And taking 20 from M in the step S100 to obtain 20 groups of oil deposit data.
According to the combination of the parameter a and the parameter gamma in the reservoir data, 20 groups of discrete fracture geometric models are generated, as shown in fig. 2, the fracture opening distribution is shown in fig. 3, and the fracture area is shown in fig. 4.
And setting the production pressure difference in the step S500 to be 2Mpa, setting the time step to be 100 days, and setting the time step to be 0.1 day. The effect of the unstructured network subdivision optimization is shown in fig. 5.
Setting the production pressure difference in the step S700 as follows in sequence: 4MPa, 6MPa, 8MPa, 10MPa, 12MPa and 14MPa.
After step S800, the following data statistics table of numerical simulation results is obtained, as shown in table 1.
Table 1 data statistics table of numerical simulation results
In step S920, determining a production pressure difference, drawing a line graph with the parameter a as an abscissa and the accumulated oil production as an ordinate, where the line graph has a plurality of broken lines with different fracture strengths γ, as shown in fig. 6, it can be known from fig. 6 that the fracture strength γ affects the productivity of the oil reservoir under the condition determined by the production pressure difference, and the greater the production pressure difference, the greater the influence of the fracture strength on the productivity.
And (3) drawing a scatter diagram with a trend line by taking a plurality of fracture strengths gamma corresponding to one value a as a horizontal coordinate and taking the accumulated oil production corresponding to each production pressure difference as a vertical coordinate, wherein the scatter diagram has a plurality of trend lines with different production pressure differences, repeating the step S910 for a plurality of times, using different values a each time, thereby obtaining a plurality of scatter diagrams, and sequentially arranging the scatter diagrams corresponding to each value a according to the value a to obtain a granite buried hill oil reservoir development rule chart, such as the graph 7.
In step S1000, the current reservoir A, B and C three blocks are explored, and parameters (a, γ) are determined as (1.5,2.5), (3,5) and (1.5,1.25), respectively, according to the reservoir data.
In the step S1100, the production pressure difference is measured to be 6Mpa according to the field experience, and the oil reservoir development potentials of A, B and the C block can be known according to the granite buried hill oil reservoir development rule chart of figure 7 as follows: b is more than A and more than C.
In step S1200, a reasonable development plan is made according to the development potential of each block.
In the test, 20 groups of oil reservoir data are selected as a basis to draw a granite buried hill oil reservoir development rule chart, in actual use, oil reservoir data of A, B, C three blocks in an oil reservoir are selected to obtain three groups (a, gamma), then, production pressure difference is taken according to production conditions, the granite buried hill oil reservoir development rule chart is referred to obtain accumulated oil production data of the three blocks, the development potentials of the three blocks can be easily obtained, and a reasonable development scheme can be formulated according to the development potentials.
The working principle of the drawing method of the granite buried hill reservoir development rule chart is as follows:
the method comprises the steps of firstly obtaining oil deposit data of a plurality of granite oil deposits, establishing a discrete fracture numerical simulation model according to the oil deposit data of each oil deposit, carrying out numerical simulation under different production pressure differences by using the discrete fracture numerical simulation model, obtaining data of accumulated oil production of different production pressure differences under different oil deposit data, drawing a development rule chart according to the data of the accumulated oil production, and predicting the productivity of the current oil deposit according to actual oil deposit data and the development rule chart before developing the specific oil deposit, so that a proper development scheme is formulated. In the invention, the oil deposit data of a plurality of oil deposits are adopted to generate different discrete fracture models for numerical simulation, so that the development rule prediction chart has stronger applicability, can be suitable for the productivity prediction of most granite buried hill oil deposits, and has convenient use and high accuracy.
Finally, although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that various changes and modifications may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (10)
1. A drawing method of a granite buried hill reservoir development rule chart is characterized by comprising the following steps: the drawing method of the chart comprises the following steps;
s100: detecting M granite buried hill oil reservoirs, coring a rock core, performing conventional imaging logging and seismic data interpretation treatment, and correspondingly obtaining M groups of oil reservoir data, wherein each group of oil reservoir data comprises oil reservoir fracture data and oil reservoir basic data;
s200: performing effective characteristic parameter analysis and statistics on fracture distribution characteristics of the oil reservoir according to oil reservoir fracture data in the ith group of oil reservoir data to obtain ith group of fracture data, wherein the ith group of fracture data comprises fracture parameters and fracture distribution characteristics, and i =1,2 and … M;
s300: generating an ith two-dimensional discrete fracture geometric model by using a matlab program by taking the ith group of fracture data as constraints;
s400: importing the ith two-dimensional discrete fracture geometric model into a comsol program, inputting oil reservoir basic data of the ith group of oil reservoir data into the comsol program for parameter supplement, and inputting a control equation of a reservoir matrix and a fracture area into the comsol program to obtain an ith discrete fracture numerical simulation model;
the input parameters of the ith discrete fracture numerical simulation model comprise production pressure difference, time step and time step length, the output parameters of the discrete fracture numerical simulation model comprise outlet flow rate, and meanwhile, the convergence of the ith discrete fracture numerical simulation model is obtained by a comsol program;
s500: setting the production pressure difference of an ith discrete fracture numerical simulation model according to actual production conditions, setting the time step and the time step length of the ith discrete fracture numerical simulation model according to the convergence of the ith discrete fracture numerical simulation model, and carrying out numerical simulation on the actual production process by using the ith discrete fracture numerical simulation model to obtain the outlet flow rate;
s600: integrating the outlet flow rate to obtain daily oil production, integrating the daily oil production with time conditions to obtain accumulated oil production, and recording the production pressure difference and the accumulated oil production;
s700: repeating the steps S500-S600 for multiple times, and setting production pressure differences with different values each time so as to obtain data of accumulated oil production under multiple different production pressure differences;
s800: repeating the steps S200-S700 for multiple times, sequentially using the oil deposit data of other groups to obtain the data of the accumulated oil production of each group of oil deposit data under different production pressure differences, and counting the data of the accumulated oil production of each group of oil deposit data under different production pressure differences together;
s900: and establishing a granite buried hill reservoir development rule chart according to the data of the accumulated oil production of each group of reservoir data under different production pressure differences.
2. The method for drawing the granite buried hill reservoir development law plate as claimed in claim 1, wherein: the crack distribution characteristics in the S200 comprise the length distribution characteristics of the cracks, and a power law length model is used when the length distribution characteristics of the cracks are counted, wherein the power law length model is shown as a formula (1-1);
n(l,L)＝αL D l -a ,l∈[l min ,l max ] (1-1)
in the formula (1-1), n (L, L) is the number of cracks with the size of L, the interval is [ L, L + dl ] and dl is far less than L, L is a modeling area, D is a fractal dimension, a is a characteristic power law index of length distribution, alpha is a constant related to fracture density, and the only intrinsic characteristic length scales of the power law length model are the minimum fracture length and the maximum fracture length, namely lmin and lmax.
3. The method for drawing the granite buried hill reservoir development law plate as claimed in claim 2, wherein: the fracture parameters in the S200 comprise fracture strength, and a formula used for analyzing the fracture strength is shown as (1-2);
gamma in the formula (1-2) is the fracture Strength, A L Is the modeled zone area, L' is the fracture length in the modeled zone, the L fracture length dimension, the dl fracture length variation dimension.
4. The method for drawing the granite buried hill reservoir development law plate as claimed in claim 3, wherein: in the step S100, each group of oil deposit basic data comprises oil deposit permeability, porosity, crude oil density and dynamic viscosity;
in the step S300, the i-th group of fracture data is used as a constraint, which means that the parameter a and the parameter γ obtained from the i-th group of fracture data are used as constraints;
in the step S400, inputting the oil reservoir basic data of the ith group of oil reservoir data into the comsol program for parameter supplementation refers to inputting the oil reservoir permeability, porosity, crude oil density and dynamic viscosity in the oil reservoir basic data of the ith group of oil reservoir data into the comsol program for parameter supplementation.
5. The method for drawing the granite buried hill reservoir development law plate as claimed in claim 4, wherein: the reservoir matrix control equation in the S400 adopts Darcy seepage flow equations, and the equations are (1-3) and (1-4);
in equations (1-3) and (1-4)
Is the matrix porosity, ρ is the fluid density,
is the velocity vector, Q is the source term, k is the matrix permeability, μ dynamic viscosity, p is the pressure, and t is the time.
6. The method for drawing the granite buried hill reservoir development law plate as claimed in claim 5, wherein: the control equation of the crack area in the S400 adopts a crack flow equation, and the equations are (1-5) and (1-6);
d in equations (1-5) and (1-6) f The f index is the crack opening and the f index is the crack parameter.
7. The method for drawing the granite buried hill reservoir development law plate as claimed in claim 6, wherein: in the step S800, the data statistics of the accumulated oil production of each group of oil deposit data under different production pressure differentials means that a value, a plurality of production pressure differentials, and the accumulated oil production under corresponding production pressure differentials in each group of oil deposit data are counted to obtain a data statistics table of numerical simulation results.
8. The method for drawing the granite buried hill reservoir development law plate as claimed in claim 7, wherein: the step S900 is detailed as follows;
s910: according to the numerical simulation result data statistical table obtained in the step S800, drawing a scatter diagram with trend lines by taking a plurality of crack strengths gamma corresponding to one value a as horizontal coordinates and taking the accumulated oil production corresponding to each production pressure difference as vertical coordinates, wherein the scatter diagram comprises a plurality of trend lines with different production pressure differences;
s920: and repeating the step S910 for multiple times, using different a values each time to obtain a plurality of scatter diagrams, and sequentially arranging the scatter diagrams corresponding to the a values according to the a values to obtain a granite buried hill reservoir development rule chart.
9. The method for drawing the granite buried hill reservoir development law plate as claimed in claim 8, wherein: the use method of the granite buried hill reservoir development rule chart comprises the following steps;
s1000: determining that the type of the oil reservoir needing to be developed at present is a granite buried hill oil reservoir, and exploring the oil reservoir needing to be developed at present to obtain current oil reservoir data;
s1100: according to the corresponding a value and gamma value in the current oil deposit data, referring to a granite buried hill oil deposit development rule plate to obtain the data of the accumulated oil production quantity of the current oil deposit under different production pressure differences;
s1200: and formulating a reasonable oil reservoir development scheme according to the data of the accumulated oil production obtained in the step S1100.
10. The method for drawing the granite buried hill reservoir development law plate as claimed in claim 1, wherein: and in the step S500, before the numerical simulation is carried out on the actual production process by using the discrete fracture numerical simulation model, the unstructured network is used for carrying out subdivision optimization on the discrete fracture numerical simulation model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010684890.8A CN111810137B (en) | 2020-07-16 | 2020-07-16 | Drawing method of granite buried hill reservoir development rule chart |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010684890.8A CN111810137B (en) | 2020-07-16 | 2020-07-16 | Drawing method of granite buried hill reservoir development rule chart |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111810137A CN111810137A (en) | 2020-10-23 |
| CN111810137B true CN111810137B (en) | 2023-03-03 |
Family
ID=72865370
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010684890.8A Active CN111810137B (en) | 2020-07-16 | 2020-07-16 | Drawing method of granite buried hill reservoir development rule chart |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111810137B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104730596A (en) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | Discrete fracture modeling method based on multiscale factor restraint |
| CN106777628A (en) * | 2016-06-29 | 2017-05-31 | 中国石油大学(华东) | Consider the oil reservoir injectivity and productivity plate method for drafting of non-Darcy flow |
| CN109116428A (en) * | 2018-07-02 | 2019-01-01 | 中国石油天然气股份有限公司 | A kind of fractured-cavernous carbonate reservoir uncertainties model method and device thereof |
| CN109558631A (en) * | 2018-10-24 | 2019-04-02 | 北京大学 | A kind of densification oil-gas reservoir multi-scale facture medium automatic history matching method |
| CN109558614A (en) * | 2017-09-27 | 2019-04-02 | 中国石油化工股份有限公司 | The analogy method and system that gas flows in shale gas reservoir multi-scale facture |
| CN109829217A (en) * | 2019-01-21 | 2019-05-31 | 中国石油大学(北京) | Pressure break Fractured Reservoir productivity simulation method and device |
| CN111222271A (en) * | 2020-01-03 | 2020-06-02 | 中国石油大学(华东) | Numerical reservoir fracture simulation method and system based on matrix-fracture unsteady state channeling |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7925482B2 (en) * | 2006-10-13 | 2011-04-12 | Object Reservoir, Inc. | Method and system for modeling and predicting hydraulic fracture performance in hydrocarbon reservoirs |
| CN102110183B (en) * | 2010-12-30 | 2011-09-28 | 中国石油化工股份有限公司胜利油田分公司地质科学研究院 | Numerical simulation method for reflecting fluid channeling of fluid along great rifts of reservoir |
| FR3041026B1 (en) * | 2015-09-15 | 2017-10-20 | Ifp Energies Now | METHOD FOR CHARACTERIZING THE NETWORK OF FRACTURES OF A FRACTURE SLOT AND METHOD FOR OPERATING IT |
-
2020
- 2020-07-16 CN CN202010684890.8A patent/CN111810137B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104730596A (en) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | Discrete fracture modeling method based on multiscale factor restraint |
| CN106777628A (en) * | 2016-06-29 | 2017-05-31 | 中国石油大学(华东) | Consider the oil reservoir injectivity and productivity plate method for drafting of non-Darcy flow |
| CN109558614A (en) * | 2017-09-27 | 2019-04-02 | 中国石油化工股份有限公司 | The analogy method and system that gas flows in shale gas reservoir multi-scale facture |
| CN109116428A (en) * | 2018-07-02 | 2019-01-01 | 中国石油天然气股份有限公司 | A kind of fractured-cavernous carbonate reservoir uncertainties model method and device thereof |
| CN109558631A (en) * | 2018-10-24 | 2019-04-02 | 北京大学 | A kind of densification oil-gas reservoir multi-scale facture medium automatic history matching method |
| CN109829217A (en) * | 2019-01-21 | 2019-05-31 | 中国石油大学(北京) | Pressure break Fractured Reservoir productivity simulation method and device |
| CN111222271A (en) * | 2020-01-03 | 2020-06-02 | 中国石油大学(华东) | Numerical reservoir fracture simulation method and system based on matrix-fracture unsteady state channeling |
Non-Patent Citations (7)
| Title |
|---|
| Joint influence of in-situ stress and fracture network geometry on heat transfer in fractured geothermal reservoirs;Zhixue Sun等;《International Journal of Heat and Mass Transfer》;20191231(第149期);全文 * |
| M凹陷太古界潜山储层评价研究;李华;《2014年环渤海、泛珠三角区域地球物理学术研讨会论文集》;20141231;全文 * |
| Numerical study of dynamic fracture aperture during production of pressure-sensitive reservoirs;Shaohua Gu等;《International Journal of Rock Mechanics and Mining Sciences》;20141231(第70期);全文 * |
| 多因素约束下的致密砂岩气藏离散裂缝特征及地质模型研究;冯建伟等;《中国石油大学学报（自然科学版）》;20161231(第01期);全文 * |
| 花岗岩潜山双重孔隙介质油藏地质建模方法;潘晓庆等;《西南石油大学学报(自然科学版)》;20191231(第04期);全文 * |
| 蚂蚁追踪技术在潜山油藏裂缝预测中的应用;张淑娟等;《断块油气田》;20111231(第001期);全文 * |
| 页岩气藏近井精细数值模拟方法研究;梁梦吟;《中国优秀硕博士学位论文全文数据库》;20200715(第07期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111810137A (en) | 2020-10-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109441422B (en) | Shale gas well spacing optimization mining method | |
| CN111441758B (en) | Shale oil gas dessert area prediction method and device | |
| CN102041995B (en) | System for monitoring complicated oil deposit flooding conditions | |
| CN107977480B (en) | Shale gas reservoir gas production performance rapid evaluation method | |
| CN109236273B (en) | Dynamic data processing method for oil field development and production | |
| CN104899411B (en) | A kind of reservoir productivity prediction model method for building up and system | |
| CN110056346B (en) | Oil reservoir three-dimensional original water saturation simulation method based on trend change function | |
| CN108150160B (en) | Method for solving under-compaction and over-pressure in stratum | |
| CN109063403B (en) | Optimal design method for slickwater fracturing | |
| CN107728227A (en) | A kind of well pattern maturation zone Quick buried channel sand body method | |
| CN105089662B (en) | The bearing calibration of carbonate formation well logging composite columnar section and system | |
| CN112049629B (en) | Fracture-cavity type oil reservoir recovery ratio prediction method based on A-type water drive characteristic curve | |
| CN106481315B (en) | Land sandstone oil reservoir individual well recoverable reserves quickly determines model and method for building up | |
| CN111810137B (en) | Drawing method of granite buried hill reservoir development rule chart | |
| CN114357750A (en) | Goaf water filling state evaluation method | |
| CN116306118A (en) | Crack modeling method | |
| CN112989257B (en) | Gas production amount measuring method for sea shale oil-gas reservoir | |
| CN115705452A (en) | Novel recovery ratio prediction method for middle and later stages of integrated sandstone reservoir development | |
| CN107369101A (en) | The movable evaluation method of target fast, economical is commented in rolling | |
| CN105888656A (en) | Liquid-measurement method for quantitatively evaluating penetration rate of natural-micro-fracture-development tight reservoir under pressure applying condition | |
| CN107942400B (en) | A kind of method and device for predicting drilling risk | |
| CN110991084A (en) | Reservoir permeability calculation method based on streamline numerical well testing | |
| CN112523748A (en) | Multi-stage and multi-dimensional evaluation method and device for tight oil volume fracturing effect | |
| CN108875109B (en) | Method and system for predicting abnormal formation pressure | |
| CN116563484B (en) | Compact sandstone reservoir property characterization method based on configuration theoretical analysis |
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 |