We aim at the development of a general modelling workflow for design and optimization of the well flowback and startup operation on hydraulically fractured wells.Fracture flowback model developed earlier by the author...We aim at the development of a general modelling workflow for design and optimization of the well flowback and startup operation on hydraulically fractured wells.Fracture flowback model developed earlier by the authors is extended to take into account several new fluid mechanics factors accompanying flowback,namely,viscoplastic rheology of unbroken cross-linked gel and coupled“fracture-reservoir”numerical submodel for influx from rock formation.We also developed models and implemented new geomechanical factors,namely,(i)fracture closure in gaps between proppant pillars and in proppant-free cavity in the vicinity of the well taking into account formation creep;(ii)propagation of plastic deformations due to tensile rock failure from the fracture face into the fluid-saturated reservoir.We carried out parametric calculations to study the dynamics of fracture conductivity during flowback and its effect on well production for the set of parameters typical of oil wells in Achimov formation of Western Siberia,Russia.The first set of calculations is carried out using the flowback model in the reservoir linear flow regime.It is obtained that the typical length of hydraulic fracture zone,in which tensile rock failure at the fracture walls occurs,is insignificant.In the range of rock permeability in between 0.01 mD and 1 D,we studied the effect of non-dimensional governing parameters as well as bottomhole pressure drop dynamics on oil production.We obtained a map of pressure drop regimes(fast,moderate or slow)leading to maximum cumulative oil production.The second set of parametric calculations is carried out using integrated well production modelling workflow,in which the flowback model acts as a missing link in between hydraulic fracturing and reservoir commercial simulators.We evaluated quantitatively effects of initial fracture aperture,proppant diameter,yield stress of fracturing fluid,pressure drop rate and proppant material type(ceramic and sand)on long-term well production beyond formation linear regime.The third set of parametric calculations is carried out using the flowback model history-matched to field data related to production of four multistage hydraulically fractured oil wells in Achimov formation of Western Siberia,Russia.On the basis of the matched model we evaluated geomechanics effects on fracture conductivity degradation.We also performed sensitivity analysis in the framework of the history-matched model to study the impact of geomechanics and fluid rheology parameters on flowback efficiency.展开更多
The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm t...The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm to fixed-point(FP) iteration method. It computes the approximation to the solutions at each iteration based on the history of vectors in extended space, which includes the vector of unknowns, the discrete form of the operator, and the equation's right-hand side. Several constraints are applied to AA algorithm, including a limitation of the time step variation during the iteration process, which allows switching to the base FP iterations to maintain convergence. Compared to the base FP algorithm, the improved version of the AA algorithm enables a reliable and rapid convergence of the iterative solution for the quasi-linear elliptic pressure equation describing the flow of particle-laden yield-stress fluids in a narrow channel during hydraulic fracturing, a key technology for stimulating hydrocarbon-bearing reservoirs. In particular, the proposed AA algorithm allows for faster computations and resolution of unyielding zones in hydraulic fractures that cannot be calculated using the FP algorithm. The quasi-linear elliptic pressure equation under consideration describes various physical processes, such as the displacement of fluids with viscoplastic rheology in a narrow cylindrical annulus during well cementing,the displacement of cross-linked gel in a proppant pack filling hydraulic fractures during the early stage of well production(fracture flowback), and multiphase filtration in a rock formation. We estimate computational complexity of the developed algorithm as compared to Jacobian-based algorithms and show that the performance of the former one is higher in modelling of flows of viscoplastic fluids. We believe that the developed algorithm is a useful numerical tool that can be implemented in commercial simulators to obtain fast and converged solutions to the non-linear problems described above.展开更多
文摘We aim at the development of a general modelling workflow for design and optimization of the well flowback and startup operation on hydraulically fractured wells.Fracture flowback model developed earlier by the authors is extended to take into account several new fluid mechanics factors accompanying flowback,namely,viscoplastic rheology of unbroken cross-linked gel and coupled“fracture-reservoir”numerical submodel for influx from rock formation.We also developed models and implemented new geomechanical factors,namely,(i)fracture closure in gaps between proppant pillars and in proppant-free cavity in the vicinity of the well taking into account formation creep;(ii)propagation of plastic deformations due to tensile rock failure from the fracture face into the fluid-saturated reservoir.We carried out parametric calculations to study the dynamics of fracture conductivity during flowback and its effect on well production for the set of parameters typical of oil wells in Achimov formation of Western Siberia,Russia.The first set of calculations is carried out using the flowback model in the reservoir linear flow regime.It is obtained that the typical length of hydraulic fracture zone,in which tensile rock failure at the fracture walls occurs,is insignificant.In the range of rock permeability in between 0.01 mD and 1 D,we studied the effect of non-dimensional governing parameters as well as bottomhole pressure drop dynamics on oil production.We obtained a map of pressure drop regimes(fast,moderate or slow)leading to maximum cumulative oil production.The second set of parametric calculations is carried out using integrated well production modelling workflow,in which the flowback model acts as a missing link in between hydraulic fracturing and reservoir commercial simulators.We evaluated quantitatively effects of initial fracture aperture,proppant diameter,yield stress of fracturing fluid,pressure drop rate and proppant material type(ceramic and sand)on long-term well production beyond formation linear regime.The third set of parametric calculations is carried out using the flowback model history-matched to field data related to production of four multistage hydraulically fractured oil wells in Achimov formation of Western Siberia,Russia.On the basis of the matched model we evaluated geomechanics effects on fracture conductivity degradation.We also performed sensitivity analysis in the framework of the history-matched model to study the impact of geomechanics and fluid rheology parameters on flowback efficiency.
基金partial financial support from Gazpromneft Science and Technology Center。
文摘The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm to fixed-point(FP) iteration method. It computes the approximation to the solutions at each iteration based on the history of vectors in extended space, which includes the vector of unknowns, the discrete form of the operator, and the equation's right-hand side. Several constraints are applied to AA algorithm, including a limitation of the time step variation during the iteration process, which allows switching to the base FP iterations to maintain convergence. Compared to the base FP algorithm, the improved version of the AA algorithm enables a reliable and rapid convergence of the iterative solution for the quasi-linear elliptic pressure equation describing the flow of particle-laden yield-stress fluids in a narrow channel during hydraulic fracturing, a key technology for stimulating hydrocarbon-bearing reservoirs. In particular, the proposed AA algorithm allows for faster computations and resolution of unyielding zones in hydraulic fractures that cannot be calculated using the FP algorithm. The quasi-linear elliptic pressure equation under consideration describes various physical processes, such as the displacement of fluids with viscoplastic rheology in a narrow cylindrical annulus during well cementing,the displacement of cross-linked gel in a proppant pack filling hydraulic fractures during the early stage of well production(fracture flowback), and multiphase filtration in a rock formation. We estimate computational complexity of the developed algorithm as compared to Jacobian-based algorithms and show that the performance of the former one is higher in modelling of flows of viscoplastic fluids. We believe that the developed algorithm is a useful numerical tool that can be implemented in commercial simulators to obtain fast and converged solutions to the non-linear problems described above.
