Numerical study on hydrate bearing sediments during gas production

The fully coupled methane hydrate model developed in Cambridge was adopted in this numerical study on gas production trial at the Eastern Nankai Trough, Japan 2013. Based on the latest experimental data of hydrate soil core samples, the clay parameters at Eastern Nankai site were successfully calibrated. With updated clay parameters and site geometry, a 50 days gas production trail was numerically simulated in FLAC2D. The geomechanical behaviour of hydrate bearing sediments under 3 different depressurization strategies were explored and discussed. The results from both axisymmetrical and plane-strain models suggest, the slope of the seabed only affects mechanical properties while no significant impact on the dissociation, temperature and pore pressure. For mechanical deformation after PT recovery, there are large settlements above the perforation zone and small uplift underneath the production zone. To validate the fully coupled model, numerical simulation with finer mesh in the hydrate production zone was carried out. The simulation results suggest good agreement between our model and JOE's results on history matching of gas and water production during trial. Parameter sensitivity of gas production is also investigated and concluded the sea water salinity is a dominant factor for gas production.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

Coupled deformation-flow analysis for methane hydrate production by depressurized wells

A formulation for coupled flow-deformation analysis of methane-hydrate extraction problems is presented. By assuming that the hydrate does not flow, a two phase flow formulation is considered, based on Darcy's law and capillary pressure relation. The formulation is implemented in the finite difference code FLAC. The code was used to investigate the stability of a methane extraction well by depressurizing the well. © 2005 Taylor & Francis Group, London.

Investigation of damping effects on statistical energy analysis of coupled structures

The effects of damping on energy sharing in coupled systems are investigated. The approach taken is to compute the forced response patterns of various idealised systems, and from these to calculate the parameters of Statistical Energy Analysis model for the systems using the matrix inversion approach [1]. It is shown that when SEA models are fitted by this procedure, the values of the coupling loss factors are significantly dependent on damping except when it is sufficiently high. For very lightly damped coupled systems, varying the damping causes the values of the coupling loss factor to vary in direct proportion to the internal loss factor. In the limit of zero damping, the coupling loss factors tend to zero. This is a view which contrasts strongly with 'classical' SEA, in which coupling loss factors are determined by the nature of the coupling between subsystems, independent of subsystem damping. One implication of the strong damping dependency is that equipartition of modal energy under low damping does not in general occur. This is contrary to the classical SEA prediction that equipartition of modal energy always occurs if the damping can be reduced to a sufficiently small value. It is demonstrated that the use of this classical assumption can lead to gross overestimates of subsystem energy ratios, especially in multi-subsystem structures. © 1996 Academic Press Limited.

Numerical analysis of the coupled circuit and cooling holes for an electromagnetic shaker

This paper presents a time-stepping shaker modeling scheme. The new method improves the accuracy of analysis of armature-position-dependent inductances and force factors, analysis of axial variation of current density in copper plates (short-circuited turns), and analysis of cooling holes in the magnetic circuit. Linear movement modeling allows armature position to be precisely included in the shaker analysis. A more accurate calculation of eddy currents in the coupled circuit is in particular crucial for the shaker analysis in a mid-or high-frequency operation range. Large currents in a shaker, including eddy currents, incur large Joule losses, which in turn require the use of a cooling system to keep temperature at bay. Sizable cooling holes have influence on the saturation state of iron poles, and hence have to be properly taken into account.

Coupled foundation-superstructure analysis and influence of building stiffness on foundation response

This paper proposes a simple method to include superstructure stiffness in foundation analyses. The method involves extracting a small "condensed structural matrix" from finite element models of the superstructure, which can then be incorporated into pile group or piled raft analyses using common approaches such as elastic continuum or load transfer methods. The matrix condensation method directly couples structural and geotechnical analyses, and eliminates the need for iterative analyses between structural and geotechnical engineers. Effectiveness of the approach is illustrated through analyses of several buildings designed with a typical floor plan but with varying heights. The parametric study illustrates that superstructure stiffness can have a significant influence on foundation settlement estimates, and the stiffening effects are dominated by the lower stories of the superstructure. The proposed method aims to bridge the gap between structural and geotechnical analyses. Also, being a computationally simple and accurate approach, it is applicable to parametric or optimization studies that would otherwise involve large amounts of analyses. © 2010 ASCE.

Coupled thermal-magnetic analysis on thermally actuated superconducting flux pump

This paper presents numerical analysis of the thermally actuated superconducting flux pump. Visualization of the behavior of the magnetic flux helps our understanding of flux injection mechanism. In addition, in order to confirm validity of the result, we conducted a preliminary flux pump experiment. This result qualitatively agrees well with the experimental one. The flux pump system utilizes a particular behavior that permeability of some materials such as Gadolinium is sensitive to the temperature. In this paper a simple heater is used to control the flux pump system. © 2010 IEEE.