CN112836860A - Method and system for determining dynamic yield of fractured well in whole period of yield decreasing stage - Google Patents
Method and system for determining dynamic yield of fractured well in whole period of yield decreasing stage Download PDFInfo
- Publication number
- CN112836860A CN112836860A CN202110038328.2A CN202110038328A CN112836860A CN 112836860 A CN112836860 A CN 112836860A CN 202110038328 A CN202110038328 A CN 202110038328A CN 112836860 A CN112836860 A CN 112836860A
- Authority
- CN
- China
- Prior art keywords
- yield
- decreasing
- rate
- stage
- well
- 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.)
- Granted
Links
- 230000003247 decreasing effect Effects 0.000 title claims abstract description 177
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000004364 calculation method Methods 0.000 claims abstract description 47
- 238000004519 manufacturing process Methods 0.000 claims description 69
- 230000007423 decrease Effects 0.000 claims description 41
- 238000010276 construction Methods 0.000 claims description 4
- 238000013480 data collection Methods 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 2
- 230000004048 modification Effects 0.000 description 6
- 238000012986 modification Methods 0.000 description 6
- 230000009467 reduction Effects 0.000 description 5
- 230000007547 defect Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Animal Husbandry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Mining & Mineral Resources (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention belongs to the technical field of unconventional reservoir hydraulic fracturing exploitation well yield prediction, and discloses a method and a system for determining dynamic yield in a whole period of a fracturing well yield decrement stage, wherein an unconventional reservoir fracturing well yield decrement whole period yield decrement rate is constructed; integrating the constructed full-period yield decreasing rate by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the yield decreasing stage of the fracturing well; collecting data on yield versus time; fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well to obtain an initial decreasing rate, a later decreasing rate, a time coefficient and a decreasing index of undetermined parameters; under the condition that the yield value of a certain oil reservoir block of the same type and a well of the same type at the initial decrement stage is determined, the yield at any moment of the decrement stage is calculated by using a yield calculation expression of the yield decrement stage of the fracturing well, and finally the dynamic yield of the unconventional reservoir hydraulic fracturing well in the whole period of the yield decrement stage is determined.
Description
Technical Field
The invention belongs to the technical field of yield prediction of unconventional reservoir hydraulic fracturing development wells, and particularly relates to a method and a system for determining dynamic yield of a fracturing well in a whole period at a yield decreasing stage.
Background
At present: unconventional resources such as dense oil gas, shale oil gas and the like in China are wide in distribution and large in reserves, and have a material basis for efficient development. Because the unconventional reservoir usually has poor physical properties, the economic and efficient development of unconventional resources is realized by adopting a hydraulic fracturing reservoir modification technology. The modification effect (such as modification area, modification degree and the like) of the hydraulic fracturing reservoir is controlled by the process technology and the lithology of the reservoir, and the modification effect is limited, so that the production process of the hydraulic fracturing well quickly enters a yield decreasing stage, and the yield decreasing usually presents an L-shaped characteristic: the yield decrease is fast in the initial stage, slow in the later stage and low. The "L" type yield decrement feature subjects the initial and later stages of yield decrement to different decrement types, which a single yield decrement type (Arps decrement, power function decrement, linear decrement, etc.) has not been fully described. The L-shaped characteristic describing the yield decrement of the unconventional reservoir fractured well is a basis for accurately predicting the yield of the fractured well in the whole period of the yield decrement stage. At present, some researches for describing the L-shaped characteristic of the yield decline of the fractured well based on data driving exist, but no description method with physical significance is given, and no yield value in the whole cycle of the yield decline is given.
Through the above analysis, the problems and defects of the prior art are as follows: at present, some researches for describing the characteristic of the yield decline of the fractured well in an L shape based on data driving exist, but no description method with physical significance is given, the decline type of the yield decline curve in the L shape is difficult to identify, and the yield value in the whole period of the yield decline stage is not given.
The difficulty in solving the above problems and defects is: how to establish the describing method of the L-shaped descending curve with physical significance has certain difficulty.
The significance of solving the problems and the defects is as follows: a description method for establishing an L-shaped descending curve with physical significance has important significance for recognizing the descending type of the yield in the whole period of the descending stage of the yield and predicting the whole period of the descending stage of the yield.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a method and a system for determining the dynamic yield of a fractured well in a whole period of a yield decreasing stage.
The invention is realized in such a way that a method for determining the dynamic yield of a fractured well in the whole period of the yield decreasing stage comprises the following steps:
constructing a full-period yield reduction rate of an unconventional reservoir fracturing well in a yield reduction stage;
integrating the constructed full-period yield decreasing rate by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the yield decreasing stage of the fracturing well;
collecting data on yield versus time;
fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well, solving the fitting problem by adopting a Levenberg-Marquardt method, and obtaining an initial decreasing rate, a later decreasing rate, a time coefficient and a decreasing index of undetermined parameters; the solving process is realized by a computer (the model of the computer is OPTIPLEX 7040, the memory is 8G, and the processor is intel (r) core (tm) i7-6700 CPU @3.4 GHz).
And calculating the yield at any moment of the decrement stage by using a yield calculation expression of the yield decrement stage of the fracturing well under the condition of determining the decrement initial yield value of a certain oil reservoir block of the same type and a well of the same type by using the initial decrement rate, the later decrement rate, the time coefficient and the decrement index of the undetermined parameters obtained by fitting, and finally determining the dynamic yield of the unconventional reservoir hydraulic fracturing well in the whole period of the yield decrement stage.
Further, fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well, solving the fitting problem by adopting a Levenberg-Marquardt method, and obtaining an initial decreasing rate D of undetermined parametersiLate stage decreasing rate D∞A time coefficient m and a decreasing index n; initial decrement rate D of parameter to be determined obtained by fittingiLate stage decreasing rate D∞Time coefficient m and decreasing index n, and determining the decreasing initial yield q of a certain oil reservoir block and a well of the same typeiUnder the condition of the value, the yield at any moment of the decreasing stage is calculated by using a yield calculation expression of the decreasing stage of the fracturing well yield, and finally the dynamic yield of the unconventional reservoir hydraulic fracturing well in the whole period of the decreasing stage of the yield is determined.
Further, the method for determining the dynamic production in the whole period of the yield decreasing stage of the fractured well constructs the yield decreasing rate of the unconventional reservoir fractured well in the whole period of the yield decreasing stage in the form of
Wherein D is the rate of decrease, 1/day; d i1/day for initial decline rate of yield decline; d∞The later decline rate of the yield decline is 1/day; t is decreasing time, days; m is a time coefficient and is dimensionless; n is a decreasing index and is dimensionless.
Further, the method for determining the dynamic production within the whole period of the production decline stage of the fractured well defines the production decline rate
q is the yield, in square/day, and is substituted into the constructed full-period yield decreasing rate, and a yield calculation expression of the yield decreasing stage of the fracturing well can be obtained by applying separation variable integration
Wherein q isiThe initial yield is decreased progressively, and the yield is decreased progressively,recipe/day; the yield calculation expression of the yield decrement stage of the fractured well comprises a time coefficient with the unit of m, and the parameters are used for segmenting yield decrement data, so that the L-shaped characteristic of the yield decrement of the unconventional reservoir fractured well can be described finally.
Further, the method for determining the dynamic yield of the fracturing well in the whole period of the yield decreasing stage collects the relation data of the yield of the fracturing A well of a certain tight oil reservoir and the time.
Further, the method for determining the dynamic production of the fractured well in the whole period of the production decline stage applies an equation
Fitting data of the output changing with time, solving the nonlinear least square problem by adopting a Levenberg-Marquardt method to obtain the initial decrement rate D of the undetermined parameteri0.001, late stage decreasing rate D∞0.016, 20.768 and 0.792.
Further, the method for determining the dynamic yield of the fractured well in the whole period of the yield decreasing stage utilizes the initial decreasing rate D obtained by fittingi0.001, late stage decreasing rate D∞When a certain oil reservoir block of the same type and a B well of the same type are determined, the initial yield q is decreased gradually, wherein the time coefficient m is 20.768 and the decreasing index n is 0.792iIn the case of 14.78, the formula is applied
And calculating to obtain the yield at any moment of the decreasing stage, and finally determining the dynamic yield of the hydraulic fracturing well in the whole period of the decreasing stage.
It is a further object of the invention to provide a computer device comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of:
constructing a yield decrement rate of the unconventional reservoir fracturing well in the whole period;
integrating the constructed full-period yield decreasing rate by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the yield decreasing stage of the fracturing well;
collecting data on yield versus time;
fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well, solving the fitting problem by adopting a Levenberg-Marquardt method, and obtaining an initial decreasing rate, a later decreasing rate, a time coefficient and a decreasing index of undetermined parameters;
and calculating the yield at any moment of the decrement stage by using a yield calculation expression of the yield decrement stage of the fracturing well under the condition of determining the decrement initial yield value of a certain oil reservoir block of the same type and a well of the same type by using the initial decrement rate, the later decrement rate, the time coefficient and the decrement index of the undetermined parameters obtained by fitting, and finally determining the dynamic yield of the unconventional reservoir hydraulic fracturing well in the whole period of the yield decrement stage.
Another object of the present invention is to provide a dynamic production system for determining a dynamic production rate of a fractured well within a whole period of a decreasing production rate phase of the fractured well, which implements the method for determining a dynamic production rate of a fractured well within a whole period of a decreasing production rate phase of the fractured well, the system comprising:
the decrement rate construction module is used for constructing the decrement rate of the production of the unconventional reservoir fracturing well in the whole period;
the yield expression obtaining module is used for integrating the constructed full-period yield decreasing rate by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the fractured well in the yield decreasing stage;
the data collection module is used for collecting data of yield and time;
the parameter acquisition module is used for fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well, solving the fitting problem by adopting a Levenberg-Marquardt method and obtaining an initial decreasing rate, a later decreasing rate, a time coefficient and a decreasing index of undetermined parameters;
and the yield result output module is used for calculating the yield at any moment in the decreasing stage by using a yield calculation expression in the yield decreasing stage of the fractured well under the condition of determining a certain oil reservoir block of the same type and the decreasing initial yield value of the same type of well by using the initial decreasing rate, the later decreasing rate, the time coefficient and the decreasing index of the undetermined parameters obtained by fitting, and finally determining the dynamic yield in the full period of the yield decreasing stage of the unconventional reservoir hydraulic fractured well.
The invention also aims to provide a yield prediction terminal for the unconventional reservoir hydraulic fracturing exploitation well, which is used for realizing the method for determining the dynamic yield of the fracturing well in the decreasing stage of the yield in the whole period. s
By combining all the technical schemes, the invention has the advantages and positive effects that: in order to further verify the innovativeness of the patent, the HP30 well yield is calculated by respectively applying the calculation result of the method provided by the invention and a method provided by a new model for analyzing the unconventional oil 2 well yield decrement law, which is published in 2016 by the paper "the university of mining in china" journal "in the 4 th, and the comparison result is shown in fig. 3 below. As can be seen from FIG. 3, the calculation results of the present invention are more accurate than those of the Qidong calculation.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed to be used in the embodiments of the present application will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present application, and it is obvious for those skilled in the art that other drawings can be obtained from the drawings without creative efforts.
FIG. 1 is a flow chart of a method for determining a dynamic production rate of a fractured well during a full cycle of a production decline phase according to an embodiment of the invention.
FIG. 2 is a schematic diagram of a system for determining the dynamic production rate of a fractured well over a full cycle of a production decline phase according to an embodiment of the present invention;
in fig. 2: 1. a decrement rate construction module; 2. a yield expression obtaining module; 3. a data collection module; 4. a parameter acquisition module; 5. and a yield result output module.
Fig. 3 is a comparison diagram of the calculation results provided by the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
In view of the problems in the prior art, the present invention provides a method and a system for determining the dynamic production of a fractured well in a whole period of a production decrement stage, and the present invention is described in detail with reference to the accompanying drawings.
As shown in fig. 1, the method for determining the dynamic production of the fractured well in the whole period of the production decline stage provided by the invention comprises the following steps:
s101: constructing a yield decrement rate of the unconventional reservoir fracturing well in the whole period;
s102: integrating the constructed full-period yield decreasing rate by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the yield decreasing stage of the fracturing well;
s103: collecting data on yield versus time;
s104: fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well, solving the fitting problem by adopting a Levenberg-Marquardt method, and finally obtaining the initial decreasing rate, the later decreasing rate, the time coefficient, the decreasing index and the like of undetermined parameters;
s105: and calculating the yield at any moment of the decreasing stage by using a yield calculation expression of the fracturing well yield decreasing stage under the condition of determining the decreasing initial yield value of a certain oil reservoir block of the same type and a well of the same type by using the initial decreasing rate, the later decreasing rate, the time coefficient, the decreasing index and the like of the undetermined parameters obtained by fitting, and finally determining the dynamic yield of the unconventional reservoir hydraulic fracturing well yield decreasing stage in the whole period.
Those skilled in the art can also implement other steps, and the method for determining the dynamic production of the fractured well within the whole period of the production decrement stage provided by the invention of fig. 1 is only one specific example.
As shown in fig. 2, the present invention provides a system for determining the dynamic production rate of a fractured well during a full cycle of a production rate decreasing phase, comprising:
the decrement rate construction module 1 is used for constructing the decrement rate of the production of the unconventional reservoir fracturing well in the whole period;
the yield expression obtaining module 2 is used for integrating the full-period yield decreasing rate constructed by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the fractured well in the yield decreasing stage;
the data collection module 3 is used for collecting data of yield and time;
the parameter obtaining module 4 is used for fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well, solving the fitting problem by adopting a Levenberg-Marquardt method, and finally obtaining an initial decreasing rate, a later decreasing rate, a time coefficient, a decreasing index and the like of undetermined parameters;
and the yield result output module 5 is used for calculating the yield at any time of the decreasing stage by using the yield calculation expression of the fracturing well yield decreasing stage under the condition of determining the oil reservoir block of a certain same type and the decreasing initial yield value of the same type of well by using the initial decreasing rate, the later decreasing rate, the time coefficient, the decreasing index and the like of the undetermined parameters obtained by fitting, and finally determining the dynamic yield of the unconventional reservoir hydraulic fracturing well in the whole period of the yield decreasing stage.
The technical solution of the present invention is further described below with reference to the accompanying drawings.
The method for determining the dynamic yield of the fractured well in the whole period of the yield decreasing stage comprises the following steps: and constructing the yield decrement rate of the unconventional reservoir fractured well in the whole period. And integrating the constructed full-period yield decreasing rate by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the fractured well in the yield decreasing stage. Yield versus time data were collected. Fitting the data of the yield changing along with the time by using a yield calculation expression in the yield decreasing stage of the fracturing well, solving the fitting problem by using a Levenberg-Marquardt method, and finally obtaining the data to be obtainedInitial rate of decline (D) with fixed parametersi) Late stage rate of decline (D)∞) Time coefficient (m), and decreasing index (n). Preliminary rate of decline (D) of parameter to be determined obtained using fittingi) Late stage rate of decline (D)∞) Time coefficient (m), decreasing index (n), etc. to determine the decreasing initial production q of a certain oil reservoir block and well of the same typeiUnder the condition of the value, the yield at any moment of the decreasing stage is calculated by using a yield calculation expression of the decreasing stage of the fracturing well yield, and finally the dynamic yield of the unconventional reservoir hydraulic fracturing well in the whole period of the decreasing stage of the yield is determined.
The method for determining the dynamic yield of the fractured well in the whole period of the yield decreasing stage provided by the embodiment of the invention specifically comprises the following steps:
the method comprises the following steps: constructing the whole-period yield reduction rate of the unconventional reservoir fracturing well in the yield reduction stage in the form of
(wherein D is the rate of decrease, 1/day; D i1/day for initial decline rate of yield decline; d∞The later decline rate of the yield decline is 1/day; t is decreasing time, days; m is a time coefficient and is dimensionless; n is a decreasing index, dimensionless).
Step two: defining the yield reduction rate
(q is yield, square/day), substituting into the constructed full-period yield decrement rate, and obtaining a yield calculation expression of the yield decrement stage of the fractured well by applying separation variable integration
(wherein q isiDecreasing initial yield, square/day). The time coefficient (m) is included in the yield calculation expression of the yield decreasing stage of the fractured well, and the parameter is used for segmenting yield decreasing data, so that the L-shaped characteristic of the yield decreasing of the unconventional reservoir fractured well can be described finally.
Step three: data were collected on the relationship of production of a tight reservoir fracturing well a to time (table 1).
TABLE 1 data of relationship between the fracturing A well yield of a tight reservoir and time
Step four: using equations
And fitting the data of the yield changing along with the time, wherein the fitting process is a nonlinear least square problem. The Levenberg-Marquardt method is adopted to solve the nonlinear least square problem, and finally the initial decrement rate (D) of the parameter to be determined can be obtainedi0.001), late stage rate of decrease (D)∞0.016), time coefficient (m 20.768), decreasing index (n 0.792), etc.
Step five: initial rate of decline (D) obtained using fittingi0.001), late stage rate of decrease (D)∞0.016), time coefficient (m 20.768), decreasing index (n 0.792), and decreasing initial yield q in determining a certain reservoir block and a B well of the same typeiIn the case of 14.78, the formula can be applied
And (4) calculating to obtain the yield at any time of the decreasing stage, and finally determining the dynamic yield of the hydraulic fracturing well in the whole period of the decreasing stage (Table 2). Actual value and meter of dynamic production in whole period of B well production decrementThe sum of the squares of the differences between the calculated values is small, 149.58; the method can accurately predict the full-period dynamic yield of the descending stage of the fractured well.
TABLE 2 comparison of calculated and actual values of dynamic production over the full period of the yield decline phase for fracturing B well
The technical effects of the present invention will be described in detail with reference to experiments.
In order to further verify the innovativeness of the patent, the HP30 well yield is calculated by respectively applying the calculation result of the method provided by the invention and a method provided by a new model for analyzing the unconventional oil 2 well yield decrement law, which is published in 2016 by the paper "the university of mining in china" journal "in the 4 th, and the comparison result is shown in fig. 3 below. As can be seen from FIG. 3, the calculation results of the present invention are more accurate than those of the Qidong calculation.
It should be noted that the embodiments of the present invention can be realized by hardware, software, or a combination of software and hardware. The hardware portion may be implemented using dedicated logic; the software portions may be stored in a memory and executed by a suitable instruction execution system, such as a microprocessor or specially designed hardware. Those skilled in the art will appreciate that the apparatus and methods described above may be implemented using computer executable instructions and/or embodied in processor control code, such code being provided on a carrier medium such as a disk, CD-or DVD-ROM, programmable memory such as read only memory (firmware), or a data carrier such as an optical or electronic signal carrier, for example. The apparatus and its modules of the present invention may be implemented by hardware circuits such as very large scale integrated circuits or gate arrays, semiconductors such as logic chips, transistors, or programmable hardware devices such as field programmable gate arrays, programmable logic devices, etc., or by software executed by various types of processors, or by a combination of hardware circuits and software, e.g., firmware.
The above description is only for the purpose of illustrating the present invention and the appended claims are not to be construed as limiting the scope of the invention, which is intended to cover all modifications, equivalents and improvements that are within the spirit and scope of the invention as defined by the appended claims.
Claims (10)
1. A method for determining the dynamic production of a fractured well in a whole period of a yield decreasing stage is characterized by comprising the following steps:
constructing a yield decrement rate of the unconventional reservoir fracturing well in the whole period;
integrating the constructed full-period yield decreasing rate by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the yield decreasing stage of the fracturing well;
collecting data on yield versus time;
fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well, solving the fitting problem by adopting a Levenberg-Marquardt method, and obtaining an initial decreasing rate, a later decreasing rate, a time coefficient and a decreasing index of undetermined parameters;
and calculating the yield at any moment of the decrement stage by using a yield calculation expression of the yield decrement stage of the fracturing well under the condition of determining the decrement initial yield value of a certain oil reservoir block of the same type and a well of the same type by using the initial decrement rate, the later decrement rate, the time coefficient and the decrement index of the undetermined parameters obtained by fitting, and finally determining the dynamic yield of the unconventional reservoir hydraulic fracturing well in the whole period of the yield decrement stage.
2. The method for determining the dynamic production within the whole period of the decreasing production phase of the fractured well as recited in claim 1, wherein the data of the production change along with the time are fitted by using a production calculation expression of the decreasing production phase of the fractured well, and the fitting problem is solved by using a Levenberg-Marquardt method to obtain the initial decreasing rate D of the undetermined parametersiLate stage decreasing rate D∞A time coefficient m and a decreasing index n; initial decrement rate D of parameter to be determined obtained by fittingiLate stage decreasing rate D∞Time coefficient m and decreasing index n, and determining the decreasing initial yield q of a certain oil reservoir block and a well of the same typeiUnder the condition of the value, the yield at any moment of the decreasing stage is calculated by using a yield calculation expression of the decreasing stage of the fracturing well yield, and finally the dynamic yield of the unconventional reservoir hydraulic fracturing well in the whole period of the decreasing stage of the yield is determined.
3. The method for determining the dynamic production rate over a full cycle of a declining phase of production from a fractured well of claim 1 wherein the method for determining the dynamic production rate over a full cycle of a declining phase of production from a fractured well establishes the rate of production decrease over a full cycle of a declining phase of production from an unconventional reservoir fractured well in the form of a rate of production decrease over a full cycle of a declining phase of production from a fractured well
Wherein D is the rate of decrease, 1/day; di1/day for initial decline rate of yield decline; d∞The later decline rate of the yield decline is 1/day; t is decreasing time, days; m is a time coefficient and is dimensionless; n is a decreasing index and is dimensionless.
4. The method for determining the dynamic production rate over a complete cycle of a decreasing production phase of a fractured well according to claim 3, wherein the method for determining the dynamic production rate over a complete cycle of a decreasing production phase of a fractured well defines a rate of production decline
q is the yield, in square/day, and is substituted into the constructed full-period yield decreasing rate, and a yield calculation expression of the yield decreasing stage of the fracturing well can be obtained by applying separation variable integration
Wherein q isiDecreasing the initial yield per day; the yield calculation expression of the yield decrement stage of the fractured well comprises a time coefficient with the unit of m, and the parameters are used for segmenting yield decrement data, so that the L-shaped characteristic of the yield decrement of the unconventional reservoir fractured well can be described finally.
5. The method for determining the dynamic production over the full cycle of the decreasing fractured well production stage of claim 4, wherein the method for determining the dynamic production over the full cycle of the decreasing fractured well production collects data on the fracturing a well production of a tight reservoir versus time.
6. The method of claim 5, wherein the method for determining the dynamic production rate of the fractured well over the period of the decreasing production phase applies an equation
Fitting data of the output changing with time, solving the nonlinear least square problem by adopting a Levenberg-Marquardt method to obtain the initial decrement rate D of the undetermined parameteri0.001, late stage decreasing rate D∞0.016, 20.768 and 0.792.
7. The method for determining the dynamic production rate of a fractured well during the whole cycle during the production decline stage of the fractured well according to claim 5, wherein the method for determining the dynamic production rate of the fractured well during the whole cycle during the production decline stage obtains the initial decline rate D by using fittingi＝0.001Late stage decreasing rate D∞When a certain oil reservoir block of the same type and a B well of the same type are determined, the initial yield q is decreased gradually, wherein the time coefficient m is 20.768 and the decreasing index n is 0.792iIn the case of 14.78, the formula is applied
And calculating to obtain the yield at any moment of the decreasing stage, and finally determining the dynamic yield of the hydraulic fracturing well in the whole period of the decreasing stage.
8. A computer device, characterized in that the computer device comprises a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to carry out the steps of:
constructing a yield decrement rate of the unconventional reservoir fracturing well in the whole period;
integrating the constructed full-period yield decreasing rate by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the yield decreasing stage of the fracturing well;
collecting data on yield versus time;
fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well, solving the fitting problem by adopting a Levenberg-Marquardt method, and obtaining an initial decreasing rate, a later decreasing rate, a time coefficient and a decreasing index of undetermined parameters;
and calculating the yield at any moment of the decrement stage by using a yield calculation expression of the yield decrement stage of the fracturing well under the condition of determining the decrement initial yield value of a certain oil reservoir block of the same type and a well of the same type by using the initial decrement rate, the later decrement rate, the time coefficient and the decrement index of the undetermined parameters obtained by fitting, and finally determining the dynamic yield of the unconventional reservoir hydraulic fracturing well in the whole period of the yield decrement stage.
9. A dynamic production system for determining the dynamic production of a fractured well within the whole period of the descending stage of the production of the fractured well, which implements the method for determining the dynamic production of the fractured well within the whole period of the descending stage of the production of the fractured well according to any one of claims 1 to 7, wherein the system for determining the dynamic production of the fractured well within the whole period of the descending stage of the production of the fractured well comprises:
the decrement rate construction module is used for constructing the decrement rate of the production of the unconventional reservoir fracturing well in the whole period;
the yield expression obtaining module is used for integrating the constructed full-period yield decreasing rate by combining the definition of the yield decreasing rate to obtain a yield calculation expression of the fractured well in the yield decreasing stage;
the data collection module is used for collecting data of yield and time;
the parameter acquisition module is used for fitting data of yield changing along with time by using a yield calculation expression in a yield decreasing stage of the fracturing well, solving the fitting problem by adopting a Levenberg-Marquardt method and obtaining an initial decreasing rate, a later decreasing rate, a time coefficient and a decreasing index of undetermined parameters;
and the yield result output module is used for calculating the yield at any moment in the yield decreasing stage by using the yield calculation expression in the yield decreasing stage of the fracturing well under the condition of determining the initial yield value of the oil reservoir block of the same type and the well of the same type by using the initial yield rate, the later yield rate, the time coefficient and the decreasing index of the undetermined parameters obtained by fitting, and finally determining the dynamic yield in the whole period of the yield decreasing stage of the hydraulic fracturing well of the unconventional reservoir.
10. The unconventional reservoir hydraulic fracturing exploitation well yield prediction terminal is used for realizing the method for determining the dynamic yield of the fracturing well in the whole period in the stage of decreasing the yield of the fracturing well according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110038328.2A CN112836860B (en) | 2021-01-12 | 2021-01-12 | Method and system for determining dynamic yield of fractured well in whole period of yield decreasing stage |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110038328.2A CN112836860B (en) | 2021-01-12 | 2021-01-12 | Method and system for determining dynamic yield of fractured well in whole period of yield decreasing stage |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112836860A true CN112836860A (en) | 2021-05-25 |
| CN112836860B CN112836860B (en) | 2022-08-02 |
Family
ID=75927832
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110038328.2A Active CN112836860B (en) | 2021-01-12 | 2021-01-12 | Method and system for determining dynamic yield of fractured well in whole period of yield decreasing stage |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112836860B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113723846A (en) * | 2021-09-06 | 2021-11-30 | 中国石油天然气股份有限公司 | Oil and gas asset depreciation processing method and device |
| CN115660294A (en) * | 2022-12-22 | 2023-01-31 | 西南石油大学 | Horizontal well full life cycle EUR tracking evaluation method, equipment and readable storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106194154A (en) * | 2016-07-11 | 2016-12-07 | 西南石油大学 | A kind of unconventionaloil pool hides medium-term and long-term PRODUCTION FORECASTING METHODS |
| US20160376885A1 (en) * | 2015-06-23 | 2016-12-29 | Petrochina Company Limited | Method and Apparatus for Performance Prediction of Multi-Layered Oil Reservoirs |
| CN107563557A (en) * | 2017-09-05 | 2018-01-09 | 中国石油天然气股份有限公司 | Method and device for determining oil well yield reduction rate |
| CN109252855A (en) * | 2018-10-15 | 2019-01-22 | 中国石油天然气股份有限公司 | Method and device for determining final cumulative yield of gas well |
| CN109409559A (en) * | 2017-08-17 | 2019-03-01 | 中国石油天然气股份有限公司 | Method and device for determining oilfield output reduction rate |
-
2021
- 2021-01-12 CN CN202110038328.2A patent/CN112836860B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160376885A1 (en) * | 2015-06-23 | 2016-12-29 | Petrochina Company Limited | Method and Apparatus for Performance Prediction of Multi-Layered Oil Reservoirs |
| CN106194154A (en) * | 2016-07-11 | 2016-12-07 | 西南石油大学 | A kind of unconventionaloil pool hides medium-term and long-term PRODUCTION FORECASTING METHODS |
| CN109409559A (en) * | 2017-08-17 | 2019-03-01 | 中国石油天然气股份有限公司 | Method and device for determining oilfield output reduction rate |
| CN107563557A (en) * | 2017-09-05 | 2018-01-09 | 中国石油天然气股份有限公司 | Method and device for determining oil well yield reduction rate |
| CN109252855A (en) * | 2018-10-15 | 2019-01-22 | 中国石油天然气股份有限公司 | Method and device for determining final cumulative yield of gas well |
Non-Patent Citations (1)
| Title |
|---|
| 吴忠维: "特高含水期产量递减规律及递减率的计算公式", 《中文科技论文》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113723846A (en) * | 2021-09-06 | 2021-11-30 | 中国石油天然气股份有限公司 | Oil and gas asset depreciation processing method and device |
| CN113723846B (en) * | 2021-09-06 | 2024-07-12 | 中国石油天然气股份有限公司 | Oil-gas asset consumption-reducing processing method and device |
| CN115660294A (en) * | 2022-12-22 | 2023-01-31 | 西南石油大学 | Horizontal well full life cycle EUR tracking evaluation method, equipment and readable storage medium |
| CN115660294B (en) * | 2022-12-22 | 2023-03-10 | 西南石油大学 | Horizontal well full life cycle EUR tracking evaluation method, equipment and readable storage medium |
| WO2024130959A1 (en) * | 2022-12-22 | 2024-06-27 | 西南石油大学 | Method for tracking and evaluating eur throughout full life cycle of horizontal well, device and readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112836860B (en) | 2022-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116186634B (en) | Intelligent management system for construction data of building engineering | |
| CN112836860B (en) | Method and system for determining dynamic yield of fractured well in whole period of yield decreasing stage | |
| CN104808587B (en) | A kind of mobility statistical method based on machining apparatus running status | |
| CN104657791B (en) | A kind of wind farm group wind speed profile prediction technique based on correlation analysis | |
| WO2020192675A1 (en) | Productivity prediction method for fractured horizontal well in tight oil reservoir | |
| CN103268069A (en) | Adaptive prediction control method based on Hammerstein model | |
| CN104865827B (en) | Oil pumping unit oil extraction optimization method based on multi-working-condition model | |
| CN103336433A (en) | Back stepping based mixed adaptive predication control system and predication control method thereof | |
| CN106779131A (en) | Grease hole evaluates phase cumulative production Forecasting Methodology | |
| CN113236228B (en) | Method and system for rapidly predicting single well yield | |
| CN117743280B (en) | Intelligent management method for highway bridge construction data | |
| CN115169089A (en) | Wind power probability prediction method and device based on kernel density estimation and copula | |
| CN107605477B (en) | Method for determining breakthrough of oil gas extracted by underground pump and stabilization time of water content | |
| CN114117654A (en) | Horizontal well perforation optimization design method and device based on machine learning | |
| CN105069574A (en) | New method for analyzing business flow behavior similarity | |
| CN117236716A (en) | Diagnosis method and system for low-yield and low-efficiency reasons of coal-bed gas well | |
| CN103971015A (en) | Gypsum slurry stirring time analysis method | |
| CN103077435A (en) | SEC (Securities and Exchange Commission) index evaluation method based on combination weighting comprehensive evaluation model | |
| CN115932704A (en) | CVT error measurement method and system based on model and data dual drive | |
| CN115169731A (en) | Smart campus energy consumption prediction method, device, equipment and medium | |
| CN103761432B (en) | Method and device for analyzing board-level test data | |
| CN116136580A (en) | Line loss grading detection method and device based on sensing equipment | |
| CN102636341B (en) | Waveform fitting based retard detection method | |
| CN110455902A (en) | A kind of method that the more standard specimens of environment measuring are quickly calibrated | |
| CN113761743B (en) | Method, device, terminal equipment and storage medium for calculating regional shale gas reserves |
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 |