CN110083852B - Method and device for determining bottom hole pressure and storage medium - Google Patents
Method and device for determining bottom hole pressure and storage medium Download PDFInfo
- Publication number
- CN110083852B CN110083852B CN201810158974.0A CN201810158974A CN110083852B CN 110083852 B CN110083852 B CN 110083852B CN 201810158974 A CN201810158974 A CN 201810158974A CN 110083852 B CN110083852 B CN 110083852B
- Authority
- CN
- China
- Prior art keywords
- oil well
- bottom hole
- source
- determining
- unstructured
- 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
Images
Classifications
-
- 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
Abstract
The invention discloses a method and a device for determining bottom hole pressure and a storage medium, and belongs to the technical field of numerical reservoir simulation. The method comprises the steps of determining the flow conductivity of each non-structural unit grid in at least one non-structural unit grid where the control area at the source and sink node of the oil well is located, and further respectively establishing a first bottom hole pressure calculation model of the source and sink node when the oil well is in a quasi-steady state and a second bottom hole pressure calculation model of the source and sink node when the oil well is in an unsteady state based on the flow conductivity of each non-structural unit grid. And then, determining the bottom hole pressure of the oil well based on the first bottom hole pressure calculation model in the quasi-steady state and the second bottom hole pressure calculation model in the non-steady state, avoiding determining the pressure of the non-structural unit grid with a certain area as the bottom hole pressure of the oil well, improving the accuracy of the bottom hole pressure of the oil well, and providing more reliable data support for the design of a subsequent development scheme of the oil well.
Description
Technical Field
The invention relates to the technical field of numerical reservoir simulation, in particular to a method and a device for determining bottom hole pressure and a storage medium.
Background
The numerical reservoir simulation technology is an important means for quantitatively describing the multiphase flow rule of a heterogeneous reservoir corresponding to an oil well at present, and particularly has irreplaceable effects on the aspects of development scheme design, deployment optimization and adjustment, recovery ratio prediction and the like of the oil well. Therefore, during the process of exploiting the oil reservoir through the oil well, the bottom hole pressure of the oil well can be determined through an oil reservoir numerical simulation technology, and then the development scheme of the oil well is designed, the deployment of the oil well is adjusted and optimized, or the recovery ratio of the oil well is predicted and the like through the bottom hole pressure of the oil well.
In the related technology, with the continuous development of computer hardware and the development of unconventional oil reservoirs, a non-structural grid is adopted to carry out fine description on complex flow boundaries at the bottom of an oil well, and numerical solution is carried out by combining numerical discrete methods such as finite elements, finite volumes and the like, so that the method becomes a main research direction in the technical field of numerical simulation of oil reservoirs. Specifically, when a complex flow boundary at the bottom of an oil well is described in detail through an unstructured grid, a plurality of unstructured unit grids are obtained by dispersing a complex flow region at the bottom of the oil well through a dispersion method, then the pressure of the unstructured unit grid at the boundary of a shaft is directly calculated through combining methods such as finite elements and finite volumes, and the pressure of the unstructured unit grid at the boundary of the shaft is determined as the pressure of a source-sink node of the oil well, namely the bottom pressure of the oil well. The source-sink node of the oil well refers to a point which is similar to a flowing area at the bottom of the oil well relative to a bottom hole area corresponding to a shaft of the oil well in a discrete process.
However, since the unstructured unit mesh located at the wellbore boundary of the oil well has a certain area, and the source-sink node of the oil well can be approximated to a point with respect to the complex flow region at the bottom of the well, the pressure of the unstructured unit mesh located at the wellbore boundary, which is directly calculated by combining finite elements, finite volumes, and the like, is not equal to the pressure of the source-sink node of the oil well, and thus an error may be caused when designing a development scheme of the oil well, adjusting and optimizing the deployment of the oil well, or predicting the recovery ratio of the oil well, etc., based on the pressure of the unstructured unit mesh located at the wellbore boundary.
Disclosure of Invention
In order to avoid the problem that the pressure of an unstructured unit grid with a certain area is determined as the bottom pressure, so that the bottom pressure of an oil well is deviated, the invention provides a bottom pressure determining method, a device and a storage medium. The technical scheme is as follows:
in a first aspect, a method of determining a bottom hole pressure is provided, the method comprising:
determining at least one unstructured cell grid in which a control area at a source-sink node of an oil well is located;
determining a flow conductivity of each of the at least one grid of unstructured cells;
establishing a first bottom hole pressure calculation model of the source and sink nodes when the oil well is in a quasi-steady state and establishing a second bottom hole pressure calculation model of the source and sink nodes when the oil well is in an unsteady state based on the flow conductivity of each unstructured unit grid;
determining a bottom hole pressure of the well based on the first bottom hole pressure computational model and the second bottom hole pressure computational model.
Optionally, the determining at least one unstructured cell grid in which the control area at the source-sink node of the oil well is located comprises:
determining a seepage zone at the bottom of the oil well;
dispersing the seepage area to obtain a plurality of unstructured unit grids included at the bottom of the oil well;
determining at least one unstructured cell grid from the plurality of unstructured cell grids in which a control area at a source-sink node of the well resides.
Optionally, the establishing a first bottom hole pressure calculation model at the well point source sink when the oil well is in a quasi-steady state based on the flow conductivity of each unstructured unit grid comprises:
acquiring the well bore radius and the bottom hole yield of the oil well, and the viscosity of fluid in the oil well, the permeability of rock and the thickness of a reservoir;
the first downhole pressure calculation model is built based on the flow conductivity of each grid of unstructured cells, the wellbore radius and downhole production of the well, and the viscosity of the fluid in the well, permeability of the rock, and thickness of the reservoir.
Optionally, the establishing a second downhole pressure calculation model of the point-source junction when the well is in an unsteady state comprises:
obtaining theoretical pressure of the source and sink node at the current moment, bottom hole yield of the oil well, viscosity of fluid in the oil well, permeability of rock, thickness of a reservoir and a reservoir pressure guiding coefficient;
and establishing the second bottom-hole pressure calculation model based on the theoretical pressure of the source-sink node at the current moment, the bottom-hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rock, the thickness of a reservoir and a reservoir pressure-guiding coefficient.
Optionally, the determining the bottom hole pressure of the well based on the first bottom hole pressure calculation model and the second bottom hole pressure calculation model comprises:
determining a first equivalent hole diameter and a second equivalent hole diameter of the oil well at the current moment;
when the first equivalent hole diameter is smaller than or equal to the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to the first bottom hole pressure calculation model based on the first equivalent hole diameter;
and when the first equivalent hole diameter is larger than the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to the second bottom hole pressure calculation model based on the second equivalent hole diameter.
In a second aspect, there is provided an apparatus for determining bottom hole pressure, the apparatus comprising:
the first determination module is used for determining at least one unstructured unit grid where the control area at the source and sink nodes of the oil well is located;
a second determining module for determining a traffic conductivity of each of the at least one grid of unstructured cells;
the establishing module is used for establishing a first bottom hole pressure calculation model of the source sink node when the oil well is in a quasi-steady state and establishing a second bottom hole pressure calculation model of the source sink node when the oil well is in an unsteady state based on the flow conductivity of each unstructured unit grid;
a third determination module to determine a bottom hole pressure of the well based on the first bottom hole pressure computational model and the second bottom hole pressure computational model.
Optionally, the first determining module includes:
a first determination unit for determining a seepage area at the bottom of the oil well;
the discrete unit is used for dispersing the seepage area to obtain a plurality of non-structural unit grids included at the bottom of the oil well;
and the second determination unit is used for determining at least one unstructured unit grid in which the control area at the source and sink nodes of the oil well is located from the plurality of unstructured unit grids.
Optionally, the establishing module includes:
the first acquisition unit is used for acquiring the shaft radius and the bottom-hole yield of the oil well, and the viscosity of fluid in the oil well, the permeability of rock and the thickness of a reservoir;
a first establishing unit for establishing the first bottom hole pressure calculation model based on the flow conductivity of each unstructured cell grid, the wellbore radius and the bottom hole production of the oil well, and the viscosity of the fluid in the oil well, the permeability of the rock and the thickness of the reservoir.
Optionally, the establishing module includes:
the second acquisition unit is used for acquiring the theoretical pressure of the source-sink node at the current moment, the bottom-hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rocks, the thickness of a reservoir and a reservoir pressure-guiding coefficient;
and the second establishing unit is used for establishing the second bottom hole pressure calculation model based on the theoretical pressure of the source-sink node at the current moment, the bottom hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rock, the thickness of a reservoir and the reservoir pressure guiding coefficient.
Optionally, the determining the bottom hole pressure of the well based on the first bottom hole pressure calculation model and the second bottom hole pressure calculation model comprises:
determining a first equivalent hole diameter and a second equivalent hole diameter of the oil well at the current moment;
when the first equivalent hole diameter is smaller than or equal to the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to the first bottom hole pressure calculation model based on the first equivalent hole diameter;
and when the first equivalent hole diameter is larger than the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to the second bottom hole pressure calculation model based on the second equivalent hole diameter.
In a third aspect, a computer-readable storage medium is provided, in which a computer program is stored, which, when executed by a processor, implements any of the methods provided in the first aspect above.
The technical scheme provided by the invention has the beneficial effects that: the method comprises the steps of determining the flow conductivity of each unstructured unit grid in at least one unstructured unit grid where the control area at the source-sink node of the oil well is located, and further respectively establishing a first bottom-hole pressure calculation model of the source-sink node when the oil well is in a quasi-steady state and a second bottom-hole pressure calculation model of the source-sink node when the oil well is in an unsteady state based on the flow conductivity of each unstructured unit grid. And then, determining the bottom hole pressure of the oil well based on the first bottom hole pressure calculation model in the quasi-steady state and the second bottom hole pressure calculation model in the non-steady state, avoiding determining the pressure of the non-structural unit grid with a certain area as the bottom hole pressure of the oil well, improving the accuracy of the bottom hole pressure of the oil well, and providing more reliable data support for the design of a subsequent development scheme of the oil well.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a schematic flow diagram of a first method of determining bottom hole pressure provided by an embodiment of the present invention;
FIG. 2A is a schematic flow chart of a second method of determining bottom hole pressure provided by an embodiment of the present invention;
FIG. 2B is a schematic structural diagram of a control area at a source/sink node and at least one grid of unstructured cells according to an embodiment of the present invention;
FIG. 2C is a schematic diagram of a comparison of a bottom hole pressure determined based on a bottom hole pressure calculation model and a bottom hole pressure determined based on an analytical model of a seepage differential equation, according to an embodiment of the present invention;
FIG. 3A is a schematic diagram of a first apparatus for determining bottom hole pressure according to an embodiment of the present invention;
FIG. 3B is a schematic diagram of a second apparatus for determining bottom hole pressure provided by an embodiment of the present invention;
FIG. 3C is a schematic diagram of a third apparatus for determining bottom hole pressure provided by an embodiment of the present invention;
FIG. 3D is a schematic diagram of a fourth apparatus for determining bottom hole pressure provided by an embodiment of the present invention;
FIG. 3E is a schematic diagram of a fifth apparatus for determining bottom hole pressure according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of a terminal according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
Fig. 1 is a schematic flow chart of a method for determining a bottom hole pressure according to an embodiment of the present invention. Referring to fig. 1, the method includes the following steps.
Step 101: at least one grid of unstructured cells is determined in which the control area at the source-sink node of the well is located.
Step 102: a flow conductivity for each of the at least one grid of unstructured cells is determined.
Step 103: and establishing a first bottom hole pressure calculation model of the source and sink node when the oil well is in a quasi-steady state and establishing a second bottom hole pressure calculation model of the source and sink node when the oil well is in an unsteady state based on the flow conductivity of each unstructured unit grid.
Step 104: a bottom hole pressure of the well is determined based on the first bottom hole pressure calculation model and the second bottom hole pressure calculation model.
In the embodiment of the invention, the flow conductivity of each unstructured unit grid in at least one unstructured unit grid where the control area at the source-sink node of the oil well is located is determined, and then a first bottom-hole pressure calculation model of the source-sink node when the oil well is in a quasi-steady state and a second bottom-hole pressure calculation model of the source-sink node when the oil well is in an unsteady state are respectively established based on the flow conductivity of each unstructured unit grid. And then, determining the bottom hole pressure of the oil well based on the first bottom hole pressure calculation model in the quasi-steady state and the second bottom hole pressure calculation model in the non-steady state, avoiding determining the pressure of the non-structural unit grid with a certain area as the bottom hole pressure of the oil well, improving the accuracy of the bottom hole pressure of the oil well, and providing more reliable data support for the design of a subsequent development scheme of the oil well.
Optionally, determining at least one grid of unstructured cells in which the control area at the source-sink node of the well resides comprises:
determining a seepage area at the bottom of the oil well;
dispersing the seepage area to obtain a plurality of unstructured unit grids included at the bottom of the oil well;
at least one grid of unstructured cells in which the control area at the source-sink node of the well is located is determined from the plurality of grids of unstructured cells.
Optionally, establishing a first bottom hole pressure calculation model at the well point source sink when the well is in a quasi-steady state based on the flow conductivity of each unstructured cell grid, comprising:
obtaining the well bore radius and the bottom hole yield of the oil well, and the viscosity of fluid in the oil well, the permeability of rock and the thickness of a reservoir stratum;
a first downhole pressure calculation model is established based on the flow conductivity of each grid of unstructured cells, the wellbore radius and the downhole production of the well, and the viscosity of the fluid in the well, the permeability of the rock, and the thickness of the reservoir.
Optionally, establishing a second downhole pressure calculation model of the well point source sink when the well is in an unsteady state comprises:
obtaining theoretical pressure of the source and sink node at the current moment, bottom hole yield of the oil well, viscosity of fluid in the oil well, permeability of rock, thickness of a reservoir and a reservoir pressure guiding coefficient;
and establishing a second bottom hole pressure calculation model based on the theoretical pressure of the source-sink node at the current moment, the bottom hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rock, the thickness of a reservoir and the reservoir pressure guiding coefficient.
Optionally, determining a bottom hole pressure of the well based on the first bottom hole pressure calculation model and the second bottom hole pressure calculation model comprises:
determining a first equivalent hole diameter and a second equivalent hole diameter of the oil well at the current moment;
when the first equivalent hole diameter is smaller than or equal to the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to a first bottom hole pressure calculation model based on the first equivalent hole diameter;
and when the first equivalent hole diameter is larger than the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to the second bottom hole pressure calculation model based on the second equivalent hole diameter.
All the above optional technical solutions can be combined arbitrarily to form an optional embodiment of the present invention, which is not described in detail herein.
Fig. 2A is a schematic flow chart of a method for determining a bottom hole pressure according to an embodiment of the present invention. Referring to fig. 2A, the method includes the following steps.
Step 201: at least one grid of unstructured cells is determined in which the control area at the source-sink node of the well is located.
The control area refers to a control area corresponding to a source and sink node of the oil well.
Since the control area at the source sink node is part of the seepage zone at the bottom of the well, the seepage zone at the bottom of the well can be determined first. And then dispersing the seepage area to obtain a plurality of unstructured unit grids included at the bottom of the oil well, and then determining at least one unstructured unit grid where the control area at the source and sink nodes of the oil well is located from the plurality of unstructured unit grids.
Specifically, when determining the seepage area at the bottom of the oil well, the coordinate point where the source/sink node of the oil well is located and the coordinate point where each boundary point of the plurality of boundary points of the seepage area is located may be obtained in advance, and the seepage area at the bottom of the oil well may be drawn based on the coordinate point of the source/sink node and the coordinate point of each boundary point. After the seepage area at the bottom of the oil well is determined, the seepage area is a continuous area, so that in order to reduce the calculation amount and improve the operation efficiency of determining the control area at the source and sink node from the seepage area, the seepage area can be discretized to obtain a plurality of unstructured unit grids included in the seepage area. And then selecting at least one non-structural unit grid where a source and sink node of the oil well is located from the plurality of non-structural unit grids, and determining the selected at least one non-structural unit grid as the at least one non-structural unit grid where the control area at the source and sink node is located.
Wherein each unstructured cell mesh may be in the shape of an acute triangle, i.e., each unstructured cell mesh has three nodes. When the seepage region is discretized, a plurality of non-structural unit meshes included in the seepage region may be determined based on open-source non-structural mesh generation software, and of course, the seepage region may also be discretized into a plurality of non-structural unit meshes by other methods, which is not limited in the embodiment of the present invention.
For example, after the seepage area at the bottom of the well is discretized, the at least one unstructured-cell grid in which the control area at the source-sink node of the well is located includes 6 unstructured-cell grids, and the 6 unstructured-cell grids are shown in fig. 2B.
Step 202: a flow conductivity for each of the at least one grid of unstructured cells is determined.
Since the flow conductivity of each of the at least one unstructured-unit mesh is the same, a detailed description will be given below by taking any one of the at least one unstructured-unit mesh as an example. Since the flow conductivity of the percolation region at the bottom of the well can be determined by the following formula (1), the flow conductivity of the unstructured cell grid determined based on the flow conductivity of the percolation region can be determined by the following formula (2).
Wherein, in the above formula (1), FSRefers to the flow conductivity of the zone, S refers to the area of the zone, k refers to the permeability of the rock in the well, μ refers to the viscosity of the fluid in the well, B refers to the permeability of the rock in the wellgRefers to the volume factor of the fluid in the well, p refers to the fluid pressure in the zone of seepage, and Ω refers to the zone extent of the seepage zone. In the above formula (2), FeMeans the flow conductivity, s, of the control area at the source-sink node in the unstructured cell grid eeRefers to the area of the control area at the source-sink node in the unstructured cell grid e, peRefers to the fluid pressure, omega, within the grid of unstructured cells eeReferring to the area range of the unstructured cell grid e, the physical meaning of the remaining parameters can be taken into the meaning of the corresponding parameters in the above formula (1).
Then, the flow conductivity of the unstructured cell grid is reduced by a gaussian-green formula, a linear interpolation function shown in the following formula (3) is introduced, and the linear interpolation function is substituted into the flow conductivity after the reduction processing and solved, so that the flow conductivity of the unstructured cell grid indicated by the following formula (4) is obtained.
Wherein, in the above formula (3), NlIs a linear interpolation function, A is the area of the unstructured unit grid e, x is the abscissa of any point in the unstructured unit grid e, y is the ordinate of the any point, and x isi、xj、xkRespectively refers to the abscissa, y, of the node i, j, k on the unstructured cell grid ei、yj、ykRespectively, the ordinate of the node i, j, k on the unstructured cell grid. In the above formula (4), NiIs the interpolation term, p, corresponding to the node i in the unstructured unit grid eiRefers to the pressure, p, of node i in the grid e of unstructured cellslRefers to the pressure, T, of a node j or k in the grid e of unstructured cellsilThe flow conductivity coefficient between a node j or k and a node i in the unstructured unit grid e is referred, and the physical meanings of the remaining parameters can be the meanings of the corresponding parameters in the formula (1) or the formula (3).
Wherein, the volume coefficient of the oil product is 1, namely Bg1. Therefore, the above formula (4) can be simplified as the following formula (5):
during the process of exploiting the oil deposit based on the oil well, the control area at the source and sink nodes of the oil well can be in a quasi-steady state or an unsteady state, and the determination method of the pressure of the source and sink nodes at the quasi-steady state and the unsteady state is different. The pressure determination method for the source/sink node when the oil well is in the quasi-steady state and the non-steady state respectively is determined next, and specifically, the method can be realized by the following steps 203 to 204. The quasi-steady state refers to that parameters such as flow velocity and pressure of fluid are constant and do not change along with the change of time in the flowing process of the fluid; unsteady state refers to the change of parameters such as flow rate and pressure of fluid with time during the fluid flow process.
Step 203: and establishing a first bottom hole pressure calculation model of the source-sink node when the oil well is in a quasi-steady state based on the flow conductivity of each unstructured unit grid.
Specifically, a wellbore radius and a bottom hole production rate of the oil well, and a viscosity of a fluid in the oil well, a permeability of rock and a thickness of a reservoir are obtained, and a first bottom hole pressure calculation model is established based on a flow conductivity of each unstructured cell grid, the wellbore radius and the bottom hole production rate of the oil well, and the viscosity of the fluid in the oil well, the permeability of rock and the thickness of the reservoir.
When the control area at the source-sink node of the oil well is in a pseudo-steady state, the fluid within the control area satisfies a flow continuity condition, and an expression of the following formula (6) can be established for the source-sink node of the oil well. And determining a pressure determination method of the source-sink node when the oil well is in a quasi-steady state based on the following formula (6) and determining the following formula (7) as a first bottom hole pressure calculation model of the oil well.
Wherein, in the above formula (6), pwRefers to the fact of the source/sink nodePressure, i.e. the bottom hole pressure of the well, poIs the theoretical pressure, r, of the source-sink nodewRefers to the wellbore radius, r, of the welloThe theoretical equivalent hole diameter of the oil well, h the thickness of a reservoir layer of the oil well, n the number of the unstructured unit grids included in the at least one unstructured unit grid, q the bottom hole yield of the oil well, pi is a constant and is usually 3.14, and the physical meanings of the rest parameters can be shown in the meanings of the corresponding parameters of the formula (1) or the formula (5). In the above formula (7), the physical meaning of each parameter may be taken into the meaning of the corresponding parameter in the above formula (6).
Step 204: and establishing a second bottom hole pressure calculation model of the source and sink node when the oil well is in an unsteady state.
Specifically, the theoretical pressure of the source sink node at the current moment and the bottom hole yield of the oil well are obtained, the viscosity of fluid in the oil well, the permeability of rock, the thickness of a reservoir and a reservoir pressure guiding coefficient are obtained, and a second bottom hole pressure calculation model is established based on the theoretical pressure of the source sink node at the current moment and the bottom hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rock, the thickness of the reservoir and the reservoir pressure guiding coefficient.
When the well is in a non-steady state, the pressure at any point in the vadose zone of the well can be indicated by the following equation (7) based on the analytical solution of the vadose differential equation for the well, at which time the bottom hole production of the well is determined:
wherein, in the above formula (8), p refers to the fluid pressure in the seepage region, r refers to the distance between any point in the seepage region and the source-sink node, t refers to the time difference between the current time and the initial time, p refers to the time difference between the current time and the initial time, andfthe initial pressure at any point in the well at the initial time, η, is the pressure coefficient of the fluid in the well, and the physical significance of the remaining parameters may be the significance of the corresponding parameters in equation (6) above.
The analytic solution of the seepage differential equation of the oil well can refer to the analytic solution of the seepage differential equation described in the formula (5) -formula (20) on page 101 of the underground oil and gas seepage mechanics, which is not described herein again in the embodiments of the present invention.
Further, an expression of the actual pressure at the source-sink node and the theoretical pressure at the source-sink node of the well may be determined based on equation (8) above. And subtracting the expression of the actual pressure at the source sink node and the expression of the theoretical pressure at the source sink node, determining a pressure determination method of the source sink node when the oil well is in an unsteady state, wherein the pressure determination method is represented by the following formula (9), and determining the following formula (9) as a second bottom hole pressure calculation model of the oil well.
In the above formula (9), the physical meaning of each parameter may be the meaning of the corresponding parameter in the above formula (6) or formula (8).
After determining the first bottom hole pressure calculation model when the well is in a quasi-steady state and the second bottom hole pressure calculation model when the well is in an unsteady state, the bottom hole pressure of the well may be determined based on the first bottom hole pressure calculation model and the second bottom hole pressure calculation model. Specifically, it can be realized as the following steps 205 to 206.
Step 205: and determining a first equivalent hole diameter and a second equivalent hole diameter of the oil well at the current moment.
The first and second equivalent calipers of the well at the current time may be determined prior to determining the bottom hole pressure of the well.
Specifically, for a first equivalent well diameter, the viscosity of a fluid in the oil well, the permeability of the rock and the thickness of the reservoir may be obtained, and then, based on the flow conductivity of each unstructured unit grid in the at least one unstructured unit grid and the coordinate points of all nodes, the average distance between the source and sink nodes and all adjacent nodes on the at least one unstructured unit grid, and the viscosity of the fluid in the oil well, the permeability of the rock and the thickness of the reservoir, the first equivalent well diameter of the oil well at the current time may be determined according to the following formula (10);
wherein, in the above formula (10), ro1Is the first equivalent hole diameter, raMeans an average distance, T, between the source and sink node and all neighboring nodes on the at least one unstructured cell gridilThe flow conductivity coefficient between the node j or k and the node i on the unstructured cell grid e is referred to, and the remaining parameters may be defined as corresponding parameters in the above formula (6) or formula (9).
For a second equivalent well diameter, acquiring the bottom hole yield of the oil well and the time difference between the current time and the initial time, the initial pressure of the source-sink node and the theoretical pressure at the current time, as well as the viscosity of fluid in the oil well, the permeability of rocks, the thickness of a reservoir and a reservoir pressure conductivity coefficient, and determining the second equivalent well diameter of the oil well at the current time by a Newton iteration method according to the following steps (1) - (2);
(1) setting j to 0, and determining the j +1 th hole diameter of the oil well by the following formula (11) based on the preset j th hole diameter;
wherein the content of the first and second substances,
refers to the j +1 th hole diameter at the source-sink node,
refers to the jth hole diameter, p, at the source-sink nodeofRefers to the initial pressure, p, of the source-sink nodeoThe theoretical pressure of the source and sink node at the time t, q the bottom hole yield of the oil well, t the time difference between the current time and the initial time, and the physical meanings of the rest parameters can be the corresponding parameters in the formula (9)The physical meaning of (1).
The preset jth hole diameter may be a wellbore radius of the oil well, and may also be other values, which is not limited in the embodiment of the present invention.
(2) And when the difference value between the j +1 th hole diameter and the j th hole diameter is larger than a preset value, making j equal to j +1, returning to the step of determining the j +1 th hole diameter of the oil well through the following formula (11) based on the j th hole diameter, and when the difference value between the j +1 th hole diameter and the j th hole diameter is smaller than or equal to the preset value, determining the j +1 th hole diameter as the second equivalent hole diameter of the oil well at the current moment.
The preset value may be preset, for example, the preset value may be 0.0001 meter, and the like, which is not limited in the embodiment of the present invention.
Step 206: and when the first equivalent hole diameter is smaller than or equal to the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to the first bottom hole pressure calculation model based on the first equivalent hole diameter.
After the first equivalent hole diameter and the second equivalent hole diameter at the current moment are determined, when the first equivalent hole diameter is smaller than or equal to the second equivalent hole diameter and indicates that the oil well is in a quasi-steady state at the current moment, the bottom hole pressure of the oil well can be determined according to a first bottom hole pressure calculation model of the oil well on the basis of the first equivalent hole diameter. That is, the first equivalent hole diameter is determined as the theoretical equivalent hole diameter of the oil well, and the bottom hole pressure of the oil well at the current time is determined according to the above formula (7).
Step 207: and when the first equivalent hole diameter is larger than the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to the second bottom hole pressure calculation model based on the second equivalent hole diameter.
After the first equivalent hole diameter and the second equivalent hole diameter at the current moment are determined, when the first equivalent hole diameter is larger than the second equivalent hole diameter, the current moment of the oil well is in an unsteady state, and then the bottom hole pressure of the oil well can be determined according to a second bottom hole pressure calculation model based on the second equivalent hole diameter. That is, the second equivalent hole diameter is determined as the theoretical equivalent hole diameter of the oil well, and the bottom hole pressure of the oil well at the current time is determined according to the above formula (9).
Further, after determining the bottom hole pressure of the well through steps 206 and 207, to verify the utility of the first pressure calculation model and the second bottom-hole pressure calculation model of the source-sink node, see figure 2C, the profile of the difference between the bottom hole pressure of the well at any one time and the initial pressure of the well at the initial time, and the variation trend curve of the difference value between the bottom hole pressure of the oil well at any moment and the initial pressure of the oil well at the initial moment is determined by the analytical model of the oil well seepage differential equation, and the change trend curve of the difference value between the bottom hole pressure of the oil well at any moment and the initial pressure of the oil well at the initial moment is the same, thereby determining the utility of the first pressure calculation model and the second bottom-hole pressure calculation model for the source-sink node. In fig. 2C, Δ p refers to the difference between the bottom hole pressure of the well at any one time and the initial pressure at the initial time, Δ p't refers to the trend of change in the difference between the bottom hole pressure of the well at any one time and the initial pressure at the initial time, and t refers to the time difference between any one time and the initial time.
In the embodiment of the invention, at least one non-structural unit grid where the control area at the source and sink node is located is determined after the seepage area of the oil well is dispersed, so that the efficiency of determining the control area is improved. And then determining the flow conductivity of each unstructured unit grid in the at least one unstructured unit grid, and respectively establishing a first bottom hole pressure calculation model of the source-sink node when the oil well is in a quasi-steady state and a second bottom hole pressure calculation model of the source-sink node when the oil well is in an unsteady state through a flow continuity condition and an analytic solution of a seepage differential equation. And then, determining the current state of the oil well through the first equivalent borehole diameter and the second equivalent borehole diameter of the oil well at the current moment, and further determining the borehole pressure of the oil well based on the first borehole pressure calculation model in the pseudo-steady state or the second borehole pressure calculation model in the non-steady state, so that the situation that the pressure of the non-structural unit grid with a certain area is determined as the borehole pressure of the oil well is avoided, the borehole pressure accuracy of the oil well is improved, and more reliable data support is provided for the development scheme design of the oil well.
Fig. 3A is a schematic structural diagram of a device for determining a bottom hole pressure according to an embodiment of the present invention. Referring to fig. 3A, the apparatus includes:
a first determining module 301, configured to determine at least one unstructured cell grid in which a control area at a source/sink node of an oil well is located;
a second determining module 302 for determining a traffic conductivity of each of the at least one grid of unstructured cells;
the establishing module 303 is configured to establish a first bottom hole pressure calculation model of the sourcing and sinking node when the oil well is in a quasi-steady state and a second bottom hole pressure calculation model of the sourcing and sinking node when the oil well is in an unsteady state, based on the flow conductivity of each unstructured unit grid;
a third determination module 304 for determining a bottom hole pressure of the well based on the first bottom hole pressure calculation model and the second bottom hole pressure calculation model.
Optionally, referring to fig. 3B, the first determining module 301 includes:
a first determination unit 3011 for determining a seepage zone at the bottom of the well;
a discrete unit 3012, configured to perform a discrete operation on the seepage area to obtain a plurality of non-structural unit grids included at the bottom of the oil well;
a second determining unit 3013, configured to determine, from the multiple unstructured-cell grids, at least one unstructured-cell grid in which the control area at the source-sink node of the oil well is located.
Optionally, referring to fig. 3C, the establishing module 303 includes:
a first obtaining unit 3031, configured to obtain the wellbore radius and bottom hole production of the oil well, and the viscosity of the fluid in the oil well, the permeability of the rock, and the thickness of the reservoir;
a first establishing unit 3032 for establishing a first bottom hole pressure calculation model based on the flow conductivity of each unstructured cell grid, the wellbore radius and bottom hole production of the oil well, and the viscosity of the fluid in the oil well, the permeability of the rock and the thickness of the reservoir.
Optionally, referring to fig. 3D, the establishing module 303 further includes:
a second obtaining unit 3033, configured to obtain the theoretical pressure of the source/sink node at the current time, the bottom-hole yield of the oil well, the viscosity of a fluid in the oil well, the permeability of rock, the thickness of a reservoir, and a reservoir pressure-guiding coefficient;
a second establishing unit 3034, configured to establish a second downhole pressure calculation model based on the theoretical pressure of the source-sink node at the current time and the downhole yield of the oil well, as well as the viscosity of the fluid in the oil well, the permeability of the rock, the thickness of the reservoir, and the reservoir pressure-leading coefficient.
Optionally, referring to fig. 3E, the third determining module 304 includes:
a third determining unit 3041, configured to determine a first equivalent hole diameter and a second equivalent hole diameter of the oil well at the current time;
a first calculating unit 3042, configured to determine, based on the first equivalent borehole diameter and when the first equivalent borehole diameter is less than or equal to the second equivalent borehole diameter, a bottom hole pressure of the oil well according to the first bottom hole pressure calculation model;
and a second calculating unit 3043, configured to determine, based on the second equivalent borehole diameter, a bottom hole pressure of the oil well according to the second bottom hole pressure calculation model when the first equivalent borehole diameter is larger than the second equivalent borehole diameter.
In the embodiment of the invention, the flow conductivity of each unstructured unit grid in at least one unstructured unit grid where the control area at the source and sink node of the oil well is located is determined, and then a first bottom hole pressure calculation model of the source and sink node when the oil well is in a quasi-steady state and a second bottom hole pressure calculation model of the source and sink node when the oil well is in an unsteady state are respectively established based on the flow conductivity of each unstructured unit grid. And then, determining the bottom hole pressure of the oil well based on the first bottom hole pressure calculation model in the quasi-steady state and the second bottom hole pressure calculation model in the non-steady state, avoiding determining the pressure of the non-structural unit grid with a certain area as the bottom hole pressure of the oil well, improving the accuracy of the bottom hole pressure of the oil well, and providing more reliable data support for the design of a development scheme of the oil well.
It should be noted that: the device for determining the bottom hole pressure provided in the above embodiment is only illustrated by dividing the above functional modules when determining the bottom hole pressure of the oil well, and in practical applications, the above function distribution may be completed by different functional modules according to needs, that is, the internal structure of the equipment is divided into different functional modules to complete all or part of the above described functions. In addition, the determining device of the bottom hole pressure and the determining method embodiment of the bottom hole pressure provided by the above embodiments belong to the same concept, and the specific implementation process is described in the method embodiment, and is not described again.
Fig. 4 shows a block diagram of a terminal 400 according to an exemplary embodiment of the present invention. Referring to fig. 4, the terminal 400 may be: a smartphone, a tablet, a laptop, or a desktop computer. The terminal 400 may also be referred to by other names such as user equipment, portable terminal, laptop terminal, desktop terminal, etc. Referring to fig. 4, the terminal 400 may include a processor 401 and a memory 402.
In some embodiments, the terminal 400 may further optionally include: a peripheral interface 403 and at least one peripheral. The processor 401, memory 402 and peripheral interface 403 may be connected by bus or signal lines. Each peripheral may be connected to the peripheral interface 403 via a bus, signal line, or circuit board. Specifically, the peripheral device includes: at least one of radio frequency circuitry 404, a display screen 405, a positioning component 406, and a power supply 407.
The peripheral interface 403 may be used to connect at least one peripheral related to I/O (Input/Output) to the processor 401 and the memory 402. In some embodiments, processor 401, memory 402, and peripheral interface 403 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 401, the memory 402 and the peripheral interface 403 may be implemented on a separate chip or circuit board, which is not limited by this embodiment.
The Radio Frequency circuit 404 is used for receiving and transmitting RF (Radio Frequency) signals, also called electromagnetic signals. The radio frequency circuitry 404 communicates with communication networks and other communication devices via electromagnetic signals. The rf circuit 404 converts an electrical signal into an electromagnetic signal to transmit, or converts a received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuit 404 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and so forth. The radio frequency circuitry 404 may communicate with other terminals via at least one wireless communication protocol. The wireless communication protocols include, but are not limited to: the world wide web, metropolitan area networks, intranets, generations of mobile communication networks (2G, 3G, 4G, and 5G), Wireless local area networks, and/or WiFi (Wireless Fidelity) networks. In some embodiments, the rf circuit 404 may further include NFC (Near Field Communication) related circuits, which are not limited in this application.
The display screen 405 is used to display a UI (User Interface). The UI may include graphics, text, icons, video, and any combination thereof. When the display screen 405 is a display screen, the display screen 405 also has the ability to capture touch signals on or over the surface of the display screen 405. The touch signal may be input to the processor 401 as a control signal for processing. At this point, the display screen 405 may also be used to provide virtual buttons and/or a virtual keyboard, also referred to as soft buttons and/or a soft keyboard. In some embodiments, the display screen 405 may be one, providing the front panel of the terminal 400; in other embodiments, the display screen 405 may be at least two, respectively disposed on different surfaces of the terminal 400 or in a folded design; in still other embodiments, the display 405 may be a flexible display disposed on a curved surface or a folded surface of the terminal 400. Even further, the display screen 405 may be arranged in a non-rectangular irregular pattern, i.e. a shaped screen. The Display screen 405 may be made of LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), and other materials.
The positioning component 406 is used to locate the current geographic Location of the terminal 400 for navigation or LBS (Location Based Service). The Positioning component 406 may be a Positioning component based on the Global Positioning System (GPS) in the united states, the beidou System in china, or the galileo System in russia.
The power supply 407 is used to supply power to the various components in the terminal 400. The power source 407 may be alternating current, direct current, disposable or rechargeable. When the power source 407 includes a rechargeable battery, the rechargeable battery may be a wired rechargeable battery or a wireless rechargeable battery. The wired rechargeable battery is a battery charged through a wired line, and the wireless rechargeable battery is a battery charged through a wireless coil. The rechargeable battery may also be used to support fast charge technology.
Those skilled in the art will appreciate that the configuration shown in fig. 4 is not intended to be limiting of terminal 400 and may include more or fewer components than those shown, or some components may be combined, or a different arrangement of components may be used.
In the above embodiments, there is also provided a non-transitory computer-readable storage medium comprising instructions for storing at least one instruction for execution by a processor to implement the method provided by the above-described embodiment shown in fig. 1 or fig. 2A.
It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program instructing relevant hardware, where the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (9)
1. A method of determining a bottom hole pressure, the method comprising:
determining at least one unstructured cell grid in which a control area at a source-sink node of an oil well is located;
determining a flow conductivity of each of the at least one unstructured cell grid based on a permeability of rock in the oil well, a viscosity of a fluid in the oil well, a volume factor of the fluid in the oil well, an area of a control area at a source-sink node in the unstructured cell grid, a fluid pressure within the unstructured cell grid, and a zone extent of the unstructured cell grid;
when the control area at the source sink node of the oil well is in a quasi-steady state, establishing a first bottom hole pressure calculation model of the source sink node when the oil well is in the quasi-steady state on the basis of the flow conductivity of each unstructured unit grid, the shaft radius and the bottom hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rocks and the thickness of a reservoir;
when the control area of the source and sink node of the oil well is in an unsteady state, establishing a second bottom hole pressure calculation model of the source and sink node when the oil well is in the unsteady state on the basis of the flow conductivity of each unstructured unit grid, the theoretical pressure of the source and sink node at the current moment, the bottom hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rocks, the thickness of a reservoir and a reservoir pressure conduction coefficient;
determining a first equivalent hole diameter and a second equivalent hole diameter of the oil well at the current moment;
when the first equivalent hole diameter is smaller than or equal to the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to the first bottom hole pressure calculation model based on the first equivalent hole diameter;
and when the first equivalent hole diameter is larger than the second equivalent hole diameter, determining the bottom hole pressure of the oil well according to the second bottom hole pressure calculation model based on the second equivalent hole diameter.
2. The method of claim 1, wherein determining at least one grid of unstructured cells in which to locate a control area at a source-sink node of a well comprises:
determining a seepage zone at the bottom of the oil well;
dispersing the seepage area to obtain a plurality of unstructured unit grids included at the bottom of the oil well;
determining at least one unstructured cell grid from the plurality of unstructured cell grids in which a control area at a source-sink node of the well resides.
3. The method of claim 1, wherein the establishing a first downhole pressure calculation model of the source-sink node while the well is in a pseudo-steady state based on the flow conductivity of each grid of unstructured cells comprises:
acquiring the well bore radius and the bottom hole yield of the oil well, and the viscosity of fluid in the oil well, the permeability of rock and the thickness of a reservoir;
the first downhole pressure calculation model is built based on the flow conductivity of each grid of unstructured cells, the wellbore radius and downhole production of the well, and the viscosity of the fluid in the well, permeability of the rock, and thickness of the reservoir.
4. The method of claim 1 or 3, wherein said establishing a second downhole pressure calculation model of the source-sink node when the well is in an unsteady state comprises:
obtaining theoretical pressure of the source and sink node at the current moment, bottom hole yield of the oil well, viscosity of fluid in the oil well, permeability of rock, thickness of a reservoir and a reservoir pressure guiding coefficient;
and establishing the second bottom-hole pressure calculation model based on the theoretical pressure of the source-sink node at the current moment, the bottom-hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rock, the thickness of a reservoir and a reservoir pressure-guiding coefficient.
5. An apparatus for determining bottom hole pressure, the apparatus comprising:
the first determination module is used for determining at least one unstructured unit grid where the control area at the source and sink nodes of the oil well is located;
a second determining module for determining a traffic conductivity of each of the at least one grid of unstructured cells;
the establishing module is used for establishing a first bottom hole pressure calculation model of the source sink node when the oil well is in a quasi-steady state and establishing a second bottom hole pressure calculation model of the source sink node when the oil well is in an unsteady state based on the flow conductivity of each unstructured unit grid;
a third determination module to determine a bottom hole pressure of the well based on the first bottom hole pressure computational model and the second bottom hole pressure computational model.
6. The apparatus of claim 5, wherein the first determining module comprises:
a first determination unit for determining a seepage area at the bottom of the oil well;
the discrete unit is used for dispersing the seepage area to obtain a plurality of non-structural unit grids included at the bottom of the oil well;
and the second determination unit is used for determining at least one unstructured unit grid in which the control area at the source and sink nodes of the oil well is located from the plurality of unstructured unit grids.
7. The apparatus of claim 5, wherein the establishing module comprises:
the first acquisition unit is used for acquiring the shaft radius and the bottom-hole yield of the oil well, and the viscosity of fluid in the oil well, the permeability of rock and the thickness of a reservoir;
a first establishing unit for establishing the first bottom hole pressure calculation model based on the flow conductivity of each unstructured cell grid, the wellbore radius and the bottom hole production of the oil well, and the viscosity of the fluid in the oil well, the permeability of the rock and the thickness of the reservoir.
8. The apparatus of claim 5 or 7, wherein the establishing module comprises:
the second acquisition unit is used for acquiring the theoretical pressure of the source-sink node at the current moment, the bottom-hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rocks, the thickness of a reservoir and a reservoir pressure-guiding coefficient;
and the second establishing unit is used for establishing the second bottom hole pressure calculation model based on the theoretical pressure of the source-sink node at the current moment, the bottom hole yield of the oil well, the viscosity of fluid in the oil well, the permeability of rock, the thickness of a reservoir and the reservoir pressure guiding coefficient.
9. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the method of any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810158974.0A CN110083852B (en) | 2018-02-26 | 2018-02-26 | Method and device for determining bottom hole pressure and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810158974.0A CN110083852B (en) | 2018-02-26 | 2018-02-26 | Method and device for determining bottom hole pressure and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110083852A CN110083852A (en) | 2019-08-02 |
| CN110083852B true CN110083852B (en) | 2022-02-01 |
Family
ID=67412697
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810158974.0A Active CN110083852B (en) | 2018-02-26 | 2018-02-26 | Method and device for determining bottom hole pressure and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110083852B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111832164B (en) * | 2020-07-03 | 2023-05-12 | 中国石油大学(北京) | Shale gas yield prediction method, device and equipment |
| CN112329355B (en) * | 2020-10-22 | 2022-10-04 | 中国石油天然气股份有限公司 | Method and device for determining single-well control area, computer equipment and storage medium |
| CN113361771B (en) * | 2021-06-04 | 2023-04-18 | 合肥工业大学 | Method and device for determining the pressure of a storage tank |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103400020B (en) * | 2013-04-16 | 2016-03-16 | 西南石油大学 | A kind of numerical reservoir simulation method calculating many crossing discrete fractures flow conditions |
| DE112013007651T5 (en) * | 2013-11-27 | 2016-09-01 | Landmark Graphics Corporation | Method and apparatus for optimized vacuum drilling |
| CN104695928B (en) * | 2015-01-27 | 2017-05-24 | 中国地质调查局油气资源调查中心 | Method for evaluating volume transformation capacity of vertical well of fractured tight oil reservoir |
| CN105019894B (en) * | 2015-07-29 | 2017-11-14 | 长江大学 | A kind of stratified reservoir inter well connectivity method for establishing model and system |
| CN107578342B (en) * | 2017-07-17 | 2020-09-08 | 中国石油大学(华东) | Model coupling exhaustion method-based method for realizing low-permeability reservoir interval working system optimization |
-
2018
- 2018-02-26 CN CN201810158974.0A patent/CN110083852B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110083852A (en) | 2019-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110083852B (en) | Method and device for determining bottom hole pressure and storage medium | |
| US20210115763A1 (en) | Algorithm for optimal icd configuration using a coupled wellbore-reservoir model | |
| EP2956618B1 (en) | Systems, methods, and computer-readable media for modeling complex wellbores in field-scale reservoir simulation | |
| CA3068260C (en) | Parallel-processing of invasion percolation for large-scale, high-resolution simulation of secondary hydrocarbon migration | |
| CN111396013B (en) | Method and device for determining shale gas well fracturing modification scheme and storage medium | |
| CN110083851B (en) | Method and device for determining bottom hole pressure of gas well and storage medium | |
| CN109978303B (en) | Method and device for determining exploitation measures of heavy oil reservoir and storage medium | |
| CN108487902B (en) | Method and device for determining bottom hole flow pressure and storage medium | |
| US20150331976A1 (en) | Forecasting oil production for multiply-fractured horizontal wells | |
| CN115984503B (en) | Geological profile generation method, system, electronic equipment and medium | |
| CN110109178B (en) | Method and device for determining dominant sand body distribution area and storage medium | |
| CN109961195B (en) | Method and device for determining energy efficiency influence factors of oilfield water injection system and storage medium | |
| CN114320236B (en) | Method and device for adjusting oil-gas well operation scheme and storage medium | |
| CN112507551B (en) | Unstructured dynamic mesh generation method and device | |
| CN113931621B (en) | Method, device and storage medium for determining gas well hydrops information | |
| CN110894782B (en) | Method and device for determining gas storage capacity of oil reservoir and storage medium | |
| CN108629463B (en) | Ground stress change prediction method and device | |
| CN103425810A (en) | Behavioral modeling methods for clock and data recovery circuit and analog circuits | |
| CN111188608B (en) | Method and device for determining shale gas well fracturing modification scheme and storage medium | |
| CN111861780A (en) | Method and device for determining production of production layer in production well and storage medium | |
| CN111075440A (en) | Unconformity oil and gas reservoir scale prediction method and device | |
| CN113496074B (en) | Sand stabilizing method and device | |
| Absalyamov et al. | Coupled static and dynamic high resolution modelling on high performance computer clusters | |
| CN113107449B (en) | Method and device for determining fracturing parameters of shale reservoir multi-stage fractured horizontal well | |
| CN114065459B (en) | Pipe fitting selection method and device and storage medium |
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 |