An analytical solution for describing the transient temperature distribution in an aquifer thermal energy storage system

A mathematical model is developed for predicting the temperature distribution in an aquifer thermal energy storage (ATES) system, which consists of a confined aquifer bounded from above and below by the rocks of different geological properties. The main transfer processes of heat include the conduction and advection in the aquifer and the conduction in the rocks. The semi‐analytical solution in dimensionless form for the model is developed by Laplace transforms and its corresponding time‐domain solution is evaluated by the modified Crump method. Field geothermal property data are used to simulate the temperature distribution in an ATES system. The results show that the heat transfer in the aquifer is fast and has a vast effect on the vicinity of the wellbore. However, the aquifer temperature decreases with increasing radial and vertical distances. The temperature in the aquifer may be overestimated when ignoring the effect of thermal conductivity. The temperature distribution in an ATES system depends on the vertical thermal conduction in the rocks and the horizontal advection and thermal conduction in the aquifer. The present solution is useful in designing and simulating the heat injection facility in the ATES systems. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical analysis of the effect of natural microcracks on the supercritical CO2 fracturing crack network of shale rock based on bonded particle models

The problem of predicting the geometric structure of induced fractures is highly complex and significant in the fracturing stimulation of rock reservoirs. In the traditional continuous fracturing models, the mechanical properties of reservoir rock are input as macroscopic quantities. These models neglect the microcracks and discontinuous characteristics of rock, which are important factors influencing the geometric structure of the induced fractures. In this paper, we simulate supercritical CO₂ fracturing based on the bonded particle model to investigate the effect of original natural microcracks on the induced‐fracture network distribution. The microcracks are simulated explicitly as broken bonds that form and coalesce into macroscopic fractures in the supercritical CO₂ fracturing process. A calculation method for the distribution uniformity index (DUI) is proposed. The influence of the total number and DUI of initial microcracks on the mechanical properties of the rock sample is studied. The DUI of the induced fractures of supercritical CO₂ fracturing and hydraulic fracturing for different DUIs of initial microcracks are compared, holding other conditions constant. The sensitivity of the DUI of the induced fractures to that of initial natural microcracks under different horizontal stress ratios is also probed. The numerical results indicate that the distribution of induced fractures of supercritical CO₂ fracturing is more uniform than that of common hydraulic fracturing when the horizontal stress ratio is small.     

Effect of land‐use changes on nonpoint source pollution in the Xizhi River watershed,Guangdong, China

The Annualized Agricultural Non‐point Source (AnnAGNPS) pollution model has been widely used to assess and predict runoff, soil erosion, sediment and nutrient loading with a geographic information system. This article presents a case study of the effect of land‐use changes on nonpoint source (NPS) pollution using the AnnAGNPS model in the Xizhi River watershed, eastern Pearl River Delta of Guangdong province, China. The land‐use changes in the Xizhi River watershed between 1998 and 2003 were examined using the multitemporal remote sensing data. The runoff, soil erosion, sediment transport and nutrient loading 1998 and 2003 were assessed using AnnAGNPS. The effects of land‐use changes on NPS were studied by comparing the simulation results of each year. Our results showed that (i) the NPS loadings increased when forest and grass land converted into paddy, orchard and farmland land, and population size and gross domestic product size as well as the usage amounts of fertilizer and pesticide in the entire watershed were firmly correlated with the NPS loadings; (ii) the land‐use change during fast urbanization in particular when other land types were converted into the development land and buildup land led to increasing of NPS pollution; and (iii) urban land expansion showed more important effects on total organic carbon (TOC) loading compared with nitrogen and phosphorus loadings. Copyright © 2012 John Wiley & Sons, Ltd.     

One‐dimensional consolidation with a threshold gradient: a Stefan problem with rate‐dependent latent heat

During the process of one‐dimensional consolidation with a threshold gradient, the seepage front moves downward gradually, and the problem is indicated as a Stefan problem. The novel feature in this Stefan problem is a latent heat that varies inversely with the rate of the moving boundary. An exact solution for the external load that increases in proportion to the square root of time is constructed using the similarity transformation technique. Computational examples concerning the effect of different parameters on the motion of the seepage front are presented. The exact solution provides a worthwhile benchmark for verifying the accuracy of numerical and approximate methods. Copyright © 2013 John Wiley & Sons, Ltd.     

Hydrodynamic characteristics and sediment transport of a tidal river under influence of wading engineering groups

