Initial Result of P Wave Tomography Model in Sunda-Banda Arc Transition using FMTOMO

Our study area is located near island Sumbawa, Sumba, Flores, West Timor, Indonesia and East Timor, popularly known as Sunda-Banda arc transition zone. The tectonic setting is mainly controlled by the movement of the oceanic lithosphere Indo-Australian plate subducting the Eurasian plate and Northward migration of Australian continental lithosphere into western Banda-arc in the region of Flores, Sumba and Timor island. We tried to image velocity structure beneath these regions using regional events and tomography inversion model. We collected 5 years of regional events from the Indonesian Agency of Meteorology, Climatology and Geophysics. In total, we reserved 3186 events recorded on 29 stations. For data processing, we used fast marching method as ray tracing between sources and receiver. We then employed subspace inversion as the tomography procedure to estimate the best velocity model representing the tectonic model in the region. Hypocenter data distribution is concentrated on shallow parts of the region and along the Benioff zone down to a maximum depth of 400 km. One of challenge of this study is that although events are abundance, the stations used are mostly located onshore and does not extend in the south-north direction that leads us to under determined problem in the inversion process. However, checker-board models show most our target area can be retrieved to its initial model with sign of smearing effects shown start from a depth of 50 km. After six iteration and optimized selection of damping and smoothing parameters, we observed low velocity anomaly under Bali, Lombok, Sumba, East Nusa Tenggara at shallow depth that may be related with volcanic activity. Deeper low anomaly can also be seen that may be related with partial melting process. A band of fast velocity is clearly seen that goes deepen to the north depicting subducting slabs own to a depth of 300 km. We also observed a possible of fast velocity in the northern part of our stations at shallow depth that we believe may represent the back arc thrust.

Development and Application of Fast Simulation Based on the PSS Pressure as a Spatial Coordinate

Recent studies have shown the utility of the Fast Marching Method and the Diffusive Time of Flight for the rapid simulation and analysis of Unconventional reservoirs, where the time scale for pressure transients are long and field developments are dominated by single well performance. We show that similar fast simulation and multi-well modeling approaches can be developed utilizing the PSS pressure as a spatial coordinate, providing an extension to both Conventional and Unconventional reservoir analysis. We reformulate the multi-dimensional multi-phase flow equations using the PSS pressure drop as a spatial coordinate. Properties are obtained by coarsening and upscaling a fine scale 3D reservoir model, and are then used to obtain fast single well simulation models. We also develop new 1D solutions to the Eikonal equation that are aligned with the PSS discretization, which better represent superposition and finite sized boundary effects than the original 3D Eikonal equation. These solutions allow the use of superposition to extend the single well results to multiple wells. The new solutions to the Eikonal equation more accurately represent multi-fracture interference for a horizontal MTFW well, the effects of strong heterogeneity, and finite reservoir extent than those obtained by the Fast Marching Method. The new methodologies are validated against a series of increasingly heterogeneous synthetic examples, with vertical and horizontal wells. We find that the results are systematically more accurate than those based upon the Diffusive Time of Flight, especially as the wells are placed closer to the reservoir boundary or as heterogeneity increases. The approach is applied to the Brugge benchmark study. We consider the history matching stage of the study and utilize the multi-well fast modeling approach to determine the rank quality of the 100+ static realizations provided in the benchmark dataset against historical data. The multi-well calculation uses superposition to obtain a direct calculation of the interaction of the rates and pressures of the wells without the need to explicitly solve flow equations within the reservoir model. The ranked realizations are then compared against full field simulation to demonstrate the significant reduction in simulation cost and the corresponding ability to explore the subsurface uncertainty more extensively. We demonstrate two completely new methods for rapid reservoir analysis, based upon the use of the PSS pressure as a spatial coordinate. The first approach demonstrates the utility of rapid single well flow simulation, with improved accuracy compared to the use of the Diffusive Time of Flight. We are also able to reformulate and solve the Eikonal equation in these coordinates, giving a rapid analytic method of transient flow analysis for both single and multi-well modeling.

Microseismic P-Wave Travel Time Computation and 3D Localization Based on a 3D High-Order Fast Marching Method

The travel time computation of microseismic waves in different directions (particularly, the diagonal direction) in three-dimensional space has been found to be inaccurate, which seriously affects the localization accuracy of three-dimensional microseismic sources. In order to solve this problem, this research study developed a method of calculating the P-wave travel time based on a 3D high-order fast marching method (3D_H_FMM). This study focused on designing a high-order finite-difference operator in order to realize the accurate calculation of the P-wave travel time in three-dimensional space. The method was validated using homogeneous velocity models and inhomogeneous layered media velocity models of different scales. The results showed that the overall mean absolute error (MAE) of the two homogenous models using 3D_H_FMM had been reduced by 88.335%, and 90.593% compared with the traditional 3D_FMM. On that basis, the three-dimensional localization of microseismic sources was carried out using a particle swarm optimization algorithm. The developed 3D_H_FMM was used to calculate the travel time, then to conduct the localization of the microseismic source in inhomogeneous models. The mean error of the localization results of the different positions in the three-dimensional space was determined to be 1.901 m, and the localization accuracy was found to be superior to that of the traditional 3D_FMM method (mean absolute localization error: 3.447 m) with the small-scaled inhomogeneous model.