A semi-analytical model for coupled flow in stress-sensitive multi-scale shale reservoirs with fractal characteristics

A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.

Numerical Investigation on the Flow and Temperature Fields in an Inductively Coupled Plasma Reactor

This paper gives a numerical study on the flow and temperature fields in an induced plasma reactor, which worked in 0.5 ATM with air as a working gas. We employed a two-dimensional mode of an inductively coupled plasma to calculate the temperature and flow field of the reactor as well as the generator. The algorithm is based on the solutions of the two-dimensional continuity, momentum, and energy equations in term of vorticity, stream function and enthalpy. An upwind finite-difference scheme was adopted to solve those equations with appropriate boundary conditions. The computed results show that there is a flat region with little parameter change in the reactor, that the diameter of the region is not much larger than that of the generator and that a deep change of parameter exists in the outer side of the region.

Influence of underground water seepage flow on surrounding rock deformation of multi-arch tunnel

Based on a typical multi-arch tunnel in a freeway, the fast Lagrangian analysis of continua in 3 dimensions(FLAC3D) was used to calculate the surrounding rock deformation of the tunnel under which the effect of underground water seepage flow was taken into account or not. The distribution of displacement field around the multi-arch tunnel, which is influenced by the seepage field, was gained. The result indicates that the settlement values of the vault derived from coupling analysis are bigger when considering the seepage flow effect than that not considering. Through the contrast of arch subsidence quantities calculated by two kinds of computation situations, and the comparison between the calculated and measured value of tunnel vault settlement, it is found that the calculated value(5.7-6.0 mm) derived from considering the seepage effect is more close to the measured value(5.8-6.8 mm). Therefore, it is quite necessary to consider the seepage flow effect of the underground water in aquiferous stratum for multi-arch tunnel design.

Two-way coupling of unsaturated-saturated flow by integrating the SWAT and MODFLOW models with application in an irrigation district in arid region of West China

This paper presents the realization of two-way coupling of the unsaturated-saturated flow interactions of the SWAT2000 and MODFLOW96 models on the basis of the integrated surface/groundwater model SWATMOD99, and its application in Hetao Irrigation District （HID）, Inner Mongolia, China. Major revisions and enhancements were made to the SWAT2000 and MODFLOW models for simulating the detailed hydrologic budget and coupled unsaturated and saturated interactions, and irrigation canal hydrology for the HID. The simulation results of seasonal groundwater recharge to and evaporate from the shallow groundwater, and the annual water budget over the district are presented and discussed. The results implied the necessity of two-way coupling of the unsaturated-saturated interactions when groundwater is shallow, and the feasibility of making comprehensive use of the information coming from both the surface water and groundwater models to make a more physically-based assessment of the coupled interactions.

Numerical modeling of influence of thermal flow coupling on flow characteristics of molten steel

Using the equation of continuity and the double equation of Navier-Stokes and k-ε, numerical modeling on a single outlet continuous casting tundish has been carried out during the process of non-thermal flow coupling. The flow field distribution inside the tundish was calculated and the viscosity response time was calculated with the mass transfer equation based on the flow field distribution. The flow characteristics of the molten steel inside the tundish were analyzed, with the results of the numerical modeling compared to the hydraulic modeling. The results showed that the Resident Time Distribution （RTD） curves in the latter anatomosed comparatively better. This certified the validity established by the mathematical model. Numerical modeling was carried out on both large and small tundishes during the processes of thermal flow coupling and also thermal non-flow coupling. The results showed that in regards to large tundishes with relatively simple flow processes, using numerical modeling for thermal flow coupling is necessary.

A well test analysis model of generalized tube flow and seepage coupling

"Generalized mobility"is used to realize the unification of tube flow and seepage in form and the unification of commonly used linear and nonlinear flow laws in form,which makes it possible to use the same form of motion equations to construct unified governing equations for reservoirs of different scales in different regions.Firstly,by defining the generalized mobility under different flow conditions,the basic equation governing fluid flow in reservoir coupling generalized tube flow and seepage is established.Secondly,two typical well test analysis models for coupling tube flow and seepage flow are given,namely,pipe-shaped composite reservoir model and partially open cylindrical reservoir model.The log-log pressure draw-down type-curve of composite pipe-shaped reservoir model can show characteristics of two sets of linear flow.The log-log pressure drawdown plot of partially opened cylindrical reservoir model can show the characteristics of spherical flow and linear flow,as well as spherical flow and radial flow.The pressure build-up derivative curves of the two models basically coincide with their respective pressure drawdown derivative curves in the early stage,pulling down features in the late stage,and the shorter the production time is,the earlier the pulling down feature appears.Finally,the practicability and reliability of the models presented in this paper are verified by three application examples.

