CN115822532B - Method, device, equipment and medium for determining minimum miscible pressure of carbon dioxide-crude oil - Google Patents
Method, device, equipment and medium for determining minimum miscible pressure of carbon dioxide-crude oil Download PDFInfo
- Publication number
- CN115822532B CN115822532B CN202211441721.7A CN202211441721A CN115822532B CN 115822532 B CN115822532 B CN 115822532B CN 202211441721 A CN202211441721 A CN 202211441721A CN 115822532 B CN115822532 B CN 115822532B
- Authority
- CN
- China
- Prior art keywords
- miscible
- data points
- piezoelectric transducer
- pressure
- fitting
- 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
- 238000000034 method Methods 0.000 title claims abstract description 65
- 239000010779 crude oil Substances 0.000 title claims abstract description 49
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 30
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 30
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims abstract description 46
- 229910002092 carbon dioxide Inorganic materials 0.000 claims abstract description 23
- 239000001569 carbon dioxide Substances 0.000 claims abstract description 20
- 238000006073 displacement reaction Methods 0.000 claims description 44
- 239000000919 ceramic Substances 0.000 claims description 25
- 230000008859 change Effects 0.000 claims description 21
- 238000002347 injection Methods 0.000 claims description 18
- 239000007924 injection Substances 0.000 claims description 18
- 230000006870 function Effects 0.000 claims description 17
- 239000012530 fluid Substances 0.000 claims description 16
- 238000002474 experimental method Methods 0.000 claims description 14
- 230000000694 effects Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 230000001965 increasing effect Effects 0.000 claims description 9
- 238000004519 manufacturing process Methods 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 4
- 239000002184 metal Substances 0.000 claims description 4
- 229910052751 metal Inorganic materials 0.000 claims description 4
- 230000007704 transition Effects 0.000 abstract description 2
- 239000002609 medium Substances 0.000 description 20
- 239000003921 oil Substances 0.000 description 18
- 230000008569 process Effects 0.000 description 10
- 230000005684 electric field Effects 0.000 description 9
- 239000007789 gas Substances 0.000 description 9
- 210000005239 tubule Anatomy 0.000 description 7
- 229920001971 elastomer Polymers 0.000 description 6
- 230000035699 permeability Effects 0.000 description 6
- 239000000243 solution Substances 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 239000008398 formation water Substances 0.000 description 5
- 239000011435 rock Substances 0.000 description 4
- UHOVQNZJYSORNB-UHFFFAOYSA-N Benzene Chemical compound C1=CC=CC=C1 UHOVQNZJYSORNB-UHFFFAOYSA-N 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000004568 cement Substances 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 239000010453 quartz Substances 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 239000006004 Quartz sand Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000003756 stirring Methods 0.000 description 2
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- PIICEJLVQHRZGT-UHFFFAOYSA-N Ethylenediamine Chemical compound NCCN PIICEJLVQHRZGT-UHFFFAOYSA-N 0.000 description 1
- 241000629264 Halicreas minimum Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 239000003822 epoxy resin Substances 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 239000003546 flue gas Substances 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 229920000647 polyepoxide Polymers 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/70—Combining sequestration of CO2 and exploitation of hydrocarbons by injecting CO2 or carbonated water in oil wells
Abstract
The invention provides a method, a device, equipment and a medium for determining the minimum miscible pressure of carbon dioxide-crude oil, which find the state transition point of the medium in a miscible state and a non-miscible state by utilizing different propagation characteristics of sound waves in the miscible medium and the non-miscible medium, thereby establishing the relationship between the sound waves and the pressure, finding a path for determining the minimum miscible pressure of the carbon dioxide and the crude oil, and solving the technical problem in the field from a new angle. The data obtained by the method is accurate, and the error range is only 1% -2% by comparing with the actual minimum miscible pressure, so that the accuracy is greatly improved. Meanwhile, the method, the device, the equipment and the medium can be popularized and applied to the determination of the mixed phase pressure of different media.
Description
Technical Field
The invention relates to the technical field of oilfield analysis, in particular to a method, a device, equipment and a medium for determining the minimum miscible pressure of carbon dioxide-crude oil.
Background
The carbon dioxide flooding technology has become an important method for improving the recovery ratio of the low permeability reservoir by the characteristics of high injectability, remarkable oil increasing effect, recycling and the like. According to the oil displacement mechanism, the technology is divided into carbon dioxide miscible flooding and non-miscible flooding technologies.
In the prior art, when determining the minimum miscible pressure, the common devices and methods for the minimum miscible pressure are a tubule experiment method, a bubble rising instrument method, an interfacial tension vanishing method and the like. While the tubule test method simulates the porous medium of the reservoir by filling glass beads or quartz sand in the long tubule, the size of the porous medium of the tubule is not consistent with the actual reservoir due to the fact that no cementing agent is used in the filling process, thus resulting in inaccuracy of the simulation. For example, see micro-nano confined space crude oil-natural gas minimum miscible pressure prediction method [1]. The mixed phase generating positions in the test process of the bubble rising instrument method and the interfacial tension method are respectively arranged in the glass tubule and the hanging drop kettle, so that the porous medium factors are not considered in the two methods. For example, see Minimum miscibility pressure determination in confined nanoporesconsidering pore size distribution of tight/scale formations [2], whereas the porous media size has a large effect on the miscible pressure of carbon dioxide and crude oil, so that the miscible pressure measured by the prior art method has a large error compared to the actual minimum miscible pressure, and the error range exceeds 15%. Therefore, it is desirable to find a new method that can accurately determine the minimum miscible pressure.
Reference document:
[1] weikrg, zhong Mengying, zhao Jinzhou, et al, micro-nano limited space crude oil-natural gas minimum miscible pressure prediction method [ J ]. Petroleum journal, 2022,9 (05): 1-10;
[2]Sun H,Li H.Minimum miscibility pressure determination in confined nanopores considering pore size distribution of tight/shale formations[J].Fuel,2021,286:119-128.。
disclosure of Invention
According to the proposal, only the related factors such as gas impurities, oil reservoir temperature, crude oil components and the like which influence the miscible pressure are considered, and the technical problem of application of a carbon dioxide flooding technology in a shale oil reservoir is not considered, so that the method, the device, the equipment and the medium for determining the minimum miscible pressure of the carbon dioxide-crude oil are provided.
The invention adopts the following technical means:
a method of determining a minimum miscible pressure of carbon dioxide-crude oil, the method comprising:
firstly, manufacturing a plurality of standby core models according to physical parameters of a target reservoir to be measured; any one of the standby core models is taken and placed into a core holder;
secondly, respectively implanting an inlet end piezoelectric transducer and an outlet end piezoelectric transducer at the centers of the left end face and the right end face of a core model placed in a core holder; the inlet end piezoelectric transducer is used for transmitting sound wave signals, and the outlet end piezoelectric transducer is used for receiving the sound wave signals transmitted by the core model arranged in the core holder and converting the sound wave signals into voltage amplitude signals; the average voltage value U of the voltage amplitude signal is collected and recorded by an industrial computer; the average voltage value U represents the difference between the highest voltage and the lowest voltage in the period from the beginning of the experiment to the breakthrough of the gas to be tested in the piezoelectric transducer at the outlet end;
thirdly, placing the core holder into a carbon dioxide-crude oil displacement detection device, performing a displacement experiment from carbon dioxide-crude oil non-miscible phase to miscible phase under a first injection end pressure P1, and recording an obtained first voltage amplitude signal U through the industrial computer; then, optionally replacing the core model in the core holder by one of the standby core models, increasing the injection end pressures P1 to Pn step by step, repeatedly executing the second step to the third step to obtain a series of voltage amplitude signals U1 to Un, and establishing a P-U coordinate system according to the pressures;
fourth step, define qn= (Un +1 -Un)/Un is the rate of change of voltage, n > 1, calculating the rate of change of voltage of the series of voltage amplitude signals U1 to Un obtained via the third step item by item, creating a voltage rate data set from the rate of change of voltage of said U1 to Un;
fifthly, when Qn is defined to be larger than a first set value, the data points in the corresponding P-U coordinate system are non-miscible data points; when Qn is defined to be smaller than a second set value, the data points in the corresponding P-U coordinate system are mixed phase data points; in a preferred embodiment, the first set value is set to 30% and the second set value is set to 5%, and it is understood that in other embodiments, the first and second set values are determined according to the Qn change range.
A sixth step of performing data fitting on the non-miscible data points and the miscible data points obtained in the fifth step in the P-U coordinate system to obtain two fitting curves; and the P value corresponding to the intersection of the two fitting curves is the minimum miscible pressure of the carbon dioxide and the crude oil.
Further, the data fitting in the sixth step is performed according to the following paths:
firstly, adopting a y=ax+b function model to respectively perform first fitting on non-miscible data points and miscible data points, and if absolute values of correlation coefficients R of a non-miscible data point curve and a miscible data point curve obtained after the first fitting are both larger than 0.9, the data change trend is linear, and the two obtained curves are effective data fitting curves; if the absolute values of the correlation coefficients R of the non-miscible data point curve and the miscible data point curve obtained after the first fitting are smaller than or equal to 0.9, the obtained two curves are invalid data fitting curves; if the absolute value of the correlation coefficient R of only one fitting curve is less than or equal to 0.9, then y=ax is used for the data point corresponding to the fitting curve 2 The second fitting is performed by a +bx+c function model, the parameters a, b, c in the function model being obtained by the least squares method.
Further, in the third step, the injection port pressures P1 to Pn are increased step by step with the third set value as a step. In a preferred embodiment, the third setting value is set to 5MPa in the present application, and in another embodiment, the third setting value is determined according to the Qn change width.
Further, the invention also comprises an apparatus for determining the minimum miscible pressure of carbon dioxide-crude oil, comprising a carbon dioxide-crude oil displacement test assembly, wherein the apparatus further comprises: the sound wave transmitting unit is used for transmitting sound waves at the displacement fluid inlet of the core model; the sound wave receiving unit is used for receiving sound waves at the displacement fluid outlet of the core model and converting the sound waves into voltage signals to be output; the central control unit has the following functions: the starting voltage signal is used for outputting a starting voltage signal for realizing sound wave emission to the sound wave emission unit; the voltage waveform signal is used for receiving the voltage waveform signal output by the sound wave receiving unit; the method comprises the steps of performing data processing on received voltage waveform signals to obtain an average voltage value U, and generating a P-U coordinate system under the control of a first built-in program after receiving input different displacement pressure values P; the method is used for distinguishing mixed phase data points and non-mixed phase data points in the P-U coordinate system under the control of a second built-in program; the method comprises the steps of respectively carrying out data fitting on screened mixed phase data points and non-mixed phase data points under the control of a third built-in program; after the effective data fitting curve is obtained, the displacement pressure value at the intersection point of the two fitting curves is determined.
Further, the acoustic wave transmitting unit includes an inlet-side piezoelectric transducer: the inlet end piezoelectric transducer is arranged at the center of one end of the core model in the core holder; the sound wave receiving unit comprises an outlet end piezoelectric transducer; the outlet end piezoelectric transducer is arranged at the center of the other end of the core model in the core holder and corresponds to the inlet end piezoelectric transducer.
Further, the inlet-side piezoelectric transducer and the outlet-side piezoelectric transducer each include: the piezoelectric transducer comprises a front cover plate, a rear cover plate, piezoelectric ceramics, an insulating layer, a high-strength stress rod, an acoustic matching layer and a pipeline, wherein the rear cover plate is arranged corresponding to the front cover plate and used for isolating external fluid, the piezoelectric ceramics is arranged at two sides of an internal cavity of the piezoelectric transducer and used for connecting an alternating current power supply, the insulating layer is connected with the piezoelectric ceramics, the high-strength stress rod is used for fixing the front cover plate and the rear cover plate and the piezoelectric ceramics, the acoustic matching layer is arranged at one end of the piezoelectric transducer and used for transmitting and absorbing sound waves, and the pipeline is arranged in an annular space between a metal shell of the piezoelectric transducer and the front cover plate and used for injecting fluid; the high-strength stress rod is fixed in the middle position of the piezoelectric transducer; the piezoelectric ceramic generates ultrasonic waves due to an inverse piezoelectric effect; the sound matching layer at one end of the piezoelectric transducer at the inlet end and the sound matching layer at one end of the piezoelectric transducer at the outlet end are placed in opposite directions.
Further, the present invention also comprises an apparatus for determining the minimum miscible pressure of carbon dioxide-crude oil, said apparatus comprising: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to enable the at least one processor to: outputting a starting voltage signal for realizing sound wave emission to the sound wave emission unit; receiving a voltage waveform signal output by the sound wave receiving unit; the method comprises the steps of performing data processing on received voltage waveform signals to obtain an average voltage value U, and generating a P-U coordinate system under the control of a first built-in program after receiving different displacement pressure values P input by an operator; under the control of a second built-in program, distinguishing mixed phase data points and non-mixed phase data points in the P-U coordinate system; under the control of a third built-in program, respectively performing data fitting on the screened mixed phase data points and the screened non-mixed phase data points; after the effective data fitting curve is obtained, the displacement pressure value at the intersection point of the two fitting curves is determined.
Further, the invention also comprises a medium for determining the minimum miscible pressure of carbon dioxide-crude oil, storing computer executable instructions configured for: outputting a starting voltage signal for realizing sound wave emission to the sound wave emission unit; receiving a voltage waveform signal output by the sound wave receiving unit; the method comprises the steps of performing data processing on received voltage waveform signals to obtain an average voltage value U, and generating a P-U coordinate system under the control of a first built-in program after receiving different displacement pressure values P input by an operator; under the control of a second built-in program, distinguishing mixed phase data points and non-mixed phase data points in the P-U coordinate system; under the control of a third built-in program, respectively performing data fitting on the screened mixed phase data points and the screened non-mixed phase data points; after the effective data fitting curve is obtained, the displacement pressure value at the intersection point of the two fitting curves is determined.
Compared with the prior art, the invention has the following advantages:
1. the invention provides a method, a device, equipment and a medium for determining the minimum miscible pressure of carbon dioxide-crude oil, creatively utilizes different propagation characteristics of sound waves in a miscible medium and a non-miscible medium, finds state transition points of the medium in a miscible state and a non-miscible state, thereby establishing the relationship between the sound waves and the pressure, finding a path for determining the minimum miscible pressure of the carbon dioxide and the crude oil, and solving the technical problem in the field from a new angle;
2. the data obtained by the method is accurate, and the error range is only 1% -2% by comparing with the actual minimum miscible pressure, so that the precision is greatly improved;
3. the method, the device, the equipment and the medium can be popularized and applied to the determination of the mixed phase pressure of different media, such as: reinjection of flue gas or recovered gas improves the field of reservoir recovery, and so forth.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings may be obtained according to the drawings without inventive effort to a person skilled in the art.
FIG. 1 is a schematic overall flow chart of the present invention.
FIG. 2 is a schematic diagram of an experimental apparatus for determining minimum miscible pressure according to the present invention.
Fig. 3 is a schematic view of the piezoelectric transducer according to the present invention.
Fig. 4 is a schematic diagram of the internal structure of the piezoelectric transducer according to the present invention.
FIG. 5 shows the deformation of piezoelectric ceramics by the inverse piezoelectric effect of the present invention, i.e., by an applied electric field. Wherein, (a) is when no electric field is applied; (b) an externally applied electric field; (c) is an externally applied reverse electric field.
FIG. 6 is a flowchart of a first internal process according to the present invention.
FIG. 7 is a flow chart of a second internal process according to the present invention.
FIG. 8 is a flow chart of a third internal process according to the present invention.
FIG. 9 is a graph showing the voltage amplitude of Liaohe oilfield according to the pressure.
Wherein 1 is an ISCO pump, 2 is a six-way valve, 3 is a steel pipeline, 4 is a CO2 piston container, 5 is a stratum water piston container, 6 is a stratum water piston container, 7 is a pressure gauge, 8 is a rubber sleeve, 9 is a core model, 10 is a piston container IV,11 is a back pressure valve, 12 is a measuring cylinder, 13 is a constant temperature box, 14 is a cable, 15 is a computer, 16 is a rear cover plate, 17 is a high-strength stress rod, 18 is an insulating layer, 19 is piezoelectric ceramic, 20 is an acoustic matching layer, 21 is an alternating current power supply, and 22 is a pipeline for injecting fluid.
Detailed Description
In order that those skilled in the art will better understand the present invention, a technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in which it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present invention without making any inventive effort, shall fall within the scope of the present invention.
It should be noted that the terms "first," "second," and the like in the description and the claims of the present invention 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 the embodiments of the invention described herein may be implemented in sequences other than those illustrated or otherwise described herein. 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.
As shown in fig. 1, the present invention provides a method of determining a minimum miscible pressure of carbon dioxide-crude oil, the method comprising the steps of:
firstly, manufacturing a plurality of standby core models according to physical parameters of a target reservoir to be measured; and (5) taking one of the standby core models and placing the standby core models into a core holder. It will be appreciated that in this embodiment, the reservoir physical parameters according to the target to be measured include mainly porosity and permeability. Porosity refers to the percentage of pore volume in the core to the volume of the rock; permeability refers to the ability of fluid to pass through the core. The main consideration of the core used in the displacement experiment is the two parameters.
Secondly, respectively implanting an inlet piezoelectric transducer and an outlet piezoelectric transducer at the centers of the left end face and the right end face of a core model placed in a core holder, wherein the inlet piezoelectric transducer is used for transmitting sound wave signals, and the outlet piezoelectric transducer is used for receiving the sound wave signals transmitted by the core model placed in the core holder and converting the sound wave signals into voltage amplitude signals as shown in fig. 3; the average voltage value U of the voltage amplitude signal is collected and recorded by an industrial computer; the average voltage value U represents the difference between the highest voltage and the lowest voltage in the period from the beginning of the experiment to the breakthrough of the gas to be tested in the piezoelectric transducer at the outlet end. It can be understood that in the application, the experiment starting time is recorded as T0, the time of the gas breaking through the piezoelectric transducer at the outlet end is recorded as T1, then the highest voltage Umax and the lowest voltage Umin are taken in the time period of T0-T1, and then the average voltage value U is expressed as:
U＝Umax-Umin。
thirdly, placing the core holder into a carbon dioxide-crude oil displacement detection device, performing a displacement experiment from carbon dioxide-crude oil non-miscible phase to miscible phase under a first injection end pressure P1, and recording an obtained first voltage amplitude signal U through the industrial computer; and then optionally replacing the core model in the core holder by one of the standby core models, increasing the injection end pressures P1 to Pn step by step, repeatedly executing the second step to the third step, obtaining a series of voltage amplitude signals U1 to Un, and establishing a P-U coordinate system according to the pressures. In the present application, 5MPa is selected as a primary stage, and the injection port pressures P1 to Pn are increased stepwise.
Fourth step, define qn= (Un +1 -Un)/Un is the rate of change of voltage, n > 1, the rate of change of voltage of the series of voltage amplitude signals U1 to Un obtained via the third step is calculated item by item, and a voltage rate data set is established from said rate of change of voltage of U1 to Un.
Fifthly, when Qn is more than 30%, the data points in the corresponding P-U coordinate system are non-miscible data points; when Qn is less than 5%, the data points in the corresponding P-U coordinate system are mixed phase data points.
A sixth step of performing data fitting on the non-miscible data points and the miscible data points obtained in the fifth step in the P-U coordinate system to obtain two fitting curves; and the P value corresponding to the intersection of the two fitting curves is the minimum miscible pressure of the carbon dioxide and the crude oil.
Specifically, in the application, first fitting is performed on non-miscible data points and miscible data points by adopting a y=ax+b function model, if absolute values of correlation coefficients R of a non-miscible data point curve and a miscible data point curve obtained after the first fitting are both greater than 0.9, a data change trend is linear, and the obtained two curves are effective data fitting curves; if the absolute values of the correlation coefficients R of the non-miscible data point curve and the miscible data point curve obtained after the first fitting are smaller than or equal to 0.9, the obtained two curves are invalid data fitting curves; if the absolute value of the correlation coefficient R of only one fitting curve is less than or equal to 0.9, then y=ax is used for the data point corresponding to the fitting curve 2 The second fitting is performed by a +bx+c function model, the parameters a, b, c in the function model being obtained by the least squares method.
Example 1
As one embodiment of the present disclosure, an apparatus for determining minimum miscible pressure of carbon dioxide-crude oil is included herein, including a carbon dioxide-crude oil displacement test assembly. Further, the apparatus further comprises: the sound wave transmitting unit is used for transmitting sound waves at the displacement fluid inlet of the core model; the sound wave receiving unit is used for receiving sound waves at the displacement fluid outlet of the core model and converting the sound waves into voltage signals to be output; the central control unit has the following functions: the starting voltage signal is used for outputting a starting voltage signal for realizing sound wave emission to the sound wave emission unit; the voltage waveform signal is used for receiving the voltage waveform signal output by the sound wave receiving unit; the method comprises the steps of performing data processing on received voltage waveform signals to obtain an average voltage value U, and generating a P-U coordinate system under the control of a first built-in program after receiving input different displacement pressure values P; the method is used for distinguishing mixed phase data points and non-mixed phase data points in the P-U coordinate system under the control of a second built-in program; the method comprises the steps of respectively carrying out data fitting on screened mixed phase data points and non-mixed phase data points under the control of a third built-in program; after the effective data fitting curve is obtained, the displacement pressure value at the intersection point of the two fitting curves is determined.
The sound wave transmitting unit comprises an inlet end piezoelectric transducer: the inlet end piezoelectric transducer is arranged at the center of one end of the core model in the core holder; the sound wave receiving unit comprises an outlet end piezoelectric transducer; the outlet end piezoelectric transducer is arranged at the center of the other end of the core model in the core holder and corresponds to the inlet end piezoelectric transducer.
The inlet-side piezoelectric transducer and the outlet-side piezoelectric transducer each include: the piezoelectric transducer comprises a front cover plate, a rear cover plate, piezoelectric ceramics, an insulating layer, a high-strength stress rod, an acoustic matching layer and a pipeline, wherein the rear cover plate is arranged corresponding to the front cover plate and used for isolating external fluid, the piezoelectric ceramics is arranged at two sides of an internal cavity of the piezoelectric transducer and used for connecting an alternating current power supply, the insulating layer is connected with the piezoelectric ceramics, the high-strength stress rod is used for fixing the front cover plate and the rear cover plate and the piezoelectric ceramics, the acoustic matching layer is arranged at one end of the piezoelectric transducer and used for transmitting and absorbing sound waves, and the pipeline is arranged in an annular space between a metal shell of the piezoelectric transducer and the front cover plate and used for injecting fluid; the high-strength stress rod is fixed in the middle position of the piezoelectric transducer; the piezoelectric ceramic generates ultrasonic waves due to an inverse piezoelectric effect; the sound matching layer at one end of the piezoelectric transducer at the inlet end and the sound matching layer at one end of the piezoelectric transducer at the outlet end are placed in opposite directions.
Example 2
As one embodiment of the present specification, an apparatus for determining a minimum miscible pressure of carbon dioxide-crude oil, the apparatus comprising: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to enable the at least one processor to: outputting a starting voltage signal for realizing sound wave emission to the sound wave emission unit; receiving a voltage waveform signal output by the sound wave receiving unit; the received voltage waveform signals are subjected to data processing to obtain an average voltage value U, and after different displacement pressure values P input by an operator are received, a P-U coordinate system is generated under the control of a first built-in program as shown in fig. 6; under the control of a second built-in program, as shown in fig. 7, the mixed phase data points and the non-mixed phase data points in the P-U coordinate system are screened out; under the control of a third built-in program, as shown in fig. 8, data fitting is performed on the screened mixed phase data points and the screened non-mixed phase data points respectively; after the effective data fitting curve is obtained, the displacement pressure value at the intersection point of the two fitting curves is determined.
Example 3:
as one embodiment of the present specification, a medium for determining a minimum miscible pressure of carbon dioxide-crude oil stores computer executable instructions configured to: outputting a starting voltage signal for realizing sound wave emission to the sound wave emission unit; receiving a voltage waveform signal output by the sound wave receiving unit; the method comprises the steps of performing data processing on received voltage waveform signals to obtain an average voltage value U, and generating a P-U coordinate system under the control of a first built-in program after receiving different displacement pressure values P input by an operator; under the control of a second built-in program, distinguishing mixed phase data points and non-mixed phase data points in the P-U coordinate system; under the control of a third built-in program, respectively performing data fitting on the screened mixed phase data points and the screened non-mixed phase data points; and after obtaining the effective data fitting curve, determining the intersection point of the two fitting curves.
Example 4:
in the embodiment, in a certain block of Liaohe oil field, in order to realize carbon dioxide miscible displacement, the oil field adopts a tubule experiment method to measure the minimum miscible pressure of carbon dioxide and crude oil, and the value is 17.65MPa. However, at this value, no miscible displacement effect is produced during the actual displacement phase. In order to determine the actual minimum miscible pressure of carbon dioxide and crude oil in a stratum, an additional 3 production wells are drilled between injection wells and production wells in Liaohe oil fields, and the minimum miscible pressure of oil and gas reaching the oil and gas miscible displacement effect is finally determined to be 20.04MPa by increasing the injection end pressure and analyzing the oil and gas components. The method takes Liaohe oil field as a target stratum to be measured, and comprises the following specific implementation processes:
firstly, manufacturing a plurality of standby core models according to physical parameters of a reservoir to be measured, and putting one standby core model into a core holder.
For Liaohe oil field, the average permeability of the reservoir of Liaohe oil field is known to be 45×10 -3 μm 2 The average porosity is 23.54%, and the spare core model of Liaohe oil field is 30×4.5×4.5cm (length×width×height) according to the prefabricated cylindrical core model size phi 2.5×30 cm. Further, according to the known average permeability and porosity of the Liaoh oilfield reservoir, the quartz mesh number and cement ratio, which is the ratio of epoxy resin to ethylenediamine, is determined, as an embodiment, in the present application, according to the average permeability and porosity of the Liaoh oilfield reservoir, the quartz mesh number is 140 mesh, and the cement ratio is 95:5. And placing the quartz meshes and the cementing agent into a square plate and uniformly stirring. After the stirring is completed, the quartz sand and the cement are placed in a mold. The mold was placed in a pressure tester and pressurized for 15min. After pressurization, the core is taken out from the mould and placed in an incubator to be baked for 6-8 hours. And after the core is naturally cooled, drilling the core into a required specification by a drilling machine. Repeating the steps to obtain a plurality of standby core models.
Secondly, an inlet end piezoelectric transducer and an outlet end piezoelectric transducer are respectively implanted in the centers of the left end face and the right end face of a core model placed in the core holder. The piezoelectric transducer mentioned in the present application is schematically shown in fig. 3 and 4, and the piezoelectric transducer 9 can mutually convert an acoustic wave signal and an electric signal, and as a preferred embodiment, an inlet piezoelectric transducer and an outlet piezoelectric transducer are provided in the present application; and the inlet-side piezoelectric transducer and the outlet-side piezoelectric transducer each include: the piezoelectric transducer comprises a front cover plate, a rear cover plate 16 which is arranged corresponding to the front cover plate and is used for isolating external fluid, piezoelectric ceramics 19 which are arranged at two sides of an internal cavity of the piezoelectric transducer and are used for connecting an alternating current power supply 21, an insulating layer 18 which is connected with the piezoelectric ceramics 19, high-strength stress rods 17 which are used for fixing the front cover plate, the rear cover plate and the piezoelectric ceramics 19, an acoustic matching layer 20 which is arranged at one end of the piezoelectric transducer and is used for transmitting and absorbing sound waves, and a pipeline 22 which is arranged in an annular space between a metal shell of the piezoelectric transducer and the front cover plate 16 and is used for injecting fluid; the high-strength stress rod 17 is fixed in the middle position of the piezoelectric transducer; the piezoelectric ceramic 19 generates ultrasonic waves due to the inverse piezoelectric effect.
The transducer used in the application is improved based on a TGM-PTZ piezoelectric transducer produced by Hunan Tiangong measurement and control technology Co. The piezoelectric transducer operates using the positive and negative piezoelectric effects of piezoelectric ceramics. The piezoelectric ceramics are stressed to generate charges, which is called positive piezoelectric effect; conversely, applying an electric field to piezoelectric ceramics produces mechanical stress and deformation, a phenomenon known as the inverse piezoelectric effect. Specifically, as shown in fig. 5, wherein (a) is when no electric field is applied; (b) an externally applied electric field; (c) is an externally applied reverse electric field. When the applied electric field is an alternating signal, a corresponding form of elastic wave can be excited in the piezoelectric ceramic. The piezoelectric transducer thus has the dual function of transmitting sound waves and receiving sound waves. As a preferred embodiment, a piezoelectric transducer is deployed at the inlet end to generate a stable acoustic wave, and the same is deployed at the outlet end to receive the acoustic wave, and the minimum miscible pressure is predicted by collecting the outlet end series voltage amplitude U data. Thus, in the present application, the acoustic matching layer at the end of the piezoelectric transducer at the inlet end is disposed opposite to the acoustic matching layer at the end of the piezoelectric transducer at the outlet end.
Further, the core holder is placed in a carbon dioxide-crude oil displacement detection device, and an average voltage value U of the voltage amplitude signal is acquired and recorded through an industrial computer. The specific operation steps are as follows:
firstly, pumping a core in a holder to a negative pressure state by using a vacuum pump, preparing simulated formation water according to the composition of reservoir formation water, and injecting the simulated formation water into the core by using a hand pump; injecting crude oil into the rock core at a constant flow rate by adopting an ISCO pump at the stratum temperature of the target block of 54.3 ℃, and completing saturated oil when the outlet end is not in water outlet; setting ISCO pump as constant pressure injection mode, injection end pressure P 1 Is an initial immiscible pressure of 5MPa; regulating the back pressure valve pressure P 2 Making the pressure lower than the injection end pressure by 0.2MPa, and opening the piezoelectric transducer; carrying out 5MPa displacement experiments under the non-miscible pressure, and receiving acoustic wave signals emitted by an inlet end by an outlet end piezoelectric transducer so as to obtain voltage amplitude U (the amplitude refers to the maximum voltage amplitude before crude oil in a porous medium flows out) after the acoustic waves are converted by the outlet end piezoelectric transducer under 5MPa; changing the core, taking 5MPa as a level, and gradually increasing the injection end pressure P 1 Sequentially carrying out displacement experiments under 10MPa, 15MPa, 20MPa, 26.35MPa, 30MPa and 35MPa to obtain a series P 1 And U data; according to the voltage amplitude U under the non-miscible and miscible pressures along with the pressure P 1 And (3) selecting a function model to fit the data. Definition qn= (Un) +1 -Un/Un is the voltage rate of change, n > 1, the voltage rate of change of the acquired voltage amplitude signals U1 to Un is calculated item by item, and a voltage rate data set is established from the voltage rate of change of said U1 to Un. When Qn is defined to be more than 30%, the data points in the corresponding P-U coordinate system are non-miscible data points; when Qn is less than 5%, the data points in the corresponding P-U coordinate system are mixed phase data points.
And fitting the non-miscible data points and the miscible data points through a y=ax+b function model, wherein the non-miscible section fitting formula is y= -5.522x+158.7R= -0.9987, and the miscible section fitting formula is y= -0.03x+49.267R= -0.9816. The absolute value of the two correlation coefficients R is larger than 0.9, so that the data change trend is linear. The function model a and b are obtained by the least square formulas (1) and (2), and the related data are shown in the tables 1 and 2.
TABLE 1 non-miscible phase section Linear fitting calculation Table
Table 2 miscible phase section linear fit calculation table
Whereby U and P under non-mixed phase and mixed phase can be obtained 1 The functional relation of (2) is the sound wave change rule under the non-miscible and miscible phases, namely the functional relation of U and P1, and when under the non-miscible phase condition, U is reduced along with the increase of P1; while under miscible conditions, U is nearly unchanged as P1 increases.
Wherein i represents a data number, n represents a data amount, and x i Representing injection port pressure Pi, y i Representing a voltage amplitude signal U i 。
As a preferred embodiment, the carbon dioxide-crude oil displacement experimental assembly in the present application is shown in fig. 2, and the device for measuring the minimum miscible pressure of the carbon dioxide-crude oil is connected with an ISCO pump 1 which is used for powering the experimental device and realizes constant-pressure driving and a six-way valve 2 for providing a plurality of channel switch controls through a plurality of steel pipelines 3; a plurality of steel lines 3 are respectively connected to the CO2 piston container 4 (internal CO by displacing the bottom piston 2 Outflow), formation water piston reservoir 5 (inner formation water outflow by displacing the bottom piston), and crude oil piston reservoir 6 (inner crude oil outflow by displacing the bottom piston); CO2 the outlet ends of the piston container 4, the stratum water piston container 5 and the crude oil piston container 6 are respectively provided with a pressure gauge 7 (for monitoring the pressure in the displacement process) and are connected with the inlet end of the clamp holder through pipelines; a core model 9 is fixed in the piston container IV10 through a rubber sleeve 8, and an annulus exists between the rubber sleeve and the whole clamp holder. White oil is injected into the annular space through a valve at the upper part of the holder until the pressure is 2MPa higher than the pressure of the injection end, so that the rubber sleeve is tightly attached to the piezoelectric transducer and the rock core, and the gas is ensured not to overflow through gaps between the rubber sleeve and the piezoelectric transducer and between the rubber sleeve and the rock core in the injection process. At the same time, the outlet end of the holder is connected with the outlet end of the back pressure valve 11 to the measuring cylinder 12. The measurement device further includes: a thermostat 13; the entire displacement system is brought to the formation temperature by placing said incubator 13 in the incubator. The plunger piston container IV10 is also connected to a computer 14 by means of a cable 13.
By comparing the obtained sound wave change rules under the condition of non-miscible phase and miscible phase, obvious difference can be found in the relation between the voltage amplitude and the pressure under the condition of non-miscible phase and miscible phase, because the sound wave transmission medium is changed from 4 types (porous medium, water, carbon dioxide and oil) to 3 types (porous medium, water and carbon dioxide and oil mixture) in the process from non-miscible phase to miscible phase. Therefore, the minimum miscible pressure is defined as P1 corresponding to the abrupt change of the voltage amplitude U, namely P1 where the voltage amplitudes U of the two fitting formulas are equal. The minimum miscible pressure of carbon dioxide and oil was calculated by the non-miscible section fitting equation and the miscible section fitting equation and was 19.93MPa (fig. 9).
The actual minimum miscible pressure of the Liaohe oilfield is 20.04MPa, and the minimum miscible pressure obtained by the method is 19.93MP which is very similar to the actual minimum miscible pressure. The accuracy of the acquisition is far higher than the minimum miscible pressure of carbon dioxide and crude oil measured by a tubule experiment method.
The foregoing embodiment numbers of the present invention are merely for the purpose of description, and do not represent the advantages or disadvantages of the embodiments.
In the foregoing embodiments of the present invention, the descriptions of the embodiments are emphasized, and for a portion of this disclosure that is not described in detail in this embodiment, reference is made to the related descriptions of other embodiments.
In the several embodiments provided in the present application, it should be understood that the disclosed technology content may be implemented in other manners. The above-described embodiments of the apparatus are merely exemplary, and the division of the units, for example, may be a logic function division, and may be implemented in another manner, for example, a plurality of units or components may be combined or may be integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be through some interfaces, units or modules, or may be in electrical or other forms.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied essentially or in part or all of the technical solution or in part in the form of a software product stored in a storage medium, including instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a removable hard disk, a magnetic disk, or an optical disk, or other various media capable of storing program codes.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention.
Claims (8)
1. A method of determining a minimum miscible pressure of carbon dioxide-crude oil, the method comprising:
firstly, manufacturing a plurality of standby core models according to physical parameters of a target reservoir to be measured; any one of the standby core models is taken and placed into a core holder;
secondly, respectively implanting an inlet end piezoelectric transducer and an outlet end piezoelectric transducer at the centers of the left end face and the right end face of a core model placed in a core holder; the inlet end piezoelectric transducer is used for transmitting sound wave signals, and the outlet end piezoelectric transducer is used for receiving the sound wave signals transmitted by the core model arranged in the core holder and converting the sound wave signals into voltage amplitude signals; the average voltage value U of the voltage amplitude signal is collected and recorded by an industrial computer; the average voltage value U represents the difference between the highest voltage and the lowest voltage in the period from the beginning of the experiment to the breakthrough of the gas to be tested in the piezoelectric transducer at the outlet end;
thirdly, placing the core holder into a carbon dioxide-crude oil displacement detection device, performing a displacement experiment from carbon dioxide-crude oil non-miscible phase to miscible phase under a first injection end pressure P1, and recording an obtained first voltage amplitude signal U through an industrial computer; then, optionally replacing the core model in the core holder by one of the standby core models, increasing the injection end pressures P1 to Pn step by step, repeatedly executing the second step to the third step to obtain a series of voltage amplitude signals U1 to Un, and establishing a P-U coordinate system according to the pressures;
fourth step, define qn= (Un +1 -Un)/Un is the rate of change of voltage, n > 1, calculating the rate of change of voltage of the series of voltage amplitude signals U1 to Un obtained via the third step item by item, creating a voltage rate data set from the rate of change of voltage of said U1 to Un;
fifthly, when Qn is defined to be larger than a first set value, the data points in the corresponding P-U coordinate system are non-miscible data points; when Qn is defined to be smaller than a second set value, the data points in the corresponding P-U coordinate system are mixed phase data points;
a sixth step of performing data fitting on the non-miscible data points and the miscible data points obtained in the fifth step in the P-U coordinate system to obtain two fitting curves; and the P value corresponding to the intersection of the two fitting curves is the minimum miscible pressure of the carbon dioxide and the crude oil.
2. A method of determining the minimum miscible pressure of a carbon dioxide crude oil as claimed in claim 1, wherein said data fitting of the sixth step is performed according to the following path:
firstly, adopting a y=ax+b function model to respectively perform first fitting on non-miscible data points and miscible data points, and if absolute values of correlation coefficients R of a non-miscible data point curve and a miscible data point curve obtained after the first fitting are both larger than 0.9, the data change trend is linear, and the two obtained curves are effective data fitting curves; if the absolute values of the correlation coefficients R of the non-miscible data point curve and the miscible data point curve obtained after the first fitting are smaller than or equal to 0.9, the obtained two curves are invalid data fitting curves; if the absolute value of the correlation coefficient R of only one fitting curve is less than or equal to 0.9, then y=ax is used for the data point corresponding to the fitting curve 2 The second fitting is performed by a +bx+c function model, the parameters a, b, c in the function model being obtained by the least squares method.
3. A method of determining the minimum miscible pressure of a carbon dioxide crude oil as claimed in claim 1, wherein: in the third step, the injection end pressures P1 to Pn are increased step by taking a third set value as a first level; the third set value is a preset pressure value.
4. An apparatus for determining a minimum miscible pressure of a carbon dioxide-crude oil comprising a carbon dioxide-crude oil displacement experiment assembly, the apparatus further comprising:
the sound wave transmitting unit is used for transmitting sound waves at the displacement fluid inlet of the core model;
the sound wave receiving unit is used for receiving sound waves at the displacement fluid outlet of the core model and converting the sound waves into voltage signals to be output;
the central control unit has the following functions:
the starting voltage signal is used for outputting a starting voltage signal for realizing sound wave emission to the sound wave emission unit;
the voltage waveform signal is used for receiving the voltage waveform signal output by the sound wave receiving unit;
the method comprises the steps of performing data processing on received voltage waveform signals to obtain an average voltage value U, and generating a P-U coordinate system under the control of a first built-in program after receiving input different displacement pressure values P;
the method is used for distinguishing mixed phase data points and non-mixed phase data points in the P-U coordinate system under the control of a second built-in program;
the method comprises the steps of respectively carrying out data fitting on screened mixed phase data points and non-mixed phase data points under the control of a third built-in program; after the effective data fitting curve is obtained, the displacement pressure value at the intersection point of the two fitting curves is determined.
5. The apparatus for determining a minimum miscible pressure of a carbon dioxide crude oil as claimed in claim 4, wherein said acoustic wave transmitting unit comprises an inlet side piezoelectric transducer: the inlet end piezoelectric transducer is arranged at the center of one end of the core model in the core holder;
the sound wave receiving unit comprises an outlet end piezoelectric transducer; the outlet end piezoelectric transducer is arranged at the center of the other end of the core model in the core holder and corresponds to the inlet end piezoelectric transducer.
6. The apparatus for determining the minimum miscible pressure of a carbon dioxide crude oil of claim 5, wherein said inlet and outlet piezoelectric transducers each comprise: the piezoelectric transducer comprises a front cover plate, a rear cover plate, piezoelectric ceramics, an insulating layer, a high-strength stress rod, an acoustic matching layer and a pipeline, wherein the rear cover plate is arranged corresponding to the front cover plate and used for isolating external fluid, the piezoelectric ceramics is arranged at two sides of an internal cavity of the piezoelectric transducer and used for connecting an alternating current power supply, the insulating layer is connected with the piezoelectric ceramics, the high-strength stress rod is used for fixing the front cover plate and the rear cover plate and the piezoelectric ceramics, the acoustic matching layer is arranged at one end of the piezoelectric transducer and used for transmitting and absorbing sound waves, and the pipeline is arranged in an annular space between a metal shell of the piezoelectric transducer and the front cover plate and used for injecting fluid; the high-strength stress rod is fixed in the middle position of the piezoelectric transducer; the piezoelectric ceramic generates ultrasonic waves due to an inverse piezoelectric effect;
the sound matching layer at one end of the piezoelectric transducer at the inlet end and the sound matching layer at one end of the piezoelectric transducer at the outlet end are placed in opposite directions.
7. An apparatus for determining a minimum miscible pressure of carbon dioxide-crude oil, the apparatus comprising:
at least one processor; the method comprises the steps of,
a memory communicatively coupled to the at least one processor; wherein,
the memory stores instructions executable by the at least one processor to enable the at least one processor to:
outputting a starting voltage signal for realizing sound wave emission to a sound wave emission unit;
receiving a voltage waveform signal output by an acoustic wave receiving unit;
the method comprises the steps of performing data processing on received voltage waveform signals to obtain an average voltage value U, and generating a P-U coordinate system under the control of a first built-in program after receiving different displacement pressure values P input by an operator;
under the control of a second built-in program, distinguishing mixed phase data points and non-mixed phase data points in the P-U coordinate system;
under the control of a third built-in program, respectively performing data fitting on the screened mixed phase data points and the screened non-mixed phase data points; after the effective data fitting curve is obtained, the displacement pressure value at the intersection point of the two fitting curves is determined.
8. A medium for determining a minimum miscible pressure of carbon dioxide-crude oil, storing computer executable instructions, the computer executable instructions configured to:
outputting a starting voltage signal for realizing sound wave emission to a sound wave emission unit;
receiving a voltage waveform signal output by an acoustic wave receiving unit;
the method comprises the steps of performing data processing on received voltage waveform signals to obtain an average voltage value U, and generating a P-U coordinate system under the control of a first built-in program after receiving different displacement pressure values P input by an operator;
under the control of a second built-in program, distinguishing mixed phase data points and non-mixed phase data points in the P-U coordinate system;
under the control of a third built-in program, respectively performing data fitting on the screened mixed phase data points and the screened non-mixed phase data points; after the effective data fitting curve is obtained, the displacement pressure value at the intersection point of the two fitting curves is determined.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211441721.7A CN115822532B (en) | 2022-11-17 | 2022-11-17 | Method, device, equipment and medium for determining minimum miscible pressure of carbon dioxide-crude oil |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211441721.7A CN115822532B (en) | 2022-11-17 | 2022-11-17 | Method, device, equipment and medium for determining minimum miscible pressure of carbon dioxide-crude oil |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115822532A CN115822532A (en) | 2023-03-21 |
| CN115822532B true CN115822532B (en) | 2024-01-30 |
Family
ID=85528834
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211441721.7A Active CN115822532B (en) | 2022-11-17 | 2022-11-17 | Method, device, equipment and medium for determining minimum miscible pressure of carbon dioxide-crude oil |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115822532B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117005853A (en) * | 2023-07-06 | 2023-11-07 | 中国石油大学(北京) | Mixed-phase pressure measurement device and method |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4455860A (en) * | 1982-12-07 | 1984-06-26 | Mobil Oil Corporation | Method and apparatus for determining CO2 minimum miscibility pressure of reservoir oil |
| CN104005741A (en) * | 2014-06-04 | 2014-08-27 | 中国石油大学(华东) | Multi-phase fluid oil displacement and pulse unblocking integrated physical simulation experiment device and method |
| CN105258840A (en) * | 2015-11-12 | 2016-01-20 | 东北石油大学 | Method for measuring optimum miscible pressure and miscible region of carbon dioxide between injection and production wells |
| CN106837269A (en) * | 2017-03-03 | 2017-06-13 | 中国石油大学(北京) | Low, the extra-low permeability oil reservoirs CO of one kind2Drive nearly miscible pressure area determination method |
| CN107542436A (en) * | 2017-10-23 | 2018-01-05 | 大庆东油睿佳石油科技有限公司 | A kind of device simulated parallel water horizontal well mixed phase and drive exploitation of gas hydrate |
| CN108222899A (en) * | 2017-12-30 | 2018-06-29 | 东北石油大学 | Compare the method and apparatus in carbon dioxide non-phase-mixing driving stage |
| CN112627783A (en) * | 2019-09-24 | 2021-04-09 | 中国石油天然气股份有限公司 | Experimental device for improving gas injection recovery ratio through low-frequency pressure transformation |
| CN113187465A (en) * | 2021-05-31 | 2021-07-30 | 中国石油大学(北京) | Fluid dynamic monitoring method and system of non-condensable gas huff and puff three-dimensional physical model |
| CN115146556A (en) * | 2022-06-09 | 2022-10-04 | 中国石油天然气股份有限公司 | Carbon dioxide-crude oil minimum miscible pressure prediction method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2441143C1 (en) * | 2010-05-19 | 2012-01-27 | Общество с ограниченной ответственностью "Нефтяная технологическая компания" НЕТЕК | Method of increasing oil reservoir recovery and device to this end |
| US9851339B2 (en) * | 2014-08-15 | 2017-12-26 | Energy And Environmental Research Center Foundation | Sampling an oil composition and determining minimum miscibility pressure of an oil composition with a fluid |
-
2022
- 2022-11-17 CN CN202211441721.7A patent/CN115822532B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4455860A (en) * | 1982-12-07 | 1984-06-26 | Mobil Oil Corporation | Method and apparatus for determining CO2 minimum miscibility pressure of reservoir oil |
| CN104005741A (en) * | 2014-06-04 | 2014-08-27 | 中国石油大学(华东) | Multi-phase fluid oil displacement and pulse unblocking integrated physical simulation experiment device and method |
| CN105258840A (en) * | 2015-11-12 | 2016-01-20 | 东北石油大学 | Method for measuring optimum miscible pressure and miscible region of carbon dioxide between injection and production wells |
| CN106837269A (en) * | 2017-03-03 | 2017-06-13 | 中国石油大学(北京) | Low, the extra-low permeability oil reservoirs CO of one kind2Drive nearly miscible pressure area determination method |
| CN107542436A (en) * | 2017-10-23 | 2018-01-05 | 大庆东油睿佳石油科技有限公司 | A kind of device simulated parallel water horizontal well mixed phase and drive exploitation of gas hydrate |
| CN108222899A (en) * | 2017-12-30 | 2018-06-29 | 东北石油大学 | Compare the method and apparatus in carbon dioxide non-phase-mixing driving stage |
| CN112627783A (en) * | 2019-09-24 | 2021-04-09 | 中国石油天然气股份有限公司 | Experimental device for improving gas injection recovery ratio through low-frequency pressure transformation |
| CN113187465A (en) * | 2021-05-31 | 2021-07-30 | 中国石油大学(北京) | Fluid dynamic monitoring method and system of non-condensable gas huff and puff three-dimensional physical model |
| CN115146556A (en) * | 2022-06-09 | 2022-10-04 | 中国石油天然气股份有限公司 | Carbon dioxide-crude oil minimum miscible pressure prediction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115822532A (en) | 2023-03-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107063963B (en) | Device and method for testing micro-crack expansion and seepage characteristics of tight reservoir | |
| RU2331057C2 (en) | Method and device for evaluation of physical parametres of undeground deposit of mineral wealth on base of study of rock fragments selected from this deposit | |
| US7472588B2 (en) | Petrophysical fluid flow property determination | |
| CN107563899B (en) | Oil-gas well productivity prediction method and device | |
| CN105738252B (en) | The measurement method of the flowable aperture boundary of viscous crude in a kind of crack | |
| CN115822532B (en) | Method, device, equipment and medium for determining minimum miscible pressure of carbon dioxide-crude oil | |
| CN109408855B (en) | Method and device for calculating fracture conductivity value of horizontal well to be fractured | |
| US7753118B2 (en) | Method and tool for evaluating fluid dynamic properties of a cement annulus surrounding a casing | |
| WO2011133885A1 (en) | Total storage capacity and total porosity of porous media | |
| CN103335928A (en) | Method and apparatus for measuring permeability of porous rock | |
| CN113484216B (en) | Method for evaluating water phase flowback rate and reasonable flowback pressure difference of tight sandstone gas reservoir | |
| CN107130955B (en) | Method for determining bottom hole flowing pressure and method for determining natural energy of reservoir | |
| CN106682384B (en) | Wellbore enlarging rate calculation method and device | |
| US11834620B2 (en) | Centrifuge energy harvesting chamber | |
| Abdulhadi et al. | Experimental study of wellbore instability in clays | |
| Guo et al. | Water invasion and remaining gas distribution in carbonate gas reservoirs using core displacement and NMR | |
| WO2004099552A3 (en) | Determining gradients using a multi-probed formation tester | |
| CN111963149A (en) | Post-fracturing formation pressure calculation method considering stagnant ground fluid volume pressurization | |
| Liu | Estimation of different coal compressibilities of coalbed methane reservoirs under replicated in situ condition | |
| CN209182206U (en) | Buried hydraulic tunnel osmotic gradient simulation experiment system | |
| CN113670769B (en) | Method for simulating gas content change in marine shale stratum lifting process | |
| CN110095584A (en) | A kind of reservoir oil water saturation bearing calibration | |
| CN106285657B (en) | For determining that the method and system of collapsing event occur for fracture-pore reservoir solution cavity | |
| CN104948175B (en) | A kind of Experimental Method in Laboratory for monitoring interlayer liquid-producing capacity difference | |
| Zhang et al. | Determining the hydraulic properties of saturated, low-permeability geological materials in the laboratory: Advances in theory and practice |
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 |