By taking the Yong River for example in this paper, based on the multiple measured data during 1957 to 2009, the change process of runoff, tide feature, tidal wave, tidal influx and sediment transport are analyzed. Then a mathematical model is used to reveal the influence mechanism on hydrodynamic characteristics and sediment transport of the wading engineering groups such as a tide gate, a breakwater, reservoirs, bridges and wharves, which were built in different periods. The results showed the hydrodynamic characteristics and sediment transport of the Yong River changed obviously due to the wading engineering groups. The tide gate induced deformation of the tidal wave, obvious reduction of the tidal influx and weakness of the tidal dynamic, decrease of the sediment yield of flood and ebb tide and channel deposition. The breakwater blocked estuarine entrances, resulting in the change of the tidal current and the reduction of the tidal influx in the estuarine area. The large-scale reservoirs gradually made the decrease of the Yong River runoff. The bridge and wharf groups took up cross-section areas, the cumulative affection of which caused the increase of tidal level in the tidal river.     

Gas generation and expulsion characteristics of Middle–Upper Triassic source rocks,eastern Kuqa Depression,Tarim Basin,China: implications for shale gas resource potential

Frontier exploration in the Kuqa Depression, western China, has identified the continuous tight-sand gas accumulation in the Lower Cretaceous and Lower Jurassic as a major unconventional gas pool. However, assessment of the shale gas resource in the Kuqa Depression is new. The shale succession in the Middle–Upper Triassic comprises the Taliqike Formation (T₃t), the Huangshanjie Formation (T₃h) and the middle–upper Karamay Formation (T_2–3k), with an average accumulated thickness of 260 m. The high-quality shale is dominated by type III kerogen with high maturity and an average original total organic carbon (TOC) of about 2.68 wt%. An improved hydrocarbon generation and expulsion model was applied to this self-contained source–reservoir system to reveal the gas generation and expulsion (intensity, efficiency and volume) characteristics of Middle–Upper Triassic source rocks. The maximum volume of shale gas in the source rocks was obtained by determining the difference between generation and expulsion volumes. The results indicate that source rocks reached the hydrocarbon expulsion threshold of 1.1% VR and the hydrocarbon generation and expulsion reached their peak at 1.0% VR and 1.28% VR, with the maximum rate of 56 mg HC/0.1% TOC and 62.8 mg HC/0.1% TOC, respectively. The volumes of gas generation and expulsion from Middle–Upper Triassic source rocks were 12.02 × 10¹² m³ and 5.98 × 10¹² m³, respectively, with the residual volume of 6.04 × 10¹² m³, giving an average gas expulsion efficiency of 44.38% and retention efficiency of 55.62%. Based on the gas generation and expulsion characteristics, the predicted shale gas potential volume is 6.04 × 10¹² m³, indicating a significant shale gas resource in the Middle–Upper Triassic in the eastern Kuqa Depression.     

Modelling combined open channel flow by artificial neural networks

Combined open channel flow is encountered in many hydraulic engineering structures and processes, such as irrigation ditches and wastewater treatment facilities. Extensive experimental studies have conducted to investigate combined flow characteristics. Nevertheless, there is no simple relationship that can fully describe the velocity profiles in a turbulent flow. The artificial neural network (ANN) has great computational capability for solving various complex problems, such as function approximation. The main objective of this study is to evaluate the applicability of the ANN for simulating velocity profiles, velocity contours and estimating the discharges accordingly. The velocity profiles measured by an acoustic doppler velocimeter in the open channel of the Chihtan purification plant, Taipei, with different discharges at fixed measuring section and different depths are presented. The total number of data sets is 640 and the data sets are split into two subsets, i.e. training and validation sets. The backpropagation algorithm is used to construct the neural network. The results demonstrate that the velocity profiles can be modelled by the ANN, and the ANN constructed can nicely fit the velocity profiles and can precisely predict the discharges for the conditions investigated. Copyright © 2005 John Wiley & Sons, Ltd.     

Dynamic analysis of a cylindrical casing–cement structure in a poroelastic stratum

The transient response of a cylindrical casing–cement structure in a poroelastic stratum under dynamic radial tractions is one of the significant issues during the analysis of downhole operations and the selection of safe material. Based on the Biot theory and general elastic mechanics, this paper gives a set of exact solutions for radial displacement, stresses for the casing–cement system and the pore pressure of the infinite surrounding poroelastic stratum in the Laplace transform space. Solutions are presented for three different types of transient radial loadings acting on the surface of casing, i.e., suddenly applied constant load, gradually applied step load and triangular pulse load. Time domain solutions are obtained using a reliable numerical method of inverse Laplace transforms. A detailed parametric study about the transient response is presented both at the casing–cement interface and the cement–stratum interface, and the distributions of the pore pressure and the effective stresses in the stratum are also examined. Copyright © 2017 John Wiley & Sons, Ltd.