Cenozoic basement structure of the South China Sea and adjacent areas by modeling and interpreting gravity data

The South China Sea and adjacent areas is a large region with diversely complicated geological conditions. In spite of investigations carried out over the past many years, the marine geological structure in many places has remained poorly understood because of the deficient data, a thick seawater layer as well as of the sensitive areas among the countries in the region. In this paper, the authors study and apply a model-based methodology of the modeling and interpretation of the newest gravity data and others, which are 3D modeling, frequency filtering, horizontal gravity gradient and maximum horizontal gravity gradient, in order to determine clearly the basement structure. The basement features such as the main faults systems, uplift-depression zones and seafloor spreading axis, oceanic boundary in the region have been defined. The achieved results are checked by the seismic data available for the region. From the studied results, the authors have brought out some initial remarks on the structure and form of the basement in the South China Sea and adjacent areas.     

Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids

We consider the exterior Neumann problem of the Laplacian with boundary condition on spheroids. We propose to use spherical radial basis functions in the solution of the boundary integral equation arising from the Dirichlet-to-Neumann map. Our meshless approach with radial basis functions is particularly suitable for handling scattered satellite data. We also propose a preconditioning technique based on an overlapping domain decomposition method to deal with ill-conditioned matrices arising from the approximation problem.     

Observations of interstellar HOCO+: abundance enhancements toward the galactic center

A survey of the 4(04)-3(03) and 1(01)-0(00) transitions of HOCO+ has been made toward several molecular clouds. The HOCO+ molecule was not observed in any sources except Sgr B2 and Sgr A. The 5(05)-4(04) and 4(14)-3(13) transitions were also detected toward Sgr B2. The results indicate that gas phase CO2 is not a major carbon reservoir in typical molecular clouds. In Sgr B2, the HOCO+ antenna temperature exhibits a peak approximately 2' north of the Sgr B2 central position (Sgr B2[M]) and the 4(04)-3(03) line emission is extended over a approximately 10' x 10' region. The column density of HOCO+ at the northern peak in Sgr B2 is approximately 3 x 10(14) cm-2, and the fractional abundance relative to H2 > or = 3 x 10(-10), which is about 2 orders of magnitude greater than recent predictions of quiescent cloud ion-molecule chemistry.     

The HNCO ring in the Sgr B2 region

The large-scale structure associated with the 2′N HNCO peak in Sgr B2 [Minh, Y.C., Haikala, L., Hjalmarson, Å., Irvine, W.M., 1998. ApJ 498, 261 (Paper I)] has been investigated. A ring-like morphology of the HNCO emission has been found; this structure may be colliding with the Principal Cloud of Sgr B2. This “HNCO Ring” appears to be centered at (l,b) = (0.7°,−0.07°), with a radius of 5 pc and a total mass of 1.0 × 10⁵ to 1.6 × 10⁶ M. The expansion velocity of the Ring is estimated to be 30–40 km s⁻¹, which gives an expansion time scale of 1.5 × 10⁵ year. The morphology suggests that collision between the Ring and the Principal Cloud may be triggering the massive star formation in the Sgr B2 cloud sequentially, with the latest star formation taking place at the 2′N position. The chemistry related to HNCO is not certain yet, but if it forms mainly via reaction with the evaporated OCN⁻ from icy grain mantles, the observed enhancement of the HNCO abundance can be understood as resulting from shocks associated with the collision between the Principal Cloud and the expanding HNCO Ring.     

Enhancement of Seismic Data Processing and Interpretation of Fracture Zones on the Upper Part of Granitic Basement in Cuu Long Basin,Vietnam

The fractured granite basement is the primary oil and gas reservoir in the Cuu Long Basin, Vietnam. Due to the complexity of this nonlayered unconventional target, combined with complicated fault and fracture systems, the seismic data quality near and within the basement section is very low. For this reason, it is important to apply improved seismic data processing workflows, filtering and migration techniques, as wells as attribute processing methods to enhance the imaging quality.Our studies show that applying different types of filters, including the f-k, Radon transform and Tau-P, improves signal to noise ratio, removing multiples, revealing basement’s top and its related fractured and fault zones. In addition, the application of multi-arrival-solution migration algorithms, such as Kirchhoff Migration and Controlled Beam Migration, provides improved imaging for identifying basement top and faults and fractures within the basement. Furthermore, the application of seismic attributes such as curvature, apparent dip, or energy gradient, is important in locating faults and fractures, whereas mapping of intensity and orientation of such structures assists the delineation of “sweet spots” and assists the planning of exploration.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

Construction and analysis of higher order finite difference schemes for the 1D wave equation

In this article, we propose to: 1. Establish most of the properties conjectured in [2] about the higher order finite difference approximation of the 1D Laplace operator. 2. Generalize to any order the fourth-order accurate scheme in space and time of Shubin and Bell [20] and Cohen [6]. For this new family of 2m–2m schemes, we establish, via elementary mathematics, various stability and dispersion results that are helpful to compare these schemes to the 2–2m schemes of Anné et al. [2]. This revised version was published online in July 2006 with corrections to the Cover Date.     

Experimental study on the mechanical behaviour of a heat exchanger pile using physical modelling

This study aims to provide knowledge on the thermo-mechanical behaviour of heat exchanger piles, through a laboratory scale model. The model pile (20 mm in external diameter) was embedded in dry sand. The behaviour of the axially loaded pile under thermal cycles was investigated. After applying the axial load on the pile head, the pile temperature was varied between 5 and 30 °C. Seven tests, corresponding to various axial loads ranging from 0 to 70 % of the pile estimated bearing capacity, were performed. The results on pile head displacement show that heating under low axial load induced heave and cooling induced settlement; the pile temperature-displacement curve was found to be reversible and compatible with the thermal expansion curve of the pile. However, at higher axial loads, irreversible settlement of the pile head was observed after a few thermal cycles. The axial load profile measured by the strain gauges evidenced that the pile head load was mainly transferred to the pile toe. Nevertheless, thermal cycles modified significantly the mobilised skin friction along the pile. The total pressure measured at various locations in the soil mass was also slightly influenced by the thermal cycles.