Vertical Heterogeneity Research Articles

A Consistent Representation of Cloud Overlap and Cloud Subgrid Vertical Heterogeneity

AbstractMany global climate models underestimate the cloud cover and overestimate the cloud albedo, especially for low‐level clouds. We determine how a correct representation of the vertical structure of clouds can fix part of this bias. We use the 1D McICA framework and focus on low‐level clouds. Using Large Eddy Simulations results as reference, we propose a method based on exponential‐random overlap that represents the cloud overlap between layers and the subgrid cloud properties over several vertical scales, with a single value of the overlap parameter. Starting from a coarse vertical grid, representative of atmospheric models, this algorithm is used to generate the vertical profile of the cloud fraction with a finer vertical resolution, or to generate it on the coarse grid but with subgrid heterogeneity and cloud overlap that ensures a correct cloud cover. Doing so we find decorrelation lengths are dependent on the vertical resolution, except if the vertical subgrid heterogeneity and interlayer overlap are taken into account coherently. We confirm that the frequently used maximum‐random overlap leads to a significant error by underestimating the low‐level cloud cover with a relative error of about 50%, that can lead to an error of SW cloud albedo as big as 70%. Not taking into account the subgrid vertical heterogeneity of clouds can cause a relative error of 20% in brightness, assuming the cloud cover is correct.

Concept and key technology of “multi-scale high-density” fracturing technology: A case study of tight sandstone gas reservoirs in the western Sichuan Basin

There are abundant tight sandstone gas resources in the Sichuan Basin, which are the important objects of reserve and production increase and large-scale development. Due to their discontinuous sandbody distribution, narrow channels, and strong horizontal and vertical heterogeneity, however, conventional fracturing technologies cannot achieve the ideal stimulation effect here. In order to address this difficulty, this paper dissects the geology engineering characteristics of tight sandstone gas reservoirs in the western Sichuan Basin. Starting from the seepage mechanics theory, the concept of “multi-scale high-density” tight gas fracturing technology is put forward by fully referring to the experience of previous multi-round reservoir stimulation in the western Sichuan Basin and the idea of unconventional volume fracturing technology. In addition, its conceptual connotation, key technologies and implementation effects are illustrated. The following research results are obtained. First, the seepage characteristics make it necessary for the efficient production of tight gas reserves to increase fracture density and stimulated reservoir volume (SRV). Second, the “multi-scale high-density” fracturing technology emphasizes the rationality of high-density hydraulic fracture creation and the matching of multi-scale fracture flow capacity, and aims at establishing a multi-level fracture system with effective and steady gas flow in tight reservoirs through fracturing. Third, the “wide, dense, support, stable, and precise” fracturing technology is applied to improve single well production and estimated ultimate recovery (EUR). Fourth, the engineering practice of “multi-scale high-density” fracturing technology in the tight reservoirs of Jurassic Shaximiao Formation and Triassic Xujiahe Formation in the ZJ Gas Field realizes the average single well production rate of 15.6 × 104 m3/d, which is 1.96 times higher than that before the stimulation. Obviously, it provides powerful support for the operation of the ZJ Gas Field into a giant gas field with the reserves of 100 billion cubic meters. In conclusion, the formation of the concept and key technologies of “multi-scale high-density” fracturing technology effectively supports the efficient development of tight gas in the western Sichuan Basin and points out the following research direction of tight gas reservoir stimulation. The research results provide reference and guidance for the large-scale benefit development of tight gas in China.

Characterization of the layered SIF distribution through hyperspectral observation and SCOPE modeling for a subtropical evergreen forest