Numerical Simulation of Bidirectional Flow and Solid Coupling of Single Pile and Single Row Pile Column and Wave

The three-dimensional model of the wave and pile structure was established by using Ansysworkbench software. The bidirectional fluid-structure interaction between the pile and wave was simulated by using Stokes wave theory and the Volume of Fluid Method. The pressure distribution diagram of the fluid-structure interface and the pile structure were observed and analyzed. The maximum equivalent stress on the pile was studied and its variation range was obtained. Compared with Von Mises yield criterion, it was concluded that the pile was stable under such wave conditions.

A New Unified Stabilized Mixed Finite Element Method of the Stokes-Darcy Coupled Problem: Isotropic Discretization

In this paper, we develop an a-priori error analysis of a new unified mixed finite element method for the coupling of fluid flow with porous media flow in RN, N ∈ {2,3}, on isotropic meshes. Flows are governed by the Stokes and Darcy equations, respectively, and the corresponding transmission conditions are given by mass conservation, balance of normal forces, and the Beavers-Joseph-Saffman law. The approach utilizes a modification of the Darcy problem which allows us to apply a variant nonconforming Crouzeix-Raviart finite element to the whole coupled Stokes-Darcy problem. The well-posedness of the finite element scheme and its convergence analysis are derived. Finally, the numerical experiments are presented, which confirm the excellent stability and accuracy of our method.

The fast method and convergence analysis of the fractional magnetohydrodynamic coupled flow and heat transfer model for the generalized second-grade fluid

In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.

Integrated wellbore-reservoir-geomechanics modeling for enhanced interpretation of distributed fiber-optic strain sensing data in hydraulicfracture analysis

Fiber-optic distributed strain sensing(FO-DSS)has been successful in monitoring strain changes along horizontal wellbores in hydraulically fractured reservoirs.However,the mechanism driving the various FO-DSS responses associated with near-wellbore hydraulic fracture properties is still unclear.To address this knowledge gap,we use coupled wellbore-reservoir-geomechanics simulations to study measured strain-change behavior and infer hydraulic fracture characteristics.The crossflow among fractures is captured through explicit modeling of the transient wellbore flow.In addition,local grid refinement is applied to accurately capture strain changes along the fiber.A Base Case model was designed with four fractures of varying properties,simulating strain change signals when the production well is shut-in for 10 d after 240 d of production and reopened for 2 d.Strain-pressure plots for different fracture clusters were used to gain insights into inferring fracture properties using DSS data.When comparing the model with and without the wellbore,distinct strain change signals were observed,emphasizing the importance of incorporating the wellbore in FO-DSS modeling.The effects of fracture spacing and matrix permeability on strain change signals were thoroughly investigated.The results of our numerical study can improve the understanding of the relation between DSS signals and fracture hydraulic properties,thus maximizing the value of the dataset for fracture diagnostics and characterization.

Optimization of water-urban-agricultural-ecological land use pattern:A case study of Guanzhong Basin in the southern Loess Plateau of Shaanxi Province,China

Extensive land use will cause many environmental problems.It is an urgent task to improve land use efficiency and optimize land use patterns.In recent years,due to the flow decrease,the Guanzhong Basin in Shaanxi Province is confronted with the problem of insufficient water resources reserve.Based on the Coupled Ground-Water and Surface-Water Flow Model(GSFLOW),this paper evaluates the response of water resources in the basin to changes in land use patterns,optimizes the land use pattern,improves the ecological and economic benefits,and the efficiency of various spatial development,providing a reference for ecological protection and high-quality development of the Yellow River Basin.The research shows that the land use pattern in the Guanzhong Basin should be further optimized.Under the condition of considering ecological and economic development,the percentage change of the optimum area of farmland,forest,grassland,water area,and urban area compared with the current land use area ratio is+2.3,+2.4,-6.1,+0.2,and+1.6,respectively.The economic and ecological value of land increases by14.1%and 3.1%,respectively,and the number of water resources can increase by 2.5%.

A novel multi-fidelity coupled simulation method for flow systems

