CN107506948B  Shale oil gas comprehensive yield analysis method based on dynamic drainage volume  Google Patents
Shale oil gas comprehensive yield analysis method based on dynamic drainage volume Download PDFInfo
 Publication number
 CN107506948B CN107506948B CN201710890029.5A CN201710890029A CN107506948B CN 107506948 B CN107506948 B CN 107506948B CN 201710890029 A CN201710890029 A CN 201710890029A CN 107506948 B CN107506948 B CN 107506948B
 Authority
 CN
 China
 Prior art keywords
 reservoir
 shale
 shale oil
 oil
 volume
 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
 239000003079 shale oil Substances 0.000 title claims abstract description 204
 238000004458 analytical method Methods 0.000 title claims abstract description 55
 238000004519 manufacturing process Methods 0.000 claims abstract description 167
 239000012530 fluid Substances 0.000 claims abstract description 44
 239000000463 material Substances 0.000 claims abstract description 26
 230000001186 cumulative Effects 0.000 claims abstract description 23
 238000009792 diffusion process Methods 0.000 claims abstract description 18
 238000011156 evaluation Methods 0.000 claims abstract description 12
 230000004048 modification Effects 0.000 claims abstract description 11
 238000006011 modification reaction Methods 0.000 claims abstract description 11
 239000010779 crude oil Substances 0.000 claims description 53
 239000003921 oil Substances 0.000 claims description 51
 238000004364 calculation method Methods 0.000 claims description 48
 230000014509 gene expression Effects 0.000 claims description 33
 230000035699 permeability Effects 0.000 claims description 27
 XLYOFNOQVPJJNPUHFFFAOYSAN water Substances data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTAwLjUwMScgeT0nMTcwJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+SDwvdGV4dD4KPHRleHQgeD0nMTI2LjExNCcgeT0nMTg2JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjI2cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+MjwvdGV4dD4KPHRleHQgeD0nMTM4JyB5PScxNzAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8L3N2Zz4K data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTMuMzEyOCcgeT0nNTMuNTkwOScgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToyM3B4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPkg8L3RleHQ+Cjx0ZXh0IHg9JzI4LjE1NjknIHk9JzYyLjg2MzYnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MTVweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID4yPC90ZXh0Pgo8dGV4dCB4PSczNS4wNDU1JyB5PSc1My41OTA5JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjIzcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPC9zdmc+Cg== O XLYOFNOQVPJJNPUHFFFAOYSAN 0.000 claims description 21
 239000011159 matrix material Substances 0.000 claims description 18
 239000007789 gas Substances 0.000 claims description 14
 230000001052 transient Effects 0.000 claims description 13
 230000015572 biosynthetic process Effects 0.000 claims description 12
 238000005755 formation reaction Methods 0.000 claims description 12
 238000000605 extraction Methods 0.000 claims description 11
 239000008398 formation water Substances 0.000 claims description 11
 230000002093 peripheral Effects 0.000 claims description 11
 238000007906 compression Methods 0.000 claims description 10
 230000000875 corresponding Effects 0.000 claims description 10
 238000000034 method Methods 0.000 claims description 9
 230000000704 physical effect Effects 0.000 claims description 8
 239000011435 rock Substances 0.000 claims description 8
 238000005457 optimization Methods 0.000 claims description 7
 238000011084 recovery Methods 0.000 claims description 7
 238000003795 desorption Methods 0.000 claims description 5
 230000003247 decreasing Effects 0.000 claims description 4
 238000010606 normalization Methods 0.000 claims description 4
 238000010224 classification analysis Methods 0.000 claims description 3
 238000005553 drilling Methods 0.000 claims description 3
 238000004090 dissolution Methods 0.000 claims 2
 239000004215 Carbon black (E152) Substances 0.000 claims 1
 150000002430 hydrocarbons Chemical class 0.000 claims 1
 230000001131 transforming Effects 0.000 abstract description 9
 238000002360 preparation method Methods 0.000 abstract 1
 239000000126 substance Substances 0.000 description 10
 238000009826 distribution Methods 0.000 description 7
 238000004422 calculation algorithm Methods 0.000 description 4
 239000011148 porous material Substances 0.000 description 4
 230000003750 conditioning Effects 0.000 description 3
 238000009795 derivation Methods 0.000 description 3
 230000001808 coupling Effects 0.000 description 2
 238000010168 coupling process Methods 0.000 description 2
 238000005859 coupling reaction Methods 0.000 description 2
 238000007405 data analysis Methods 0.000 description 2
 230000001360 synchronised Effects 0.000 description 2
 125000004122 cyclic group Chemical group 0.000 description 1
 238000011157 data evaluation Methods 0.000 description 1
 238000010586 diagram Methods 0.000 description 1
 230000000694 effects Effects 0.000 description 1
 238000000197 pyrolysis Methods 0.000 description 1
 230000036632 reaction speed Effects 0.000 description 1
 238000000611 regression analysis Methods 0.000 description 1
 238000006467 substitution reaction Methods 0.000 description 1
Classifications

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06Q—DATA PROCESSING SYSTEMS OR METHODS, SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL, SUPERVISORY OR FORECASTING PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL, SUPERVISORY OR FORECASTING PURPOSES, NOT OTHERWISE PROVIDED FOR
 G06Q10/00—Administration; Management
 G06Q10/06—Resources, workflows, human or project management, e.g. organising, planning, scheduling or allocating time, human or machine resources; Enterprise planning; Organisational models
 G06Q10/063—Operations research or analysis
 G06Q10/0631—Resource planning, allocation or scheduling for a business operation
 G06Q10/06315—Needsbased resource requirements planning or analysis

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06Q—DATA PROCESSING SYSTEMS OR METHODS, SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL, SUPERVISORY OR FORECASTING PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL, SUPERVISORY OR FORECASTING PURPOSES, NOT OTHERWISE PROVIDED FOR
 G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
 G06Q50/02—Agriculture; Fishing; Mining