Solar-induced chlorophyll fluorescence (SIF) has been widely used as a prospective proxy of plant photosynthesis to accurately detect the response of vegetation to climate change and study terrestrial ecosystem carbon cycle. However, the challenge of estimating vegetation gross primary product (GPP) based on SIF observations remains unresolved due to the confounding effects of leaf physiology and canopy structure, especially the vertical heterogeneity of plant attributes. Radiative transfer models are ideally suited to investigate the variability in radiance signals including fluorescence, but previous studies lack vertically layered spectra to validate the model. In this study, we developed a vertically layered hyperspectral system and proposed an innovative framework for assessing vertical characteristics of SIF and improving the estimation of GPP based on SIF in a subtropical evergreen forest. The global sensitivity analysis on Soil-Canopy-Observation of Photosynthesis and Energy fluxes (SCOPE) model with vertical profiles showed that chlorophyll content (Cab), leaf inclination distribution (LIDFa), leaf area index (LAI), and senescent material (Cs) dominate the vertical variations in reflectance (Ref) and SIF. We found that the vertical characteristics of SIF were mainly impacted by these parameters of the adjacent vertical layer, except for the effects of LAI across the vertical layers on the observed understory SIFU, which were approximately equivalent. 7.8 ± 1.7% of the SIFU was transferred to the top of canopy, contributing up to 17% of the SIF observed on the top of canopy (SIFTOC) during the year. Furthermore, our findings suggest that substituting the observed SIFTOC with the total emitted SIF, particularly from the simulated overstory and understory SIF, could enhance the correlation between GPP and SIF with an increased R2 (ΔR2 = 0.09). This study highlights the importance of accounting for the contribution of the layered SIF in the quantification of the total SIF emission and can benefit the GPP estimation based on SIF signals in subtropical evergreen forests with high canopy density.

Prediction of Permian karst reservoirs in the Yuanba gas field, northern Sichuan Basin, China

Karst reservoir has great hydrocarbon potential, however karst reservoir prediction is inhibited by the strong lateral and vertical heterogeneity of karstification which in turn results in recognition difficulty when geological and geophysical methods are used in isolation. Combined geological and geophysical methods, including core observation, thin section analysis, well log interpretation, seismic attributes and seismic inversion are applied to understand the depositional environment, types of karstification, the geophysical response and distribution of karst reservoir and the controls over and evolution model of karst reservoir in Permian Yuanba gas field, northern Sichuan Basin. The results show that there are four types of major seismic facies recognized in the study area, which correspond to platform margin reef, carbonate platform, platform margin inter-bay and platform margin slope. Karstification can be divided into three zones: the supergene karst zone, the vertical seepage zone and the horizontal underflow zone, among which the supergene karst zone have the strongest karstification. Karstification is highlighted with low seismic impedance, low Poisson's ratio, and high seismic attenuation on seismic inversion as well as micro scale paleo geomorphology recognized by trend surface analysis. By comparing the distribution of karstification with paleo geomorphology, fault distribution and gas contents, it can be observed that the karstification have a good matching relationship with these factors. The dominant type of karstification is epigenic karst, which is strongly influenced by depositional facies and paleo geomorphology. Fracture networks have contributed to karstification but are not the dominant factor of karst formation. The integration of geological and geophysical methods can predict karst reservoir with high accuracy within large area and can be applied to karst reservoir hydrocarbon exploration of similar geological setting.

Heterogeneity of hydrological connectivity in coral limestone groundwater pool from vertical, spatial and temporal tracing of groundwater chemicals and isotopes

Carbonate aquifers are critical for the health and well-being of human communities. Despite considerable efforts on the structure of karst aquifers, the variable and invisible flow paths or the vertical and spatial heterogeneity of karst systems remain too complex to be fully understood. An artificially constructed small limestone aquifer with an extensive network of observation wells in the Miyako Island, Japan, provides a unique opportunity to illustrate the vertical and spatial hydrological connectivity for water movement and storage. Water chemical data at two different depths over six years, including electrical conductivity (EC), dissolved oxygen (DO), oxidation–reduction potential (ORP), major ions (Na+, K+, Ca2+, Mg2+, SO42−, and HCO3−), δ18O and δ2H capture water mixing processes during groundwater recharge in the aquifer influenced by irrigation pumping and recharge by precipitations throughout drying and wetting cycles. Results suggest that the intensive weathering of coral limestone under the impact of agricultural fertilization, as indicated by high EC values and high concentrations of Ca2+ and HCO3−. The high DO and low ORP values imply a lack of organic decomposition activity, which is typical in the karst groundwater condition. More importantly, groundwater recharge in the small aquifer was found to be influenced by following three processes. Spatially, both Cl− and NO3−−N indicate the large heterogeneity of water mixing as rainfall or drainage water percolates differently to various parts of the aquifer. Vertically, there are apparent differences in the concentration of Cl− and NO3−−N between the top and bottom layers in several wells, suggesting that, even within a single well, the vertical profile can be weakly connected between the top and the bottom layers. The bottom layer should have been recharged by lateral water inflow from fractures or conduits. And thirdly, the even distribution of Cl− over the long precipitation season from October 2016 and March 2021 and another nearly even distribution of NO3−−N in the low fertilization season in March 2021, demonstrate that the aquifer system is basically connected and porous. In summary, the insights gained from our study reveal that: 1) coral limestone can be treated as a well-developed porous aquifer and 2) the hydrological connectivity is highly heterogeneous and influenced by rock properties, including formation structure, degree of weathering, and clay content. Such findings are essential for understanding the hydrogeological properties of coral limestone aquifers globally.