For the numerical simulation of flow systems with various complex components, the traditional one-dimensional （1D） network method has its comparative advantage in time consuming and the CFD method has its absolute advantage in the detailed flow capturing. The proper coupling of the advantages of different dimensional methods can strike balance well between time cost and accuracy and then significantly decrease the whole design cycle for the flow systems in modern machines. A novel multi-fidelity coupled simulation method with numerical zooming is developed for flow systems. This method focuses on the integration of one-, two-and three-dimensional codes for various components. Coupled iterative process for the different dimensional simulation cycles of sub-systems is performed until the concerned flow variables of the whole system achieve convergence. Numerical zooming is employed to update boundary data of components with different dimen-sionalities. Based on this method, a highly automatic, multi-discipline computing environment with integrated zooming is developed. The numerical results of Y-Junction and the air system of a jet engine are presented to verify the solution method. They indicate that this type of multi-fidelity simulationmethod can greatly improve the prediction capability for the flow systems.

Power flow transmission in a coupled flexible system with active executive elements

Based on its prototype of machine-isolator-foundation systems, a theoretical model for dynamic coupled linear system is established, in which both the passive and active control factors are considered. Power flow is used as the cost function to evaluate the isolation effectiveness. And the transmission of vibratory power flow from a vibrating rigid body into a simply supported thin panel through passive isolators and actuators is investigated numerically. The active control strategy is summarized in the conclusion.

Cast-rolling force model in solid-liquid cast-rolling bonding(SLCRB) process for fabricating bimetal clad strips

Based on twin-roll casting, a cast-rolling force model was proposed to predict the rolling force in the bimetal solid-liquid cast-rolling bonding(SLCRB) process. The solid-liquid bonding zone was assumed to be below the kiss point(KP). The deformation resistance of the liquid zone was ignored. Then, the calculation model was derived. A 2D thermal-flow coupled simulation was established to provide a basis for the parameters in the model, and then the rolling forces of the Cu/Al clad strip at different rolling speeds were calculated. Meanwhile, through measurement experiments, the accuracy of the model was verified. The influence of the rolling speed, the substrate strip thickness, and the material on the rolling force was obtained. The results indicate that the rolling force decreases with the increase of the rolling speed and increases with the increase of the thickness and thermal conductivity of the substrate strip. The rolling force is closely related to the KP height. Therefore, the formulation of reasonable process parameters to control the KP height is of great significance to the stability of cast-rolling forming.

NUMERICAL SIMULATION OF CASTING PROCESS

In this paper, the development status of casting numerical simulation technology is introduced. In additional, mathematical model, solution method, initial condition, boundary condition and defect predicting scheme of foundry process are also analyzed, which include the mold filling process, solidification process and the process coupling fluid flow with heat transfer. Finally, a practical casting is taken out to show how to predict defects and optimize foundry process with numerical simulation technology.

4D-stress evolution of tight sandstone reservoir during horizontal wells injection and production: A case study of Yuan 284 block,Ordos Basin,NW China

To investigate the 4D stress change during injection and production in tight sandstone reservoirs, a multi-physical fields modeling method is proposed considering the reservoir heterogeneity, hydraulic fracture and complex injection-production system. The 4D stress evolution of tight sandstone reservoir in Yuan 284 block of Huaqing oilfield, Ordos Basin,during injection-production in horizontal well network is investigated by modeling coupled flow and geomechanics. Results show:(1) Induced by injection and production, the 3D stress increases near the injectors but decreases near the producers, and the horizontal stresses are distributed in obvious strips along their respective stress directions.(2) The horizontal stress difference is the highest at the horizontal wellbore beside injectors during injection and production, while it is the lowest in undeveloped zone between the injectors, and the orientation of maximum horizontal principal stress changes the most near the injectors, which is distributed radially.(3) The hydraulic fracture in re-fracturing well was observed to be asymmetrical in geometry and deflected as the stress changed. The results provide theoretical guidance for horizantal well network modification and re-fracturing optimization design in tight sandstone reservoir.

Temperature influence on macro-mechanics parameter of intact coal sample containing original gas from Baijiao Coal Mine in China