Abstract
The invention discloses a shale oil gas comprehensive yield analysis method based on dynamic drainage volume, which comprises the steps of arranging preparation data; initially estimating the limit leakage flow volume and the effective crack half length, and calculating the dynamic leakage flow volume, the average reservoir pressure change and the pressure correction factor; drawing a yield normalized simulated pressure curve, calculating the fluid diffusion coefficient in a fracturing transformation area, drawing a square root time characteristic curve, calculating the effective fracture halflength, and determining the ultimate discharge volume and the final output degree according to a modified Duong method; if the difference between the calculated values of the effective crack half length and the limit leakage volume and the input value is large, substituting the calculated values again, and circulating until the convergence is consistent to obtain the accurate effective crack half length and the limit leakage volume; substituting the dynamic drainage volume formula into a macroscopic material balance model, and comparing to obtain the cumulative contribution of the size of the fracturing modification area to the shale oil and gas production; and determining the optimal fracture interval and the optimal well spacing of the fractured horizontal well. The invention improves the accuracy of shale oil gas yield dynamic analysis and evaluation.
Description
Technical Field
The invention relates to the technical field of oilfield development and optimization, in particular to a shale oil gas comprehensive yield analysis method based on a dynamic drainage volume.
Background
In the shale oil and gas development process, the historical yield analysis of a production well is an important means for acquiring formation information, evaluating the production quality and optimizing the oil and gas production. However, complex reservoir characteristics of shale oil and gas, dynamic desorption of oil and gas, micronano multiscale nonDarcy seepage, fracture closure and shale shrinkage in the production process and other complex factors cause the conventional oil and gas dynamic analysis method to be inaccurate and not applicable any more, thereby bringing difficult challenges for evaluating and optimizing shale oil and gas development by utilizing production data.
The main reasons are as follows: (1) the existing method is suitable for simple well types and conventional fracturing wells, and the production dynamics of the volume fracturing horizontal well adopted in shale oil and gas development is difficult to accurately represent. (2) The existing method does not have complex distribution attributes of coupling shale oil gas nanometerscale matrix pores and a multiscale fracture network. (3) The existing method does not couple the influence of shale oil pyrolysis, adsorption and nonDarcy seepage characteristics on the production well yield. (4) Existing methods are based on the assumption of steadystate flow of reservoir fluids; the shale permeability is extremely low, the steadystate flow stage is difficult to achieve in a short time, and the whole production life is in the transient flow stage, so that the existing method is not suitable any more. (5) The prior method does not consider the influence of dynamic damage change of reservoir rock and cracks on yield in the development process; the shale oilgas well yield decreases rapidly and the pressure failure is rapid, so that the shale oilgas reservoir has large rock shrinkage and crack closure damage which are not negligible. (6) The existing method is difficult to synchronously optimize and predict based on historical data, and achieves the coupling analysis of realtime evaluation and prediction. Therefore, the conventional oil and gas dynamic analysis method is not suitable for shale oil and gas, and even brings large errors, thereby bringing difficult challenges for effectively and optimally developing shale oil and gas.
Disclosure of Invention
Aiming at the technical defects of the existing analysis method, the invention provides a shale oil and gas comprehensive yield analysis method based on dynamic drainage volume, which fully considers the characteristics of shale oil and gas development complex reservoirs and fluid characteristics, and comprises the influences of extremely low permeability of shale nanopore matrix, complex fracture network distribution of a fractured horizontal well, oil and gas dynamic desorption, limited fracture transformation range, shale shrinkage deformation, fracture closure and the like on the shale oil and gas development, so that the shale oil and gas comprehensive yield analysis method is a comprehensive yield analysis method integrating historical evaluation, development optimization and yield prediction, and the shale oil and gas yield dynamic analysis and evaluation accuracy is improved.
In order to achieve the purpose, the invention adopts the following technical scheme:
a shale oil and gas comprehensive yield analysis method based on a dynamic leakage flow volume comprises the following steps:
(I): preparing data arrangement; the data includes: shale reservoir petrophysical properties, fluid physical properties, well parameters, and production history of the target well;
the shale reservoir rock physical properties comprise compressibility, shale permeability and oil gas desorption parameters, the fluid physical properties comprise viscosity, volume coefficient and fluid compressibility, the well parameters comprise drilling and fracturing completion parameters, and the production history of the target well comprises oil gas water daily yield and cumulative yield.
(II): and initially estimating the limit leakage flow volume and the effective fracture halflength to obtain a calculation expression of the saturation and the dynamic leakage flow volume of the shale oilgas well, substituting the calculation expression and the production history of the target well into a macroscopic material balance model, and calculating the dynamic leakage flow volume, the average reservoir pressure change and a pressure correction factor in the production process.
(III): drawing a relation curve of the yield normalized simulated pressure and time, identifying the flowing stage of the shale fractured horizontal well, and calculating the fluid diffusion coefficient in the fractured modified area according to the ending time of the early transient linear flow and the known fracture interval; then drawing a square root time characteristic curve to obtain the slope of the square root time characteristic curve, and calculating the effective crack halflength; and carrying out production data decreasing analysis according to the modified Duong method, and determining the limit leakage flow volume and the final output degree;
(IV): and (3) if the difference between the effective crack half length and the limit leakage volume calculated in the step (three) and the input value in the step (two) is larger, substituting the value calculated in the step (three) into the step (two), repeating the step (two) and the step (three) in sequence until the input value in the step (two) is consistent with the convergence of the effective crack half length and the limit leakage volume calculated in the step (three), and stopping calculating to obtain the accurate effective crack half length and the limit leakage volume.
(V): the shale oil and gas flow area comprises two parts: and (3) calculating and expressing the dynamic drainage volume in the shale oil and gas production process of the volume fractured horizontal well in the fracturing modified area and the nonmodified matrix area, substituting the dynamic drainage volume into a macroscopic material balance model, and comparing to obtain the cumulative contribution of the fracturing modified area and the nonmodified matrix area to the shale oil and gas production.
(VI): and (4) performing optimization determination on fracture parameters and well spacing of the shale oil and gas fracturing horizontal well after analysis, evaluation and prediction based on the shale oil and gas production data.
Preferably, in the step (ii), the limit drainage volume and the effective fracture halflength are initially estimated to obtain a shale oil and gas well saturation calculation expression (1) and a dynamic drainage volume calculation expression (2) in the shale oil and gas production process applicable to the volume fracture horizontal well, then the calculation expressions and the production history of the target well are substituted into the macroscopic material balance model (3), and the dynamic drainage volume, the reservoir average pressure change and the pressure correction factor in the production process are calculated according to the production history of the shale oil and gas well.
A shale oilgas well saturation calculation expression:
wherein, subscript n represents parameter values at different production moments; subscript i represents the values of the parameters at the initial reservoir conditions; RF is shale oil production degree and fraction at different production moments; b is_{oi}Is the initial reservoir pressure p_{i}Volume coefficient of shale oil of m^{3}/m^{3}；B_{ti}Is the initial reservoir pressure p_{i}Fluid volume coefficient of shale oil reservoir, m^{3}/m^{3}；S_{wi}The original water saturation and fraction of the shale oil reservoir; v_{p(n)}Is production of t_{n}Dynamic drainage volume m of shale oil well at any moment^{3}；B_{o}Is the average pressure condition of the reservoirLower crude oil volume coefficient, m^{3}/m^{3}；S_{w(n)}，S_{o(n)}，S_{g(n)}At a dynamic discharge volume V_{p(n)}Formation crude oil, water and gas saturations in; (ii) a B is_{t(n)}Is t_{n}Fluid volume coefficient, m, of shale reservoir at that moment^{3}/m^{3}；c_{t(n)}Is t_{n}Integral compression coefficient of shale oil reservoir at any moment, Pa^{1}；V_{p},_{max}Is the ultimate drainage volume, m, of a shale oil and gas production well^{3}。
The dynamic drainage volume calculation expression is applicable to the shale oil and gas production process of the volume fracturing horizontal well:
a. and (3) keeping the shale gas well production bottom flow pressure constant:
wherein，x_{f}Is the effective fracture halflength; t different production moments of the fractured horizontal well, day; t is t_{e}Is the time of the end of the early linear flow of the fractured horizontal well, h is the shale reservoir thickness, m, η_{o}Is the diffusion coefficient of the peripheral unmodified area of the shale reservoir, m^{2}/s； η_{I}Is the diffusion coefficient of the fracturereformed region, m^{2}/s；y_{e}Is the crack half spacing, m; v_{p}Is the dynamic drainage volume, m, of shale oil and gas production wells^{3}。
Macroscopic material balance model:
wherein, V_{n}Is t_{n}Dynamic drainage volume m of fractured horizontal well of shale oil reservoir at moment^{3}(ii) a The subscript i represents the original reservoir pressure p_{i}Each parameter value under the condition; subscript n represents the shale reservoir mean pressure p_{n}Conditioning each parameter value; phi is a_{i}Is the average porosity, fraction, of the shale reservoir under initial conditions; c. C_{f}Shale oil reservoir rock compressibility factor, Pa^{1}；B_{o}，B_{w}，B_{g}Is the volume coefficient of crude oil, formation water and gas of the shale oil reservoir, m^{3}/m^{3}；R_{v}Is the condensate gasoil ratio of shale reservoir, m^{3}/m^{3}；R_{s}Shale oil reservoir crude oil dissolved gasoil ratio, m^{3}/m^{3}；S_{w}，S_{o}，S_{g}Is the formation crude oil, water and gas saturations; n is a radical of_{op},N_{wp},N_{gp}Is the cumulative oil production, water production and gas production of the shale oil well at different moments m^{3}。
According to the invention, a shale oilgas well saturation degree calculation expression (1) is substituted into a material balance equation (3), production prediction and historical analysis are coupled into a whole, and an iterative algorithm is used for synchronous calling calculation, so that the integration of 'historical evaluationscheme optimizationproduction prediction' is realized.
The method also considers that the permeability of the shale oil and gas reservoir nanoscale pores is extremely low, and the transient flow of the volume fracture horizontal well can penetrate through the whole production stage, so that the traditional leakage flow volume calculation method based on steadystate linear flow or radial flow is not applicable any more. The method establishes the shale oilgas fracturing horizontal well linear flow model, deduces the formula of the dynamic leakage flow volume model (2) in the shale oilgas production process which is accurately suitable for the volume fracturing horizontal well, combines the formula (2) and the formula (3), and can represent the contribution of the dynamic leakage flow volume change and the fracturing transformation range to the shale oilgas production based on historical data.
Preferably, in the step (three), because the shale oil and gas production process has rapid pressure failure, obvious matrix shrinkage and crack closure and obvious influence of stronger compressibility of shale fluid on the shale oil and gas yield, the method is different from the traditional yield analysis parameter definition method, the pseudo pressure m (p) considering permeability timevarying property and fluid compressibility, namely the formula (4a), is introduced, and the pressure correction factor f is introduced_{cp}I.e., (4b), improves the desired parameters in the yield analysis curve: yield normalized pseudo pressure (4c) equation and production material equilibrium time (4d) equation:
wherein m (p) is a pseudo pressure, Pa, introduced in consideration of timevarying permeability and compressibility of the fluid; m (RNP) is the shale reservoir pseudo pressure after yield normalization, Pa/(m)^{3}/d)；m(P_{wf}) Shale oil reservoir production bottom hole pressure p_{wf}A corresponding pseudo pressure value, Pa; m (p)_{i}) Initial pressure p of shale oil reservoir_{i}A corresponding pseudo pressure value, Pa; k is a radical of_{i}The matrix permeability mD under the condition of the average pressure p of the shale oil reservoir; k is a radical of_{ro}Is the average pressure p and the average saturation S of the shale oil reservoir_{o}Relative permeability of the oil phase under the condition, and fraction; b is_{o}Is the volume coefficient of crude oil m under the condition of the average pressure p of a shale oil reservoir^{3}/m^{3}；μ_{o}The method comprises the steps of (1) determining the viscosity of crude oil in Pa.s under the condition of shale oil reservoir average pressure p; f. of_{cp}Is a pressure correction factor, fraction, gamma_{k}Is the reservoir permeability modulus, Pa^{1}；c_{f}Is a reservoir under the condition of average pressure p of a shale oil reservoirCoefficient of compression, Pa^{1}；c_{μ}Is a crude oil viscosity factor Pa under the condition of shale oil reservoir average pressure p^{1}；q_{o,sc}For crude oil production at surface conditions, m^{3}A day; n is a radical of_{op}Is the cumulative oil production m of the shale oil well at different times t^{3}；t_{ma}Considering the substance balance time corresponding to different moments t after the influence of daily output change of the shale oil well is considered, and measuring the substance balance time;
plotting normalized pseudopressure m (RNP) and time t_{ma}Performing instantaneous yield Analysis (RNP), identifying the flowing stage of the shale fractured horizontal well according to different slopes of a fitted regression line, and calculating the fluid diffusion coefficient in the fractured and transformed region by using a formula (2) according to the ending time te of the early Transient linear flow (the slope is 1/2) and the known fracture spacing; then drawing square root time characteristic curve (5) to obtain its slope mL, then according to Wattenbarger (6) formula calculating effective crack halflength x_{f}And carrying out production data descending analysis according to the modified Duong method, and determining the limit leakage flow volume and the final output degree.
Wherein m is_{L}Is a transient yield analysis (m (RNP) vs. t_{ma}) The slope of the relationship; mu.s_{o}Is crude oil viscosity, pa.s; h is_{f}Is the shale oil well fracture height, m; x is the number of_{f}Is the crack length, m; q. q.s_{o,sc}Is the crude oil yield m under the surface condition of the shale oil well^{3}A day; k is shale reservoir matrix permeability, mD; phi is shale reservoir porosity, fraction; c. C_{t}Is the overall compression coefficient, Pa, of the shale oil reservoir^{1}。
Preferably, in the step (III), the limit drainage volume V is determined first due to the calculation of the crude oil recovery factor in the formula (1)_{pmax}Which can be obtained by modifying the formula of Duong method (7)(from Duong, 2010). The specific calculation procedure includes first determining coefficients a and m using equation (7a), and then plotting 7 to determine the initial production q based on equation (7b)_{1}。
Wherein q is_{o}Is the daily oil production m of the shale oil well at different times t^{3}A day; n is a radical of_{op}Is the cumulative oil production m of the shale oil well at different times t^{3}；q_{o，1}Is the initial oil production of the shale oil well, m^{3}A day; and a and m are yield curve fitting coefficients.
Preferably, in the step (six), based on analysis, evaluation and prediction of shale oil and gas production data, in order to eliminate interference caused by length difference of horizontal sections of different production wells, a horizontal section is introduced to normalize the final production degree, namely the ratio of the RF to the length of the horizontal section, namely the final production degree of a unit horizontal section; determining the optimal fracture interval of the fractured horizontal well by normalizing the horizontal segment lengths of all target production wells and classifying and analyzing the final extraction degree and the fracture interval; and meanwhile, determining the optimal fractured horizontal well spacing by dividing the predicted ultimate drainage volume by the known horizontal section length.
The beneficial effect of the invention is that,
1. comprehensive 'historical evaluationdevelopment optimizationyield prediction' integration
Different from the existing dynamic analysis method, the shale oilgas well saturation degree calculation expression (1), namely a relational expression between fluid saturation degreereservoir pressureextraction degree, is introduced into the macroscopic material balance model (3) for the first time, and the integration of 'historical evaluationscheme optimizationproduction prediction' is realized by using an iterative algorithm for synchronous adjustment and calculation. Wherein, the history evaluation part: drawing a production flow stage identification curve and a square root time characteristic curve; and a future prediction part: predicting the final recovery rate and the ultimate discharge volume of shale oil and gas by adopting a corrected Duong curve; and an optimization design part: and obtaining the optimal fracture spacing and the optimal fractured horizontal well spacing by summarizing and analyzing all production well evaluation results of the target block. The invention mutually calls the input and output results of each part of the evaluationoptimizationprediction process, and realizes the comprehensive dynamic analysis of shale oil and gas production.
2. Introducing a dynamic drainage volume concept to represent the cumulative contribution of the fracturing modification range to the shale oil production
The permeability of the nanoscale pores of the shale oil and gas reservoir is extremely low, and the volume fracturing horizontal well with complex fracture distribution causes transient flow to possibly penetrate through the whole production stage, so that the traditional leakage flow volume calculation method based on steadystate linear flow or radial flow is not applicable any more. According to the method, the shale oilgas fracturing horizontal well linear flow model is established, the dynamic drainage volume model which is accurately suitable for the shale oilgas production process of the volume fracturing horizontal well is deduced, and the contribution of the fracturing modification range and the peripheral nonmodification range to the shale oilgas production can be represented respectively.
3. Shale matrix shrinkage and fracture closure influence are coupled, and iterative algorithm accurately represents fracture parameters
In order to consider the influences of matrix shrinkage and obvious fracture closure caused by fast pressure failure in the shale oilgas production process and stronger compressibility of shale fluid on the shale oilgas yield, the invention defines the simulated pressure introduced by considering permeability timevarying property and fluid compressibility, introduces the pressure correction primer, improves the yield normalization simulated pressure and the production material balance time in the shale oil yield analysis curve, and can obtain more accurate effective fracture length of the shale oil volume fracturing horizontal well.
4. Based on historical data analysis, shale oil and gas well arrangement and fracturing optimization are synchronously realized, and yield prediction is realized
Because the limit leakage flow volume is obtained firstly when calculating the shale oil dynamic production degree, and the limit leakage flow volume is obtained by predicting through yield decrement analysis. Therefore, the method couples shale oil and gas historical production data analysis with future production prediction by introducing an iterative algorithm. By adopting a modified Duong method, the final recovery rate and the ultimate drainage volume of shale oil and gas are predicted, then the final recovery rate and the ultimate drainage volume are substituted into a dynamic drainage volume model in the shale oil and gas production process to carry out yield dynamic analysis and cyclic iteration, the accurate ultimate drainage volume and the recovery degree are finally obtained, the shale oil and gas yield dynamic analysis and evaluation accuracy is improved, and technical support is provided for optimizing and developing shale oil and gas.
Drawings
FIG. 1 is a detailed flow chart of shale oil reservoir production data evaluationoptimization designfuture prediction comprehensive analysis according to the invention;
FIG. 2 is a graph of results from determining different flow stages of a shale oil volume fractured horizontal well of a real block using RNP curve analysis according to an embodiment of the present invention;
FIG. 3 is a flow chart of decreasing analysis and prediction of a shale oil volume fractured horizontal well of a real block by using a modified Duong method according to an embodiment of the invention
FIG. 4 is a graph illustrating the contribution of the unmodified area and the fracture modified area of a shale oil volume fractured horizontal well of a certain practical block to the production of the horizontal well according to the embodiment of the invention;
FIG. 5 is a graph of the outcome of determining the optimal fracture spacing for a classification analysis of a real block of shale oil target production wells according to an embodiment of the present invention;
FIG. 6 is a schematic of the trilinear flow of a shale reservoir volume fractured horizontal well of the present invention;
FIG. 7 is a graph of the results of comparison of the dynamic drainage volume calculation analysis model and the CMG calculation results under the constant shale oil well production borehole bottom flow pressure according to the present invention;
FIG. 8 is a graph of a comparison result between a dynamic drainage volume calculation analysis model and a CMG calculation result under a constant shale oil well production speed condition according to the present invention;
FIG. 9 is a schematic diagram of a dynamic drainage volume change process during a longterm transient flow of a shale oil fractured horizontal well according to the present invention.
Detailed Description
The present invention will be further described with reference to the following examples.
A shale oil and gas comprehensive yield analysis method based on a dynamic leakage flow volume comprises the following steps:
(I): preparing data arrangement; the required data includes: shale reservoir petrophysical properties, fluid physical properties, well parameters, and production history of the target well;
the shale reservoir rock physical properties comprise compressibility, shale permeability and oil gas desorption parameters, the fluid physical properties comprise viscosity, volume coefficient and fluid compressibility, the well parameters comprise drilling and fracturing completion parameters, and the production history of the target well comprises oil gas water daily yield and cumulative yield.
Taking an example reservoir of Niobrara shale oil in the United states as an example, the data is collated as shown in Table 1.
TABLE 1 Niobrara shale oil reservoir and fluid parameters for a certain example of U.S.A
(II): initially estimating the limit leakage flow volume and the effective fracture halflength, respectively substituting into a shale oilgas well saturation calculation expression (1) and a dynamic leakage flow volume calculation expression (2) suitable for a volume fracture horizontal well in a shale oilgas production process, then substituting the calculation expressions and the production history of the target well into a macroscopic material balance model (3), calculating and determining the dynamic leakage flow volume of the shale oil well and the average reservoir pressure change in the production process according to the production history of the shale oilgas well, and determining a pressure correction factor (4b) formula.
A shale oilgas well saturation calculation expression:
wherein, subscript n represents parameter values at different production moments; subscript i represents the values of the parameters at the initial reservoir conditions; RF is shale oil production degree and fraction at different production moments; b is_{oi}Is the initial reservoir pressure p_{i}Volume coefficient of shale oil of m^{3}/m^{3}；B_{ti}Is the initial reservoir pressure p_{i}Shale oil reservoirVolume coefficient of fluid, m^{3}/m^{3}；S_{wi}The original water saturation and fraction of the shale oil reservoir; v_{p(n)}Is t_{n}Dynamic drainage volume m of shale oil well at any moment^{3}；B_{o}Is the average pressure condition of the reservoirLower crude oil volume coefficient, m^{3}/m^{3}；S_{w(n)}，S_{o(n)}，S_{g(n)}At a dynamic discharge volume V_{p(n)}Formation crude oil, water and gas saturations in; v_{p,max}Ultimate discharge volume, m^{3}；B_{t(n)}Is t_{n}Fluid volume coefficient, m, of shale reservoir at that moment^{3}/m^{3}；c_{t(n)}Is t_{n}Integral compression coefficient, Pa, of shale reservoir at time^{1}；V_{p,max}Is the ultimate drainage volume, m, of a shale oil and gas production well^{3}。
The dynamic drainage volume calculation expression is applicable to the shale oil and gas production process of the volume fracturing horizontal well:
a. and (3) keeping the shale gas well production bottom flow pressure constant:
wherein x is_{f}Is the effective fracture halflength of the shale well; t is the different production times of the fractured horizontal well, day; t is t_{e}The time of the completion of the early linear flow of the shale fractured horizontal well, i.e., the time of the completion of the early linear flow of the shale fractured horizontal well, h is the thickness of a shale reservoir, m, η_{o}Is the diffusion coefficient of the peripheral unmodified area of the shale reservoir, m^{2}/s；η_{I}Is the diffusion coefficient of the fracturereformed region, m^{2}/s；y_{e}Is the crack half spacing, m; v_{p}Is the dynamic drainage volume, m, of the shale oil and gas well^{3}。
Macroscopic material balance model:
wherein, V_{n}Is t_{n}Dynamic drainage volume m of fractured horizontal well of shale oil reservoir at moment^{3}(ii) a The subscript i represents the original reservoir pressure p_{i}Each parameter value under the condition; subscript n represents the shale reservoir mean pressure p_{n}Conditioning each parameter value; phi is a_{i}Is the average porosity, fraction, of the shale reservoir under initial conditions; c. C_{f}Shale oil reservoir rock compressibility factor, Pa^{1}；B_{o}，B_{w}，B_{g}Is the volume coefficient of crude oil, formation water and gas of the shale oil reservoir, m^{3}/m^{3}；R_{v}Is the condensate gasoil ratio of shale reservoir, m^{3}/m^{3}；R_{s}Shale oil reservoir crude oil dissolved gasoil ratio, m^{3}/m^{3}；S_{w}，S_{o}，S_{g}Is the formation crude oil, water and gas saturations; n is a radical of_{op},N_{wp},N_{gp}Is the cumulative oil production, water production and gas production of the shale oil well at different moments m^{3}。
(III): and (3) drawing a yield normalized simulated pressure and production time relation curve, as shown in figure 2, identifying the flow stage of the shale fractured horizontal well according to different slope sizes of the fitted regression line, and according to the ending time t of the early transient linear flow (the slope is 1/2)_{e}And knowing the crack spacing, and calculating the fluid diffusion coefficient in the fracture transformation area by using the formula (2); then, the square root time characteristic curve (5) is drawn to obtain the slope m_{L}Then, according to the formula of Wattenbarger (6), the effective halflength of the crack is calculated, and the production data is subjected to descending analysis according to the modified Duong method, as shown in figure 3(c), so as to determine the ultimate leakage flow volume and the final output degree. Different from the definition method of parameters used by the traditional oil deposit yield analysis method, due to the fast pressure failure, obvious matrix shrinkage and crack closure in the shale oil and gas production process and the influence of stronger compressibility of shale fluid on the shale oil and gas yield, pseudo pressure m (p) considering permeability timevarying property and fluid compressibility is introduced, as shown in formula (4a), and a pressure correction factor f is introduced_{cp}I.e. (4b)) Formula (xxxvii), parameters required in improving RNP yield analysis curves: yield normalized pseudo pressure (4c) equation and production material equilibrium time (4d) equation:
wherein m (p) is a pseudo pressure, Pa, introduced in consideration of timevarying permeability and compressibility of the fluid; m (RNP) is the shale reservoir pseudo pressure after yield normalization, Pa/(m)^{3}/d)；m(P_{wf}) Shale oil reservoir production bottom hole pressure p_{wf}A corresponding pseudo pressure value, Pa; m (p)_{i}) Initial pressure p of shale oil reservoir_{i}A corresponding pseudo pressure value, Pa; k is a radical of_{i}The matrix permeability mD under the condition of the average pressure p of the shale oil reservoir; k is a radical of_{ro}Is the average pressure p and the average saturation S of the shale oil reservoir_{o}Relative permeability of the oil phase under the condition, and fraction; b is_{o}Is the volume coefficient m of crude oil under the condition of average pressure p of shale oil reservoir^{3}/m^{3}；μ_{o}The method comprises the steps of (1) determining the viscosity of crude oil in Pa.s under the condition of shale oil reservoir average pressure p; f. of_{cp}Is a pressure correction factor, fraction, gamma_{k}Is the reservoir permeability modulus, Pa^{1}；c_{f}Is the reservoir compression coefficient Pa under the condition of the average pressure p of the shale oil reservoir^{1}；c_{μ}Is a crude oil viscosity factor Pa under the condition of shale oil reservoir average pressure p^{1}；q_{o,sc}Is the crude oil yield under surface conditions, m^{3}A day; n is a radical of_{op}Is the cumulative oil production m of the shale oil well at different times t^{3}；t_{ma}Considering the substance balance time corresponding to different moments t after the influence of daily output change of the shale oil well is considered, and measuring the substance balance time;
drawing a relation curve of yield normalized simulated pressure and production time, performing instantaneous yield Analysis (RateTransient Analysis, RNP), identifying the flowing stage of the shale fractured horizontal well according to different slopes of a fitted regression line, identifying the flowing stage of the shale fractured horizontal well according to different slopes of the fitted regression line, and calculating a fluid diffusion coefficient in a fracturing modification area (SRV) by using a formula (2) according to the ending time te of early Transient linear flow (the slope is 1/2) and the known fracture interval; then theDrawing square root time characteristic curve (5) to obtain its slope m_{L}Then, according to the formula of Wattenbarger (6), calculating the effective halflength x of the crack_{f}And carrying out production data descending analysis according to the modified Duong method to determine the limit leakage flow volume and the final output degree.
Wherein m is_{L}Is a transient yield analysis (m (RNP) vs. t_{ma}) The slope of the relationship; mu.s_{o}Is crude oil viscosity, pa.s; h is_{f}Is the shale oil well fracture height, m; x is the number of_{f}Is the effective fracture length of the shale production well, m; q. q.s_{o,sc}Is the crude oil yield m under the surface condition of the shale oil well^{3}A day; k is shale reservoir matrix permeability, mD; phi is shale reservoir porosity, fraction; c. C_{t}Is the comprehensive compression coefficient, Pa, of shale oil reservoir^{1}。
Since the calculation of the crude oil recovery in equation (1) requires first determining the ultimate drainage volume V_{pmax}It can be determined by modifying the Duong method (7) equation, (from Duong, 2010). The specific calculation process comprises the steps of firstly determining coefficients a and m by using the formula (7a), fitting production data to obtain values of parameters a and m as shown in the figure 3(a), and then drawing and determining an initial yield q according to the formula (7b)_{1}Fitting the production data acquisition parameters q as in FIG. 3(b)_{1}，。
Wherein q is_{o(n)}Is the daily oil production m of the shale oil well at different times t^{3}A day; n is a radical of_{op(n)}Is the cumulative oil production m at different times t of the shale oil well^{3}；q_{o1}Is the initial oil production of the shale oil well, m^{3}A day; and a and m are yield curve fitting coefficients.
(IV): and (3) if the difference between the effective crack half length and the limit leakage volume calculated in the step (III) and the input value in the step (II) is larger, substituting the value calculated in the step (III) into the step (II), repeating the step (II) and the step (III) in sequence until the input value in the step (II) is consistent with the convergence of the effective crack half length and the limit leakage volume calculated in the step (III), stopping calculation, and obtaining the accurate effective crack half length and the limit leakage volume.
(V): according to the formulas shown in fig. 6 and (2), the shale oil and gas flow area comprises two parts: a fracture engineered zone and an unmodified matrix zone. Substituting the formula (2) into a macroscopic material balance model (4), and comparing to obtain the cumulative contribution of the size of the fracturing modification area to the shale oil and gas production. As shown in fig. 4, in a certain practical shale oil calculation example, the contribution of the unmodified area to the 1000day cumulative yield of the volume horizontal well is not more than 5%, that is, the size of the range of the practical fracture modified area basically determines most of the yield of the volume fracture horizontal well.
(VI): by integrating the production data evaluation and prediction results of all shale wells, as shown in table 2, in order to eliminate the interference caused by the length difference of the horizontal sections of different production wells, the horizontal section is introduced to normalize the final extraction degree (the ratio of the RF to the length of the horizontal section, i.e., the final extraction degree of the unit horizontal section). By normalizing the horizontal segment lengths of all target production wells to obtain the final extraction degree and the fracture spacing, classification analysis can be carried out, and the optimal fracture spacing of the fractured horizontal well can be determined to be about 20m, as shown in fig. 5. And then dividing the predicted ultimate drainage volume by the known horizontal segment length and the reservoir height of 15m to obtain the optimal fractured horizontal well spacing of about 180270m by estimation.
TABLE 2. Niobrara shale oil certain practical block volume fracturing horizontal well production history evaluation and prediction result summarization
The invention relates to a determination of a dynamic leakage flow volume calculation expression in the shale oil and gas production process applicable to a volume fracture horizontal well:
the model derivation hypothesis is that ① shale oil reservoir outer boundary is closed, a homogeneous box shape is formed, physical properties do not change along with time, a single infinite largeflowconductivity vertical fracture is fractured on a horizontal well section, a ② finiteflowconductivity vertical artificial fracture completely penetrates through the reservoir, the fracture height is equal to the oil reservoir thickness, ③ reservoir fluid can only flow into a borehole from a fracture of a perforation section, pressure loss of the ④ horizontal well section is ignored, ⑤ hydraulic support fractures are symmetric doublewing fractures, are perpendicular to a horizontal shaft, fracture interference is considered, noflow boundaries exist among fractures, and the model deduced by ⑥ is suitable for singlephase fluid flow.
As shown in fig. 6, the subscript "O" represents the fracture horizontal well peripheral unmodified zone, "I" represents the fracture modified zone, and "f" represents the fracture modified primary fracture. To facilitate the derivation of the formula for calculating the leakage area, the following dimensionless variables are defined:
dimensionless pressure:
dimensionless time:
dimensionless distance:
in the shale fracturing modified area and the peripheral nonmodified area, the fluid flow diffusion coefficient is as follows:
the fluid in the shale oil reservoir has dimensionless flow capacity:
wherein k is_{I},k_{O}Respectively the permeability, mD, of a shale oil reservoir fracturing modification area and a peripheral nonmodification area; h is shale reservoir height, m; q is the daily yield of crude oil from shale reservoir, m^{3}A day; c. C_{t}Comprehensive compression coefficient, Pa, of shale oil reservoir^{1}；p_{i}Is the original reservoir pressure, Pa; x is the number of_{f}Shale oil reservoir fracturing horizontal well fracture length m η_{o}Is the diffusion coefficient of the peripheral unmodified area of the shale reservoir, m^{2}/s；η_{I}Is the diffusion coefficient of the fracturereformed region, m^{2}/s；y_{e}Is the crack half spacing, m.
From the above model assumptions, the following set of seepage equations can be written for FIG. 6:
in the peripheral fracture nonmodified area, the shale reservoir fluid flows into the internal fracture modified area linearly, the flow equation is as follows,
through Laplace transformation and combining boundary conditions, a mathematical model can be obtained:
the solution of the equation is that,
similar to the above, the mathematical model of seepage in the Laplace space of the internal fracture transformation area is listed as,
wherein the content of the first and second substances,
1. under the production condition of fixed bottom hole flowing pressure:
based on the above assumptions, only after the fluid flow pressure boundary in the fracture reformation region reaches the noflow boundary between the primary fractures, i.e., time t_{eD}The flow is started only in the outer unmodified range, while the bottom hole flow pressure is determined according to the obtained conditions of (Behmanesh et al, 2015) at t_{eD}At the moment, the pressure distribution in the fracture transformation range is,
at this time, the pressure conditions at the nonflowing boundary between fractures and along the periphery of the main fracture are,
p＝p_{wf}→p_{D}＝1，aty_{D}＝0 (20)
combining the expressions (19) to (20) of the boundary conditions and the expression (18) of the initial condition, the solution of the expression of the nonhomogeneous partial differential equation (14) is,
wherein the content of the first and second substances,
b_{0}＝b_{2}＝0
C_{2}＝1b_{0}C_{1}
then, according to the equation (15), the pressure distribution in the outer unmodified region is,
2. under the condition of definite production yield
Similarly, the solution obtained according to (Behmanesh et al, 2015) at t_{eD}At the moment, the pressure distribution in the fracture transformation range is,
the pressure conditions at the nonflowing boundaries between fractures and along the main fracture perimeter are,
at y_{D}＝0，t_{D}>0 (24)
At y_{D}＝y_{eD}，t_{D}>0 (25)
Combining the expressions (24) to (25) of the boundary conditions and the expression (23) of the initial condition, the solution of the expression of the nonhomogeneous partial differential equation (14) is,
wherein the content of the first and second substances,
b_{3}＝b_{5}＝0
according to equation (15), the pressure distribution of the dynamic drainage volume of the outer unmodified region is expressed as,
3. determining the dynamic drainage volume of the volume fractured horizontal well:
according to fig. 6, the change of the leakage flow area in the shale oil and gas production process can be divided into two processes: in a fracturing modification area, the drainage area is continuously increased along the direction (y direction) vertical to the main crack until a nonflowing boundary in the middle of the crack is reached; then, the fluid in the peripheral nonmodified area starts to flow into the fracture modified area along the x direction, and the drainage area starts to continuously expand and increase along the x direction. Thus, the blowoff volume of a volume fractured horizontal well may be expressed as:
V_{p}fracture transformation range + outer unmodified zone dynamic drainage volume (28)
According to the method for calculating the maximum reaction speed of the pressure wave adopted by the leakage flow boundary by Kuchuk (2009), the invention also determines the dynamic leakage flow volume expansion speed by calculating the leakage flow volume at the moment when the dimensionless pressure and the time second derivative are equal to zero, and the method comprises the following steps:
in consideration of the complexity of the quadratic derivation calculation of the expressions (22) and (27), it is difficult to directly obtain an analytical expression. Therefore, the invention adopts a multivariate regression analysis method to obtain a calculation method of the dynamic leakage flow volume by fitting.
a. Constant shale oil well production borehole bottom flow pressure:
wherein x is_{f}Is the effective fracture halflength of the shale fractured horizontal well; t is the different production times of the fractured horizontal well, day; t is t_{e}Is the time of the end of the early linear flow of the fractured horizontal well, h is the shale reservoir thickness, m, η_{o}Is the diffusion coefficient of the peripheral unmodified area of the shale reservoir, m^{2}/s；η_{I}Is the diffusion coefficient of the fracturereformed region, m^{2}/s；y_{e}Is the crack half spacing, m; v_{p}Is the dynamic drainage volume of the shale oil and gas production well, m^{3}。
In order to verify the accuracy of the formula (2) of the dynamic leakage flow volume calculation, a volume fracturing horizontal well calculation example is set by using a commercial numerical simulator CMG black oil model. As shown in fig. 7 and 8, the dynamic leakage flow volume obtained by the formula (2) calculation has good fitting consistency with the CMG calculation result, which shows that the dynamic leakage flow volume calculation expression newly derived in the shale oil and gas production process applicable to the volume fractured horizontal well is accurate enough.
A macroscopic material balance model for a volume fractured horizontal well based on dynamic drainage volume:
in the shale oil and gas production process, the dynamic drainage volume of the volume fracturing horizontal well changes, as shown in fig. 9, at the production time t_{n}Volume of the relief volume of the volume fractured horizontal well is V_{n}At the next production time t_{n+1}Dynamic drainage volume increase to V_{n+1}And the change of reservoir fracturing and fluid saturation in the dynamic process can be represented by using a macroscopic material balance relation.
In the shale oil and gas well production process, the material balance equation of the formation water is as follows:
wherein the content of the first and second substances,is t_{n}At time of dynamic drainage volume V_{(n)}The amount of formation water inventory within,is t_{n+1}t_{n}The volume of formation water increases over a period of time due to the expansion of the blowoff volume,time dynamic discharge volume V_{(n+1)}The amount of formation water therein. To the right of the equation at t_{n+1}t_{n}Cumulative water production from the well over the time period.
In the shale oilgas well production process, the material balance equation of shale oil is as follows:
wherein the content of the first and second substances,is t_{n}At time of dynamic drainage volume V_{(n)}The inventory of crude oil in the formation within,is t_{n+1}Time dynamic discharge volume V_{(n+1)}The reservoir of formation oil in the right side of the equation is at t_{n+1}t_{n}Cumulative oil production from the well over the time period. If the shale reservoir is not a condensate reservoir, the condensate oil content in the equation can be ignored
In the shale oilgas well production process, the material balance equation of shale gas is as follows:
wherein the content of the first and second substances,represents the amount of dissolved gas released from the crude oil as the reservoir pressure drops, which is zero if the reservoir pressure is above the crude oil saturation pressure.
Adding and summing the material balance relation equations of oil, gas and water to obtain a macroscopic material balance equation in the shale oilgas production process:
wherein, V_{n}Is t_{n}Dynamic drainage volume m of fractured horizontal well of shale oil reservoir at moment^{3}(ii) a The subscript i represents the original reservoir pressure p_{i}Each parameter value under the condition; subscript n represents the shale reservoir mean pressure p_{n}Conditioning each parameter value; phi is a_{i}Is the average porosity, fraction, of the shale reservoir under initial conditions; c. C_{f}Shale oil reservoir rock compressibility factor, Pa^{1}；B_{o}，B_{w}，B_{g}Is the volume coefficient of crude oil, formation water and gas of the shale oil reservoir, m^{3}/m^{3}；R_{v}Is the condensate gasoil ratio of shale reservoir, m^{3}/m^{3}；R_{s}Shale oil reservoir crude oil dissolved gasoil ratio, m^{3}/m^{3}；S_{w}，S_{o}，S_{g}Is the formation crude oil, water and gas saturations; n is a radical of_{op},N_{wp},N_{gp}Is the cumulative oil production, water production and gas production of the shale oil well at different moments m^{3}。
In the production process of shale oil gas, the production data of crude oil, water and gas, the change of the dynamic drainage volume of the shale oil gas well can be obtained by calculating according to the formula (30) and substituted into the formula (34), and the formula is combined with S_{w}+S_{O}+S_{g}The average formation pressure of the shale reservoir and the corresponding crude oil, gas and formation water saturations can be found as 1. In addition, according to the definition of the extent of production of a generic reservoir (Craft et al, 1991), the extent of dynamic production of shale reservoirs can be characterized as:
the shale oil reservoir extraction degree calculation comprises two parts: crude oil production efficiency and drainage volume sweep efficiency. Different from general oil reservoirs, the method introduces the drainage volume sweep efficiency, namely the ratio of the dynamic drainage volume to the limit drainage volume, and quantifies the longterm transient flow influence of the shale oil well. Wherein the limiting leakage volume is determined by a yielddecreasing analysis method. Meanwhile, the oil reservoir crude oil, formation water and gas saturation can be linked with the shale oil reservoir dynamic production degree by using the formula (36), as follows:
finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention, and not for limiting the scope of the present invention, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions can be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention.
Claims (5)
1. A shale oil gas comprehensive yield analysis method based on dynamic leakage flow volume is characterized by comprising the following steps:
firstly, preparing data arrangement; the data includes: shale reservoir petrophysical properties including compressibility, shale permeability and hydrocarbon desorption parameters, fluid physical properties including viscosity, volume coefficient and fluid compressibility, well parameters including drilling and fracturing completion parameters, and production history of a target well including oil, gas and water daily yield and cumulative yield;
initially estimating the limit drainage volume and the effective fracture halflength to obtain a shale oilgas well saturation and dynamic drainage volume calculation expression, substituting the calculation expression and the production history of the target well into a macroscopic material balance model, and calculating the dynamic drainage volume, the reservoir average pressure change and a pressure correction factor in the production process;
drawing a relation curve of the yield normalized simulated pressure and time, identifying the flowing stage of the shale fractured horizontal well, and calculating the fluid diffusion coefficient in the fractured modified area according to the ending time of the early transient linear flow and the known fracture interval; then drawing a square root time characteristic curve to obtain the slope of the square root time characteristic curve, and calculating the effective crack halflength; and carrying out production data decreasing analysis according to the modified Duong method, and determining the limit leakage flow volume and the final output degree;
if the difference between the effective crack half length and the limit leakage volume calculated in the step three and the input value in the step two is large, the value calculated in the step three needs to be substituted into the step two again, the step two and the step three are repeated in sequence until the effective crack half length and the limit leakage volume calculated in the step two are converged and consistent, the calculation is stopped, and the accurate effective crack half length and the accurate limit leakage volume are obtained;
fifthly, the shale oil and gas flow area comprises two parts: a fracture engineered zone and an unmodified matrix zone; substituting a dynamic leakage flow volume calculation expression in the shale oil and gas production process suitable for the volume fracturing horizontal well into a macroscopic material balance model, and comparing to obtain the cumulative contribution of a fracturing modification area and an unmodified matrix area to shale oil and gas production respectively;
and sixthly, performing optimization determination on fracture parameters and well spacing of the shale oil and gas fracturing horizontal well based on analysis, evaluation and prediction of shale oil and gas production data.
2. The shale oil and gas comprehensive yield analysis method based on the dynamic leakage flow volume as claimed in claim 1, wherein in the second step, firstly, the limit leakage flow volume and the effective fracture halflength are initially estimated to obtain a shale oil and gas well saturation degree calculation expression (1) and a dynamic leakage flow volume calculation expression (2) in the shale oil and gas production process suitable for the volume fracture horizontal well, then the shale oil and gas saturation degree calculation expression and the shale oil and gas production process suitable for the volume fracture horizontal well are substituted into a macroscopic material balance model (3), and the dynamic leakage flow volume, the reservoir average pressure change and the pressure correction factor in the production process are calculated according to the production history of the shale oil and gas well;
a shale oilgas well saturation calculation expression:
wherein, subscript n represents parameter values at different production moments; subscript i represents the values of the parameters at the initial reservoir conditions; RF (radio frequency)_{(n)}Shale oil extraction degrees and fractions at different production moments; b is_{oi}Is the initial reservoir pressure p_{i}Volume coefficient of shale oil of m^{3}/m^{3}；B_{ti}Is the initial reservoir pressure p_{i}Fluid volume coefficient of shale oil reservoir, m^{3}/m^{3}；S_{wi}The original water saturation and fraction of the shale oil reservoir; v_{p(n)}Is production of t_{n}Dynamic drainage volume m of shale oil well at any moment^{3}；B_{o}Is the reservoir mean pressure condition P_{n}Lower crude oil volume coefficient, m^{3}/m^{3}；S_{w(n)}，S_{o(n)}，S_{g(n)}At a dynamic discharge volume V_{p(n)}Formation crude oil, water and gas saturations in; b is_{t(n)}Is t_{n}Fluid volume coefficient, m, of shale reservoir at that moment^{3}/m^{3}；c_{t(n)}Is t_{n}Integral compression coefficient, Pa, of shale oil reservoir at any moment^{1}；V_{p,max}Is the ultimate drainage volume, m, of the shale production well^{3}，B_{o}(n) is the shale reservoir mean pressure p_{n}Volume coefficient of shale reservoir crude oil, m^{3}/m^{3}；
Shale suitable for volume fracturing horizontal well; the dynamic leakage flow volume calculation expression in the oil and gas production process is as follows:
a. and (3) keeping the shale gas well production bottom flow pressure constant:
wherein x is_{f}Is the effective fracture halflength, m, of the shale fractured horizontal well; t is the different production times of the fractured horizontal well, day; t is t_{e}Is the time of the end of the early linear flow of the fractured horizontal well, h is the shale reservoir thickness, m, η_{o}Is the diffusion coefficient of the peripheral unmodified area of the shale reservoir, m^{2}/s；η_{I}Is the diffusion coefficient of the fracturereformed region, m^{2}/s；y_{e}Is the crack half spacing, m; v_{p}Is the dynamic drainage volume, m, of a shale oil and gas production well^{3}；
Macroscopic material balance model:
wherein, V_{(n)}Is t_{n}Dynamic drainage volume m of fractured horizontal well of shale oil reservoir at moment^{3}；φ_{(i)}Is the average porosity, fraction, of the shale reservoir under initial conditions; c. C_{(f)}Shale oil reservoir rock compressibility factor, Pa^{1}；P_{n}Namely the average pressure of the oil reservoir, Pa; p_{i}To the original reservoir pressureForce, Pa;the formation crude oil saturation under the average pressure of the shale oil reservoir;is the volume coefficient of formation water, m, at the average pressure of the reservoir^{3}/m^{3}；Gas saturation under the average pressure of the shale oil reservoir;is the volume coefficient of gas, m, at the average pressure of the reservoir^{3}/m^{3}；R_{v(n)}Is the condensate oilgas ratio, m, of the shale oilgas reservoir under the average pressure of the shale oil reservoir^{3}/m^{3}；Water saturation under the average pressure of the shale oil reservoir;is the volume coefficient m of the crude oil of the shale oil reservoir under the average pressure of the shale oil reservoir^{3}/m^{3}；R_{s(n)}The shale oil reservoir crude oil dissolution gasoil ratio m under the average pressure of the shale oil reservoir^{3}/m^{3}；S_{w(i)}The formation crude oil saturation under the original reservoir pressure of the shale; b is_{w(i)}Is the volume coefficient of formation water, m, at the original reservoir pressure^{3}/m^{3}；S_{o(i)}Water saturation at original reservoir pressure; b is_{o(i)}Is the volume coefficient of shale oil reservoir crude oil under the original reservoir pressure, m^{3}/m^{3}；R_{s(i)}Is the shale oil reservoir crude oil dissolution gasoil ratio m under the original oil reservoir pressure^{3}/m^{3}；S_{g(i)}Gas saturation at the original reservoir pressure; b is_{g(i)}Is the volume factor of the gas at the original reservoir pressure,m^{3}/m^{3}；R_{v(i)}is the condensate gasoil ratio, m, of the shale oilgas reservoir at the original reservoir pressure^{3}/m^{3}；N_{wp(n)}The cumulative water yield m of the shale oil well at different moments under the average pressure of the shale oil reservoir^{3}；N_{op(n)}The cumulative oil production m of the shale oil well at different moments under the average pressure of the shale oil reservoir^{3}；N_{gp(n)}The cumulative oil production m of the shale oil well at different moments under the average pressure of the shale oil reservoir^{3}。
3. The method for analyzing shale oil and gas comprehensive yield based on dynamic drainage volume as claimed in claim 2, wherein in the third step, because the shale oil and gas production process has rapid pressure failure, obvious matrix shrinkage and crack closure, and strong compressibility of shale fluid has obvious influence on shale oil and gas yield, the method is different from the traditional yield analysis parameter definition method, the method introduces pseudo pressure m (p) considering permeability timevarying property and fluid compressibility, namely formula (4a), and introduces a pressure correction factor f_{cp}I.e., (4b), improves the desired parameters in the yield analysis curve: yield normalized pseudo pressure (4c) equation and production material equilibrium time (4d) equation:
wherein m (p) is a pseudo pressure, Pa, introduced in consideration of timevarying permeability and compressibility of the fluid; m (RNP) is the shale reservoir pseudo pressure after yield normalization, Pa/(m)^{3}/d)；m(P_{wf}) Shale oil reservoir production bottom hole pressure p_{wf}A corresponding pseudo pressure value, Pa; m (p)_{i}) Initial pressure p of shale oil reservoir_{i}A corresponding pseudo pressure value, Pa; k is a radical of_{i}The matrix permeability mD under the average pressure condition of the shale oil reservoir; k is a radical of_{ro}Is the average pressure and average saturation S of shale oil reservoir_{o}Relative permeability of the oil phase under the condition, and fraction; b is_{o}Is the volume coefficient m of crude oil under the average pressure condition of the shale oil reservoir^{3}/m^{3}；μ_{o}The viscosity of crude oil is Pa.s under the condition of average pressure of the shale oil reservoir; f. of_{cp}Is pressureForce correction factor, fraction, gamma_{k}Is the reservoir permeability modulus, Pa^{1}；c_{f}Is the reservoir compression coefficient Pa of the shale oil reservoir under the condition of average pressure^{1}；c_{μ}Is a crude oil viscosity factor Pa under the condition of average pressure of shale oil reservoir^{1}；q_{o,sc}For crude oil production at surface conditions, m^{3}A day; n is a radical of_{op}Is the cumulative oil production m of the shale oil well at different times t^{3}；t_{ma}Considering the influence of daily output change of the shale oil well, the material balance time, day and mu corresponding to different moments t_{oi}Is the crude oil viscosity, Pa.s, under the original reservoir pressure condition; b is_{oi}Means volume coefficient of shale reservoir crude oil under original reservoir pressure condition, m^{3}/m^{3}；
Plotting normalized pseudopressure m (RNP) and time t_{ma}Performing instantaneous yield Analysis (RateTransient Analysis, namely RNP), identifying the flow stage of the shale fractured horizontal well according to different slopes of a fitted regression line, and calculating the fluid diffusion coefficient in the fractured reconstruction region by using a formula (2) according to early Transient linear flow, namely the end time te with the slope of 1/2 and the known fracture spacing; then drawing square root time characteristic curve (5) to obtain its slope mL, then according to Wattenbarger (6) formula calculating effective crack halflength x_{f}Carrying out production data decreasing analysis according to a modified Duong method, and determining the limit leakage flow volume and the final output degree;
where the index i represents the value of each parameter under the original reservoir pressure conditions, m_{L}Is the slope of the instantaneous yield analysis relationship curve; mu.s_{o}Is crude oil viscosity, pa.s; h is_{f}Is the shale oil well fracture height, m; x is the number of_{f}Is the crack length, m; q. q.s_{o,sc}Is the crude oil yield under the surface condition of the shale oil well，m^{3}A day; k is shale reservoir matrix permeability, mD; phi is shale reservoir porosity, fraction; c. C_{t}Is the overall compression coefficient, Pa, of the shale oil reservoir^{1}。
4. The method for analyzing shale oil and gas comprehensive yield based on dynamic drainage volume as claimed in claim 3, wherein in the third step, the limit drainage volume V is determined first due to the calculation of the crude oil recovery factor in the formula (1)_{pmax}It can be determined by modifying the Duong method (7) by a calculation process comprising first determining the coefficients a and m using equation (7a) and then plotting and determining the initial yield q according to equation (7b)_{o,1}：
Wherein q is_{o}Is the daily oil production m of the shale oil well at different times t^{3}A day; n is a radical of_{op}Is the cumulative oil production m of the shale oil well at different times t^{3}；q_{o,1}Is the initial oil production of the shale oil well, m^{3}A day; and a and m are yield curve fitting coefficients.
5. The method according to claim 1, wherein in the sixth step, based on analysis, evaluation and prediction of shale oil and gas production data, in order to eliminate interference caused by length difference of horizontal sections of different production wells, a horizontal section is introduced to normalize the final extraction degree, i.e. the ratio of the RF to the length of the horizontal section, i.e. the ratio of the shale oil extraction degree to the length of the horizontal section at different production moments, and the fraction is the final extraction degree of the unit horizontal section; determining the optimal fractured horizontal well fracture interval by normalizing the horizontal segment lengths of all target production wells to obtain the final extraction degree and the fracture interval in a classification analysis manner; and meanwhile, determining the optimal fractured horizontal well spacing by dividing the predicted ultimate drainage volume by the known horizontal segment length.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201710890029.5A CN107506948B (en)  20170927  20170927  Shale oil gas comprehensive yield analysis method based on dynamic drainage volume 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201710890029.5A CN107506948B (en)  20170927  20170927  Shale oil gas comprehensive yield analysis method based on dynamic drainage volume 
Publications (2)
Publication Number  Publication Date 