Vertical heterogeneity of hydrocarbon-degrading bacteria in a core sediment sample from the Central Indian Ridge

Hydrocarbons are ubiquitous in marine environments and might fuel hydrocarbon-metabolizing microbes in the ocean. Numerous studies have documented microbial hydrocarbon degradation in water columns and deep-sea surface sediment. However, the degradation potential and biogeochemical cycling of hydrocarbons in subsurface sediments remain largely unknown. In this study, we used two different hydrocarbons, n-hexadecane (HEX) and methylcyclohexane (MCH), to investigate the distribution and diversity of hydrocarbon-consuming bacteria in a core sediment sample from the Central Indian Ridge (CIR), which is adjacent to mid-ridge hydrothermal vents in the Indian Ocean. We observed different vertical profiles of HEX- and MCH-degrading bacteria in the core sediments. Specifically, HEX-degrading bacteria were universally distributed, while MCH-degrading bacteria were found only in the intermediate layers of the core sediments. Changing factors including dissolved oxygen might affect the natural distribution of different hydrocarbon consumers. We found that a novel species of the genus C1-B045 might play a pivotal role in metabolizing MCH in the CIR deep biosphere. Through amino acid identity comparison with published sequences, we determined that C1-B045 harbors two novel classes of cyclohexanone monooxygenases involved in MCH metabolism. This study sheds light on the structure and function of hydrocarbon-consuming microbes in deep biospheres.

Observing the Near‐Surface Properties of the Greenland Ice Sheet

AbstractSpaceborne radar altimetry over ice sheets has exclusively focused on assessing volume changes through changes in surface elevation. For use in mass balance calculations, these measurements are supplemented with surface property information derived from regional climate models with limited large‐scale observational validation. Simultaneously, the strength at which a radar signal is reflected from the surface contains information on these same near‐surface properties. Here we show that a quantitative interpretation of European (ESA) CryoSat‐2 and French/Indian (CNES/ISRO) SARAL surface echo powers yields timeseries of pan‐Greenland wavelength‐scale roughness and surface density. Individually, the CryoSat‐2 and SARAL results strive toward providing an observational lens to enhance our current understanding of the surface processes affecting ice sheet mass balance and validate their representation by computational models. Taken together, they highlight how diversity in currently operational and future radar satellite altimetry missions can shed light on near‐surface vertical heterogeneity.

Effects of Astronomical Cycles on Laminated Shales of the Paleogene Shahejie Formation in the Dongying Sag, Bohai Bay Basin, China

