- Home
- Conferences
- Conference Proceedings
- Conferences
The 11th International Symposium on Recent Advances in Exploration Geophysics (RAEG 2007)
- Conference date: 15 Apr 2007 - 16 Apr 2007
- Location: Kyoto, Japan
- Published: 15 April 2007
-
-
2D and 3D Simulation of Elastic Wave Propagation by Using Element Free Galerkin's Method
Authors M. Katou, T. Matsuoka, H. Mikada, Y. Sanada and Y. AshidaThe element free Galerkin’s (EFG) technique is proposed for solving 2D and 3D elastic wave equation. In this paper, we compare the waveforms calculated by three methodologies; Finite Difference Method with 4th order accuracy in space (FDM4), traditional FEM and EFG. Traditional FEM and EFG are tested on two individual interpolating functions (IFUNC); 1st order (plane) IFUNC and 2nd order (curved) IFUNC. From the comparison to the exact waveforms of Lamb’s problem, the proposed algorism gives accurate waveforms. By two-dimensional implementation, FEM4 demonstrates more accurate result than traditional FEM with plane IFUNC and EFG with plane IFUNC demonstrates more accurate one than FEM4. However, by employing curved IFUNC, traditional FEM and EFG demonstrate almost same accuracy. Moreover, we expand 2D program to 3D one. On the case of 1st order (plane) interpolating function, EFG improves numerical accuracy of traditional FEM without the increase of computational memory.
-
-
-
Numerical Simulation of Rock Test for Engineering Properties Using Mesh-Free Method
Authors J. Takekawa, Y. Yamada, H. Mikada and Y. AshidaAn element free Galerkin method which does not require connectivity between elements and nodes is examined its efficiency for rock experiment. In this study, uni-axial tensile and compressive tests and collapse of cavity are examined. In this process, the rock specimen is assumed to heterogeneous assigned by Weibull distribution. Acoustic Emission which is microscopic behavior of the specimen agrees with laboratory tests, and failure mode of the specimen also agrees with general laboratory test results.
-
-
-
Seismic Wave Simulation based on Method of Characteristics
Authors K. Shiraishi and T. MatsuokaIn this paper, we apply the CIP (Cubic Interpolated Profile) method to the simulation of elastic waves as a highly accurate and stable algorithm to solve first-order wave equations. The key idea of the CIP is that not only the physical value itself but also its first spatial derivative obeys the same equations. Using this property the solution is interpolated by cubic polynomials and interpolation coefficients can be evaluated arithmetically. We implemented this idea to the elastic wave simulation by derivation of the first-order wave equations from the basic equations of motion. The derived equations can be interpreted as the combined first-order wave equations for each mode of wave. Then, we define boundary conditions by using the merit of one-side propagation; free surface, solid-fluid boundary, irregular topography. From simulation study and the stability evaluation, we recognize the method of characteristics with the CIP is a very powerful simulation technique for the elastic wave propagation. Numerical dispersion is negligible, requiring about half the number of grid cells per wavelength than other solvers. This allows accurate, high-frequency, full-wavefield simulation in models with highly variable, random elastic contrasts with fluid-solid mixed media and complex topographic media.
-
-
-
Feasibility Study of Marine Controlled-Source Electromagnetic for Gas Hydrate
Authors K. Yamane, T. Saeki and T. InamoriMarine controlled source electromagnetic method may be one of the most significant new technologies for hydrocarbon exploration. The promise for the technology lies in its ability to differentiate resistive, potentially oil-bearing intervals from surrounding, more conductive water-bearing units. The principle is the same as that used in well logging devices to identify hydrocarbon zones in well bores. The technique is not new but the capability to resolve relatively thin resistive intervals in the depth offers new promise to lower the risk through direct hydrocarbon indicators in conjunction with modern seismic methods. In this paper, we discuss the capability of marine CSEM method towards gas hydrate exploration by numerical studies of 1-D modeling and 2.5-D inversion.
-
-
-
Temporal Change in Reflected Waves from Deep Crust Possibly Synchronized with the Activity of DLFTs Detected Using ACROSS
Authors T. Watanabe, T. Soma, R. Ikuta and K. YamaokaACROSS stands for Accurately Controlled, Routinely Operated Signal System. Repeated transmission of a controlled seismic signal and long-term continuous observation of the seismic waves enable precise monitoring of time-evolving physical properties of the Earth’s interior. An observation using the ACROSS source was conducted in the Tokai distinct, Japan to establish a method to monitor the state of coupling of the plate boundary and to understand the time-evolving phenomena such as slow-slip and deep low-frequency tremors (DLFTs). A seismic signal emitted from the ACROSS source was observed by a short-span surface seismometer array and a Hi-net station (N.HOUH) located at the distance of 57 km from the source for 10 months. The transfer function clearly shows the P-wave, S-wave and their later-phases. Coherent later phases were detected using the semblance analysis of the array of transfer functions. They were interpreted as the reflections from the Moho and the plate boundary. The temporal changes of these phases were examined using Hi-net data which were free from environmental effects. The temporal changes of these phases show correlation with the episodic activity of the DLFTs in 2005. This implies the possibility of seismic monitoring of seismogenic plate-subduction zone using ACROSS.
-
-
-
Evaluation of AVO Attributes Derived from Angle Sub Stacks
More LessAmplitude-versus-offset (AVO) attribute analysis becomes a defact-standard technique for the contemporary hydrocarbon exploration. A basic product for this analysis is angle sub stacks, such as so-called near and far, and sometime middle in order to evaluate an amplitude change along offset or incident angle. Detail study of AVO interpretation is facilitated by crossplotting AVO intercept (P0) and gradient (G) under Shuey’s two term approximation. Further elastic attributes are derived basically from those P0 and G. Estimation of the P0 and G is normally derived via least-square fit of a prestack-migrated seismic gather, which is transformed from offset into incident angle. A partial stack technique is commonly introduced to enhance signal-to-noise ratio (S/N) anticipating better estimation. This ultimate case of the partial stack is the angle sub stack utilizing for amplitude study at an interpretation. In this paper, we are using synthetic and real seismic data contaminated by Gaussian noise to demonstrate that P0 and G estimation based on angle sub stacks is as reliable as computed on a full gather. This implies our standard interpretation data in our hands will be enough for in-house calculation of P0 and G, which can be done by an interpretation workstation, and this contributes budget saving for further P0 and G based processing.
-
-
-
Seismic Imaging for 3-D RVSP Data Using Image Point Transform
Authors C. Lee, T. Matsuoka, K. Ishigaki, M. Tanoue, S. Yamaguchi, T. Tsuchiya and T. MatsuokaImage Points are defined through geometry of sources, receivers and reflectors. IP (Image Point) transform is one type of Radon transform which change common shot section to IP domain (r-x). By the IP transform, reflection signals in time domain can be accumulated and enhanced in IP domain. We applied the IP transform to the field data of RVSP that had been carried out at the Mizunami Underground Research Laboratory (MIU) site. We could get the fault image around the MIU site using IP transform.
-
-
-
A New Method for Stress Orientation Measurements by Using Borehole Deformation and its Applications to the Active Faults
Authors T. Kiguchi and Y. KuwaharaA new sensor tool for borehole deformation using a laser displacement sensor has been developed for stress orientation measurement at shallow depths. The new method for estimation of stress field using this tool is based on following principles: 1) The borehole drilled in a stress field is quasistatically deformed due to viscous property of rock around a borehole. 2) The orientation of maximum horizontal compressive stress (SHmax) is determined from a minor axis of an ellipse fitted to a deformed shape of borehole cross section. One of the advantages of this method to estimate stress orientation is a remarkably lower cost than existing methods, such as stress release and hydraulic fracturing methods. Radial displacements of the borehole wall under creep deformation are measured with the laser displacement sensor rotating continuously around a borehole axis. A total resolution of 0.7μm is achieved by stacking of several thousands data. A field performance test confirms that the new tool can detect the creep deformation, about ten μm increasing of the borehole radius in a measurement of half a day just after drilling. After the performance test, we applied this new method to active faults in order to evaluate stress state of the faults. In the measurement for the Kego fault in Fukuoka Prefecture, western Japan, the average orientation of SHmax is about N60°E for the two measurement point around the southern edge of the fault. The estimated orientation is found to be unfavorable for a strike slip on the Kego fault. In the application to the Atotsugawa fault, central Japan, the orientation of SHmax is estimated to be ENE-WSW by the measurements in a mine cavity near the fault. This orientation is consistent with the trend of the directions of P-axes of the focal mechanism solutions of ultra-micro earthquakes along the fault.
-
-
-
Shallow Exploration of Bangkok Clay Using Multiple Square Arrays
More LessThis study presents on implementation of multiple square arrays, including alpha, beta, gamma in investigating microanisotropy and heterogeneity of Bangkok soft clays. Fusion of these multiple arrays was automatically performed in the field by a special set up and could help identify very well the engineering facilities with a preferential direction, e.g., a concrete water pipe buried in Bangkok clay. As the amount of data generated by automatic multiple square arrays is significant an emphasis in this study was to find out a good way to present and interpret the resistivity data measured by different square array sizes at different azimuthal angles.
-