CN107506948A CN107506948A (en)  20171222 
CN107506948B true CN107506948B (en)  20200421 
Family
ID=60699605
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201710890029.5A Active CN107506948B (en)  20170927  20170927  Shale oil gas comprehensive yield analysis method based on dynamic drainage volume 
Country Status (1)
Country  Link 

CN (1)  CN107506948B (en) 
Families Citing this family (6)
Publication number  Priority date  Publication date  Assignee  Title 

CN108547612A (en)  20180302  20180918  中国石油天然气股份有限公司  Shale oil converted insitu output oil mass, the determination method and apparatus of output tolerance 
CN108468540B (en) *  20180329  20210409  西南石油大学  Calculation method for single well drainage area of lowpermeabilityultralowpermeability oil and gas reservoir 
CN110130884A (en) *  20190428  20190816  中国石油天然气股份有限公司  A kind of gas well residual production potential determines method 
CN110084431B (en) *  20190430  20201211  四川洁能锐思石油技术有限公司  Shale gas well yield analysis and prediction method and system 
CN112308264B (en) *  20190726  20210831  中国石油天然气股份有限公司  Method and device for predicting optimal development mode of shale oil insitu conversion 
CN111502652B (en) *  20200701  20210423  西南石油大学  Yield decreasing and production dynamic prediction method for threehole medium gas reservoir horizontal well 
Citations (2)
Publication number  Priority date  Publication date  Assignee  Title 

