Search NASA⌕ Search

SEARCH · Search NASA

Results for “Oil Reservoirs”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

A comparative investigation of the effect of gas type on foam strength and flow behavior in tight carbonates

Foam generation technique has been practically applied to overcome gas mobility issue during enhanced oil recovery (EOR) processes. Incorporated with the utilized gas, the capability of foaming agent to generate persistent foam determines the flow performance under reservoir conditions. Here, in this study, we first identified the surfactants with excellent phase behavior and foaming properties and, second, evaluated their bulk foam performance with different gases, including N 2 , CO 2 , and hydrocarbon gas (CH 4 ) under high-pressure and high-temperature conditions. In bulk, the foam generated with the high-density gas, i.e., CO 2 , exhibited the greatest foamability but the lowest foam stability. With amidoamine oxide-based surfactant, CH 4 -foam was found to be more stable than N 2 -foam, while the opposite was observed for sulfobetaine-based surfactant. The most effective surfactant was later used to perform flow experiments on two tight carbonate rocks, including Minnesota Norther Cream Buff and Edwards limestone, with permeabilities of ~0.8 and 16 mD, respectively. The core-flooding experiments showed that the optimum quality of CO 2 -foam was lower than that of N 2 - and CH 4 -foams. For a given flow rate and fixed foam quality, the steady-state strength followed the descending order of CH 4 -, N 2 -, and CO 2 -foam. However, in general, the difference in foam strength between N 2 - and CH 4 -foams was insignificant when the total flow rate was varied. We also found that the variation in rock permeability had a greater effect on the strength of the CO 2 -foam compared to those generated with N 2 and CH 4 . In the presence of oil, the shear-thinning behavior of all foams was found more prominent. The increase in oil saturation had a more adverse effect on the strength of N 2 - and CH 4 -foams, particularly in the higher permeability rock. The results indicated that the foam strength, which is key in the effectiveness of foam EOR applications, is greatly influenced by several factors such as gas density, gas solubility in water, the presence of oil, and rock permeability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Newly acquired core enhances geologic understanding of the northern Paradox Basin Cane Creek play, southeastern Utah

The Cane Creek interval is an emerging unconventional heterolithic play within the Pennsylvanian Paradox Formation of the northern Paradox Basin, southeastern Utah. Structural controls (basement uplifts, salt movement, faults, fractures) likely control production from the Cane Creek, where total oil production of ~10 MMBO is just a fraction of current estimates of the undiscovered resource (up to 1.2 BBOE). Oil production has been most successful from unstimulated horizontal wells in the Cane Creek Unit in the central play area, whereas areas to the north (Greentown and Gunnison Valley) have to date proven to be largely unproductive from vertical and horizontal wells drilled without the benefit of 3D seismic or the deployment of modern hydraulic fracturing techniques. However, the limited historical drilling has shown significant promise (e.g., significant shows). The lack of core data in the northern part of the play has impeded thorough assessment of reservoir quality, character, and lithological facies heterogeneity. The Cane Creek is informally divided into three distinct zones. The upper A and basal C zones are composed of fabric destructive and wavy bedded anhydrite, dolomitic mudstone, and organic-rich algal laminated mudstones (source rock). The middle B zone is generally wave rippled and burrowed, low-permeability sandstone-siltstone (reservoir) with subordinate displacive anhydrite. Recent potash cores collected from the Salt Valley anticline near Crescent Junction, Utah (but east of the main play area), show thicker reservoir sandstone/siltstone packages (~40 ft) in the B zone, and parts of the C zone, in the northern part of the Paradox Basin in comparison to the central (~35 ft) and southern (~30 ft) areas. Source rock analyses from northern core/cuttings also indicate that it is well positioned within the oil window (VRo ~0.8), with thin organic-rich mudstone beds containing up to 20 wt% TOC. Because the siliciclastic reservoir rocks have low permeabilities (0.009–0.202 mD) and variable porosities (6–17%) due to dolomite-anhydrite cement, naturally occurring and possibly stimulated fractures may be essential for hydrocarbon recovery. To better understand the northern part of the Cane Creek play, a new core was obtained from the research stratigraphic well State 16-2 (January 2021). Details from this new core will help characterize reservoir quality via core analyses, petrography, fracture orientations, geomechanics, and correlation of facies and depositional trends. Additionally, an evaluation of matrix porosity/permeability will be conducted to determine its contribution to productivity.

02 PETROLEUM↗

