- Home
- Conferences
- Conference Proceedings
- Conferences
The 14th International Symposium on Recent Advances in Exploration Geophysics (RAEG 2010)
- Conference date: 04 Nov 2010 - 05 Nov 2010
- Location: Kyoto, Japan
- Published: 04 November 2010
-
-
Availability of Fresnel Volume Migration to 3C-Component Seismic Reflection Data Using Tau-P Transform
Authors T. Kawabayashi, H. Mikada, T. Goto and J. TakekawaIn seismic reflection surveys, the superposition of P- and S- waves often influences the resolution of obtained images. We, therefore, needs to first decompose multi-component data before imaging procedure even in Fresnel volume migration. For decomposition processing of refracted P- and S-waves, we need incident angles of received elastic waves. In this study, we calculated the incident angle for each receiver using tangential angles, one of which using amplitudes of acquired three component data and the other using tau-P transform. After the comparison, we used the latter considering that both incident angle and medium wave slowness could be estimated in the acquisition. Also, we use incident angles to the Fresnel volume migration. We found that the slant stack method showed better results for the data in which the secondary P-P reflected wave and first P-S reflected waves are overlapped.
-
-
-
Spiking Deconvolution by Means of Adaptive Simulated Annealing
Authors E.J. HondoriConvolution of seismic wavelet with reflectivity of the earth yields seismic trace which also involves additive noise. Sparse spike deconvolution can remove wavelet effect from seismic noisy data and recover impulse response of the earth. The problem of deconvolution can be viewed as an inversion problem which consists of two main steps. The first stage is to find locations of spikes and the second one is to calculate amplitudes of them. Indeed, desired solution is a set of spikes with appropriate amplitudes such that when convolved with known seismic wavelet fit the data within a predetermined misfit error. Detection of spikes time lags is a nonlinear optimization problem that can be solved using Adaptive Simulated Annealing (ASA). ASA is one of the global stochastic optimization methods which can find the global minimum of complex cost functions with good accuracy. Here for a given seismic trace we locate spikes one by one and after achieving the optimized set of spikes, amplitudes are evaluated using linear least squares.
-
-
-
Crs-Stack Method for Improved Seismic Imaging: A Case Study for Median Tectonic Line
Authors S. Minato, T. Tsuji, T. Matsuoka, Y. Ishikawa, N. Nishizaka, Y. Ohno and M. IkedaWe processed a marine multi-channel seismic (MCS) data acquired across the Median Tectonic Line (MTL) in southwest Japan using CRS-stack method. MTL is an arc-parallel strike-slip fault related to oblique subduction of the Philippine Sea plate beneath the Eurasian plate, and it separates the Ryoke belt from the Sambagawa belt throughout Southwest Japan. Because the subsurface fault geometry of the MTL have not been cleared sufficiently in seismic profiles, the high-quality seismic images are required. A common-reflection-surface (CRS) stack method is a multi-parameter stacking technique which has been proposed as an alternative to the conventional CMP-stack method. The CRS-stack method does not require the macro-velocity model and increases the stacking-fold number. Therefore the CRS-stack method has a potential to enhance reflection signals over the conventional CMP-stack method. Here we applied CRS-stack method to these low-fold seismic data acquired across the MTL and compared the results from conventional CMP-stack method. We show that the CRS-stack method enhanced the reflections signals from faults at depth as well as lithological boundaries.
-
-
-
A Trial About the Evaluation of Compound Geophysical Explorations by Self-Organizing Maps
Authors T. Yamamoto, H. Kusumi, M. Nakamura, T. Yamamoto and T. TsujiIn Japan, in the high economic growth period in 1960’s, a great number of slopes were formed to construct many roads. Now, the slopes have been aging, it is important to estimate the health of the aging slope and maintain slopes effectually. So, in situs, we usually carry out seismic wave method, surface wave method, electric method, electromagnetic wave method, frequency domain electromagnetic method and so on. However, there is not the technique to compound and interpret the result of each geophysical exploration in a numerical formula of the engineering now. Therefore, we notice to self-organizing maps (SOM) used widely in a field of the information processing engineering, and tried to interpret multidimensional data by integrating. In this paper, we classified the ground property by self-organizing maps. The classification result is relatively conformal with boring data. Therefore, it is recognized that it can be used to improve the interpretative accuracy of compound geophysical explorations. And, it can be shown that this technique is effective to estimate of the ground property of the aging slope.
-
-
-
Development of Hard- and Software of Microseismic for Geothermal Exploration and Monitoring: First Results
Authors R. Sule, D. P. Sahara, Y. Kusnadi, A. W. Kusumo, A. D. Nugraha, T. Setiawan, S. Widiyantoro and T. MatsuokaThe Ministry of Research and Technology (RISTEK) of Republic Indonesia launched a so called “Riset Insentif RISTEK” since 2005. This allows Indonesian researchers to be active in research activities, which its funding is higher than other research funding available in Indonesia for small group of researchers. We received this kind of funding for the year of 2010 with its theme is development of hard- and software of microseismic for geothermal exploration and monitoring. This paper describes the first achievements after the funding is granted to us 8 (eight) months ago. This research theme is grouped into 3 main activities, namely development of borehole seismometer by using MEMS, development of precise microseismic relocation technique, and development of first arrival travel time tomography software. The overall current achievements are beyond the target of first year, so that the funding will be extended for the 2nd year.
-
-
-
Application of Self-potential Measurements to Reservoir Engineering
Authors T. IshidoA computational tool so called the “EKP-postprocessor” has been developed to calculate space/time distributions of self-potentials caused by electrokinetic coupling resulting from histories of underground conditions (pressure, temperature, vapor saturation, concentration of dissolved species, etc.) computed by unsteady multi-dimensional geothermal reservoir simulations (Ishido and Pritchett, 1999). This forward calculation technique enables us to incorporate self-potential monitoring data into history-matching studies to improve mathematical reservoir models. In addition to history-matching studies, the EKP-postprocessor can be applied to various problems related to subsurface fluid flow in volcanic areas, sedimentary basins, etc. and characterization of the transport properties of reservoir-forming rocks. In this paper, the results of numerical simulation studies using the EKP-postprocessor to appraise the utility of a combination of pressure and self-potential transient data to characterize fractured reservoirs are presented.
-
-
-
Simulation of Tsunami Propagation with Space-varying Seafloor Topography
Authors T. Ohata, H. Mikada, T. Goto and J. TakekawaTsunami is generated by a sudden deformation of the seafloor, such as uplift and subsidence, caused by fault motion of an earthquake below the seafloor. Numerical simulation of tsunami propagation is frequently used to predict the arrival time and the order of magnitude of the inundation for disaster mitigation purposes. In the propagation process, reflected waves are generated by the change in water depths and influence the tsunami height estimation, in particular in the later phases. In this study, we try to simulate tsunami propagation to accommodate the 2-D varying seafloor topography. In our simulation code, we assume water as a non-viscous fluid. A finite difference method (FDM) is employed using three equations; the equations of continuity, motion, and barotropy. In this study, we simulate the tsunami generation by a sudden change in the water depth and the propagation, using the Pearson approximation to accommodate the spatially varying water depth. We impose the seafloor topography on the basis of the 500m-mesh bathymetry data that JODC (Japan Oceanographic Data Center) provides. We assume a domain included in the data region and simulate tsunami. By using this method, we are able to calculate not only the propagation velocity due to the change in the water depth, but also reflected waves at the same time.
-
-
-
Relationship between Coda-Q and Loaded Stress
Authors K. Okamoto, H. Mikada, T. Goto and J. TakekawaA coda-wave is formed by superposed signals caused by scatterers such as cracks and medium inhomogeneities in the underground materials. For an inhomogeneous medium, it is natural to deal with stochastic methodologies to interpret seismic data. In the past, it is revealed that coda-Q, which is derived from the coda-wave has relationship with frequency of earthquakes1). If the frequency of earthquakes is reflected the change of loaded stress, there is the relationship between the loaded stress and coda-Q. Our purpose is to reveal that relationship. We employ 2-D finite difference method to simulate scatter waves which is generated by cracks. When the stress is loaded to a simulation model, anisotropy of velocity is occurred along direction of the stress, and formation of the scatter waves is changed. Then value of coda-Q is changed. In this study, we find that coda-Q increases constantly according to the increasing stress and behavior of coda-Q has a tendency according to direction of the stress respectively (i.e. stress is loaded vertically or horizontally). In conclusion, we are able to estimate relative magnitude and the direction of the stress. It might lead a new method in underground monitoring.
-
-
-
Monitoring the Plate Boundary in Tokai Area Using Across -Role of Toyohashi Across Source-
Authors T. Watanabe, K. Yamaoka, K. Yamazaki and Y. SuzukiWe have been studying long-term continuous monitoring of seismological processes occurring along and around subducting plate boundary using ACROSS (Accurately-Controlled, Routinely-Operated Signal System). In the Tokai region, three ACROSS sources at Toki, Morimachi and Toyohashi are in operation and concentrated efforts are made for the active subsurface monitoring. Using the continuous seismic record observed at Horai seismic array, it is proved that the signal from each ACROSS source is separated without crosstalk, although the sources are operated simultaneously and continuously. The Toyohashi source is located at preferable place to detect and monitor Tokai slow slip event. Using the signals from the Toyohashi source observed by the Hi-net stations, stacking time and S/N ratio of selected stations are evaluated.
-
-
-
Numerical Simulation of Magma Plumbing System Associated with the Eruption at the Showa Crater of Sakurajima Inferred from Ground Deformation
Authors S. Minami, M. Iguchi, H. Mikada, T. Goto and J. Takekawaground deformation associated with the eruption at the Showa crater of Sakurajima, which has been active since 2006. Using the Mogi’s spherical pressure model, a volume change of magma chambers can be estimated from the displacement, tilt, or strain observations near the ground surface. The tilt and strain data observed 36 hours before an eruption in April 9, 2009, are analyzed using the Mogi’s spherical model to reveal the behavior of magma leading to eruption. From these data, there seems to be a time lag in the inflation between the two magma chambers at a depth of 4km and 0.1km, respectively. Moreover, the volume change of the shallow source is about one order less than that of the deep one. A system which consists of shallow and deep magma chambers and several vertical conduits connecting them is numerically modeled to investigate the mechanism of the time lag and why the difference in the magnitude of the volumetric changes in the two chambers appears as described above. The initial values of magma properties in the deep magma chamber are assumed from the volcanic ejecta of Sakurajima volcano. We assumed that magma is supplied with a constant rate to the deep magma chamber. The pressure limit that the magma starts to ascend to the shallow chamber is assigned to the deep one. In a one-dimensional steady flow model of a magma conduit, we consider the vesiculation of volatile-bearing magma, gas escape and overpressure in the bubble due to the viscous resistance, which largely influences the physical properties of magma. We confirmed that our hypothetical model could explain the time lag of the inflation and the difference in the volume change although our simulation results cannot exactly describe the data. We would like to propose that the numerical simulation could be a powerful tool for understanding the behavior of magma before eruption once ground deformations are well observed.
-
-
-
Marine Controlled-source Electromagnetic Sounding on Submarine Massive Sulphides Using 2.5-D Simulation
Authors N. Imamura, T. Goto, J. Takekawa and H. MikadaRecently, controlled-source electromagnetic method is widely used for shallow sub-seafloor explorations. In this paper, we propose a new controlled-source electromagnetic method using two autonomous underwater vehicles (AUVs) for the exploration of submarine massive sulphides. A numerical simulation code for 2.5-D electromagnetic field in the frequency domain is developed in order to estimate electromagnetic responses on possible conductivity structures. As a result, we confirmed that the magnetic field on the seafloor is very sensitive to the horizontal distribution of buried conductive zone. We also confirmed that the intensity of received electric field decreases near the seafloor in an area wider than the horizontal extent of the buried conductive zone, when the current source dipole is located above the conductive zone. This is because electric field is attenuated due to the return current in the conductive zone. Based on our result, we think that the new controlled-source electromagnetic method is realizable.
-
-
-
Effect of Heterogeneity of Hydraulic Conductivity on Streaming Potential
Authors Y. Ozaki, H. Mikada, T. Goto, J. Takekawa, M. Tsujimura and F. HachaniSelf-potential (SP) is electrical potential mainly generated by thermoelectric, chemical and streaming potentials in the subsurface. The flow of groundwater is often recognized as a bigger source of SP. Using this feature, there are many attempts to localize and quantify flows of liquid in the soil, including groundwater. In case that underground structure is homogeneous, electrical current density according to ground-water flow becomes uniform. Therefore SP on the surface increases monotonically from upstream to downstream of groundwater flow. As a basic interpretation of SP, the direction of increase in SP corresponds with the direction of groundwater flow. However, the anomalous fluctuations of SP start to appear for subsurface inhomogeneous groundwater flow due to the non-uniform conveyance of changes. Extra charge occurs on the boundary of these parameters. As a result, local minimum or maximum in the profile of SP generate just above the boundaries. Our simulation shows that the anomalies of permeability and coupling-coefficient in the subsurface are predominant parameters to effectively estimate the distribution of surface SP in the existence of inhomogeneous underground flow. The effect of coupling coefficient on the SP is often simulated on the underground flow scale. Some anomalies of SP are explained by the inhomogeneous of coupling coefficient. Simple model simulation, for example the well-pump model shows the effect of hydraulic conductivity on the SP. However, SP anomalies that are generated by inhomogeneity of hydraulic conductivity are not simulated on the underground flow scale. We compare the difference of these anomalies and study the different feature of SP anomaly that is generated by hydraulic conductivity from it by coupling coefficient.
-