CN112069620B - Method and device for determining migration degree of natural gas hydrate reservoir particles - Google Patents
Method and device for determining migration degree of natural gas hydrate reservoir particles Download PDFInfo
- Publication number
- CN112069620B CN112069620B CN202010926449.6A CN202010926449A CN112069620B CN 112069620 B CN112069620 B CN 112069620B CN 202010926449 A CN202010926449 A CN 202010926449A CN 112069620 B CN112069620 B CN 112069620B
- Authority
- CN
- China
- Prior art keywords
- throat
- reservoir
- particles
- migration
- reservoir particles
- 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/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application provides a method and a device for determining the migration degree of natural gas hydrate reservoir particles, which comprises the following steps: and establishing a pore network model to simulate the flow of reservoir particles in the network model under the phase equilibrium condition. And obtaining capillary force according to the transport condition parameters of the reservoir particles, and determining the throat through which the reservoir particles flow according to the capillary force. And when the concentration of the reservoir particles is less than the preset concentration, if the diameter of a single reservoir particle is less than the diameter of the throat or the volume of all the reservoir particles entering the throat is less than the volume of the throat, updating the radius size of the throat to be the first throat radius. And when the concentration of the reservoir particles is greater than or equal to the preset concentration and the diameter of the reservoir particles is smaller than one third of the diameter of the throat, if the volume of all the reservoir particles entering the throat is smaller than the volume of the throat, updating the radius size of the throat to be the first throat radius. And judging the migration degree of the reservoir particles according to a first preset rule according to the first throat radius. It is possible to understand the effect of particle transport on fluid flow during hydrate dissociation.
Description
Technical Field
The application belongs to the technical field of oil exploitation, and particularly relates to a method and a device for determining migration degree of natural gas hydrate reservoir particles.
Background
Natural gas hydrate is receiving increasing attention from the world as a potentially vast amount of unconventional energy. The seabed natural gas hydrate is mainly stored in a shallow layer of deep water, the cementation between the seabed natural gas hydrate and a reservoir stratum is poor in strength in the decomposition process, and part of reservoir stratum particles fall off and flow along with fluid, so that the sand production phenomenon is caused in the process of exploiting the natural gas hydrate.
In the process of exploitation, the natural gas hydrate loses the original stability due to the decomposition of the natural gas hydrate and the sudden and large generation of gas. It is therefore important to know the phase equilibrium of natural gas hydrates. Moreover, natural gas hydrates are sensitive to changes in temperature and pressure, and slight changes can cause the natural gas hydrates to decompose. Different dissociation equilibrium conditions of natural gas hydrate were measured by a low temperature sapphire cell, and the equilibrium temperature of the gas hydrate was found to increase with increasing carbon dioxide concentration and decrease with increasing nitrogen concentration. Sunrey et al propose a thermodynamic model that can predict the phase equilibrium condition of methane hydrate in porous media, and predict the solubility of methane in the state of coexistence of methane hydrate and water balance.
However, in the simulation process, the influence factor of reservoir particle flow is not considered, so that the model cannot fit the actual situation.
Disclosure of Invention
In order to overcome the above-mentioned defects of the prior art, the technical problem to be solved by the present invention is to provide a method and an apparatus for determining the degree of migration of particles in a gas hydrate reservoir, which can understand the influence of the migration of particles on the fluid flow during the hydrate decomposition process.
The specific technical scheme of the invention is as follows:
the invention provides a method for determining the migration degree of natural gas hydrate reservoir particles, which comprises the following steps:
establishing a three-dimensional network model comprising pores and throats, and setting condition parameters and initial conditions of the reservoir particle migration;
setting the reservoir particles in the network model, and simulating the flow condition of the reservoir particles in the network model under the phase equilibrium condition;
obtaining capillary force according to the condition parameters of the migration of the reservoir particles, and determining a throat through which the reservoir particles flow according to the capillary force;
when the concentration of the reservoir particles is less than the preset concentration, if the diameter of a single reservoir particle is less than the diameter of the throat or the volume of all the reservoir particles entering the throat is less than the volume of the throat, updating the radius size of the throat to be a first throat radius;
when the concentration of the reservoir particles is greater than or equal to a preset concentration and the diameter of the reservoir particles is smaller than one third of the diameter of the throat, if the volume of all the reservoir particles entering the throat is smaller than the volume of the throat, updating the radius size of the throat to be a first throat radius;
and judging the migration degree of the reservoir particles according to a first preset rule and the first throat radius.
In a preferred embodiment, the first predetermined rule is:
if the first throat radius =0, the throat is blocked;
if the first throat radius is larger than 0, the throat is smooth.
In a preferred embodiment, when the reservoir particle concentration is less than a predetermined concentration,
if the diameter of a single reservoir particle is larger than or equal to the diameter of the throat or the volume of all reservoir particles entering the throat is larger than or equal to the volume of the throat, judging the migration degree of the reservoir particles according to a second preset rule;
and when the concentration of the reservoir particles is greater than or equal to the preset concentration and the diameter of the reservoir particles is greater than or equal to one third of the diameter of the throat, judging the migration degree of the reservoir particles according to a second preset rule.
In a preferred embodiment, the second predetermined rule is: whether other throats are present or not is judged,
if other throats exist, continuing to simulate the flowing situation of the reservoir particles in the network model comprising the other throats under the phase equilibrium condition;
if no other throat is present, the throat is blocked.
In a preferred embodiment, the initial condition is that a throat through which the reservoir particulates do not flow exists in the network model, or the number of times the reservoir particulates flow through the throat and the pores is calculated to be less than a specified maximum number of displacements.
In a preferred embodiment, the throat radius size is updated to a first throat radius according to the following equation:
wherein, V pt Volume sum of existing reservoir particles in the throat, unit: cubic meter, V t Initial throat volume, unit: cubic meter, d t_new First throat diameter, unit: meter, L is the original throat length, unit: and (4) rice.
In a preferred embodiment, the condition parameters for reservoir particle migration include: spacing between adjacent pore centerpoints of the network model, inlet pressure, outlet pressure, wetting angle, and interfacial tension.
In a preferred embodiment, the capillary force is calculated according to the following formula:
wherein F is capillary force in N; σ is interfacial tension; the units N/m, theta are the wetting angle, the units DEG, r i Is the throat radius in m.
In a preferred embodiment, the predetermined concentration is 5%.
In a preferred embodiment, the pore radius is greater than the throat radius, and if all of the throats are plugged, the pores are also plugged.
In addition, the application also provides a device for determining the migration degree of the natural gas hydrate reservoir particles, which comprises the following modules:
the establishing module is configured to establish a three-dimensional network model comprising pores and throats, and set condition parameters and initial conditions of the reservoir particle migration;
a flow module configured to set the reservoir particle in the network model and simulate a flow of the reservoir particle in the network model under a phase equilibrium condition;
the selection module is configured to obtain capillary force according to the condition parameters of the migration of the reservoir particles and determine a throat through which the reservoir particles flow according to the capillary force;
a first update module configured to update the throat radius size to a first throat radius when the concentration of reservoir particulates is less than a predetermined concentration, if a single reservoir particulate diameter is less than the throat diameter, or if the volume of all reservoir particulates entering the throat is less than the throat volume;
a second updating module configured to update the throat radius size to a first throat radius if all reservoir particle volumes entering the throat are less than the throat volume when the reservoir particle concentration is greater than or equal to a preset concentration and the reservoir particle diameter is less than one-third of the throat diameter;
a determination module configured to determine the extent of reservoir particle migration according to a first predetermined rule based on the first throat radius.
In addition, the present application also provides an apparatus for determining a degree of migration of natural gas hydrate reservoir particles, comprising a memory and a processor, the memory having stored therein a computer program which, when executed by the processor, performs the steps of: a method of determining the extent of migration of reservoir particles as described above.
Borrow by above technical scheme, the beneficial effect of this application lies in:
the method and the device for determining the migration degree of the natural gas hydrate reservoir particles have important reference significance for the sand production phenomenon caused by the hydrate exploitation process. Under the condition of phase equilibrium, the fluid flow condition in a natural gas hydrate reservoir containing particle flow is simulated, and the influence of the continuous increase of particles on the flow in the decomposition process of the hydrate can be known. Namely the influence of the sand production phenomenon on the flow of the gas-water two-phase flow.
Specific embodiments of the present application are disclosed in detail with reference to the following description and drawings, indicating the manner in which the principles of the application may be employed. It should be understood that the embodiments of the present application are not so limited in scope. The embodiments of the application include many variations, modifications and equivalents within the spirit and scope of the appended claims.
Features that are described and/or illustrated with respect to one embodiment may be used in the same way or in a similar way in one or more other embodiments, in combination with or instead of the features of the other embodiments.
It should be emphasized that the term "comprises/comprising" when used herein, is taken to specify the presence of stated features, integers, steps or components but does not preclude the presence or addition of one or more other features, integers, steps or components.
Drawings
The drawings described herein are for illustration purposes only and are not intended to limit the scope of the present disclosure in any way. In addition, the shapes, the proportional sizes, and the like of the respective members in the drawings are merely schematic for assisting understanding of the present application, and do not specifically limit the shapes, the proportional sizes, and the like of the respective members in the present application. Those skilled in the art, having the benefit of the teachings of this application, may select various possible shapes and proportional sizes to implement the present application, depending on the particular situation. In the drawings:
fig. 1 is a flow chart of a method of determining a degree of migration of natural gas hydrate reservoir particulates in accordance with an embodiment of the present disclosure;
fig. 2 is a block diagram of an apparatus for determining a migration degree of natural gas hydrate reservoir particles according to an embodiment of the present disclosure.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only some embodiments of the present application, not all embodiments. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments in the present application without making creative efforts shall fall within the protection scope of the present application.
It will be understood that when an element is referred to as being "disposed on" another element, it can be directly on the other element or intervening elements may also be present. When an element is referred to as being "connected" to another element, it can be directly connected to the other element or intervening elements may also be present. The terms "vertical," "horizontal," "left," "right," and the like as used herein are for illustrative purposes only and do not denote a single embodiment.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein in the description of the present application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
As shown in fig. 1, the present application provides a method of determining a degree of migration of natural gas hydrate reservoir particulates, the method comprising the steps of:
s1: and establishing a three-dimensional network model comprising pores and throats, and setting condition parameters and initial conditions for the migration of reservoir particles.
S2: and setting the reservoir particles in the network model, and simulating the flow condition of the reservoir particles in the network model under the phase equilibrium condition.
S3: and obtaining capillary force according to the condition parameters of the migration of the reservoir particles, and determining a throat through which the reservoir particles flow according to the capillary force.
S4: and when the concentration of the reservoir particles is less than the preset concentration, if the diameter of a single reservoir particle is less than the diameter of the throat or the volume of all reservoir particles entering the throat is less than the volume of the throat, updating the radius size of the throat to be a first throat radius.
S5: and when the concentration of the reservoir particles is greater than or equal to a preset concentration and the diameter of the reservoir particles is smaller than one third of the diameter of the throat, if the volume of all the reservoir particles entering the throat is smaller than the volume of the throat, updating the radius size of the throat to be a first throat radius.
S6: and judging the migration degree of the reservoir particles according to a first preset rule and the first throat radius.
In this embodiment, a network model is first established containing pores and throats for simulating the flow of fluids in the reservoir, the pore space in the reservoir being simulated with pores and throats, the pores being approximately spherical, the throats being used to connect adjacent pores while being used to study the flow of fluids. This model is a regular cubic three-dimensional structure (10 (voids) × 10 (voids)), assuming that each aperture connects 6 throats, i.e., the number of bits is 6. Pore and throat diameters of 10 -5 -10 -4 Are randomly distributed among rice, andthe diameter of the pore is slightly larger than that of the throat, so that the influence of the particle flow on the whole seepage process can be more prominent.
And then setting condition parameters and initial conditions of reservoir particle migration. The condition parameters of the reservoir particle migration comprise: spacing between adjacent pore centroids, inlet pressure, outlet pressure, wetting angle, and interfacial tension of the network model. Typically, the spacing between the centers of adjacent pores (the spacing from one pore center to the next pore center) is 10 -4 And m is selected. The inlet pressure is the pressure value corresponding to different hydrate saturation degrees under the phase equilibrium condition, and the outlet pressure is 0MPa. Wetting angle theta is 100 DEG, interfacial tension sigma is 34.33x10 -3 mN/m。
The initial condition is that a throat through which the reservoir particles do not flow exists in the network model, or the times of flowing the reservoir particles through the throat and the pores are calculated, wherein the times are less than the specified maximum displacement times. And setting the reservoir particles in the network model under the condition that the initial condition is met, and simulating the flowing condition of the reservoir particles in the network model under the phase equilibrium condition.
In particular, for better study on the particle flow law in natural gas hydrate (methane hydrate), all particles in the network model are set as spheres. In the initial state, i.e. the first experiment in table 1 (temperature 277K, pressure difference (difference between inlet and outlet pressure) 3.8MPa, hydrate saturation 58%), 5 particles were placed in each pore, 30 particles in each throat, and the particle sizes were randomly distributed over a range. In the process of decomposing the hydrate, the number of reservoir particles in each throat is increased by 6 when the hydrate is decomposed by 1%, and the number of natural gas hydrate particles in each pore is increased by 3. The number of particles present in the pores and throat corresponding to different hydrate saturations before the start of the test is shown in table 1. In the network model, the change of reservoir particles in pores is not considered, and the reservoir particles in the throat are assumed to be deposited only in the original pores, and the reservoir particles in the throat pass through the pores but are not deposited in the pores and are deposited only in the throat. After the first test is finished, 2-7 sets of tests can be sequentially performed.
TABLE 1 values of the Condition parameters under phase equilibrium conditions under different test conditions
Then, capillary force can be obtained according to the condition parameters of the migration of the reservoir particles, and the throat through which the reservoir particles flow is determined according to the capillary force. As water is used as a displacement phase, every time water enters one pore and the throat adjacent to the pore, 80% of the number of the existing particles in the throat are deposited, and 20% of the particles flow into the next throat along with the displacement phase. Since each pore has at least two connected throats, it is necessary to select which throat to enter during particle flow, hydrate decomposition. The order of entry into the throat is judged according to the magnitude of capillary force, and usually the particles preferentially enter the throat with smaller capillary force. And finishing the program operation after all the particles are moved. Specifically, the capillary force may be calculated according to the following formula:
wherein F is capillary force in N; σ is interfacial tension; the units N/m, theta are the wetting angle, the units DEG, r i Is the throat radius in m.
When the concentration of the reservoir particles is less than a preset concentration (5%), if the diameter of a single reservoir particle is less than the diameter of the throat, or the volume of all reservoir particles entering the throat is less than the volume of the throat, updating the radius size of the throat to be a first throat radius:
wherein, V pt Volume sum of existing reservoir particles in the throat, unit: cubic meter, V t Is the initial throatVolume, unit: cubic meter, d t_new First throat diameter, unit: meter, L is the original throat length, unit: and (4) rice.
When the concentration of the reservoir particles is greater than or equal to a preset concentration (5%) and the diameter of the reservoir particles is smaller than one third of the diameter of the throat, if the volume of all the reservoir particles entering the throat is smaller than the volume of the throat, updating the radius size of the throat to be a first throat radius;
wherein, V pt Volume sum of existing reservoir particles in the throat, unit: cubic meter, V t Initial throat volume, unit: cubic meter, d t_new First throat diameter, unit: meter, L is the original throat length in units: and (4) rice.
And then, according to the radius of the first throat, judging the migration degree of the reservoir particles according to a first preset rule.
The first predetermined rule is specifically: if the first throat radius =0, the throat is blocked. If the first throat radius is larger than 0, the throat is unblocked.
However, when the concentration of the reservoir particles is less than the preset concentration (5%), if the diameter of a single reservoir particle is greater than or equal to the diameter of the throat, or the volume of all reservoir particles entering the throat is greater than or equal to the volume of the throat, the migration degree of the reservoir particles is judged according to a second preset rule. Or when the concentration of the reservoir particles is greater than or equal to a preset concentration (5%) and the diameter of the reservoir particles is greater than or equal to one third of the diameter of the throat, judging the migration degree of the reservoir particles according to a second preset rule.
The second predetermined rule is specifically: and judging whether other throats exist or not, if so, continuing to simulate the flow condition of the reservoir particles in the network model comprising the other throats under the phase balance condition. And if no other throat is present, determining that the throat is blocked.
It is noted that when the pore radius is greater than the throat radius, if all of the throats are plugged, the pores are also plugged.
Therefore, by using the method for determining the migration degree of the particles in the natural gas hydrate reservoir, the influence of the deposition, blockage and flow of the particles on the change of the pore space in the reservoir and the flow of the fluid in the pore space in the hydrate decomposition process can be known.
Based on the same inventive concept, the embodiment of the present invention also provides an apparatus for determining the migration degree of natural gas hydrate reservoir particles, as described in the following embodiment. Because the principle of solving the problems of the device for determining the migration degree of the particles of the natural gas hydrate reservoir is similar to that of the method for determining the migration degree of the particles of the natural gas hydrate reservoir, the implementation of the device for determining the migration degree of the particles of the natural gas hydrate reservoir can be referred to the implementation of the method for determining the migration degree of the particles of the natural gas hydrate reservoir, and repeated parts are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.
As shown in fig. 2, the present invention also provides an apparatus for determining a degree of migration of natural gas hydrate reservoir particles, the apparatus comprising the following modules:
a building module 101 configured to build a three-dimensional network model including pores and throats, set condition parameters and initial conditions for the reservoir particle migration.
A flow module 102 configured to set the reservoir particle in the network model and simulate a flow of the reservoir particle in the network model under phase equilibrium conditions.
A selection module 103 configured to obtain a capillary force according to the condition parameter of the reservoir particle migration, and determine a throat through which the reservoir particle flows according to the capillary force.
A first update module 104 configured to update the throat radius size to a first throat radius when the reservoir particle concentration is less than a preset concentration, if a single reservoir particle diameter is less than the throat diameter, or all reservoir particle volumes entering the throat are less than the throat volume.
A second updating module 105 configured to update the throat radius size to the first throat radius if all the volume of the reservoir particles entering the throat is less than the throat volume when the reservoir particle concentration is greater than or equal to a predetermined concentration and the reservoir particle diameter is less than one-third of the throat diameter.
A determining module 106 configured to determine the reservoir particle migration degree according to a first predetermined rule based on the first throat radius.
In a preferred embodiment, the first predetermined rule is:
if the first throat radius =0, the throat is blocked;
if the first throat radius is larger than 0, the throat is unblocked.
In a preferred embodiment, when the reservoir particle concentration is less than a predetermined concentration,
if the diameter of a single reservoir particle is larger than or equal to the diameter of the throat or the volume of all reservoir particles entering the throat is larger than or equal to the volume of the throat, judging the migration degree of the reservoir particles according to a second preset rule;
and when the concentration of the reservoir particles is greater than or equal to the preset concentration and the diameter of the reservoir particles is greater than or equal to one third of the diameter of the throat, judging the migration degree of the reservoir particles according to a second preset rule.
In a preferred embodiment, the second predetermined rule is: judging whether other throats exist or not, if so, continuing to simulate the flowing condition of the reservoir particles in the network model comprising the other throats under the phase balance condition; if no other throat is present, the throat is blocked.
In a preferred embodiment, the initial condition is that a throat through which the reservoir particulates do not flow exists in the network model, or that the number of times the reservoir particulates flow through the throat and the pores is calculated to be less than a specified maximum number of displacements.
In a preferred embodiment, the throat radius size is updated to a first throat radius according to the following equation:
wherein, V pt Volume sum of existing reservoir particles in the throat, unit: cubic meter, V t Initial throat volume, unit: cubic meter, d t_new First throat diameter, unit: meter, L is the original throat length, unit: and (4) rice.
In a preferred embodiment, the condition parameters for reservoir particle migration include: spacing between adjacent pore center points of the network model, inlet pressure, outlet pressure, wetting angle, and interfacial tension.
In a preferred embodiment, the capillary force is calculated according to the following formula:
wherein F is capillary force in N; σ is interfacial tension; the unit N/m, theta is the wetting angle, the unit DEG, r i Is the throat radius in mm. (Unit)
In a preferred embodiment, the predetermined concentration is 5%.
In a preferred embodiment, the aperture radius is greater than the throat radius, and if all of the throats are plugged, the apertures are also plugged.
In addition, the invention also provides a device for determining the migration degree of natural gas hydrate reservoir particles, which comprises a memory and a processor, wherein the memory stores a computer program, and the computer program realizes the steps of the method for determining the migration degree of reservoir particles.
In this embodiment, the memory may include a physical device for storing information, and typically, the information is digitized and stored in a medium using an electrical, magnetic, or optical method. The memory according to this embodiment may further include: devices that store information using electrical energy, such as RAM, ROM, etc.; devices that store information by magnetic energy, such as hard disks, floppy disks, tapes, core memories, bubble memories, usb disks; devices for storing information optically, such as CDs or DVDs. Of course, there are other ways of memory, such as quantum memory, graphene memory, and so forth.
In this embodiment, the processor may be implemented in any suitable manner. For example, the processor may take the form of, for example, a microprocessor or processor and a computer-readable medium that stores computer-readable program code (e.g., software or firmware) executable by the (micro) processor, logic gates, switches, an Application Specific Integrated Circuit (ASIC), a programmable logic controller and embedded microcontroller, and so forth.
The specific functions of the server, the processor and the memory thereof implemented by the embodiments of the present specification can be explained in comparison with the foregoing embodiments of the present specification.
In another embodiment, a software is provided, and the software is used for executing the technical solutions described in the above embodiments and preferred embodiments.
In another embodiment, a storage medium is provided, where the software is stored, and the storage medium includes but is not limited to: optical disks, floppy disks, hard disks, erasable memory, etc.
As can be seen from the above description, the embodiments of the present invention achieve the following technical effects: the method and the device for determining the migration degree of the natural gas hydrate reservoir particles have important reference significance for the sand production phenomenon caused by the hydrate exploitation process. Under the condition of phase equilibrium, the fluid flow condition in a natural gas hydrate reservoir containing particle flow is simulated, and the influence of the continuous increase of particles on the flow in the decomposition process of the hydrate can be known. Namely the influence of the sand production phenomenon on the flow of the gas-water two-phase flow.
It will be apparent to those skilled in the art that the modules or steps of the embodiments of the invention described above may be implemented by a general purpose computing device, they may be centralized on a single computing device or distributed across a network of multiple computing devices, and alternatively, they may be implemented by program code executable by a computing device, such that they may be stored in a storage device and executed by a computing device, and in some cases, the steps shown or described may be performed in an order different than that described herein, or they may be separately fabricated into individual integrated circuit modules, or multiple ones of them may be fabricated into a single integrated circuit module. Thus, embodiments of the invention are not limited to any specific combination of hardware and software.
A plurality of elements, components, parts or steps can be provided by a single integrated element, component, part or step. Alternatively, a single integrated element, component, part or step may be divided into separate plural elements, components, parts or steps. The disclosure of "a" or "an" to describe an element, ingredient, component or step is not intended to foreclose other elements, ingredients, components or steps.
It is to be understood that the above description is intended to be illustrative, and not restrictive. Many embodiments and many applications other than the examples provided will be apparent to those of skill in the art upon reading the above description. The scope of the present teachings should, therefore, be determined not with reference to the above description, but should instead be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. The disclosures of all articles and references, including patent applications and publications, are hereby incorporated by reference for all purposes. The omission in the foregoing claims of any aspect of subject matter that is disclosed herein is not a disclaimer of such subject matter, nor is it to be construed that the applicant does not consider such subject matter to be part of the disclosed subject matter.
Claims (12)
1. A method of determining the extent of migration of gas hydrate reservoir particulates comprising the steps of:
establishing a three-dimensional network model comprising pores and throats, and setting condition parameters and initial conditions for the migration of reservoir particles;
setting the reservoir particles in the network model, and simulating the flowing condition of the reservoir particles in the network model under the phase equilibrium condition;
obtaining capillary force according to the condition parameters of the migration of the reservoir particles, and determining a throat through which the reservoir particles flow according to the capillary force;
when the concentration of the reservoir particles is less than the preset concentration, if the diameter of a single reservoir particle is less than the diameter of the throat or the volume of all the reservoir particles entering the throat is less than the volume of the throat, updating the radius size of the throat to be a first throat radius;
when the concentration of the reservoir particles is greater than or equal to a preset concentration and the diameter of the reservoir particles is smaller than one third of the diameter of the throat, if the volume of all the reservoir particles entering the throat is smaller than the volume of the throat, updating the radius size of the throat to be a first throat radius;
and judging the migration degree of the reservoir particles according to a first preset rule and the first throat radius.
2. A method of determining a degree of natural gas hydrate reservoir particle migration according to claim 1, wherein the first predetermined rule is:
if the first throat radius =0, the throat is blocked;
if the first throat radius is larger than 0, the throat is smooth.
3. The method of determining a degree of natural gas hydrate reservoir particle migration according to claim 1, wherein when the reservoir particle concentration is less than a preset concentration,
if the diameter of a single reservoir particle is larger than or equal to the diameter of the throat or the volume of all reservoir particles entering the throat is larger than or equal to the volume of the throat, judging the migration degree of the reservoir particles according to a second preset rule;
and when the concentration of the reservoir particles is greater than or equal to the preset concentration and the diameter of the reservoir particles is greater than or equal to one third of the diameter of the throat, judging the migration degree of the reservoir particles according to a second preset rule.
4. A method of determining a degree of natural gas hydrate reservoir particle migration according to claim 1, wherein the second predetermined rule is: whether other throats are present or not is judged,
if other throats exist, continuing to simulate the flowing condition of the reservoir particles in the network model containing the other throats under the condition of phase equilibrium;
if no other throat is present, the throat is blocked.
5. A method of determining the extent of migration of gas hydrate reservoir particles as claimed in claim 1 wherein the initial condition is the presence of a throat in the network model through which the reservoir particles do not flow or the number of times the reservoir particles flow through the throat and the pores is calculated to be less than the specified maximum number of displacements.
6. A method of determining a degree of migration of natural gas hydrate reservoir particulates according to claim 1, wherein the throat radius size is updated to a first throat radius according to the following formula:
wherein, V pt Volume sum of existing reservoir particles in the throat, unit: cubic meter, V t Initial throat volume, unit: cubic meter, d t_new First throat diameter, unit: meter, L is the original throat length, unit: and (4) rice.
7. A method of determining a degree of gas hydrate reservoir particle migration as claimed in claim 1, wherein the condition parameters of reservoir particle migration include: spacing between adjacent pore center points of the network model, inlet pressure, outlet pressure, wetting angle, and interfacial tension.
8. A method of determining a degree of migration of natural gas hydrate reservoir particulates according to claim 1, characterized in that the capillary force is calculated according to the following formula:
wherein F is capillary force in N; σ is interfacial tension; the unit N/m, theta is the wetting angle, the unit DEG, r i Is the throat radius in m.
9. The method of determining the extent of migration of natural gas hydrate reservoir particles according to claim 1, wherein the predetermined concentration is 5%.
10. A method of determining the extent of migration of natural gas hydrate reservoir particulates according to claim 1, wherein the pore radius is greater than the throat radius, and if all of the throats are plugged, the pores are also plugged.
11. An apparatus for determining a degree of migration of gas hydrate reservoir particles, comprising:
the establishing module is configured to establish a three-dimensional network model comprising pores and throats, and set condition parameters and initial conditions of the reservoir particle migration;
a flow module configured to set the reservoir particles in the network model and simulate flow of the reservoir particles in the network model under phase equilibrium conditions;
the selection module is configured to obtain capillary force according to the condition parameters of the migration of the reservoir particles and determine a throat through which the reservoir particles flow according to the capillary force;
a first update module configured to update the throat radius size to a first throat radius when the reservoir particle concentration is less than a preset concentration, if a single reservoir particle diameter is less than the throat diameter, or all reservoir particle volumes entering the throat are less than the throat volume;
a second updating module configured to update the throat radius size to a first throat radius if all the volume of the reservoir particles entering the throat is smaller than the throat volume when the reservoir particle concentration is greater than or equal to a preset concentration and the reservoir particle diameter is smaller than one-third of the throat diameter;
a determination module configured to determine the extent of reservoir particle migration according to a first predetermined rule based on the first throat radius.
12. An apparatus for determining the extent of migration of gas hydrate reservoir particles, comprising a memory and a processor, the memory having stored therein a computer program which, when executed by the processor, performs the steps of: a method of determining the extent of reservoir particle migration as claimed in any one of claims 1 to 10.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010926449.6A CN112069620B (en) | 2020-09-07 | 2020-09-07 | Method and device for determining migration degree of natural gas hydrate reservoir particles |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010926449.6A CN112069620B (en) | 2020-09-07 | 2020-09-07 | Method and device for determining migration degree of natural gas hydrate reservoir particles |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112069620A CN112069620A (en) | 2020-12-11 |
CN112069620B true CN112069620B (en) | 2023-01-13 |
Family
ID=73662762
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010926449.6A Active CN112069620B (en) | 2020-09-07 | 2020-09-07 | Method and device for determining migration degree of natural gas hydrate reservoir particles |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112069620B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113673183B (en) * | 2021-08-20 | 2022-09-23 | 中国石油大学(北京) | Reservoir development parameter determination method and device combining particle migration influence |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106703742A (en) * | 2015-07-20 | 2017-05-24 | 中国石油化工股份有限公司 | Method for matching particle size of elastic microspheres used for profile controlling and flooding of oil reservoir deep part with pore throat |
CN108798658A (en) * | 2018-06-05 | 2018-11-13 | 中国石油天然气股份有限公司 | Method and device for determining petroleum migration parameters |
-
2020
- 2020-09-07 CN CN202010926449.6A patent/CN112069620B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106703742A (en) * | 2015-07-20 | 2017-05-24 | 中国石油化工股份有限公司 | Method for matching particle size of elastic microspheres used for profile controlling and flooding of oil reservoir deep part with pore throat |
CN108798658A (en) * | 2018-06-05 | 2018-11-13 | 中国石油天然气股份有限公司 | Method and device for determining petroleum migration parameters |
Non-Patent Citations (2)
Title |
---|
基于网络模型的油藏优势通道形成微观机制;肖康等;《大庆石油地质与开发》;20170601(第03期);全文 * |
水驱砂岩油藏储层参数变化三维网络模拟;冯其红等;《油气地质与采收率》;20130725(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112069620A (en) | 2020-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wennberg et al. | Determination of the filtration coefficient and the transition time for water injection wells | |
CN105160134B (en) | Mixed medium simulation method and device for oil-gas flow in multiple media of tight reservoir | |
CN111236931B (en) | Method and system for generating unsteady well testing chart of gas reservoir highly-deviated well | |
Kang et al. | Numerical modeling of pore-scale phenomena during CO2 sequestration in oceanic sediments | |
CN112347675B (en) | Method for cooperatively regulating and controlling reservoir natural gas hydrate phase state by drilling fluid additive and temperature and pressure field | |
CN105205273B (en) | Method and device for simulating oil-gas relay flow in multiple media of tight reservoir | |
CN105260543A (en) | Multi-medium oil-gas flow simulation method and device based on double-hole model | |
CN112069620B (en) | Method and device for determining migration degree of natural gas hydrate reservoir particles | |
CN107621350A (en) | One kind simulation supercritical CO2The method of displacement natural gas flow | |
Persova et al. | Flow balancing in FEM modelling of multi-phase flow in porous media | |
CN112052634B (en) | Method and device for analyzing particle migration influence in water injection process of low-permeability reservoir | |
Johannsen | On the validity of the Boussinesq approximation for the Elder problem | |
CN108131127B (en) | Method and device for obtaining gas-oil ratio of production of foam oil type extra heavy oil field | |
Neale | Degrees of anisotropy for fluid flow and diffusion (electrical conduction) through anisotropic porous media | |
CN113673111B (en) | Polymer adsorption damage oil-gas layer modeling method, damage degree spatial-temporal evolution 4D quantification and intelligent diagnosis technology | |
CN113887067A (en) | Compact low-permeability reservoir fracturing well pattern gas flooding effect prediction method based on LSTM | |
CN115522909A (en) | Single well model size determination method and device based on numerical simulation | |
CN112685970A (en) | Quantitative characterization method and system for seepage interface of flow unit of water-drive reservoir | |
CN110135122B (en) | Low-permeability natural gas reservoir yield prediction method and system | |
Li et al. | Pore-Scale Modellings on the Impacts of Hydrate Distribution Morphology on Gas and Water Transport in Hydrate-Bearing Sediments | |
CN110863825A (en) | Method and device for distinguishing characteristic region | |
CN112014293A (en) | Method and device for representing hydrate reservoir seepage capability | |
CN105205282B (en) | The computational methods of THREE DIMENSIONAL GROUNDWATER FLOW amount | |
Yin et al. | GPU‐Based Computation of Formation Pressure for Multistage Hydraulically Fractured Horizontal Wells in Tight Oil and Gas Reservoirs | |
If et al. | Estimation of shape factors in fractured reservoirs |
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 |