Laminated shales are widely developed in the Dongying Sag and have attracted much attention as an oil reservoir. Macroscopically, these shales generally have multi-scale cyclicity, which is closely related to the development of laminae. Therefore, analyzing the origin of their cyclicity is helpful to understanding the formation mechanism of laminated shales and the vertical heterogeneity of shale reservoirs, which are of great significance for continental shale oil exploration and development. In this study, a gamma ray (GR) logging series, high-resolution elemental geochemical data, high-resolution core scanning photos and grayscale data, and mineralogical data were used to characterize the cyclicity of shale at different scales, and their relationship with different astronomical cycles was discussed. The results show that the Es3L and Es4U shale in the Dongying Sag has cyclicity from the meter-scale to the ten-meter scale and then to the hundred-meter scale, which is mainly manifested by periodic changes in organic matter abundance, mineral composition, element abundance, and grayscale. These cycles of different scales coincide with different astronomical periods. Specifically, the hundred-meter scale cyclicity is mainly controlled by the very long orbital period; the ten-meter scale cyclicity is mainly related to the eccentricity cycle; while the precession period is the main driver of the meter-scale cyclicity. Finally, we propose a simplified model for illustrating the formation of rhythmic organic-rich shale. This study is helpful to understanding the origin of continental organic-rich shale and predicting shale reservoir properties.

High-resolution geologic modeling and CO2 flow simulation of a realistic clastic deltaic 3D model derived from a laboratory flume tank experiment

Carbon sequestration in geologic reservoirs is a proven method for reducing carbon dioxide (CO2) emissions from industrial point sources. Subsurface deltaic deposits are attractive candidates for CO2 storage projects due to their common occurrence as host to significant hydrocarbon and hydrologic resources, and are highly prospective for subsurface storage of large volumes of CO2. This research informs subsurface passive margin deltaic reservoir characterization and fluid flow performance for geologic carbon storage. The objective of this research is to create a high-resolution geocellular digital model that honors both the large-scale stratigraphic architecture and local vertical and lateral facies heterogeneity represented in deltaic environments of deposition. We present a novel approach to reservoir model generation by utilizing a previously generated meter-scale physical flume sand tank analog experiment to represent realistic deltaic stratigraphic architectures and facies distributions. The result is an extremely high-resolution geologic model that can be used for a variety of applications, including single and multiphase fluid flow simulation. An application is presented using the example model to simulate buoyant fluid (CO2) migration and resulting saturation distribution. Results emphasize the importance of recognizing sequence stratigraphic architectures in deltaic depositional settings for predicting CO2 migration and retention (storage capacity).

Effects of vegetation structure and environmental characteristics on pollinator diversity in urban green spaces

With accelerating urbanization, insect pollinators in urban ecosystems face challenges such as reduced pollen sources, habitat fragmentation, and damage to the nesting environment. Urban green spaces (UGS) are essential for the stability of pollinator communities. However, little is known about the relationship between vertical layer heterogeneity and horizontal layer complexity of vegetation structure in UGS and pollinator communities. The present study aimed to assess how vegetation structure and environmental characteristics shape the insect pollinator community in UGS. To this end, this study was conducted with seven typical vegetation types which were selected according to the biotope mapping classification system (BMCS) in the ring parks around Hefei City, in Anhui province, China. A total of 11,401 pollinators belonging to 6 orders and 34 families were identified during the eight-month survey. Among the seven habitats under the BMCS, mainly successional short-cut shrub and partly open green space, trees two- or multi-layered broad-leaved mixed forest and partly closed green space, and mainly successional tall grass and partly open green space were identified as high-quality insect pollinator habitats. According to the results of the generalized linear regression, the explanatory power of the four best-fitting generalised linear models is relatively high (over 77%). In four optimal models, the effect of vegetation structure on pollinator community was greater than that of environmental characteristics. The redundancy discriminant analysis showed that the flowering abundance of nectar plants, herb richness, and shrub coverage rate were the three most important factors influencing insect pollinator communities, with a cumulative explanatory power of up to 78.8%. Pollinator abundance was positively influenced by spontaneous herbs and low-intensity management. However, high-intensity management, low diversity of plants, low nectar plant richness, ignoring seasonal nectar plant configuration, and dense tree distribution could limit pollinator reproduction and population growth. These results reflect the status of insect pollinator community in UGS in Hefei city and present a possible direction for improving urban green habitats and plant configurations.

A new method to evaluate the effects of mechanical heterogeneity on fault architecture in sedimentary sequences