A spatial-statistical investigation of surface expressions associated with cyclic steaming in the Midway-Sunset Oil Field, California

In the Midway Sunset Oil Field in Central California, operators inject steam into the shallow diatomite formation to enhance heavy oil recovery through imbibition, wettability alteration, and viscosity reduction, among other mechanisms. The injected steam, however, does not always remain in the reservoir or return through the wells. In two zones in the study area, the steam comes out at the surface, creating sinkholes, seeps, and steam outlets. These phenomena, called “surface expressions,” pose safety and environmental hazards. Even though these surface expressions are a widespread problem in Central California, they are not well documented and understood. Possible causes of the surface expressions include: high injection pressure, structurally controlled flow patterns, leakage of steam through old improperly abandoned wells, high injection volumes, or flow along naturally occurring faults, among other possible factors. This work examines attributes of the zones with surface expressions in order to determine factors that may contribute to their occurrence. Spatial statistical analysis using logistic regression, random forests, and classification trees is used to explore the relationship between the surface expressions and geological and production-related attributes. The results point to a significant spatial correlation between the surface expressions and two predictors: concentration of plugged wells and geologic seal thickness. The results guide follow-up studies to further investigate the role of well abandonment and seal thickness in the occurrence of surface expressions.

02 PETROLEUM↗

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗

Efficient baffle prevents oil backstreaming in diffusion pumps

High vacuum diffusion pump baffle, positioned immediately above pump stack, prevents vaporized pumping oil from entering and contaminating test chamber. Pumping vapors are collected and condensed on the baffle and flow back to pump base reservoir.

Meneghetti, J. R.↗

Tuscaloosa Marine Shale: Seal or Source? Petrophysical Comparative Study of Wells in SE Louisiana and SW Mississippi

The Tuscaloosa Marine Shale (TMS) is a versatile Late Cretaceous shale formation present in central and SE Louisiana and SW Mississippi, which drew attention because of the various roles played within the Tuscaloosa Group. In this paper, it is debated whether the Tuscaloosa Marine Shale can act as a source, reservoir, or seal all throughout the shale play or only in certain areas. Well log and core data from Adams County, Mississippi, are compared to data from East Feliciana Parish in Louisiana. Conclusions were drawn based on the results of well log analysis, X-ray Diffraction (XRD), porosity–permeability measurements, programmed pyrolysis, and fracture analysis. It was shown that the Tuscaloosa Marine Shale interval in SE Louisiana consists of important amounts of calcite, exhibits multiple natural fractures, has porosity values as high as 9.3%, and shows a TOC content of up to 2.8 wt%. On the other hand, samples from a well at the Cranfield field, MS, are characterized by considerably lower TOC values of around 0.88 wt%, porosities between 0.33% and 4%, and no serious fracturing. The formation demonstrates better reservoir and source potential in SE Louisiana and reliable CO2 sealing capacity in SW Mississippi. The analysis presented in this paper represents a holistic approach to the characterization of shale formations, is applicable to other plays around the world, and can be used as an integral part of CO2 sequestration or hydraulic fracturing programs.

04 OIL SHALES AND TAR SANDS↗

Paragenetic controls on CO 2 -fluid-rock interaction and weakening in a macroporous-dominated sandstone

The injection and storage of anthropogenic CO 2 in the subsurface is being deployed as a climate change mitigation tool; however, diagenetic-paragenetic heterogeneity in sandstone reservoirs often contributes to interval specific chemomechanical changes that affect injection and can increase leakage risk.In this paper we address reservoir heterogeneities’ impact on chemomechanical changes in a macroporous-dominated lithofacies of Morrow B sandstone, a formation containing several diagenetically-distinct hydraulic facies while undergoing enhanced oil recovery (EOR) and carbon dioxide (CO 2 ) sequestration. We performed three flow-through experiments using a CO 2 -charged or uncharged formation water combined with four indirect tensile strength tests per post-test sample. We then used the microstructure and paragenetic sequence to understand chemomechanical weakening with key observations as follows: dissolution of carbonates and feldspars changed porosity; increased permeability led to reclassifying each sample in a different hydraulic flow unit; decreased ultrasonic velocity; and did not lead to a loss of tensile strength. Tensile strength maintenance occurred due to the low abundance and minor dissolution of siderite, the stability of quartz, and the relative position of diagenetic ankerite within feldspar. This macroporous-dominated lithofacies is the primary reservoir for the Morrow B Sandstone, and is analogous to other porous sandstone reservoirs. It represents an end-member of a chemomechanically low-risk siliceous CO 2 sequestration and CO 2 -EOR reservoir.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrocoagulation of high-salinity produced water: lessons learned from its early applications in unconventional reservoir plays

