Adsorbed gas may account for a significant part of the gas resources in shale gas and coalbed methane plays. Understanding gas sorption behaviors and integrating gas desorption into analytical reservoir modeling and an associated transient performance analysis are important for evaluating a system’s gas desorption ability and further analyzing its CO₂ injectability, utilization, and storage capacity. However, gas desorption, along with other pressure-dominated gas properties, increases a system’s non-linearity in theoretical studies. Few studies on analytical modeling have integrated the gas desorption feature into a non-linear system and validated the model’s accuracy. In this study, the desorbed gas due to pressure decay was treated as an additional source/sink term in the source-and-sink function methods. This method was combined with the integral image method in a semi-analytical manner to determine the amount of gas desorption. Fundamental reservoir and gas properties from the Horn River Basin shale gas play were chosen to evaluate the methodology and the performance of the associated production well. The results were compared with the commercial fine-gridding numerical simulation software, and good matches were achieved. The results showed that the desorbed gas released from rock will supply free-gas flow when the pressure significantly decreases due to gas production. The production wellbore pressure can be maintained at a higher level, and the production rate was higher than in cases where gas desorption was not considered, depending on the operating conditions.

The objective of this work is to innovatively apply the boundary element method (BEM) as a general modeling strategy to deal with complicated reservoir modeling problems, especially those related to reservoir heterogeneity and fracture systems, which are common challenges encountered in the practice of reservoir engineering. The transient flow behaviors of reservoirs containing multi-scale heterogeneities enclosed by arbitrarily shaped boundaries are modeled by applying BEM. We demonstrate that a BEM-based simulation strategy is capable of modeling complex heterogeneous reservoirs with robust solutions. The technology is beneficial in making the best use of geological modeling information. The governing differential operator of fluid flow within any locally homogeneous domain is solved along its boundary. The discretization of a reservoir system is only made on the corresponding boundaries, which is advantageous in closely conforming to the reservoir’s geological description and in facilitating the numerical simulation and computational efforts because no gridding within the flow domain is needed. Theoretical solutions, in terms of pressure and flow rate responses, are validated and exemplified for various reservoir–well systems, including naturally fractured reservoirs with either non-crossing fractures or crossing fractures; fully compartmentalized reservoirs; and multi-stage, fractured, horizontal wells with locally stimulated reservoir volumes (SRVs) around each stage of the fracture, etc. A challenging case study for a complicated fracture network system is examined. This work demonstrates the significance of adapting the BEM strategy for reservoir simulation due to its flexibility in modeling reservoir heterogeneity, analytical solution accuracy, and high computing efficiency, in reducing the technical gap between reservoir engineering practice and simulation capacity.

This study presents a comprehensive analytical modeling technology to model transient behaviors of multilayered reservoirs with inter-layer pure-planar crossflow induced by multi-stage hydraulically fractured horizontal well (MHFHW). The objective of this study is to develop an analytical model for multilayered reservoirs in conjunction with complex MHFHW and to achieve not only accurate and efficient computation, but also well-organized solutions expressed in a systematically integrated manner. The consideration of inter-layer crossflow across adjacent layers sets up the foundation for successful modeling of multilayered reservoirs. Source/sink function method (SSFM) is applied to describe fluid flow. Unsteady-state pressure or production rate solutions of MHFHW with the advantages of fast computation, accurate, and stable solutions are achieved. Comparative and consistent outcomes generated by this work and widely applied industry software have largely enhanced our technical confidence. More importantly, innovatively defined modified dimensionless terms that integrate systematic well-reservoir geometry information, as well as rock/fluid properties of each layer, have been newly applied to regulate the new modified dimensionless rate decline curve. This new technique sheds light on the reservoir characterization practice for complicated reservoir systems. Theoretical results in terms of transient pressure and rate were generated by the proposed multilayered model (SSFM-ML) for five scenarios of general concern, under various reservoir and well parameters, which were examined and discussed to demonstrate technical robustness. Not only does this study give solutions to the targeted multiple layered reservoirs, but it also provides insights into modeling three-dimensional fluid flow in heterogeneous reservoir with complex well configurations. It is recommended that future research should be conducted for more complicated two- and three-dimensional reservoirs, using the similar strategy of developing new type curves through adopting other new forms of modified dimensionless rate and time terms.

Using data from 56 tight gas wells from the study field (Y field) in the Ordos basin of China, this paper presents performance-based reservoir characterization of the study field from production data and geophysical data. Post-fracturing evaluation is realized by applying our new modified production decline type curves for fractured wells. Compared to traditional type curves, our newly proposed modified dimensionless type curves help identify field data diagnostics for various flow regimes of fractured wells and also facilitate the curve matching process with real data to obtain fruitful and crucial reservoir and fractured well information, including key parameters such as reservoir flowing capacity (kh), well productivity, fracture length, drainage area and original gas in place. This paper intends to promote the extensive application of this new technique. With the support of the reservoir information provided by production decline analysis using our modified type curves, the commercial flow units are delineated in terms of interrelated porosity-permeability of sandstone based on pore throat aperture crossplotting and corresponding flow unit productivity. Furthermore, two crossplots of well logging interpreted porosity versus resistivity are constructed, suggesting their good correlated relationships with relative flow unit productivity and initial gas abundance in place, respectively. The two crossplots enable qualitative evaluation of formation penetrated by well, which makes them very useful and practical as wireline logging is basically available for every well. The well production routine is also analyzed systematically by considering a well’s inflow performance, tubing performance and minimal liquid-carrying gas flow rate to investigate if a gas well is producing at optimal conditions or if a measure should be taken to improve the well’s production. Through analysis of the Y field, this study introduces an integrated workflow with the support of the new modified type curves to effectively help understand the reservoir characteristics and the flow behaviors of the tight gas field. The key takeaway from this study is that the new modified dimensionless production decline curves in terms of $q_{DM} vs. t_{DM}$ can be applied in field practice to achieve a systematically comparable understanding of the performance of MHFHWs globally.