Geologists have long realized the important control of mechanical stratigraphy on fault nucleation, fault linkage and fault compartmentalization. However, uncertainties still exist concerning the mechanical heterogeneity of the stratigraphy because of difficulties in quantifying the mechanical stratigraphy. In this study, the average strength Save of a section is calculated considering the strength and thickness of each individual layer. The strength difference SΔ and normalized mechanical contrast SΔ/Save are defined for each individual layer in a section. A diagram of the normalized mechanical contrast SΔ/Save versus the normalized height H/∑i=1nhi can then be employed to effectively evaluate the vertical mechanical heterogeneity of the stratigraphy in a section. High-frequency fluctuations between negative and positive SΔ/Save values represent high mechanical heterogeneity, and vice versa. The scaled mechanical heterogeneity MHscaled can be used as an index to quantify the mechanical heterogeneity of the stratigraphy in sections with different scales. However, it is also important to realize that fault throw may promote or restrict the control of mechanical heterogeneity on faulting deformation. The mechanical heterogeneity may present limited control on the fault architecture, particularly when the fault throw is several times larger than the scale of the section. Increasing deformation rate may also minimise the control of the mechanical heterogeneity on the fault architecture.

The Benefits of Using State‐Of‐The‐Art Digital Soil Properties Maps to Improve the Modeling of Soil Moisture in Land Surface Models

AbstractThis study assesses the added value of using emerging maps of soil properties to improve surface soil moisture simulations using the HydroBlocks land surface model with different soil hydraulic parameterization schemes. Simulations were run at an hourly 30‐m resolution between 2012 and 2019 and evaluated against U.S. Climate Reference Network measurements. The results show that state‐of‐the‐art soil properties maps (POLARIS and SoilGrids250m V2.0) improve the accuracy of simulated surface soil moisture when compared to the STATSGO‐derived CONUS‐SOIL map. Contemporary pedotransfer functions (multi‐linear regression and Artificial Neural Networks‐based) also improve model performance in comparison to the lookup table‐derived soil parameterization schemes. The addition of vertical heterogeneity to the soil properties further improves the mean Kling‐Gupta efficiency by 0.04 and lowers the mean Root mean square error by 0.003 over the CONUS. This study demonstrates that land surface modeling can be improved by using state‐of‐the‐art maps of soil properties, accounting for the vertical heterogeneity of soils, and advancing the use of contemporary pedotransfer functions.

Automated Reservoir Model Calibration Applied to a Complex Multizone Heavy Oil Field

_ This article, written by JPT Technology Editor Chris Carpenter, contains highlights of paper SPE 210185, “Automated Reservoir Model Calibration for Field-Development-Plan Evaluation Under Subsurface Uncertainty Applied to a Complex Multizone Heavy Oil Field,” by Luis D. Mendoza, Atahualpa J. Villarroel, SPE, and Maria F. Hurtado, Schlumberger, et al. The paper has not been peer reviewed. _ Fully integrated reservoir modeling for field-development optimization under subsurface uncertainty has been a major challenge so far for Rubiales, a major heavy oil field in Colombia with original oil in place of greater than 5 billion STB. An automated reservoir characterization work flow was developed to generate multiple history-matched models on the field and well level. The developed methodology and work flow successfully delivered field-development evaluation under subsurface uncertainty. The work-flow design is applicable for other fields with similar characteristics and delivery objectives. Introduction The Rubiales field is in the southeast of Puerto Gaitán, approximately 310 km from Bogotá. It is the most important oil field in Colombia in terms of extension, original volumes, and production but also is one of the country’s most complex fields, with a variety of technical challenges. Lithologically, the field is an unconsolidated sandstone reservoir with stratigraphic complexity that includes a high degree of vertical and lateral heterogeneity. Additionally, it has an infinite aquifer with tilted oil/water contact (OWC) in combination with high permeability. The produced fluid is a heavy oil (15 °API) with high water cut (approximately 95%). The current field-development plan (FDP) is based on horizontal wells as infill drilling to maintain oil production levels. A solution was designed to be tested in a specific area of the field implementing digital integration technologies. This solution included the capture of geological heterogeneities in terms of lateral and vertical continuity of the sand bodies. This digital methodology allows understanding seal and OWC effects on water production using a high-resolution model combined with an automated calibration process based on actual production as evidence. Work-Flow Design History Matching. In this project, a typical assisted history-matching process was run initially that included sensitivity and uncertainty analyses followed by an optimization algorithm. Once the sensitivity analysis was complete, the most-influential parameters were identified. Then, some variables were disabled and the uncertainty process was run with a Monte Carlo sampling considering only those parameters considered significant for the objective response. Through this process, a case ensemble with a wide-spectrum result was obtained. The next step was to focus on the optimization that was designed to improve the objective function value that should be minimized by tuning the set of input parameters remaining as active influencers. An evolutionary algorithm was implemented for the optimization that operated with three phases. In the first phase, simulation cases were distributed randomly in the search region, providing an initial ensemble. During the second phase, the better simulation cases were retained and the average objective function value was reduced. In the final phase, only the best simulation cases survived. Even when a modern algorithm was applied, the matching results were not sufficient at the well level. However, this work flow provided a good base case and insights regarding reservoir response, which are the inputs for the next stage of history matching. The further development of the history-matching methodology is detailed in the complete paper.