Ample research has demonstrated that electrocoagulation is a versatile technology capable of facilitating the removal of a wide range of physical, biological, organic, and inorganic constituents such as suspended solids, turbidity, bacteria, viruses, oil and grease, organic carbon (including chemical and biological oxygen demand), silica, fluoride, and so on. Despite its purported advantages and extensive lab-scale evaluations, it has not yet been widely adopted for water and wastewater treatment and reuse at commercial scale. This is symptomatic of a chasm between its technological capabilities on one hand and field performance and reputation on the other. Herein, we opine as to why this might be the case using our collective experiences with the oil and gas industry as an exemplar. In this study, we highlight scientific, technological, engineering, and business issues that need to be addressed to realize the full potential of this promising technology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prediction and uncertainty quantification of shale well performance using multifidelity Monte Carlo

Uncertainty quantification is an integral component of reservoir management, especially considering the inherent uncertainty in subsurface systems. While a standard practice to estimate the uncertainty, Monte Carlo (MC) simulation is computationally intense when the sampling population comprises high-fidelity simulations. Alternatively, the Multi-fidelity Monte Carlo (MFMC) simulation overcomes this computational intensity by integrating low- and high-fidelity simulations. Our goal is to minimize the number of expensive high-fidelity simulations while maintaining accuracy and using numerous fast and cheap low-fidelity simulations to efficiently sample to input parameter space of interest. We selected gas production from unconventional wells to demonstrate the potential speedups and accuracy of the MFMC approach. The model fidelity usually determines the trade-off between accuracy and efficiency. While the high-fidelity model is more accurate, the low-fidelity model is more efficient. Our high-fidelity simulation includes reservoir simulations of a hydraulically fractured well. On the other hand, our low-fidelity model comprises the parallel-plate flow model. We used differential programming to efficiently solve the 1D flow model, where automatic differentiation is used to efficiently compute the gradients. We matched the production profile of high-fidelity simulations with our low-fidelity simulations. Then, we used a support vector regression to map the high- and low-fidelity input parameters. The mapping function is essential to tune the low-dimensional parameter space of the low-fidelity model to the high-dimensional parameter space of the high-fidelity model. We found that we can use a combination of 9 high fidelity and 10,000 low fidelity simulations to efficiently and accurately simulate pressure management. This method is at least two orders of magnitude faster than only using high-fidelity simulations. Finally, from a broader perspective, MFMC could efficiently estimate the uncertainty of various systems and models, integrating low- and high-fidelity models.

04 OIL SHALES AND TAR SANDS↗

Identifying Abandoned Well Sites Using Database Records and Aeromagnetic Surveys

Oil and natural gas are primary sources of energy in the United States. Improved drilling and extracting techniques have led to a renewed interest in historic oil and gas fields, but limited records of legacy wells make new drilling efforts more difficult, as abandoned wells may provide conduits for liquids and gases to migrate into groundwater reservoirs or the atmosphere. Well finding using aeromagnetic surveys pinpoints the location of steel-cased wells, detecting both active and abandoned wells, including buried casings lacking aboveground markers. Here, we present six aeromagnetic surveys conducted in Pennsylvania and Wyoming as case studies, comparing the magnetic points to locations known in databases. In all study sites, more magnetic points were detected than recorded in databases. Based on differences between theoretical database well counts and the actual number of wells detected in surveys, we estimated the total number of wells in Pennsylvania to be 395 000–466 000 and 181 000–182 000 in Wyoming. Extrapolating to the national level, in this work, we estimate the average number of wells in the continental United States is 6.04 ± 19.97 million wells with 1.16 ± 3.84 million of those designated as abandoned wells, within the range of previous abandoned well count estimations. Although aeromagnetic surveys are limited to detecting steel-cased wells and do not differentiate sites based on well status, this study nevertheless demonstrates the utility of aeromagnetic surveys in well finding efforts across the country and shows limitations in database records of oil and natural gas wells.

54 ENVIRONMENTAL SCIENCES↗

Self-Sealing Mafic Sills for Carbon and Hydrogen Storage

