CN114060013B - Interference well testing method for volcanic gas reservoir horizontal well - Google Patents
Interference well testing method for volcanic gas reservoir horizontal well Download PDFInfo
- Publication number
- CN114060013B CN114060013B CN202010732921.2A CN202010732921A CN114060013B CN 114060013 B CN114060013 B CN 114060013B CN 202010732921 A CN202010732921 A CN 202010732921A CN 114060013 B CN114060013 B CN 114060013B
- Authority
- CN
- China
- Prior art keywords
- well
- horizontal well
- pressure
- gas
- dimensionless
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 44
- 230000008859 change Effects 0.000 claims abstract description 33
- 238000013178 mathematical model Methods 0.000 claims abstract description 19
- 238000010206 sensitivity analysis Methods 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims description 35
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000004519 manufacturing process Methods 0.000 claims description 16
- 230000035699 permeability Effects 0.000 claims description 13
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 239000011435 rock Substances 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 4
- 230000005465 channeling Effects 0.000 claims description 3
- 238000004146 energy storage Methods 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 2
- 230000002452 interceptive effect Effects 0.000 claims 4
- 238000010998 test method Methods 0.000 abstract description 16
- 239000000243 solution Substances 0.000 description 13
- 238000010586 diagram Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measuring Fluid Pressure (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The application provides an interference well testing method for a volcanic gas reservoir horizontal well, which comprises the following steps: measuring basic parameters in a well group of vertical wells around a horizontal well serving as a center; b, changing the working system of the horizontal well, and measuring the pressure change data of the bottom of the horizontal well and the bottom hole of the vertical well around the horizontal well; c, establishing a bottom hole pressure analysis chart of the vertical well under the interference of the horizontal well according to a mathematical model; and D, establishing a bottom hole pressure parameter sensitivity analysis chart of the vertical well under the interference of the horizontal well. The technical scheme of the application effectively solves the problems of more recorded data and large workload of the interference well test method for the volcanic gas reservoir horizontal well in the prior art; and connectivity and parameter problems between the horizontal well and the vertical well of the volcanic gas reservoir can be identified and explained.
Description
Technical Field
The application relates to the technical field of a volcanic gas reservoir horizontal well exploitation method, in particular to an interference well testing method for a volcanic gas reservoir horizontal well.
Background
The well test method is a method for researching the oil reservoir characteristics through the dynamic reaction of the oil reservoir pressure. Well testing is performed by giving a series of "signals" to the reservoir, which are typically generated by changing the well's operating regime, and testing the reservoir's pressure for dynamic response under a given signal. The well test method mainly comprises two main types, namely an unstable well test method, a well test method for continuously testing the change process of bottom hole pressure along with time caused by the change of a well working system, and a well test method for researching the characteristic parameters of a well and a reservoir by the characteristics of the pressure change process, wherein the well test method mainly comprises pressure recovery well test, pressure drop well test, variable flow well test, interference well test, pulse well test and the like; the other type is a well testing method for testing the production (injection) capability of a well by testing the well under different stable working systems, namely a production capacity well testing method, and the well testing method mainly comprises back pressure well testing, isochronal well testing, one-point well testing and the like. Of the two types of methods, the interference well test and the pulse well test are effective methods for judging the connectivity between wells.
The technical condition of the interference test is that one exciting well is used for generating interference signals, and the exciting well needs to have enough or as large as possible yield (water injection quantity) to generate the interference signals with enough intensity; the observation well has enough stable production or shut-in time before testing to make its bottom hole pressure stable. The basic operation is that a high-precision pressure gauge is put into the bottom of an observation well before testing, a pressure trend line is measured for a period of time, then the working system of an exciting well is changed according to the well test design, and the pressure change of the observation well is recorded.
Pulsed well testing is a special form of disturbing well testing, the technical conditions and basic operation of which are the same as those of disturbing well testing. Unlike the interference test, the observation well data obtained by one test is required to generate at least 3 pressure extreme points, and all open and close cycles of the activated well must be the same. So that the technical requirements are more strict than those of the interference well test.
Volcanic gas reservoirs have the characteristics of thick production layer, high yield, large reserves and the like, and are one of important exploration targets at present. However, the lithology and lithology of the volcanic gas reservoir change rapidly, the reservoir space is complex, natural cracks develop, and the reservoir is generally a medium-low pore and low permeability reservoir, so that the effective development of the volcanic gas reservoir is difficult.
At present, domestic and foreign scholars have studied single well testing methods of a plurality of volcanic gas reservoirs and have achieved some results, however, due to the fact that reservoirs of the volcanic gas reservoirs are complex, inter-well relations are complex, inter-well connectivity is difficult to determine, as interference well testing is an effective method for determining inter-well connectivity, research on the interference well testing methods of the volcanic gas reservoirs does not have a better method, the existing interference well testing is a design method for observing pressure changes of a central well by sequentially and periodically changing working systems of adjacent wells, and therefore the workload is large, and standardization of field operation and effectiveness of data logging are difficult to guarantee; and connectivity and parameter problems between the horizontal well and the vertical well of the volcanic gas reservoir cannot be judged and interpreted.
Disclosure of Invention
The application mainly aims to provide an interference well testing method for a volcanic gas reservoir horizontal well, so as to solve the problems of more recorded data and large workload of the interference well testing method for the volcanic gas reservoir horizontal well in the prior art; and connectivity and parameter problems between the horizontal well and the vertical well of the volcanic gas reservoir can be identified and explained.
In order to achieve the above object, the present application provides an interference well testing method for a volcanic gas reservoir horizontal well, comprising: measuring basic parameters in a well group of vertical wells around a horizontal well serving as a center; b, changing the working system of the horizontal well, and measuring the bottom hole pressure change data of the horizontal well and the vertical wells around the horizontal well; c, establishing a bottom hole pressure analysis chart of the vertical well under the interference of the horizontal well according to a mathematical model; and D, establishing a bottom hole pressure parameter sensitivity analysis chart of the vertical well under the interference of the horizontal well.
Further, the step C further includes: setting initial parameters of the mathematical model according to geological data of logging and well completion, calculating a change value of bottom hole pressure of the vertical well along with production time according to the initial parameters, and drawing a bottom hole pressure analysis chart of the vertical well.
Further, the initial parameters are directly acquired through the geological data, including coordinates of the horizontal well and coordinates of surrounding vertical wells.
Further, the mathematical model includes:
in the method, in the process of the application,
wherein, the length of the dimensionless horizontal well
Dimensionless distance in a fracture
Storage capacity ratio
Coefficient of cross flow
Dimensionless gas reservoir thickness
Crack z-direction dimensionless distance
Dimensionless coordinates of crack horizontal midpoint and vertical well intersection point
Dimensionless coordinates
L is the length of the horizontal well, m; l is the effective well radius, m; k (K) 0 A Bessel function for the 0 th order second class imaginary volume; r is distance, m; u is Laplace variable; omega is the elastic energy storage ratio; phi (phi) f Is the porosity in the fracture system, dimensionless; mu (mu) fi In a fracture systemIs a raw gas viscosity of mPas; c (C) tfi Is the comprehensive compression coefficient in a fracture system, MPa -1 ；φ m Is the porosity in the matrix system, dimensionless; mu (mu) mi mPas, the viscosity of the original gas in the matrix system; c (C) tmi Is the comprehensive compression coefficient in the matrix system, MPa -1 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the cross-flow coefficient; k (k) m Is the matrix system permeability, mD; r is R m A circular boundary radius, m; k (k) f Is the permeability of the fracture system, mD; i 0 、I 1 The virtual volume Bessel functions of the first type are respectively 0 th order and 1 st order; h is the reservoir thickness, m; z is a coordinate value in the z direction of the crack; z w Coordinate values of the horizontal midpoint of the crack and the intersection point of the vertical well; y is 1 Is the coordinate in the y direction of the horizontal well.
Further, the mathematical model further comprises establishing:
basic rock mass system model
ψ mD (r mD ,0)＝0
ψ mD |r mD ＝1＝ψ fD
Point source model of natural fracture system
ψ fD (r D ,z D ,0)＝0
Wherein, the pseudo-pressure and dimensionless variable define:
pseudo pressure
Dimensionless matrix pseudo-pressure
Dimensionless fracture pressure simulation
Dimensionless distance in matrix
Dimensionless time
Dimensionless opening section height
Based on t D Laplace transformation is carried out on the bedrock block system model and the natural fracture system point source model, and then a separation variable method is utilized to solve the obtained point source solution into the model
Wherein n is a discrete number, n=1, 2, infinity;
the point source solution is integrated along the horizontal section of the horizontal well and multiplied by the source intensity to obtain the point source solution (x D ,y D ,z D ) The pressure of (2) is:
wherein, psi is the true gas pseudo pressure, MPa; p is the gas pressure, MPa; mu (mu) sc The viscosity of the gas in a standard state is mPas; z is Z sc Is a gas deviation factor in a standard state, and is dimensionless; p is p sc Is ground standard pressure, MPa; mu is the viscosity of the gas, mPas; z is a gas deviation factor, dimensionless; k (k) fr Is the formation permeability, mD; t (T) sc Is the standard temperature of the stratum, K; q sc For gas well production under standard conditions, 10 4 m 3 /d; t is the temperature of the gas layer, K; psi phi type i The pressure is the original pseudo pressure of the real gas, and is MPa; psi phi type m Simulating pressure for real gas in a matrix system, and MPa; psi phi type f Simulating pressure for real gas in a crack system, and simulating MPa; r is (r) m Distance, m; r is R m The radius of the columnar matrix is m; k (k) fz The permeability of the fracture in the z direction, mD; t is time, h; epsilon is the open segment height, m.
Further, in the step a, the basic parameters are directly obtained through geological data, including horizontal well length, effective well radius, comprehensive compression coefficient, gas viscosity, gas deviation factor, gas layer temperature and inter-well distance.
Further, in the step D, a vertical well bottom hole pressure parameter sensitivity analysis chart under horizontal well interference is established, where the parameter sensitivity analysis chart includes a curve change chart when the fracture storage ratio is changed, a curve change chart when the channeling coefficient is changed, and a curve change chart when the vertical well abscissa position is changed.
Further, in the step B, changing the working regime of the horizontal well specifically includes changing the production of the horizontal well.
Further, altering the production of the horizontal well includes: 1. closing the horizontal well and surrounding vertical wells; 2. and after the pressure is restored to be stable, the horizontal well is opened first and then closed, the surrounding vertical wells are closed all the time, and the bottom hole pressure of the horizontal well and the bottom hole pressure of the surrounding vertical wells are monitored.
Further, in the step B, the bottom hole pressure change data of the horizontal well is measured by a pressure gauge, and the running position of the pressure gauge is a horizontal well deflecting position.
Further, in the step B, the bottom hole pressure change data of the vertical well is measured by a pressure gauge, and the position of the pressure gauge is a deep position in the liquid surface of the vertical well.
By applying the technical scheme of the application, the working system of the horizontal well is changed, the pressure change data of the bottom hole of the well group of the horizontal well and the vertical well around the horizontal well which are centered by the horizontal well are measured, a mathematical model is built, and a bottom hole pressure analysis chart of the vertical well under the interference of the horizontal well and a bottom hole pressure parameter sensitivity analysis chart of the vertical well under the interference of the horizontal well are built. The technical scheme of the application effectively solves the problems of more recorded data and large workload of the interference well test method for the volcanic gas reservoir horizontal well in the prior art; and connectivity and parameter problems between the horizontal well and the vertical well of the volcanic gas reservoir can be identified and explained.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application. In the drawings:
FIG. 1 shows a schematic diagram of a horizontal well manometer run-in position for an embodiment of an interference test method for a volcanic gas reservoir horizontal well according to the present application;
FIG. 2 shows a schematic view of an observation well pressure gauge run-in position of an embodiment of an interference well test method for a volcanic gas reservoir horizontal well in accordance with the present application;
FIG. 3 shows a schematic diagram of a physical model of an embodiment of an interference well test method for a volcanic gas reservoir horizontal well according to the present application;
FIG. 4 shows a schematic analysis of the bottom hole pressure of a vertical well in accordance with an embodiment of the disturbance test method for a volcanic gas reservoir horizontal well according to the present application;
FIG. 5 shows a graphical representation of the change in fracture capacitance ratio as it changes for an embodiment of an interference well test method for a volcanic gas reservoir horizontal well in accordance with the present application;
FIG. 6 shows a schematic graph of a change in a cross-flow coefficient as it changes for an embodiment of a disturbance test method for a volcanic gas reservoir horizontal well according to the present application;
fig. 7 shows a schematic diagram of a change in curve when the vertical axis position of an embodiment of the disturbance testing method for a volcanic gas reservoir horizontal well according to the present application is changed.
Wherein the above figures include the following reference numerals:
10. a horizontal well; 20. a pressure gauge; 30. a vertical well; 40. the liquid level of the vertical well.
Detailed Description
It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other. The application will be described in detail below with reference to the drawings in connection with embodiments.
It should be noted that the following detailed description is illustrative and is intended to provide further explanation of the application. 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.
Spatially relative terms, such as "above … …," "above … …," "upper surface at … …," "above," and the like, may be used herein for ease of description to describe one device or feature's spatial location relative to another device or feature as illustrated in the figures. It will be understood that the spatially relative terms are intended to encompass different orientations in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements described as "above" or "over" other devices or structures would then be oriented "below" or "beneath" the other devices or structures. Thus, the exemplary term "above … …" may include both orientations of "above … …" and "below … …". The device may also be positioned in other different ways (rotated 90 degrees or at other orientations) and the spatially relative descriptors used herein interpreted accordingly.
Exemplary embodiments according to the present application will now be described in more detail with reference to the accompanying drawings. These exemplary embodiments may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. It should be appreciated that these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of these exemplary embodiments to those skilled in the art, that in the drawings, thicknesses of layers and regions are exaggerated for clarity, and identical reference numerals are used to denote identical devices, and thus descriptions thereof will be omitted.
As shown in fig. 1 to 7, the interference well testing method for the volcanic gas reservoir horizontal well of the present embodiment includes: a measures a base parameter within a well group of vertical wells centered around the horizontal well 10. And B, changing the working system of the horizontal well, and measuring the bottom hole pressure change data of the horizontal well and the vertical wells around the horizontal well. And C, establishing a bottom hole pressure analysis chart of the vertical well under the interference of the horizontal well according to the mathematical model. And D, establishing a bottom hole pressure parameter sensitivity analysis chart of the vertical well under the interference of the horizontal well.
By applying the technical scheme of the embodiment, the working system of the horizontal well is changed, the pressure change data of the bottom hole of the well group of the horizontal well and the vertical well around the horizontal well which are centered by the horizontal well are measured, a mathematical model is built, and a bottom hole pressure analysis chart of the vertical well under the interference of the horizontal well and a bottom hole pressure parameter sensitivity analysis chart of the vertical well under the interference of the horizontal well are built. The technical scheme of the embodiment effectively solves the problems of more recorded data and large workload of the interference well test method for the volcanic gas reservoir horizontal well in the prior art; and connectivity and parameter problems between the horizontal well and the vertical well of the volcanic gas reservoir can be identified and explained.
The well group includes a horizontal well and a vertical well surrounding the horizontal well, and the above-mentioned parameters are easily valued, for example, the measurement of pressure is relatively easy.
As shown in fig. 1 and fig. 4, in the technical solution of this embodiment, step C further includes: setting initial parameters of a mathematical model according to geological data of logging and well completion, calculating a change value of bottom hole pressure of the vertical well along with production time according to the initial parameters, and drawing a bottom hole pressure analysis chart of the vertical well. The steps provide a numerical basis for the subsequent mathematical model, so that the accuracy of the interference well test method is ensured. The change value of the vertical well pressure along with the production is the bottom-hole pressure drop of the vertical well, including the bottom-hole pressure drop of the vertical well caused during the production of the horizontal well.
In the technical solution of this embodiment, the initial parameters are directly obtained through geological data, including coordinates of the horizontal well and coordinates of surrounding vertical wells. The acquisition of the geological data is carried out according to the requirement of the mathematical model, and the acquisition of the initial parameters is easier.
In the technical solution of the present embodiment, the mathematical model includes:
in the method, in the process of the application,
wherein, the length of the dimensionless horizontal well
Dimensionless distance in a fracture
Storage capacity ratio
Coefficient of cross flow
Dimensionless gas reservoir thickness
Crack z-direction dimensionless distance
Dimensionless coordinates of crack horizontal midpoint and vertical well intersection point
Dimensionless coordinates
L is the length of the horizontal well, m; l is the effective well radius, m; k (K) 0 A Bessel function for the 0 th order second class imaginary volume; r is distance, m; u is Laplace variable; omega is the elastic energy storage ratio; phi (phi) f Is the porosity in the fracture system, dimensionless; mu (mu) fi The viscosity of the original gas in the fracture system is mPa.s; c (C) tfi Is the comprehensive compression coefficient in a fracture system, MPa -1 ；φ m Is the porosity in the matrix system, dimensionless; mu (mu) mi mPas, the viscosity of the original gas in the matrix system; c (C) tmi Is the comprehensive compression coefficient in the matrix system, MPa -1 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the cross-flow coefficient; k (k) m Is the matrix system permeability, mD; r is R m A circular boundary radius, m; k (k) f Is the permeability of the fracture system, mD; i 0 、I 1 The virtual volume Bessel functions of the first type are respectively 0 th order and 1 st order; h is the reservoir thicknessM; z is a coordinate value in the z direction of the crack; z w Coordinate values of the horizontal midpoint of the crack and the intersection point of the vertical well; y is 1 Is the coordinate in the y direction of the horizontal well. The mathematical model is established by using Laplace transformation method and separation variable method.
In the technical solution of this embodiment, the mathematical model further includes establishing: a basic rock mass system model and a natural fracture system point source model.
Basic rock mass system model
ψ mD (r mD ,0)＝0
ψ mD |r mD ＝1＝ψ fD
Point source model of natural fracture system
ψ fD (r D ,z D ,0)＝0
Wherein, the pseudo-pressure and dimensionless variable define:
pseudo pressure
Dimensionless matrix pseudo-pressure
Dimensionless fracture pressure simulation
Dimensionless distance in matrix
Dimensionless time
Dimensionless opening section height
Based on t D Laplace transformation is carried out on the bedrock block system model and the natural fracture system point source model, and then a separation variable method is utilized to solve the obtained point source solution into the model
Wherein n is a discrete number, n=1, 2, infinity;
the point source solution is integrated along the horizontal section of the horizontal well and multiplied by the source intensity to obtain the point source solution (x D ,y D ,z D ) The pressure of (2) is:
wherein, psi is the true gas pseudo pressure, MPa; p is the gas pressure, MPa; mu (mu) sc The viscosity of the gas in a standard state is mPas; z is Z sc Is a gas deviation factor in a standard state, and is dimensionless; p is p sc Is ground standard pressure, MPa; mu is the viscosity of the gas, mPas; z is a gas deviation factor, dimensionless; k (k) fr Is the formation permeability, mD; t (T) sc Is the standard temperature of the stratum, K; q sc For gas well production under standard conditions, 10 4 m 3 /d; t is the temperature of the gas layer, K; psi phi type i The pressure is the original pseudo pressure of the real gas, and is MPa; psi phi type m Simulating pressure for real gas in a matrix system, and MPa; psi phi type f Simulating pressure for real gas in a crack system, and simulating MPa; r is (r) m Distance, m; r is R m The radius of the columnar matrix is m; k (k) fz The permeability of the fracture in the z direction, mD; t is time, h; epsilon is the open segment height, m.
In the technical scheme of the embodiment, the basic parameters are directly acquired through geological data, and include horizontal well length, effective well radius, comprehensive compression coefficient, gas viscosity, gas deviation factor, gas layer temperature and inter-well distance. The acquisition of the geological data is carried out according to the requirement of the mathematical model, and the acquisition of the basic parameters is easier.
As shown in fig. 5 to 6, in the technical solution of the present embodiment, in step C, a vertical well bottom hole pressure parameter sensitivity analysis chart under horizontal well interference is established, where the parameter sensitivity analysis chart includes a curve change chart when the fracture storage ratio is changed, a curve change chart when the fluid channeling coefficient is changed, and a curve change chart when the vertical well abscissa position is changed.
In the technical solution of this embodiment, in step B, changing the working regimen of the horizontal well specifically includes changing the output of the horizontal well. The working system is easy to change, and the requirement of experiments can be met.
In the technical solution of this embodiment, changing the output of the horizontal well includes: 1. firstly, closing the horizontal well and surrounding vertical wells; 2. after the pressure is restored to be stable, the horizontal well is started first and then closed, the surrounding vertical wells are closed all the time, and the bottom hole pressure of the horizontal well and the bottom hole pressure of the surrounding vertical wells are monitored. The operation is convenient and the implementation is easy.
As shown in fig. 1, in the embodiment, in the step B, the bottom hole pressure change data of the horizontal well 10 is measured by the pressure gauge 20, and the gauge run-in position is the horizontal well deflecting position.
As shown in fig. 2, in the embodiment, in step B, the bottom hole pressure change data of the vertical well 30 is measured by a pressure gauge, and the position where the pressure gauge is lowered is a deep position in the vertical well surface 40.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present application. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
It should be noted that the terms "first," "second," and the like in the description and the claims of the present application and the above figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that embodiments of the application described herein may be capable of being practiced otherwise than as specifically illustrated and described. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
The above description is only of the preferred embodiments of the present application and is not intended to limit the present application, but various modifications and variations can be made to the present application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.
Claims (10)
1. An interference well testing method for a volcanic gas reservoir horizontal well, comprising:
measuring basic parameters in a well group of vertical wells around a horizontal well serving as a center;
b, changing the working system of the horizontal well, and measuring the bottom hole pressure change data of the horizontal well and the vertical wells around the horizontal well;
c, establishing a bottom hole pressure analysis chart of the vertical well under the interference of the horizontal well according to a mathematical model;
d, establishing a bottom hole pressure parameter sensitivity analysis chart of the vertical well under the interference of the horizontal well;
the mathematical model comprises:
in the method, in the process of the application,
wherein, the length of the dimensionless horizontal well
Dimensionless distance in a fracture
Storage capacity ratio
Coefficient of cross flow
Dimensionless gas reservoir thickness
Crack z-direction dimensionless distance
Dimensionless coordinates of crack horizontal midpoint and vertical well intersection point
Dimensionless coordinates
L is the length of the horizontal well, m; l is the effective well radius, m; k (K) 0 A Bessel function for the 0 th order second class imaginary volume; r is distance, m; u is Laplace variable; omega is the elastic energy storage ratio; phi (phi) f Is the porosity in the fracture system, dimensionless; mu (mu) fi The viscosity of the original gas in the fracture system is mPa.s; c (C) tfi Is the comprehensive compression coefficient in a fracture system, MPa -1 ；φ m Is the porosity in the matrix system, dimensionless; mu (mu) mi mPas, the viscosity of the original gas in the matrix system; c (C) tmi Is the comprehensive compression coefficient in the matrix system, MPa -1 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the cross-flow coefficient; k (k) m Is the matrix system permeability, mD; r is R m A circular boundary radius, m; k (k) f Is the permeability of the fracture system, mD; i 0 、I 1 The virtual volume Bessel functions of the first type are respectively 0 th order and 1 st order; h is the reservoir thickness, m; z is a coordinate value in the z direction of the crack; z w Coordinate values of the horizontal midpoint of the crack and the intersection point of the vertical well; y is 1 Is the coordinate in the y direction of the horizontal well.
2. The method of interfering with a horizontal well of a volcanic gas reservoir as set forth in claim 1, wherein said step C further comprises: setting initial parameters of the mathematical model according to geological data of logging and well completion, calculating a change value of bottom hole pressure of the vertical well along with production time according to the initial parameters, and drawing a bottom hole pressure analysis chart of the vertical well.
3. The method of claim 2, wherein the initial parameters are obtained directly from the geological data, including the coordinates of the horizontal well and the coordinates of the surrounding vertical wells.
4. The method of interfering with a well test for a horizontal well of a volcanic gas reservoir of claim 2, wherein the mathematical model further comprises establishing:
basic rock mass system model
ψ mD (r mD ,0)＝0
Point source model of natural fracture system
ψ fD (r D ,z D ,0)＝0
Wherein, the pseudo-pressure and dimensionless variable define:
pseudo pressure
Dimensionless matrix pseudo-pressure
Dimensionless fracture pressure simulation
Dimensionless distance in matrix
Dimensionless time
Dimensionless opening section height
Based on t D Laplace transformation is carried out on the bedrock block system model and the natural fracture system point source model, and then a separation variable method is utilized to solve the obtained point source solution into the model
Wherein n is a discrete number, n=1, 2, infinity;
the point source solution is integrated along the horizontal section of the horizontal well and multiplied by the source intensity to obtain the point source solution (x D ,y D ,z D ) The pressure of (2) is:
wherein, psi is the true gas pseudo pressure, MPa; p is the gas pressure, MPa; mu (mu) sc The viscosity of the gas in a standard state is mPas; z is Z sc Is a gas deviation factor in a standard state, and is dimensionless; p is p sc Is ground standard pressure, MPa; mu is the viscosity of the gas, mPas; z is a gas deviation factor, dimensionless; k (k) fr Is the formation permeability, mD; t (T) sc Is the standard temperature of the stratum, K; q sc For gas well production under standard conditions, 10 4 m 3 /d; t is the temperature of the gas layer, K; psi phi type i The pressure is the original pseudo pressure of the real gas, and is MPa; psi phi type m Simulating pressure for real gas in a matrix system, and MPa; psi phi type f Simulating pressure for real gas in a crack system, and simulating MPa; r is (r) m Distance, m; r is R m The radius of the columnar matrix is m; k (k) fz The permeability of the fracture in the z direction, mD; t is time, h; epsilon is the open segment height, m.
5. The method of claim 1, wherein in step a, the base parameters are obtained directly from geological data, including horizontal well length, effective well radius, integrated compression factor, gas viscosity, gas deviation factor, gas formation temperature, and inter-well distance.
6. The method of claim 1, wherein in step D, a vertical well bottom hole pressure parameter sensitivity analysis map is created for the horizontal well, the parameter sensitivity analysis map including a graph of change in fracture storage ratio, a graph of change in fluid channeling coefficient, and a graph of change in vertical well abscissa position.
7. The method of claim 1, wherein in step B, changing the operating regime of the horizontal well comprises changing the production of the horizontal well.
8. The method of interfering with a well test for a horizontal well of a volcanic gas reservoir of claim 7, wherein altering the production of the horizontal well comprises: 1. closing the horizontal well and surrounding vertical wells; 2. and after the pressure is restored to be stable, the horizontal well is opened first and then closed, the surrounding vertical wells are closed all the time, and the bottom hole pressure of the horizontal well and the bottom hole pressure of the surrounding vertical wells are monitored.
9. The method of claim 1, wherein in step B, the bottom hole pressure change data of the horizontal well is measured by a pressure gauge, and the pressure gauge run-in position is a horizontal well kick-off position.
10. The method for interfering with the well logging of a horizontal well in a volcanic gas reservoir according to claim 1, wherein in the step B, the bottom hole pressure variation data of the vertical well is measured by a pressure gauge, and the gauge run-in position is a deep position in the liquid surface of the vertical well.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010732921.2A CN114060013B (en) | 2020-07-27 | 2020-07-27 | Interference well testing method for volcanic gas reservoir horizontal well |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010732921.2A CN114060013B (en) | 2020-07-27 | 2020-07-27 | Interference well testing method for volcanic gas reservoir horizontal well |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114060013A CN114060013A (en) | 2022-02-18 |
| CN114060013B true CN114060013B (en) | 2023-09-26 |
Family
ID=80226538
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010732921.2A Active CN114060013B (en) | 2020-07-27 | 2020-07-27 | Interference well testing method for volcanic gas reservoir horizontal well |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114060013B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102146788A (en) * | 2011-03-30 | 2011-08-10 | 中国石油化工股份有限公司 | Method for enhancing recovery ratio through three-dimensional streamline regulation of water drive reservoir |
| CN105019894A (en) * | 2015-07-29 | 2015-11-04 | 长江大学 | Method and system for establishing multilayer oil pool inter-well connectivity model |
| CN105626023A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Well test determination method for vertical fracturing fracture azimuth of low-permeability oil reservoir |
| CN108875122A (en) * | 2018-04-25 | 2018-11-23 | 杭州迅美科技有限公司 | The artificial intelligence approach and system of geologic parameter are calculated using well logging data |
| CN111255442A (en) * | 2020-01-14 | 2020-06-09 | 大庆油田有限责任公司 | Method for evaluating fracturing fracture by using interference well testing theory |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA3065684A1 (en) * | 2018-12-22 | 2020-06-22 | Schlumberger Canada Limited | Dynamic reservoir characterization |
-
2020
- 2020-07-27 CN CN202010732921.2A patent/CN114060013B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102146788A (en) * | 2011-03-30 | 2011-08-10 | 中国石油化工股份有限公司 | Method for enhancing recovery ratio through three-dimensional streamline regulation of water drive reservoir |
| CN105626023A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Well test determination method for vertical fracturing fracture azimuth of low-permeability oil reservoir |
| CN105019894A (en) * | 2015-07-29 | 2015-11-04 | 长江大学 | Method and system for establishing multilayer oil pool inter-well connectivity model |
| CN108875122A (en) * | 2018-04-25 | 2018-11-23 | 杭州迅美科技有限公司 | The artificial intelligence approach and system of geologic parameter are calculated using well logging data |
| CN111255442A (en) * | 2020-01-14 | 2020-06-09 | 大庆油田有限责任公司 | Method for evaluating fracturing fracture by using interference well testing theory |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114060013A (en) | 2022-02-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10436027B2 (en) | Method of geometric evaluation of hydraulic fractures | |
| CN103982179B (en) | A kind of paleopressure quantitative inversion detection method of oil reservoir | |
| CN106991236B (en) | Repeated fracturing well selection layer selection method based on four-dimensional ground stress dynamic change | |
| CN103323352B (en) | Natural gas hydrate deposit dynamic triaxial mechanic-acoustic-electrical synchronous test experimental device and method | |
| US4665984A (en) | Method of measuring crustal stress by hydraulic fracture based on analysis of crack growth in rock | |
| CA3221947A1 (en) | Mapping of fracture geometries in a multi-well stimulation process | |
| CN105653815A (en) | Reservoir fluid distribution quantitative interpretation method based on rock physical model theory | |
| CN103257151B (en) | In a kind of quantitative evaluation oil and gas secondary migration process, pore throat employs the method for rule | |
| Zhang et al. | Pressure–dependent fracture permeability of marine shales in the Northeast Yunnan area, Southern China | |
| WO2015095557A1 (en) | Method for determining hydraulic fracture orientation and dimension | |
| CN101446196A (en) | Well test analysis method and device of treble medium oil pool branch horizontal well | |
| CN112746837A (en) | Shale oil reservoir exploration data acquisition system and method based on distributed optical fiber sensing | |
| Gong et al. | Quantitative prediction of sub-seismic faults and their impact on waterflood performance: Bozhong 34 oilfield case study | |
| CN114166726B (en) | Core permeability tensor sensitivity damage testing equipment and evaluation method | |
| CN111751885A (en) | Shale gas volume fracturing micro-seismic monitoring method | |
| CN106769684B (en) | Shale gas diffusivity test macro | |
| CN113821956B (en) | Evaluation method for disturbance quantity of current geostress structure of deep shale reservoir | |
| CN114000872B (en) | Natural gas hydrate horizontal well stratified mining process soil layer deformation testing device | |
| CN114060013B (en) | Interference well testing method for volcanic gas reservoir horizontal well | |
| CN111706318B (en) | Method for determining distribution condition of residual oil of hypotonic reservoir | |
| CN115270536A (en) | Method for evaluating reconstruction gas storage reservoir containing water traps in fractured basin | |
| CN112068197A (en) | Fracture cracking scale description method based on seismic wave dynamic parameters | |
| Akin | Design and analysis of multi-well interference tests | |
| Cheng et al. | A New Numerical Simulation technology of Multistage Fracturing in Horizontal Well | |
| Agrawal | Pore-network Modeling and Analysis of Laboratory Induced Hydraulic Fractures |
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 |