Long-term impacts of ecosystem restoration on saturated hydraulic conductivity in the Loess Plateau

Since 1999, large-scale ecological restoration activities have been launched across China’s Loess Plateau, known as the ‘Grain for Green Program’, which has substantially changed the landscape as well as the physicochemical properties of soil, including saturated hydraulic conductivity (Ks), an essential parameter in hydrological cycle. This could also contribute to the reduction of water and sediment in the Yellow River in this century. However, there is little study about the changes in root zone (upper one meter) Ks profiles at large scale. Quantitative description of the revegetation impacts on Ks is needed to assess the impacts of revegetation on local water cycle. In this research, 2665 soil samples at six depths within root zone (0–10, 10–20, 20–30, 30–50, 50–70, 70–100 cm) from four land use types (forestland, shrubland, grassland, and sparse grassland) at 25 different sites across the central Loess Plateau were collected and analyzed. Our results suggest that revegetation has significantly improved Ks values with the largest increment in the near surface layers, enlarging the vertical heterogeneity of Ks. The parameters of vertical profiles of Ks were correlated to vegetation index (NDVI), especially for Ks and its decay rate at surface, and the correlations under well vegetated land use types collapsed on the same line. The change of Ks due to revegetation was also correlated to the change of NDVI. This contrasting Ks amplified by revegetation could lead to a shift of runoff generation mechanism in the Loess Plateau. Our results could provide support for estimation of local water usage and availability changes caused by revegetation, shedding light on the vegetation restoration implementation in the Loess Plateau.

Effect of Partially Melting Droplets on Polarimetric and Bi-Spectral Retrieval of Water Cloud Particle Size

Water droplet size, or effective radius, of water clouds can be retrieved using measurements of the total reflectance in near-infrared channels (bi-spectral method) or the linearly polarized radiance of the cloud rainbow (polarimetric method). The retrieved droplet sizes from the bi-spectral method can be significantly larger than those from the polarimetric method. Droplet size vertical heterogeneity is considered as a likely cause of this difference. In this radiative transfer theoretical study, we find that partially melting droplets in mixed-phase clouds can be another cause of this difference. The theoretical study suggests that, when the clouds are dominated by large particles, the existence of partially melting droplets could cause the polarimetric method to underestimate the particle size; but for moderate-sized particles, the polarimetric method should be able to retrieve accurate particle size no matter if the clouds consist of partially melting droplets or not. However, the droplet effective radius of partially melting droplets retrieved by the bi-spectral method can be overestimated if the infrared channel of the bi-spectral method is at 1.6 μm, but is relatively insensitive to the presence of partially melted spherical droplets when the infrared channel is at 2.1 μm. The retrieved droplet size difference between droplet sizes from 1.6 μm and 2.1 μm channels can be used for detecting partially melting droplets.

Network Analysis Indicates Microbial Assemblage Differences in Life Stages of Cladophora.