Tabular igneous intrusions (sills) are common features in sedimentary basins and have the potential to be useful seals for fluids (e.g., CO 2 or H 2 ) in geological storage scenarios, and may be important as the need for geologic carbon sequestration and alternative fuel storage increases. This is advantageous in regions without ready access to large-volume reservoirs in depleted hydrocarbon plays and saline aquifers, like the northeastern United States. Moreover, geological H 2 storage requires more demanding conditions than typically associated with oil and gas. An enhanced seal will have the properties of a traditional seal - competent, low permeability, and laterally extensive - with the addition of being able to self-seal pre-existing and induced fractures. Self-sealing will occur along fluid pathways like fractures where CO 2 , water, and minerals like plagioclase, olivine, and pyroxene react together. Dolerite sills from the Gettysburg Basin, Pennsylvania, have remarkably low permeability and homogeneous compositions that include minerals that will readily react with CO 2 dissolved in water. Here, we characterize the physical properties and chemical gradients within several mafic sills cored in five boreholes. In addition, preliminary CO 2 -reaction experiments on dolerite samples demonstrated rapid carbonate mineralization.

58 GEOSCIENCES↗

Novel Geochemistry Determined from High P-T Simulation Experiments of HFTS 2

A standard method in unconventional oil and gas production is the process of hydraulic fracturing followed by a shut-in period, during which the fracture fluid remains pressurized in the reservoir for up to three weeks before production begins. Despite this widely used process, very little is known about what occurs in the reservoir during this shut-in process. In order to properly delineate potential reservoir reactions that may occur during shut-in that would lead to corrosion and scaling events, experiments were conducted in high pressure, high temperature reactors to simulate conditions of the Wolfcamp Formation in the Delaware Basin. Experimental design allowed the assessment of the effect of proppant, microbiology, and time on the mineralogy and fluid chemistry in the reservoir. Results suggest the biggest impact on fluid chemistry and shale mineralogy during shut-in is time. Analyses demonstrate dissolution of the shale material, with maximum dissolved ions occurring after 7 days of the shut-in period. After 21 days, results suggest precipitation occurs. The Delaware Basin is demonstrated to be high in sulfate content, which further increases in the fluid due to dissolution reactions during shut-in. The pressure vessel experiments suggest there was no significant contribution to reactions from microbiology during the shut-in period. However, early production samples demonstrate a significant selection of the microorganism Caminicella, a genus of which has previously been correlated to corrosion and sulfide production. Results suggest shut-in conditions may provide high sulfate concentrations that could later be utilized by a shifted microbial community to drive potential well infrastructure failure. This is the first study to incorporate microbiology, mineralogy, and fluid chemistry to investigate the fundamental geochemical reactions that occur during shut-in and early phase production of the Delaware Basin. Results from this study can complement observations from the Hydraulic Fracture Test Site 2 observations.

Gulliver, Djuna↗

Analysis of Transient Pressure and Rate Data in a Complex of Enhanced Oil Recovery Fields in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. A key activity of the MRCSP monitoring program was to continuously monitor the CO 2 injection and bottom-hole reservoir pressure in the injection wells and selected monitoring wells in the study reefs from 2013 through 2019. During the monitoring period, which varied from reef to reef, the reefs were undergoing one or more phases of development, such as CO 2 injection without hydrocarbon production, CO 2 injection with production Enhanced Oil Recovery (EOR), production without CO 2 injection, and idle (no injection or production) in pressure depleted condition. A multi-year record of reservoir pressure and injection rate data are available for one or more wells for each reef for analysis. This report discusses the methodology and results of various types of analyses performed with this data, with the primary objective of characterizing the in-situ permeability of the formation of interest.

01 COAL, LIGNITE, AND PEAT↗

Positive commandable oiler for satellite bearing lubrication

The results of a feasibility study showed that on-orbit commandable lubrication of ball bearings can be accomplished by direct oil application to the moving ball surfaces. Test results for the lubricant applicator portion of the system are presented, in conjunction with a design approach for the reservoir and metering components.

James, G. E.↗

Gas Hydrate Saturation Estimation Form Acoustic Log Data in the 2018 Alaska North Slope Hydrate-01 Stratigraphic Test