Investigation of temperature effect on mechanical parameters of coal is very important for understanding the mechanical response of coal bed at high temperature.It is especially benefcial for mitigating the thermal-induced disasters occurred in those coal mines suffering from heat hazard.In this work,coal samples,obtained from the No.2442 working face of Baijiao Coal Mine,were subjected to uniaxial compression ranging from 20 to 40℃ with an interval of 5℃.The apparatus used was designed to obtain deformation of a stressed sample,as well as the emission of gases desorbing from coal matrix.The adsorbed gas desorption caused by heating is measured during the entire testing.It is evident that the concentrations of releasing gas(containing methane,carbon dioxide and ethane)slightly rise with increasing temperature.Gas movement observed is closely related to the deformation of coal sample.Both uniaxial compressive strength and elastic modulus of coal samples tend to reduce with temperature.It reveals that increasing temperature can not only result in thermal expansion of coal,but also lead to desorption of preexisting gas in coal which can in turns harden coal due to shrinks of the coal matrix.Even though desorption of adsorbed gas can contribute to the hardening effect for the heated coal,by comparison to the results,it could be inferred that the softening of coal resulted from thermal expansion still predominates changes in mechanical characters of coal sample with temperature at the range from20 to 40℃.

NUMERICAL SIMULATION OF MOLTEN POOL AND CONTROL STRATEGY OF KISS POINT IN A TWIN-ROLL STRIP CASTING PROCESS

Coupled turbulent flow, temperature fields of the twin-roll casting strip process were simulated by three-dimensional finite element method. Based on the heat balance calculation and using inverse methods between the simulations and real experiments, the relational models among casting speed, location, and coefficient of heat transfer between molten metal and rolls in different regions are given. In the simulation, the calculated surface temperatures are in good agreement with the measured values. An on-line model of kiss point is derived by simulations and the geometry of molten pool, corresponding control strategy is also proposed.

A Feasible Zone Analysis Method with Global Partial Load Scanning for Solving Power Flow Coupling Models of CCHP Systems

Heat exchanger systems(HXSs)or heat recovery steam generators(HRSGs)are commonly used in 100 kW to 50 MW combined cooling,heating,and power(CCHP)systems.Power flow coupling(PFC)is found in HXSs and is complex for researchers to quantify.This could possibly mislead the dispatch schedule and result in the inaccurate dispatch.PFC is caused by the inlet and outlet temperatures of each component,gas flow pressure variation,conductive medium flow rate,and atmosphere condition variation.In this paper,the expression of PFC is built by using quadratic functions to fit the non!inearit>of thermal dynamics.While fitting the model,the environmental condition needs prediction,which is calculated using phase space reconstruction(PSR)Kalman filter.In order to solve the complex quadratic dispatch model,a hybrid following electricity load(FEL)and following thermal load(FTL)mode for reducing the dimension of dispatch model,and a feasible zone analysis(FZA)method are proposed.As a result,the PFC problem of CCHP system is solved,and the dispatch cost,investment cost,and the maximum power requirements are optimized.In this paper,a case in Jinan,China is studied.The PFC model is proven to be more precise and accurate compared with traditional models.

Numerical analysis of thermal impact on hydro-mechanical properties of clay

As is known, high-level radioactive waste （HLW） is commonly heat-emitting. Heat output from HLWwilldissipate through the surrounding rocks and induce complex thermo-hydro-mechanical-chemical（THMC） processes. In highly consolidated clayey rocks, thermal effects are particularly significantbecause of their very low permeability and water-saturated state. Thermal impact on the integrity of thegeological barriers is of most importance with regard to the long-term safety of repositories. This studyfocuses on numerical analysis of thermal effects on hydro-mechanical properties of clayey rock using acoupled thermo-mechanical multiphase flow （TH2M） model which is implemented in the finite elementprogramme OpenGeoSys （OGS）. The material properties of the numerical model are characterised by atransversal isotropic elastic model based on Hooke＇s law, a non-isothermal multiphase flow model basedon van Genuchten function and Darcy＇s law, and a transversal isotropic heat transport model based onFourier＇s law. In the numerical approaches, special attention has been paid to the thermal expansion ofthree different phases： gas, fluid and solid, which could induce changes in pore pressure and porosity.Furthermore, the strong swelling and shrinkage behaviours of clayey material are also considered in thepresent model. The model has been applied to simulate a laboratory heating experiment on claystone.The numerical model gives a satisfactory representation of the observed material behaviour in thelaboratory experiment. The comparison of the calculated results with the laboratory findings verifies thatthe simulation with the present numerical model could provide a deeper understanding of the observedeffects. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.