Cladophora represents a microscopic forest that provides many ecological niches and fosters a diverse microbiota. However, the microbial community on Cladophora in brackish lakes is still poorly understood. In this study, the epiphytic bacterial communities of Cladophora in Qinghai Lake were investigated at three life stages (attached, floating, and decomposing). We found that in the attached stage, Cladophora was enriched with chemoheterotrophic and aerobic microorganisms, including Yoonia-Loktanella and Granulosicoccus. The proportion of phototrophic bacteria was higher in the floating stage, especially Cyanobacteria. The decomposing stage fostered an abundance of bacteria that showed vertical heterogeneity from the surface to the bottom. The surface layer of Cladophora contained mainly stress-tolerant chemoheterotrophic and photoheterotrophic bacteria, including Porphyrobacter and Nonlabens. The microbial community in the middle layer was similar to that of floating-stage Cladophora. Purple oxidizing bacteria were enriched in the bottom layer, with Candidatus Chloroploca, Allochromatium, and Thiocapsa as the dominant genera. The Shannon and Chao1 indices of epibiotic bacterial communities increased monotonically from the attached stage to the decomposing stage. Microbial community composition and functional predictions indicate that a large number of sulfur cycle-associated bacteria play an important role in the development of Cladophora. These results suggest that the microbial assemblage on Cladophora in a brackish lake is complex and contributes to the cycling of materials. IMPORTANCE Cladophora represents a microscopic forest that provides many ecological niches fostering a diverse microbiota, with a complex and intimate relationship between Cladophora and bacteria. Many studies have focused on the microbiology of freshwater Cladophora, but the composition and succession of microorganisms in different life stages of Cladophora, especially in brackish water, have not been explored. In this study, we investigated the microbial assemblages in the life stages of Cladophora in the brackish Qinghai Lake. We show that heterotrophic and photosynthetic autotrophic bacteria are enriched in attached and floating Cladophora, respectively, whereas the epiphytic bacterial community shows vertical heterogeneity in decomposing mats.

Numerical simulation of failure properties of interbedded hydrate-bearing sediment and their implications on field exploitation

Marine natural gas hydrate is widely accepted as a promising substitute energy resource, whereas geotechnical responses of hydrate-bearing sediments (HBS) are crucial for the safe development of natural gas hydrate. The hydrate spatial distribution in sediments is multitype and present different coalescence degrees, its distribution pattern significantly influences the failure behavior of HBS, but the mechanics behind this influence is unclear. Numerical simulations are of great importance in revealing such mechanics. In this paper, a series of numerical simulations, using a modified Mohr–Coulomb (MC) model, were performed to analyze the deformation mechanisms of HBS when the hydrate was laminarly distributed in its host sediment. Considering the vertical heterogeneity of the hydrate distribution, the vertical-laminar hydrate distribution samples were simplified as interbedded samples. The interbedded HBS specimens consisted of three layers: the middle layer with relatively lower hydrate saturation; and the upper and lower layers with relatively high hydrate saturation. Interbedded samples with different hydrate saturations, depressurization ratios, and thickness ratios were designed to analyze their mechanical behaviors. The results indicated that in both the hydrate saturation and thickness ratio of the interbedded samples, the plastic yield area first occurred in the lower hydrate saturation layer; the failure mechanism was mainly determined by the strength of the lower hydrate saturation layer. In addition, the shear failure behavior of the HBS was controlled by the hydrate saturation difference between the middle and upper layers. With a smaller hydrate saturation difference between the middle and upper layers, the samples were more likely to experience shear failure. In the thickness ratio samples, the HBS changed from brittle to ductile failure with an increase in the thickness of the lower hydrate saturation layer. Under different depressurization amplitudes, the larger compression and shear yield zones were mainly confined in the upper layer. This study provides a theoretical framework for controlling the strength-weakening mechanism during natural gas development, which can be used to guide reservoirs characterization, the selection of depressurization schemes, and decomposition zones monitoring.

Multi-scale characterization of organic matter pore space in deep marine shale combined with mathematical morphology

