180 likes | 337 Views
A New Approach for Location of Induced Earthquakes in Petroleum Fields. Sudipta Sarkar, M. Nafi Toksöz , William L. Rodi Massachusetts Institute of Technology. Outline. Overall Seismicity in the Region. Induced Seismicity and the Petroleum Field. Induced Earthquake Data.
E N D
A New Approach for Location of Induced Earthquakes in Petroleum Fields Sudipta Sarkar, M. NafiToksöz, William L. Rodi Massachusetts Institute of Technology
Importance of Velocity Model in Induced Earthquake Location • In petroleum fields, subsurface velocity models could change abruptly with depth due to layered sedimentary rocks. • Alternating High/Low velocity layers give rise to complex ray paths, thereby affecting travel time calculation. • An accurate velocity model is important for hypocenter determination. • Focal depth is more strongly affected by velocity model because of depth/origin time trade-off.
Depth Distribution of Events located using a “Smooth” Velocity Model
Depth Distribution of Events located using our Improved Velocity Model
Challenges Faced in Location of Reservoir-Induced Seismicity • Very high accuracy requirement • Large uncertainty in S picks • Poor station coverage • An earth structure with alternating high and low velocities in sedimentary sections • Conventional algorithms based on increasing layer velocities/smooth velocity gradients cannot handle many thin layers with high and low velocities.
Our Approach to Location of Reservoir-Induced Seismicity Use highly accurate, detailed velocity model (e.g. derived from well-logs) Pre-calculate 3D travel times using finite difference [Podvin and Lecomte (1991)] and store on disk Compute travel time differences between two stations for each event Minimize residuals of travel-time differences for all station-pairs Use exhaustive grid search method
Our Approach to Location of Reservoir-Induced Seismicity Arrival Time: Event to Station Arrival Time Difference: Event to Station Pairs Travel Time Difference Grid between Station Pairs Find (X,Y,Z) for which errors are “minimum” Data: Model: Objective Function:
Comparison Epicenter locations from an advanced non-linear, 3D location method
Comparison Epicenter locations from our location method
Comparison Depth locations from an advanced non-linear, 3D location method Depth locations from our location method