CN105243182A (en) *  20150908  20160113  中国石油天然气股份有限公司  Method for calculating dynamic reserve of dense oil fractured horizontal well 
CN105625990A (en) *  20141030  20160601  中国石油化工股份有限公司  Method for mining shale oil by virtue of horizontal well through performing volume fracturing on reservoirs 

2017
 20170927 CN CN201710890029.5A patent/CN107506948B/en active Active
Patent Citations (2)
Publication number  Priority date  Publication date  Assignee  Title 

CN105625990A (en) *  20141030  20160601  中国石油化工股份有限公司  Method for mining shale oil by virtue of horizontal well through performing volume fracturing on reservoirs 
CN105243182A (en) *  20150908  20160113  中国石油天然气股份有限公司  Method for calculating dynamic reserve of dense oil fractured horizontal well 
NonPatent Citations (1)
Title 

页岩气藏多重介质耦合流动模型;苏玉亮;《天气工业》;20160228;全文 * 
Also Published As
Publication number  Publication date 

CN107506948A (en)  20171222 
Similar Documents
Publication  Publication Date  Title 

CN107506948B (en)  Shale oil gas comprehensive yield analysis method based on dynamic drainage volume  
CN107577831B (en)  Method for calculating scale of karst cave of fracturecavity carbonate oilgas reservoir  
CN105386751A (en)  Well logging and productivity prediction method of horizontal well based on oil reservoir seepage flow model  
CN105569646A (en)  Oil and gas well technical recoverable reserves prediction method  
CN105893679B (en)  Well Test Data Analysis Method is corrected in low yield horizontal well afterflow  
RU2723769C1 (en)  Method of calculating volume of reverse flow of fluid for hydraulic fracturing of formation during hydraulic fracturing in horizontal wells in gas deposits of fractured sandstones  
CN105631184B (en)  A kind of oil/gas well Production Decline Analysis method and system  
CN111062129A (en)  Shale oil complex seam network discrete fracture continuous medium mixed numerical simulation method  
CN111502652A (en)  Yield decreasing and production dynamic prediction method for threehole medium gas reservoir horizontal well  
Chen  Production from a fractured well with finite fracture conductivity in a closed reservoir: an exact analytical solution for pseudosteadystate flow  
CN108894777B (en)  Method for determining physical properties and fracture characteristic parameters of reservoir of splitlayer fractured multilayer commingled production hydrocarbon reservoir  
Zhang et al.  An evaluation method of volume fracturing effects for vertical wells in low permeability reservoirs  
CN106932324B (en)  A kind of method of determining High water cut sandstone oil reservoir reservoir permeability changing rule  
CN107247816B (en)  Method for judging farthest construction radius of cold recovery chemical viscosity reduction for heavy oil reservoir  
CN107313759B (en)  Hypotonic heavy crude reservoir straight well thermal recovery pressure distribution forecasting method and system  
CN108518218B (en)  Unconventional oil and gas reservoir multisection fracturing horizontal well single well dynamic reserve determination method  
CN111581786B (en)  Well test interpretation method of well test interpretation model for analyzing fracturecave seriesmode doublehole composite reservoir  
CN107355200A (en)  One kind receives micron particles dispersion improving ecology well choosing method  
CN108664677A (en)  A kind of oil/gas well creation data analysis method  
RU2737043C1 (en)  Method for development of oil reservoir of multilayer oil and gas condensate deposit  
CN107704646B (en)  Modeling method after compact reservoir volume modification  
CN110656915B (en)  Shale gas multisection fracturing horizontal well multiworkingsystem productivity prediction method  
CN108468540B (en)  Calculation method for single well drainage area of lowpermeabilityultralowpermeability oil and gas reservoir  
CN113586044A (en)  Optimization method and system for gas testing working system of selfblowing shale gas  
CN109441415B (en)  Well testing interpretation method of polymer flooding oil reservoir test well based on adjacent well interference 
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  
TA01  Transfer of patent application right 
Effective date of registration: 20191227 Address after: 266580 No. 66 Changjiang West Road, Qingdao, Shandong Applicant after: China University of Petroleum (East China) Address before: 266000 internal entrepreneurship center, No. 66, No. 66 oil university, Changjiang West Road, Shandong Applicant before: As the Joint Research Institute of energy technology (Qingdao) Co. Ltd. Applicant before: China University of Petroleum (East China) 

TA01  Transfer of patent application right  
GR01  Patent grant  
GR01  Patent grant 