Completed in December 2018, the Alaska North Slope Hydrate 01 stratigraphic test well provides a wealth of logging-while-drilling (LWD) data for strata to below the base of gas hydrate stability (BGHS). This well is intended to be the first of three wells drilled for a comprehensive long-term gas hydrate production test conducted by the National Energy Technology Laboratory, the Japan Oil, Gas and Metals National Corporation, and the U.S. Geological Survey (USGS). The Hydrate 01 stratigraphic test well confirmed the presence of gas hydrate in two sand reservoirs within the hydrate stability zone, indicating the suitability of this location for a long-term gas hydrate production test.<p>The USGS, using an effective-medium-theory rock-physics approach, has estimated gas hydrate saturations from compressional (P) and shear (S) wave log data acquired in the Hydrate 01 well. We assume that gas hydrate occurs as pore-filling load-bearing material (i.e., part of the grain matrix). For Unit D, approximately 500 feet above the BGHS, both P-wave and S-wave acoustic logs indicate moderate gas hydrate saturations with S-wave results slightly lower than those for P-waves. For the Unit B, located just above the BGHS, we obtain moderate to high gas hydrate saturation estimates from both sonic logs. Our P-wave saturation estimates agree well with results from electrical-resistivity-based estimates, whereas estimates from nuclear magnetic resonance LWD data generally suggest 5 to 10 percent higher saturations; our S-wave results suggest lower saturations. These differences likely indicate complexities in the form of gas hydrate occurrence within the sediment pore space, potentially including differences between hydrate occurrence in Units D and B.</p>

Haines, Seth↗

Upscaling Methods Applied to a Fine-Scale Reservoir Model

This study was conducted as part of the Southwest Regional Partnership on Carbon Sequestration (SWP) project to evaluate how upscaling fine-scale simulation models to coarse-scale simulation models impacted the results. The focus was on the Farnsworth Unit (FWU) and its Morrow' B' Sandstone reservoir, specifically the west half of the field. Due to data limitations and the geologic characteristics of the surrounding area, the upscaling was limited to the west half of the FWU rather than a broader basinscale model. The primary aim was to explore how upscaling impacts numerical simulation models, particularly regarding CO 2 -enhanced oil recovery (EOR) and storage capacity predictions. Upscaling was necessary to reduce computational demands when transitioning from high-resolution geological models to coarser grids, as large-scale simulations with finer grids can be computationally prohibitive. This study expands on previous work by the SWP to understand how additional upscaling, applied to already fine-scale numerical simulation models, affects reservoir performance simulations (Ampomah, Balch, & Grigg, 2015). This is key to understanding how loss of resolution can affect coarsescale model results that may be used for large sensitivity analyses, uncertainty quantifications, and training data for machine learning applications.

02 PETROLEUM↗

Effects of Interactions Between Produced Formation Fluid and Rock Matrix on Pore Structure of Caney Shale, Southern Oklahoma

ABSTRACT: Rock-fluid interactions change properties of shales during exploitation. To investigate effects of rock-fluid interactions on pore structure of shales matrix after hydraulic fracturing, powder samples from two late Mississippian Caney Shale cores in the Ardmore Basin, southern Oklahoma, were used to react with formation produced fluid from the field in the batch reactor analysis. X-ray diffraction for mineralogy and Low-pressure nitrogen adsorption isotherms for pore structure were measured for original, after-7days, and after-30days samples. Results show that the samples consist mainly of quartz, followed by clay minerals, carbonates, and feldspar. The pore sizes of micropore (<2 nm) and mesopore (2-50 nm) increase 14%-233% due to dissolution of pyrite, feldspar, and carbonates after 7 days. Due to the transformation from smectite to illite and the increase of pore size, the specific surface area (SSA) decreases after 7-days interactions. After 30-days interactions, the micropore volume slightly increases and the mesopore and macropore volume decreases. Due to the decrease of pore size, the SSA of 30-days reacted samples increases correspondingly and is lower (for the clay-rich sample) or higher (for the calcareous sample) than that of the unreacted samples. Findings improve our understanding of dynamic alteration of shale properties during production. 1. INTRODUCTION Energy demand will continuously grow owing to the increasing global population as well as energy consumption (EIA, 2023). On the other hand, shale gas and oil reshaped the energy market in the United States, enabling the United States to become a net-export of natural gas country in 2017 (EIA, 2023). However, shale reservoirs are challenging tight formations that are still poorly understood in the extraction and production of hydrocarbons (Ross and Bustin, 2009; Curtis et al., 2012; Xiong et al., 2015, 2021a; Li Y. et al., 2016; Gong et al., 2019a; Benge et al., 2021; Awejori et al., 2022; Huang et al., 2022). One of the most challenging topics is the rock-fluid interactions post hydraulic fracturing and its subsequent impacts on the pore structures of fractured formation matrix.

Xiong, Fengyang↗