Organic matter (OM)-hosted pores are the primary pore type in the deep shale reservoirs of the Longmaxi Formation, in the Sichuan Basin, China, as well as the main locations for shale gas enrichment. However, the quantitative characteristics of OM pore space are not clear, and there is a lack of characterization methods. In this study, scanning electron microscopy (SEM), mercury intrusion porosimetry, low temperature gas adsorption, and logging data are all used to carry out multi-scale quantitative characterization of OM-hosted pores. The shape factor (F) model, fractal dimension model, and the OM porosity model are improved by combining them with mathematical morphology. The study primarily examines the effect of heterogeneity and the development of full scale pores in the OM pore space of the Luzhou area. The results show that (1) the ideal shape factor (F*) reflects the compaction degree of OM-hosted pores. From macropores to micropores, F* ranges from 0.61 to 0.93 as the shape gradually changes from flattened to round. (2) Shape factor and fractal dimension analyses revealed that the heterogeneity strength of the nanopore structure is mainly controlled by macropores. At the micro level and macro levels, the OM space shows lateral homogeneity and vertical heterogeneity. (3) Compared to the established method, the relative error of the improved model used in this study for predicting OM porosity was reduced by 19%. (4) Compaction and pore structure heterogeneity influence the formation of nanoscale pores in the Luzhou region and these factors also promote the development of mesopores and micropores, respectively. Furthermore, the OM porosities of the dominant shale reservoirs in Luzhou account for more than 50% of the overall porosity. This understanding will be of considerable benefit for evaluating deep shale reservoirs.

Study of the Influence of Dynamic and Static Capillary Forces on Production in Low-Permeability Reservoirs

Low-permeability reservoirs have strong heterogeneity, and the production prediction based on traditional seepage model is not accurate enough. The dynamic capillary-force seepage model can characterize the dynamic heterogeneity of seepage and more accurately describe the oil–water flow process. In this paper, the calculation formula of the dynamic capillary force is obtained through a real low-permeability core experiment, and the seepage model of dynamic capillary force is established. Based on the model, the authors quantitatively study the effects of formation pressure, heterogeneity and production speed on dynamic capillary force through numerical solutions. It is found that compared with the traditional static capillary-force seepage model, the dynamic capillary-force seepage model makes the predicted water cut increase and the recovery factor decrease. With the increase in development time, formation pressure and production rate will make the effect of dynamic capillary force more obvious. According to the comparison of heterogeneous reservoir models, results show that the horizontal heterogeneity will strengthen the dynamic capillary-force effect, while the vertical heterogeneity will weaken the dynamic capillary-force effect. In the range of research parameters, the recovery ratio predicted by the dynamic capillary-force seepage model can be reduced by 4.7%. A new oil–water seepage model is proposed, which can characterize the spatial difference and dynamic change of low-permeability reservoirs with time. It is of great significance for describing the remaining oil distribution of low-permeability reservoirs in detail and making decisions on efficient EOR measures.

How Accurately Can Warm Rain Realistically Be Retrieved with Satellite Sensors? Part II: Horizontal and Vertical Heterogeneities

Abstract There are many sources of uncertainty in satellite precipitation retrievals of warm rain. In this paper, the second of a two-part study, we focus on uncertainties related to spatial heterogeneity and surface clutter. A cloud-resolving model simulation of warm, shallow clouds is used to simulate satellite observations from three theoretical satellite architectures—one similar to the Global Precipitation Measurement Core Observatory, one similar to CloudSat, and one similar to the planned Atmosphere Observing System (AOS). Rain rates are then retrieved using a common optimal estimation framework. For this case, retrieval biases due to nonuniform beamfilling are very large, with retrieved rain rates negatively (low) biased by as much as 40%–50% (depending on satellite architecture) at 5 km horizontal resolution. Surface clutter also acts to negatively bias retrieved rain rates. Combining all sources of uncertainty, the theoretical AOS satellite is found to outperform CloudSat in terms of retrieved surface rain rate, with a bias of −19% as compared with −28%, a reduced spread of retrieval errors, and an additional 17.5% of cases falling within desired uncertainty limits. The results speak to the need for additional high-resolution modeling simulations of warm rain so as to better characterize the uncertainties in satellite precipitation retrievals.