To quantitatively evaluate the complexities and heterogeneities of pore structures in sandstone reservoirs, we apply single fractal theory and multifractal theory to explore the fractal characteristics of pore size distributions based on mercury intrusion porosimetry. The fractal parameters were calculated and the relationships between the petrophysical parameters (permeability and entry pressure) and the fractal parameters were investigated. The results show that the single fractal curves exhibit two-stage characteristics and the corresponding fractal dimensions D₁ and D₂ can characterize the complexity of pore structure in different sizes. Favorable linear relationships between log(ε) and log(μ,(ε)) indicate that the samples satisfy multifractal characteristics and ε is the sub-intervals with size ε = J × 2^−k. The multifractal singularity curves used in this study exhibit a right shape, indicating that the heterogeneity of the reservoir is mainly affected by pore size distributions in sparse regions. Multifractal parameters, D(0), D(1), and Δf, are positively correlated with permeability and entry pressure, while D(0), D(1), and Δf are negatively correlated with permeability and entry pressure. The ratio of larger pores volumes to total pore volumes acts as a control on the fractal dimension over a specific pore size range, while the range of the pore size distribution has a definite impact on the multifractal parameters. Results indicate that fractal analysis and multifractal analysis are feasible methods for characterizing the heterogeneity of pore structures in a reservoir. However, the single fractal models ignore the influence of microfractures, which could result in abnormal values for calculated fractal dimension. Compared to single fractal analysis, multifractal theory can better quantitatively characterize the heterogeneity of pore structure and establish favorable relationships with reservoir physical property parameters.

Based on the property discontinuity in the radial direction, this paper develops a new composite model to simulate the productivity of volume fractured wells in naturally fractured reservoirs. The analytical solution of this model is derived in detail and its accuracy is verified by the same model’s numerical solution. Detailed analyses of the traditional transient and cumulative rate are provided for the composite model. Results show that volume fracturing mainly contributes to the early-flow period’s production rate. As interregional diffusivity ratio increases or interregional conductivity ratio decreases, the transient rate at the same wellbore pressure increases. Three characteristic decline stages may be observed on transient rate curves and the shape of traditional rate curves in the middle- and late-flow periods depends on naturally fractured medium and boundary condition. By introducing a new pseudo-steady constant, new Blasingame type curves are also established and their features are more salient than those of traditional rate curves. Five typical flow regimes can be observed on these new type curves. Sensitivity analysis indicates that new Blasingame type curves for varied interregional diffusivity ratio, interregional conductivity ratio, interporosity coefficient and dimensionless reservoir radius, except storativity ratio, will normalize in the late-flow period.

Multilayer commingled production is the most efficient development technique of coalbed methane under the condition of multiple coal seams. However, due to the differences in physical properties between multilayer superimposed gas-bearing systems, interlayer interference severely limits coalbed methane development in commingled production. To achieve multilayer-commingled production, interlayer interference must be reduced and the combination of production layers must be optimized. Physical simulations are an effective measure to achieve this goal. According to the characteristics of multiple thin interbeds, strong reservoir heterogeneity and interlayer pressure difference in the Surat basin, a physical model is established to simulate the multilayer-commingled production process of coalbed methane reservoirs and the gas production contribution, and a pressure change of each layer is analyzed. The greater the interlayer pressure difference, the more obvious the early backflow phenomenon of the low-pressure layer, the more obvious the difference of layered production contribution in the later stage, the lower the degree of commingled production and recovery, and the stronger the interlayer interference. In view of these, this study proposes a new experimental method named the succession production. The novelty of this method is to control the commingled production time, that is, the high-pressure layer is produced first, and the low-pressure layer is combined when the interlayer pressure is consistent. The results show that this method can eliminate the early backflow phenomenon of the low-pressure layer and reduce interlayer interference. Furthermore, the characteristics of interlayer interference and the change law of multilayer-commingled production capacity of succession and commingled production are clarified, providing theoretical and technical support for reducing interlayer interference and optimizing production layer combination to promote the efficient development of multiple thin interbedded coalbed methane reservoirs.

In order to provide paleofluid evidence of hydrocarbon accumulation periods in the Amu Darya Right Bank Block, microexperiments and simulations related to the Middle-Upper Jurassic Callovian-Oxfordian carbonate reservoirs were performed. On the basis of petrographic observation, the diagenetic stages were divided by cathodoluminescence, and the entrapment stages of fluid inclusions were divided by laser Raman experiment and UV epifluorescence. The hydrocarbon generation (expulsion) curve and burial (thermal) history curve of source rocks were simulated by using real drilling data coupled with geochemical parameters of source rocks, such as total organic carbon (TOC) and vitrinite reflectance (). The above results were integrated with microthermometry of fluid inclusions by inference the timing of hydrocarbon migration into the carbonate reservoirs. The horizon-flattening technique was used to process the measured seismic profile and restore the structural evolution profile. Four diagenetic periods and three hydrocarbon accumulation periods were identified. (i) For Syntaxial stage, the fluid captured by the overgrowing cement around particles is mainly seawater; (ii) for (Early) Mesogenetic burial stage, the calcite cements began to capture hydrocarbon fluids and show yellow fluorescence under UV illumination; (iii) for (Late) Mesogenetic burial stage, two sets of cleavage fissures developed in massive calcite cements, and oil inclusions with green fluorescence were entrapped in the crystal; (iv) for Telogenetic burial stage, blue fluorescent inclusions along with hydrocarbon gas inclusions developed in dully luminescent calcite veins. Based on the accurate division of hydrocarbon migration and charging stages, combined with the structural evolution history of the traps, the hydrocarbon accumulation model was established. Because two of the three sets of source rocks are of marine origin, resulting in the lack of vitrinite in the kerogen of those source rocks, there may be some deviation between the measured value of and the real value. Some systematic errors may occur in the thermal history and hydrocarbon generation (expulsion) history of the two sets of source rocks. Due to the limitations of seismic horizon-flattening technique—such as the inability to accurately recover the inclined strata thickness and horizontal expansion of strata—the final shape of the evolution process of structural profile may also deviate from the real state in geological history. The accumulation model established in this study was based upon the fluid inclusion experiments, which can effectively characterize the forming process of large condensate gas reservoirs in the Amu Darya Right Bank Block and quantify the timing of hydrocarbon charging. However, the hydrocarbon migration and accumulation model does not take the oil-source correlation into account, but only the relationship between the mature state of source rocks and the timing of hydrocarbon charging into the reservoirs. Subsequent research needs to conduct refined oil-source correlation to reveal the relationship between gas, condensate, source rocks, and recently discovered crude oil and more strictly constrain and modify the accumulation model, so as to finally disclose the origin of the crude oil and oil reservoir forming process in the Amu Darya Right Bank Block, evaluate the future exploration potential, and point out the direction of various hydrocarbon resources (condensate gas and crude oil).