WorldWideScience

Sample records for 3-d seismic survey

  1. Marine biota sightings during 3D marine seismic surveys

    Oliveira, Joao Luiz Martinez de; Uller, George A. [CGG do Brasil, Rio de Janeiro, RJ (Brazil); Derntl, Jose Renato; Ribeiro, Camila Castroviejo da Silva; Pereira, Edisio [GEOCOOP Cooperativa de Trabalho, Rio de Janeiro, RJ (Brazil); Miranda, Cristina Maschio de [Nautilus Cooperativa de Trabalho (Brazil); Ferraz, Alexandre Almeida; Costa, Leandro Damiao Soares da [Okeanos Consultoria e Meio Ambiente Ltda. (Brazil)

    2004-07-01

    This work intends to make a correlation between the presence of the marine biota and the seismic source activity (air guns) during seismic surveys, in Campos (BM-C-25 and BM-C-16) and Santos (BM-S-3) Basin, since July 2003 until March 2004. Environmental data were acquired onboard of the Seismic Vessel CGG Harmattan by a team of four oceanographers (environmental technicians), working on the highest place of the Vessel to record and identify the animals whenever was possible. The data were recorded in forms where fields about the biotic and environmental aspects were filled. In 212 days of observations, 2580,1 hours of sighting's effort were recorded; the air guns worked during 37,6% of the time of the effort. These efforts were made during the daylight reaching an average value of 11,35 hours/day. Sightings were divided into the suborders Odontocetes and Mysticetes, and others (fishes, turtles and non identified mammals). 175 sightings were recorded, being 54% when the air gun was off (24% Mysticetes, 56% Odontocetes, 20% others). Similarly, when the air gun was working, 46% of the records were made (24% Mysticetes, 61% Odontocetes, 6% others); the major concentration (58%) of individuals was inside the 1000 m radius around the ship, followed by 14% of the individuals occurring between 3001-4000 m radius away from the ship. The analysis of the data suggests a non-evasive behavior related to the working of the seismic source, corroborating the results reached by other publications using the data collected onboard CGG Vessels. (author)

  2. Near-surface 3D reflections seismic survey; Sanjigen senso hanshaho jishin tansa

    Nakahigashi, H.; Mitsui, H.; Nakano, O.; Kobayashi, T. [DIA Consultants Co. Ltd., Tokyo (Japan)

    1997-05-27

    Faults are being actively investigated across Japan since the Great Hanshin-Awaji Earthquake. Discussed in this report is the application of the 3D near-surface reflection seismic survey in big cities. Data from trenching and drilling is used for the geological interpretation of the surroundings of a fault, and the reflection seismic survey is used to identify the position, etc., of the fault. In this article, when the results obtained from the experimental field are examined, it is found that the conventional 2D imaging reflection survey betrays the limit of its capability when the geological structure is complicated, that the 3D reflection seismic survey, on the contrary, is capable of high-precision imaging and, when augmented by drilling, etc., becomes capable of a more detailed interpretation, and that it also contributes effectively to the improvement of local disaster prevention in big cities. Using as the model the Tachikawa fault that runs near JR Tachikawa Station, embodiment of the 3D reflection seismic survey is reviewed. For the acquisition of data excellent in quality in a 3D reflection seismic survey conducted utilizing the roads in the sector chosen for experiment in the urban area, the shock generating points and receiving points should be positioned by taking into account the parameters in the bin arranging process so that the mid-points will be regularly distributed on the surface. 3 refs., 11 figs., 1 tab.

  3. Case History of 3D Seismic Survey in Aershan Oil Field

    Lu Yousheng

    1996-01-01

    @@ Aershan Oil Field is located at the eastern end of the Manite depression .in the northeastern part of Erlian basin of North China (Fig. 1). Three oilfields -Anan, Abei and Hanan - were discovered in the Erlian basin in the early 1980s using 2D seismic survey with a very dense 0.5 km ×0.5 km grid. However. difficulties were encountered in structural delineation and lateral reservoir prediction.partly due to low signal to noise ratio and resolution of 2D section, partly due to extensive lateral variations. In order to solve these problems, large- area 3D seismic survey was conducted in the late 1980s. 3D survey area is 80 km2.This survey accurately identified the details of reservoir structures, as well as located low-relief structures. Good results were also obtained in predicting lateral reservoir variations and ascertaining oil/water contacts and oilbearing areas.

  4. Large-scale 3-D seismic surveys in Gulf of Mexico - economic approach to exploration

    Graebner, R.J.

    1987-05-01

    Three-dimensional seismic surveys, designed to meet exploration requirements, were conducted in 465 OCS blocks in the Gulf of Mexico, primarily in the deep-water Flexure trend. The geology is very complex, with both salt-induced structures and structures similar to the Rocky Mountain Overthrust belt. To avoid the prohibitive costs in collecting 3-D seismic data at a spatial sample rate required to image the wide range of dips, interpolation methods, based on pattern recognition techniques, were used with a relatively sparsely sampled data collection grid. Within limits, the method is capable of interpolating both crossing and aliased seismic events. The economic impact of this approach is significant because it essentially reduces the cost of 3-D acquisition to the level normally associated with detailed 2-D surveying. However, the familiar benefits of 3-D migration are retained. Using interactive interpretation systems, the results are presented as a comprehensive series of time and event amplitude maps for structural and preliminary stratigraphic mapping. Regional amplitude maps illustrate a combination of structural, stratigraphic, and other geologic phenomena. In addition to conventional structural mapping, the 3-D data may be used to assist in understanding the salt growth history, especially its effect on sedimentation and hydrocarbon potential. It may provide a significant link between regional geology work and the detailed analysis required to develop drillable prospects.

  5. Faults survey by 3D reflection seismics; Sanjigen hanshaho jishin tansa ni yoru danso chosa

    Tsuchiya, T.; Ejiri, T.; Yamada, N.; Narita, N.; Aso, H.; Takano, H.; Matsumura, M. [Dia Consultants Company, Tokyo (Japan)

    1996-10-01

    This paper describes fault survey by 3D seismic reflection exploration. Survey has been conducted mainly at flat land area without pavement not in urban area in Japan. Subsurface structure is complicated with intersecting multiple faults. In this area, a lot of geological investigations have been done prior to the seismic reflection exploration. Fairly certain images of faults have been obtained. However, there were still unknown structures. Survey was conducted at an area of 170m{times}280m in the CDP range. Measurements were carried out by using 100 g of dynamite per seismic generation point combined with 40 Hz velocity geophones. Fixed distribution consisting of lattice points of 12{times}12 was adopted as an observation method. In and around the lattice, a great number of explosions were carried out. The CDP stacking method and the method of migration after stacking were used for the data processing. The 3D structures of six horizons and five faults could be interpreted. Interpreted horizons were well agreed with the logging results. 3 figs.

  6. Quantitative elastic migration. Applications to 3D borehole seismic surveys; Migration elastique quantitative. Applications a la sismique de puits 3D

    Clochard, V.

    1998-12-02

    3D VSP imaging is nowadays a strategic requirement by petroleum companies. It is used to precise in details the geology close to the well. Because of the lack of redundancy and limited coverage in the data. this kind of technology is more restrictive than surface seismic which allows an investigation at a higher scale. Our contribution was to develop an elastic quantitative imagine (GRT migration) which can be applied to 3 components borehole dataset. The method is similar to the Kirchhoff migration using sophistical weighting of the seismic amplitudes. In reality. GRT migration uses pre-calculated Green functions (travel time. amplitude. polarization). The maps are obtained by 3D ray tracing (wavefront construction) in the velocity model. The migration algorithm works with elementary and independent tasks. which is useful to process different kind of dataset (fixed or moving geophone antenna). The study has been followed with validations using asymptotic analytical solution. The ability of reconstruction in 3D borehole survey has been tested in the Overthrust synthetic model. The application to a real circular 3D VSP shows various problems like velocity model building, anisotropy factor and the preprocessing (deconvolution. wave mode separation) which can destroy seismic amplitudes. An isotropic 3 components preprocessing of the whole dataset allows a better lateral reconstruction. The choice of a big migration aperture can help the reconstruction of strong geological dip in spite of migration smiles. Finally, the methodology can be applied to PS converted waves. (author)

  7. Extending field life in offshore Gulf of Mexico using 3-D seismic survey

    Bulling, T.P.; Olsen, R.S. (ARCO Oil and Gas Co., Houston, TX (USA))

    1990-05-01

    Discovered by ARCO in 1967, the High Island 24L field (lower Miocene) is located in the Texas state waters of the Gulf of Mexico. By 1986, the field had produced 320 billion ft{sup 3} of gas and 3.0 million bbl of oil. An engineering field study completed in 1986 showed the field was declining and would be unprofitable within 3 yr. Study of reservoir maps revealed three basin problems: volumetric reserve calculations were less than reserves produced, hydrocarbon-water contacts were inconsistent between wells thought to be in communication, and ultimate recoveries could not be accurately calculated. Attempts to remap the field with the existing two-dimensional seismic data base and well data proved unsuccessful. In 1986, a three-dimensional seismic survey was acquired in an effort to evaluate the true present worth and potential of the field. Remapping of 30 reservoir horizons began in 1987. The integration of detailed well log correlations tied to the dense grid of quality three dimensional seismic data improved the reservoir maps. These maps helped resolve engineering problems by defining the configuration of the reservoirs more accurately. Reservoir maps now closely match volumetrics, fluid contacts within reservoir units are consistent, and a better definition of extension well opportunities exists. The authors study resulted in six additional wells. These wells along with engineering modifications and operations cost containment resulted in the extension of the economic life of the High Island 24-L field by at least 8 yr.

  8. Estimation of 3D positions and orientations of reflectors identified in the reflection seismic survey at the Forsmark area

    The reinterpretation of the seismic data collected in the spring of 2002 in the Forsmark area confirms that the bedrock has been imaged down to depths of at least 3 km. The data quality in general and the frequency content in particular seem to allow a good resolution in identifying and discriminating reflection events in the time domain. The main difficulty resides nonetheless in determining reliably the actual 3D positions and orientations of the reflector imaged with the relatively sparse investigation layout used. If a dense survey grid on the surface were used, it would be possible to apply 3D migration procedures, which would in principle provide the locations and orientations of outstanding reflectors directly by stacking the data. The work volume to acquire and process a full 3D seismic data set would however be in a completely different class of magnitude, compared with the investigation strategy followed so far, and even so 3D seismics would not solve all problems related to the imaging of crystalline rock, e.g. imaging steeply inclined features. With the sparse array at hand, the summation from the 3D migration approach is replaced by a voting principle. That implies that a given reflector is recognized as one and the same in several profiles and a geometrical match is sought amongst the classes of solutions fitting all events identified as belonging to the same reflector. Difficulties can arise when several combinations of events produce viable solutions. This is where the voting comes in. Let it also be noted that selecting an event to form a geometrical solution makes it ineligible for alternative solutions, which raises the problem to the level of a complex combinatory approach. Our objective with the present exercise has been to investigate viable geometrical solutions for 3D reflector positions in the vicinity of the initial guesses proposed in the SKB Report R-02-43. We ventured further from the original orientation proposed only when no group of

  9. A radial 3-D seismic survey at Exxon's Grand Isle 16 field, Gulf of Mexico

    Beebe, R.R.; Harrison, W.K. (Exxon Co., New Orleans, LA (USA))

    1990-05-01

    Grand Isle 16 field is one of several large salt diapirs that occur along the Bay Marchand salt ridge located along the southern boundary of the Terrebonne trough. Discovered in 1948, Exxon's oldest field in the Gulf of Mexico has produced 268 MB0 and 417 GCFG from upper Miocene sands. Over 400 wells have been drilled for hydrocarbons found in three major plays: (1) updip truncation of steeply dipping reservoirs by salt, (2) subunconformity accumulations, and, (3) off-structure stratigraphic traps. Despite extensive drilling, significant potential remains to be evaluated. The authors decided that a 3-D (three-dimensional) seismic survey is the best way to delineate the field's remaining potential. A rectangular survey with lines shot using streamers was eliminated because of surface obstacles, high cost, and inefficient imaging where lines interesect the salt face and steeply dipping beds at an oblique angle. As an alternative, a radial survey was designed whereby lines radiate from a central point over the dome with infill lines added when predetermined separation between lines is reached. Following extensive 2-D (two-dimensional) testing, shooting with bay cables commenced in February 1988, and was completed in October 1988. The processed data set consists of radial lines, circumferential lines, and time slices migrated using a two-pass technique. A regridded cube including inlines migrated with a one-pass (true 3-D) method is also included. This data became available for interpretation in July 1989, and has accomplished the goals of reducing acquisition cost and defining steeply dipping beds near the salt face.

  10. Land 3D-Seismic Data: Preprocessing Quality Control Utilizing Survey Design Specifications, Noise Properties, Normal Moveout, First Breaks, and Offset

    Abdelmoneam Raef

    2009-01-01

    The recent proliferation of the 3D reflection seismic method into the near-surface area of geophysical applications, especially in response to the emergence of the need to comprehensively characterize and monitor near-surface carbon dioxide sequestration in shallow saline aquifers around the world, Justifies the emphasis on cost-effective and robust quality control and assurance (QC/QA) workflow of 3D seismic data preprocessing that is suitable for near-surface applications. The main purpose of our seismic data preprocessing QC is to enable the use of appropriate header information, data that are free of noise-dominated traces, and/or flawed vertical stacking in subsequent processing steps. In this article, I provide an account of utilizing survey design specifications, noise properties, first breaks, and normal moveout for rapid and thorough graphical QC/QA diagnostics, which are easy to apply and efficient in the diagnosis of inconsistencies. A correlated vibroseis time-lapse 3D-seismic data set from n CO2-flood monitoring survey is used for demonstrating QC dlagnostles. An Important by-product of the QC workflow is establishing the number of layers for n refraction statics model in a data-driven graphical manner that capitalizes on the spatial coverage of the 3D seismic data.

  11. Land 3D-seismic data: Preprocessing quality control utilizing survey design specifications, noise properties, normal moveout, first breaks, and offset

    Raef, A.

    2009-01-01

    The recent proliferation of the 3D reflection seismic method into the near-surface area of geophysical applications, especially in response to the emergence of the need to comprehensively characterize and monitor near-surface carbon dioxide sequestration in shallow saline aquifers around the world, justifies the emphasis on cost-effective and robust quality control and assurance (QC/QA) workflow of 3D seismic data preprocessing that is suitable for near-surface applications. The main purpose of our seismic data preprocessing QC is to enable the use of appropriate header information, data that are free of noise-dominated traces, and/or flawed vertical stacking in subsequent processing steps. In this article, I provide an account of utilizing survey design specifications, noise properties, first breaks, and normal moveout for rapid and thorough graphical QC/QA diagnostics, which are easy to apply and efficient in the diagnosis of inconsistencies. A correlated vibroseis time-lapse 3D-seismic data set from a CO2-flood monitoring survey is used for demonstrating QC diagnostics. An important by-product of the QC workflow is establishing the number of layers for a refraction statics model in a data-driven graphical manner that capitalizes on the spatial coverage of the 3D seismic data. ?? China University of Geosciences (Wuhan) and Springer-Verlag GmbH 2009.

  12. Down-hole seismic survey system with fiber-optic accelerometer sensor array for 3-dimensions vertical seismic profile (3D-VSP)

    Zou, Qilin; Wang, Liwei; Pang, Meng; Tu, Dongsheng; Zhang, Min; Liao, Yanbiao

    2006-08-01

    We demonstrated a down-hole seismic survey system that can be applied in three dimensions vertical seismic profile (VSP) detection in petroleum exploration. The results of experiments show that the system has a dynamic measurement range of 80db (ratio of maximum to minimum value) and the total delay for signal collection, process and communication is less than 200ms @ 2k bit sample rates. An array consisting of six fiber-optic accelerometers (receivers) is applied in this system. Each receiver is comprised of three fiber-optic Michelson interferometers. In order to meet the requirements of high precision and real-time measurement, the high-speed DSP chips are employed to realize the algorithms of signal filters and Phase Generated Carrier (PGC) demodulation to obtain the seismic information. Multi-ARM CPUs are introduced into the system to design the fiber-optic accelerometer array controller and the receiver array local bus that are used for real-time data communication between the multi-level receivers and controller. The system interface for traditional ELIS Down-hole Instrument Bus (EDIB) is designed by the use of FPGA so that our system can attach to EDIB and cooperate with other instruments. The design and experiments of the system are given in this paper in detail.

  13. NORTH HILL CREEK 3-D SEISMIC EXPLORATION PROJECT

    Marc T. Eckels; David H. Suek; Denise H. Harrison; Paul J. Harrison

    2004-05-06

    Wind River Resources Corporation (WRRC) received a DOE grant in support of its proposal to acquire, process and interpret fifteen square miles of high-quality 3-D seismic data on non-allotted trust lands of the Uintah and Ouray (Ute) Indian Reservation, northeastern Utah, in 2000. Subsequent to receiving notice that its proposal would be funded, WRRC was able to add ten square miles of adjacent state and federal mineral acreage underlying tribal surface lands by arrangement with the operator of the Flat Rock Field. The twenty-five square mile 3-D seismic survey was conducted during the fall of 2000. The data were processed through the winter of 2000-2001, and initial interpretation took place during the spring of 2001. The initial interpretation identified multiple attractive drilling prospects, two of which were staked and permitted during the summer of 2001. The two initial wells were drilled in September and October of 2001. A deeper test was drilled in June of 2002. Subsequently a ten-well deep drilling evaluation program was conducted from October of 2002 through March 2004. The present report discusses the background of the project; design and execution of the 3-D seismic survey; processing and interpretation of the data; and drilling, completion and production results of a sample of the wells drilled on the basis of the interpreted survey. Fifteen wells have been drilled to test targets identified on the North Hill Creek 3-D Seismic Survey. None of these wildcat exploratory wells has been a dry hole, and several are among the best gas producers in Utah. The quality of the data produced by this first significant exploratory 3-D survey in the Uinta Basin has encouraged other operators to employ this technology. At least two additional 3-D seismic surveys have been completed in the vicinity of the North Hill Creek Survey, and five additional surveys are being planned for the 2004 field season. This project was successful in finding commercial oil, natural gas

  14. Anisotropic 3D Modeling for Long Offset VSP Survey Design

    Today's seismic techniques allow the geoscientist to do the interpretation more quantitatively. AVO and anisotropy measurements are the examples of DHI (Direct Hydrocarbon Indication). These measurements can be done accurately using long offset borehole seismic survey such as walk away VSP, having the geophones located down hole close to the target formation. This paper will show the importance 3D seismic modeling prior to the survey, by simulating the seismic wave propagation in three-dimensional volume filled with continuous material properties. This pre-survey modeling can help us suppressing the uncertainties and narrowing the error bars on the real survey. Some examples from offshore Nigeria showed dramatic geometrical differences between ordinary 2D compared to 3D observations Assumption that the seismic wave travels in 2D plane is not always acceptable for survey design. The examples also demonstrated the ability to observe some critical information such as the limit of incidence angle, compromise between resolution and image coverage, effects of velocity anomalies, anisotropy and dipping formations on lateral coverage. Fluid effect in 3D modeling will also be discussed here. Amplitude anomalies are predicted by replacing different type of fluids effect in the target reservoirs, as well as various types of AVO classes. A well-prepared long offset VSP survey is very critical to provide us high quality and high accuracy information that can be used to calibrate and optimise the full 3D seismic processing and interpretation in the area. This process is known as Well Driven Seismic (WDS)

  15. Additional geological insight brought by 3-D seismic data

    Liu, Y.; Milton, J.

    2002-12-01

    3-D seismic data allows geo-scientists to study the earth at a level that is much more detailed than ever. It is shown in this case study that 3-D seismic can not only be used to identify different types of oil/gas reservoirs, significantly reduce the risk of misinterpretation, but also help to understand geological history and establish paleogeography at different geological times. In the study of Southeast Maricopa Seismic Survey in southern San Joaquin Valley, two types of potential hydrocarbon traps are interpreted: stratigraphical traps due to turbidite channels, and structural traps due to faulting. The distinctive characteristics of two types of channels indicate different depositional environments. With 3-D visualization tools, it is found that localized faults had been leaking during certain geological times, resulting in structural traps of oil/gas. A geological history of the local area can be estimated by building a series of pseudo-paleogeographic maps using 3-D seismic data, which further reconfirms the existence of different depositional systems indicated by two distinctive types of channels.

  16. The business impact and value of 3-D seismic

    3-D Seismic has had a profound and lasting impact on the hydrocarbon exploration and production industry. It is a technology which is often excellent at assessing the risk associated with trap definition, seal and reservoir distribution, the very parameters which are the most critical to optimizing the economics associated with E and P projects. This paper discusses Amoco Corporation's experience with 3-D Seismic when used for new field rate acceleration, older field extension, and wildcat exploration. Its emphasis is on assessing the value added by 3-D by reviewing recent E and P experiences in a post-appraisal mode and then in applying the lessons learned from these analyses and case histories to potential new projects. This work is significant because it first assesses the impact 3-D has had on a large number of business situations at Amoco; that is, it is based on data collected on159 3-D surveys acquired at Amoco between 1991--1994. Second, it uses the data collected from these surveys and applies the business improvements observed in the data to typical international business opportunities to quantify, in expected value $ terms, the value that the technology brings to an average project. Finally, it looks at project economics not only from an oil company perspective, but from the standpoint of a host government, with a discussion of insights and implications of the data, economics and techniques utilized

  17. Deep reflection structure imaged by the 2008 3D seismic reflection Survey at the RIDGE- 2000 East Pacific Rise Integrated Studies Site

    Nedimović, M. R.; Carbotte, S. M.; Mutter, J. C.; Canales, P. J.; Carton, H.; Aghaei, O.; Marjanović, M.; Newman, K. R.; Hu, M.; Stowe, L.

    2008-12-01

    The first multi-source and multi-streamer 3D seismic reflection experiment carried out using academic resources was done aboard the R/V Marcus G. Langseth in Summer 2008 during cruise MGL0812. The targeted area was the RIDGE-2000 Integrated Studies Site at the East Pacific Rise. Our primary 3D survey grid extends from about 9° 57'N to 9° 42'N, with a smaller grid just to the south covering approximately from 9° 40'N to 9° 37.5'N. Additionally, about 1 and 0.5 km wide across-ridge-axis swaths of data were collected at 9° 36'N and 9°30'N respectively, as well as an along-ridge-axis swath about 1 km wide and extending from 10° 05'N to 9° 40'N. We here focus on a preliminary analysis of the reflection structure imaged within the lower crust and uppermost mantle. Moho reflection arrivals are imaged through much of the investigated area. The character of Moho reflection events varies from simple, single reflection wavelet to more complex arrivals indicating spatial changes in structure within the Moho transition zone. Particularly strong Moho reflections are observed in the southern half of the main 3D grid. In places, Moho reflection event appears to extend across the ridge axis potentially suggesting "zero-age" Moho development. Weak Moho arrivals are found at the north end of the main 3D box and within the smaller box to the south. Most notable place lacking Moho reflections is the Lamont seamount area where Moho is not observed on either side of the ridge axis, although the area lacking Moho reflections is wider on the western ridge flank. Further south, along the across-ridge-axis swaths, Moho reflections again become more pronounced. A suit of what mostly appear to be reflection events is recognized between the AMC and Moho. Many of them do not appear to be multiples, and their origin is not well understood. Possible origins for these events include: lower boundary of the AMC, S-converted waves, and lower crustal melt lenses. Along sections of the two 3D

  18. Surface 3-D reflection seismics - implementation at the Olkiluoto site

    Posiva Oy takes care of the final disposal of spent nuclear fuel in Finland. In year 2001 Olkiluoto was selected for the site of final disposal. Construction of the underground research facility, ONKALO, is going on at the Olkiluoto site. The aim of this work was to study the possibilities for surface 3-D seismics and to review experiences for design before field work. The physical parameters and geometric properties of the site, as well as efficient survey layout and source arrangements, were considered in this work. Reflection seismics is most used geophysical investigation method in oil exploration and earth studies in sedimentary environment. Recently method has also been applied in crystalline bedrock for ore exploration and nuclear waste disposal site investigations. The advantage of the method is high accuracy combined with large depth of investigation. The principles of seismic 2-D and 3-D soundings are well known and advanced. 3-D sounding is a straightforward expansion of 2-D line based surveying. In investigation of crystalline bedrock, the high frequency wave sources and receivers, their right use in measurements and careful processing procedure (refraction static corrections in particular) are important. Using the site parameters in 2-D numerical modeling, two cases of faulted thin layer at depths of 200, 400 and 600 meters were studied. The first case was a layer with vertical dislocation (a ramp) and the other a layer having limited width of dislocated part. Central frequencies were 100, 200, 400 and 700 Hz. Results indicate that 10 - 20 m dislocation is recognizable, but for depths greater than 600 m, over 20 meters is required. Width of the dislocated part will affect the detectability of vertical displacement. At depths of 200 m and 400 m 10 - 50 m wide parts appear as point-like scatterers, wider areas have more continuity. Dislocations larger than 20 m can be seen. From depth of 600 m over 100 m wide parts are discernible, narrower are visible

  19. ONKALO 3D tunnel seismic investigations at Olkiluoto in 2009

    Cosma, C.; Enescu, N.; Balu, L.; Jacome, M. (Vibrometric Oy, Perttuola (Finland))

    2011-02-15

    POSIVA Oy conducts bedrock investigations at the spent nuclear fuel final disposal site at Olkiluoto, in western Finland. The excavation of the access tunnel to the repository hosts the ONKALO underground rock characterization facility. The investigations carried out at ONKALO focus on the bedrock and groundwater conditions prevailing on the final disposal site and how construction work affects them. Tunnel seismic investigations were carried out in July 2009, as an extension of similar work performed in December 2007. The main objective of the tunnel seismic investigations have been to demonstrate the possibility to detect, locate and image cost effectively steeply and gently dipping fractures, at the side and/or below the tunnel and to characterize the volume of rock surrounding a 250 m long segment of the ONKALO tunnel. The survey was conducted at a depth of 350 m, over a 240 m long line of 3-components receivers, spaced at 3m intervals. Seismic signals were produced along two lines, on the tunnel wall and floor, with source points spaced at 1m. A timedistributed swept-impact, the Vibsist-250 hydraulic source, was used. The source was hosted on a mini excavator. Receiver holes approximately 0.4 m deep were drilled prior to the survey, horizontally into the tunnel wall. One of the procedures used for data stacking and migration is based on a proprietary method combining the DMO (Dip Move Out) correction and an expression of the Radon Transform. Horizontal and vertical migrated profiles were computed both for the P wave and S wave reflected wave fields. A true 3D migration technique (Image Point migration) was used to create 3D migrated sections oriented to incremental azimuths around the tunnel, the result being a cylindrical imaging volume. A general conclusion is that seismic surveys along the tunnel can economically be used for rock mass characterization. High quality results can be obtained by operations in tunnel working conditions, provided that due

  20. ONKALO 3D tunnel seismic investigations at Olkiluoto in 2009

    POSIVA Oy conducts bedrock investigations at the spent nuclear fuel final disposal site at Olkiluoto, in western Finland. The excavation of the access tunnel to the repository hosts the ONKALO underground rock characterization facility. The investigations carried out at ONKALO focus on the bedrock and groundwater conditions prevailing on the final disposal site and how construction work affects them. Tunnel seismic investigations were carried out in July 2009, as an extension of similar work performed in December 2007. The main objective of the tunnel seismic investigations have been to demonstrate the possibility to detect, locate and image cost effectively steeply and gently dipping fractures, at the side and/or below the tunnel and to characterize the volume of rock surrounding a 250 m long segment of the ONKALO tunnel. The survey was conducted at a depth of 350 m, over a 240 m long line of 3-components receivers, spaced at 3m intervals. Seismic signals were produced along two lines, on the tunnel wall and floor, with source points spaced at 1m. A timedistributed swept-impact, the Vibsist-250 hydraulic source, was used. The source was hosted on a mini excavator. Receiver holes approximately 0.4 m deep were drilled prior to the survey, horizontally into the tunnel wall. One of the procedures used for data stacking and migration is based on a proprietary method combining the DMO (Dip Move Out) correction and an expression of the Radon Transform. Horizontal and vertical migrated profiles were computed both for the P wave and S wave reflected wave fields. A true 3D migration technique (Image Point migration) was used to create 3D migrated sections oriented to incremental azimuths around the tunnel, the result being a cylindrical imaging volume. A general conclusion is that seismic surveys along the tunnel can economically be used for rock mass characterization. High quality results can be obtained by operations in tunnel working conditions, provided that due

  1. Joint-inversion of Seismic and Gravity Survey Data in Sedimentary Basins for 3-D Underground Structures in Strong Motion Simulation

    Koketsu, K.; Afnimar, A.

    2001-12-01

    Modeling of underground structures is one of important tasks for strong motion predition. Seismic surveys such as refraction and reflection explorations can give precise estimates of seismic velocities, but they are too expensive to cover the whole area of a sedimentary basin with a fine spacing of survey points. On the other hand, gravity surveys can be carried out densely and homogeneously with reasonable cost, but they cannot measure seismic velocities directly. In order to compromise this inconsistency, we combine data from the seismic and gravity surveys, and jointly invert them assuming a relation between densities and seismic velocities. An interface separating sediments and basement rocks is parameterized by the Lagrange interpolation with knot points distributed for representing a steep part related to a fault zone. The linear and absolute constraints are imposed to avoid oscillatory artifacts in the solution and to minimize the non-uniqueness of the solution. The joint inversion also overcome the ray coverage problem of a refraction survey in a complex underground structure. The validity of the formulation is verified using synthetic data for a valley-like basin with pseudo-random noises. The proposed method is successfully applied to the actual data obtained in and around the Osaka basin, Japan. The result shows a good agreement with the geological setting in this region. The instability due to the lack of data beneath the Osaka bay is improved by the introduction of the absolute constraints. We then apply the method to the Kanto basin around Tokyo assuming two layers in the sediments. This layering increases the number of variables introducing some instabilities into the inversion, and so measurements at boreholes and reflection survey lines are used as fixed points in the model.

  2. Pseudo-random data acquisition geometry in 3D seismic survey; Sanjigen jishin tansa ni okeru giji random data shutoku reiauto ni tsuite

    Minegishi, M.; Tsuburaya, Y. [Japan National Oil Corp., Tokyo (Japan). Technology Research Center

    1996-10-01

    Influence of pseudo-random geometry on the imaging for 3D seismic exploration data acquisition has been investigate using a simple model by comparing with the regular geometry. When constituting wave front by the interference of elemental waves, pseudo-random geometry data did not always provide good results. In the case of a point diffractor, the imaging operation, where the constituted wave front was returned to the point diffractor by the interference of elemental waves for the spatial alias records, did not always give clear images. In the case of multi point diffractor, good images were obtained with less noise generation in spite of alias records. There are a lot of diffractors in the actual geological structures, which corresponds to the case of multi point diffractors. Finally, better images could be obtained by inputting records acquired using the pseudo-random geometry rather than by inputting spatial alias records acquired using the regular geometry. 7 refs., 6 figs.

  3. 3D Seismic Investigations at Olkiluoto 2007 Factual Report

    This document describes the activities related to 3D seismic investigations done at Olkiluoto, Finland, in 2007. The investigations were done in order to infer the position and to assess the continuity of gently dipping features in an area at the South West of the 3D survey acquired in 2006. The receiver template neither extends over the ONKALO area, nor is it adjacent to the 3D patch acquired during 2006. Site preparation, mobilisation, data acquisition and demobilisation have been done from April 15th until June 10th. 581 shots have been recorded by 240 channels. The area of investigations is centered on East = 1526530, North = 6791770. Posiva performed staking, according to the survey plans designed by Vibrometric. True coordinates were measured with +/-0.5 m accuracy for all stations. The source used was a VIBSIST-1000 mounted on a 24-tone Liebherr tracked excavator, with a Krupp HM 720-4000 hydraulic breaker. A Summit II Plus 24-bit seismic data recording system, with 186 channels and a Summit II system with 64 channels were used for data acquisition. The geophones used were 28 Hz OYO Geospace (SNG20DM/PC-710/KCL-4F/33'). The data have been processed in November 2007 by Vibrometric Oy and Uppsala University. The main processing platform has been Claritas 3D Seismic Processing Software. A thorough noise coherency analysis was performed, because of the relatively large number of raw shot gathers on which the arrival times could not be picked due to noise contamination. Electric, weather and cultural noise patterns were identified and removed. The data were binned into 8.33 m x 12.5 m CDP bins. Statics, both refraction and residual, and choice of temporal filter had the greatest influence on the processed stack. The stack was relatively insensitive to the velocity used. DMO helped somewhat in imaging dipping reflections. Data from the 3D cube has been migrated using Stolt and FD migration algorithms. There are three main zones of reflectivity, the first at 300

  4. 3D Seismic Investigations at Olkiluoto 2007 Factual Report

    Cosma, C.; Cozma, M.; Enescu, N. (Vibrometric Oy, Vantaa (Finland))

    2008-08-15

    This document describes the activities related to 3D seismic investigations done at Olkiluoto, Finland, in 2007. The investigations were done in order to infer the position and to assess the continuity of gently dipping features in an area at the South West of the 3D survey acquired in 2006. The receiver template neither extends over the ONKALO area, nor is it adjacent to the 3D patch acquired during 2006. Site preparation, mobilisation, data acquisition and demobilisation have been done from April 15th until June 10th. 581 shots have been recorded by 240 channels. The area of investigations is centered on East = 1526530, North = 6791770. Posiva performed staking, according to the survey plans designed by Vibrometric. True coordinates were measured with +/-0.5 m accuracy for all stations. The source used was a VIBSIST-1000 mounted on a 24-tone Liebherr tracked excavator, with a Krupp HM 720-4000 hydraulic breaker. A Summit II Plus 24-bit seismic data recording system, with 186 channels and a Summit II system with 64 channels were used for data acquisition. The geophones used were 28 Hz OYO Geospace (SNG20DM/PC-710/KCL-4F/33'). The data have been processed in November 2007 by Vibrometric Oy and Uppsala University. The main processing platform has been Claritas 3D Seismic Processing Software. A thorough noise coherency analysis was performed, because of the relatively large number of raw shot gathers on which the arrival times could not be picked due to noise contamination. Electric, weather and cultural noise patterns were identified and removed. The data were binned into 8.33 m x 12.5 m CDP bins. Statics, both refraction and residual, and choice of temporal filter had the greatest influence on the processed stack. The stack was relatively insensitive to the velocity used. DMO helped somewhat in imaging dipping reflections. Data from the 3D cube has been migrated using Stolt and FD migration algorithms. There are three main zones of reflectivity, the first at

  5. ONKALO 3D tunnel seismic investigations, Olkiluoto 2013

    Enescu, N.; Cosma, C.; Crawford, J. [Vibrometric Oy Cosma, Perttula (Finland)

    2014-10-15

    POSIVA Oy conducts bedrock investigations at the spent nuclear fuel final disposal site at Olkiluoto, in western Finland. The purpose of these efforts, which include a significant R and D component, is to ensure compliance with the requirements set forth for the long-term safety of final disposal. The excavation of the access tunnel to the repository hosts the ONKALO underground rock characterization facility. The investigations carried out in ONKALO focus on the bedrock and groundwater conditions prevailing on the final disposal site and how construction work affects them. Deformation zones and hydraulically conductive zones can limit the areas which are suitable for a deposition hole to be placed. The main objective of the tunnel seismic investigations presented here is to develop a seismic investigation technique for deposition area characterization. The field work consisted of 3 receiver lines using 3-component geophones and sources locations using the Vibsist-500 for 6 profiles, the Vibsist-20 for 1 profile (sparse), and a mechanical borehole hammer, MH-70, for 5 profiles for a total of 12 profiles. This investigation took place in the ONKALO demonstration area. Design of the seismic field work is based on previous tunnel seismic works in ONKALO (Cosma et al. 2008 and 2011) and pre-survey numerical modeling for field data and processing (Heinonen et al. 2013). The fieldwork was carried out in June 2013. Tomographic inversion of P-wave first arrival times was done on combined data sets, containing data from several shot-receiver configurations along tunnels and boreholes. The velocity distributions derived by tomographic inversion were verified by comparison with the forward modeled profiles. The Image Point (IP) migration method is characterized by is its ability to accumulate reflection events measured in the time distance data sets into points in the IP domain, which permits the enhancement of coherent backscattered events. The result of the 3D IP migration

  6. ONKALO 3D tunnel seismic investigations, Olkiluoto 2013

    POSIVA Oy conducts bedrock investigations at the spent nuclear fuel final disposal site at Olkiluoto, in western Finland. The purpose of these efforts, which include a significant R and D component, is to ensure compliance with the requirements set forth for the long-term safety of final disposal. The excavation of the access tunnel to the repository hosts the ONKALO underground rock characterization facility. The investigations carried out in ONKALO focus on the bedrock and groundwater conditions prevailing on the final disposal site and how construction work affects them. Deformation zones and hydraulically conductive zones can limit the areas which are suitable for a deposition hole to be placed. The main objective of the tunnel seismic investigations presented here is to develop a seismic investigation technique for deposition area characterization. The field work consisted of 3 receiver lines using 3-component geophones and sources locations using the Vibsist-500 for 6 profiles, the Vibsist-20 for 1 profile (sparse), and a mechanical borehole hammer, MH-70, for 5 profiles for a total of 12 profiles. This investigation took place in the ONKALO demonstration area. Design of the seismic field work is based on previous tunnel seismic works in ONKALO (Cosma et al. 2008 and 2011) and pre-survey numerical modeling for field data and processing (Heinonen et al. 2013). The fieldwork was carried out in June 2013. Tomographic inversion of P-wave first arrival times was done on combined data sets, containing data from several shot-receiver configurations along tunnels and boreholes. The velocity distributions derived by tomographic inversion were verified by comparison with the forward modeled profiles. The Image Point (IP) migration method is characterized by is its ability to accumulate reflection events measured in the time distance data sets into points in the IP domain, which permits the enhancement of coherent backscattered events. The result of the 3D IP migration

  7. AutoCAD Civil 3D - Survey

    Luketić, Antonio; Padovan, Ivan

    2011-01-01

    AutoCAD Civil 3D je vrlo kompleksan i napredan program za geodeziju i građevinarstvo. Survey izbornik predstavlja mali dio tog softvera. U programu je sadržana mogućnost kreiranja vlastitih kodova kao i naredbi kojima možemo ubrzati vizualizaciju izmjerenih objekata. Njegovim korištenjem se znatno skraćuje vrijeme obrade podataka. Članak je napisan na temelju jednostavnog, ali korisnog primjera kojime smo željeli upoznati čitatelje s prilično nepoznatom aplikacijom koju nudi AutoCAD Civil 3D....

  8. 3-D seismic mapping and amplitude analysis: a Gulf of Mexico case history

    Kidney, R.L.; Silver, R.S.; Hussein, H.A. (Oryx Energy Co. (United States))

    1992-01-01

    Utilization of 3-D seismic data and Direct Hydrocarbon Indicators led to the successful drilling of appraisal and development wells in the Gulf of Mexico block South Timbalier 198 (ST 198), significantly reducing time and cost. Based on 2-D seismic mapping, a Pliocene Lower Buliminella (L BUL) prospect was drilled in ST 198. An Upper Buliminella (U BUL) gas sandstone was encountered instead. An appraisal well of the U BUL interval confirmed this discovery. It became apparent that the structural complexities and the seismic amplitude anomalies of the area could not be adequately resolved using the 2-D seismic grid. A 3-D seismic survey was shot to investigate the remaining potential of ST 198. Direct Hydrocarbon Indicators (DHIs), which are seismic anomalies resulting from the hydrocarbon effect on rock properties, are generally expected from these age sands. While the 3-D survey shows a seismic amplitude anomaly associated with the U BUL reservoir, the areal extent of the anomaly did not match the findings of the two wells. A DHI study was performed to determine if this inconsistency could be explained and if the anomaly could be used in the well planning. The two key steps which confirmed that this amplitude anomaly is a DHI were properly calibrating the seismic data to the well control and determining the theoretical seismic response of the gas sandstones. The DHI study along with the 3-D mapping led to the successful development of the ST 198 U BUL reservoir and to setting up a successful adjacent fault block play. Finally, 3-D mapping also identified a L BUL trap updip from the original L BUL prospect which resulted in a successful drilling effort. (author).

  9. Application of seismic attributes at the 3D seismic interpretation

    Bibiana Kytková

    2009-01-01

    Full Text Available This article brings a short overview of commonly used seismic attributes with their basic characterization. It also discusses the possibility of their application for the hydrocarbon prospection in two different areas. One of them is the Vienna Basin and the other one is the contact zone of the Bohemian Massif and the Western Carpathians (Vranovice and Nesvačilka Trough.

  10. 3D Modelling of Seismically Active Parts of Underground Faults via Seismic Data Mining

    Frantzeskakis, Theofanis; Konstantaras, Anthony

    2015-04-01

    During the last few years rapid steps have been taken towards drilling for oil in the western Mediterranean sea. Since most of the countries in the region benefit mainly from tourism and considering that the Mediterranean is a closed sea only replenishing its water once every ninety years careful measures are being taken to ensure safe drilling. In that concept this research work attempts to derive a three dimensional model of the seismically active parts of the underlying underground faults in areas of petroleum interest. For that purpose seismic spatio-temporal clustering has been applied to seismic data to identify potential distinct seismic regions in the area of interest. Results have been coalesced with two dimensional maps of underground faults from past surveys and seismic epicentres, having followed careful reallocation processing, have been used to provide information regarding the vertical extent of multiple underground faults in the region of interest. The end product is a three dimensional map of the possible underground location and extent of the seismically active parts of underground faults. Indexing terms: underground faults modelling, seismic data mining, 3D visualisation, active seismic source mapping, seismic hazard evaluation, dangerous phenomena modelling Acknowledgment This research work is supported by the ESPA Operational Programme, Education and Life Long Learning, Students Practical Placement Initiative. References [1] Alves, T.M., Kokinou, E. and Zodiatis, G.: 'A three-step model to assess shoreline and offshore susceptibility to oil spills: The South Aegean (Crete) as an analogue for confined marine basins', Marine Pollution Bulletin, In Press, 2014 [2] Ciappa, A., Costabile, S.: 'Oil spill hazard assessment using a reverse trajectory method for the Egadi marine protected area (Central Mediterranean Sea)', Marine Pollution Bulletin, vol. 84 (1-2), pp. 44-55, 2014 [3] Ganas, A., Karastathis, V., Moshou, A., Valkaniotis, S., Mouzakiotis

  11. 3D seismics for geothermal reservoir characterization - a case study from Schneeberg (Germany)

    Hlousek, F.; Hellwig, O.; Buske, S.

    2013-12-01

    We present the results of a 3D seismic survey acquired near Schneeberg in the western Erzgebirge (Germany). The aim of the project is to use seismic exploration methods to image and to characterize a major fault zone in crystalline rock which could be used as a geothermal reservoir at a target depth of about 5-6 km with expected temperatures between 160-180°C. For this purpose a high resolution 3D Vibroseis survey with more than 5300 source and approximately 8000 receiver locations was performed at the end of 2012 and covered an area of approximately 10 km x 13 km. The 3D survey was complemented by an additional wide-angle seismic survey using explosives along eleven profile lines radially centered at the target area. The region itself is dominated by the NW-SE striking Gera-Jáchymov fault system. The main geological features in the survey area are well known from intensive mining activities down to a depth of about 2 km. The seismic investigations aimed at imaging the partly steeply dipping fault branches at greater depth, in particular a dominant steeply NE dipping fault in the central part of the survey area. Beside this main structure, the Gera-Jáchymov fault zone consists of a couple of steeply SW dipping conjugated faults. Advanced processing and imaging methods have been applied to the data set. 3D Kirchhoff prestack depth migration delivered a clear image of the structure of the various fault branches at depths of around 2-5 km. Furthermore, focusing migration methods (e.g. coherency migration) have been applied and even sharpened the image such that the 3D seismic result allows for a profound characterization of this potential geothermal reservoir in crystalline rock.

  12. Response of 3D Free Rigid Objects under Seismic Excitations

    Yanheng, Li

    2008-01-01

    Previous studies of precariously balanced structures in seismically active regions to provide important information for aseismatic engineering and theoretical seismology are almost found on an oversimplified assumption. According to that, any 3-dimensional practical structure with special symmetry could be regarded as a 2-dimensional finite object in light of the corresponding symmetry. Thus the complex and troublesome problem of 3D rotation, in mathematics, can be reduced to a tractable one of 1D rotation but a distorted description of the real motion in physics. To gain an actual evolution of precariously balanced structures bearing various levels of ground accelerations, we should address ourselves to a 3D calculation. In this study, the responses of a cylinder under a set of half- and full-sine-wave excitations with different frequencies related to seismic ground motion are investigated in virtue of some reasonable works from a number of mechanicians. A computer program is also developed possibly to study...

  13. Advanced computational tools for 3-D seismic analysis

    Barhen, J.; Glover, C.W.; Protopopescu, V.A. [Oak Ridge National Lab., TN (United States)] [and others

    1996-06-01

    The global objective of this effort is to develop advanced computational tools for 3-D seismic analysis, and test the products using a model dataset developed under the joint aegis of the United States` Society of Exploration Geophysicists (SEG) and the European Association of Exploration Geophysicists (EAEG). The goal is to enhance the value to the oil industry of the SEG/EAEG modeling project, carried out with US Department of Energy (DOE) funding in FY` 93-95. The primary objective of the ORNL Center for Engineering Systems Advanced Research (CESAR) is to spearhead the computational innovations techniques that would enable a revolutionary advance in 3-D seismic analysis. The CESAR effort is carried out in collaboration with world-class domain experts from leading universities, and in close coordination with other national laboratories and oil industry partners.

  14. Offshore 3D seismic, geochemical data integration, Main Pass project, Gulf of Mexico

    Belt, J.Q. Jr.; Rice, G.K. [GeoFrontiers Corp., Dallas, TX (United States)

    1996-04-01

    A bottom cable, 3D seismic, and shallow piston-core geochemical survey was conducted in summer 1992 on Main Pass Blocks 41 and 58 in the Gulf of Mexico. The 15.5 sq mile study area is in 11--24 m (36--79 ft) of water. Two hundred seventy-six shallow, 2 m piston-core sediment samples were collected using a uniform grid pattern on 400 m (1,320 ft) spaced intervals. Retrieved sediment samples were immediately sealed in metal 1 pt cans containing biocide solution. The purpose of the 3D seismic survey and geochemical offshore data integration project was three-fold: determine if near-shore, low-cost shallow piston-core sediment samples would be affected by fluvial contamination; evaluate the efficiency of a shallow-core, dense-grid sample design program in detecting thermogenic hydrocarbons at depth; evaluate the benefits of integrating offshore, shallow sediment geochemistry with 3D seismic data in developing a petroleum geological model. All 3D seismic survey data and block boundaries, as described in illustrations and figures, are only generally located within the Main Pass Area.

  15. Using decision analysis to estimate 3-D seismic value-Minas field, Sumatra Indonesia

    Mangold, K.M.; Whitacre, T.P.; Seffibudianti (Caltex Pacific Indonesia, Sumatra (Indonesia))

    1996-01-01

    Decision Analysis has been used to estimate the value added from a 3-D seismic survey recorded over Minas field, Central Sumatra. The method involves comparing the expected values which result from the various decision options, such as acquiring 3-D or not. Probabilities must be assigned to the various branches of the decision tree. These include for example, the expected reliability of the 3-D data as well as the subsequent interpretation. Anticipated drilling results with and without 3-D are assessed after reviewing historical data and interviewing experts to obtain 10th, 50th and 90th percentile results for various scenarios. In this way the expected value, or cumulative distribution of the expected value of the 3-D can be computed and risk can be assessed. The Minas 3-D survey is the largest (450 square kilometers) of over 25 development 3-D surveys recorded by Caltex Pacific Indonesia (CPI) over its fields in Central Sumatra. This survey was conducted after nearly 50 years of production from more than 750 wells. CPI's 3-D experience has shown that increasing the subsurface resolution within complex high angle faulted areas results in new drilling locations in older mature fields such as Minas. Better knowledge of the oil producing reservoirs can also be used to optimize pattern waterflood locations, horizontal drilling and other tertiary recovery studies.

  16. Using decision analysis to estimate 3-D seismic value-Minas field, Sumatra Indonesia

    Mangold, K.M.; Whitacre, T.P.; Seffibudianti [Caltex Pacific Indonesia, Sumatra (Indonesia)

    1996-12-31

    Decision Analysis has been used to estimate the value added from a 3-D seismic survey recorded over Minas field, Central Sumatra. The method involves comparing the expected values which result from the various decision options, such as acquiring 3-D or not. Probabilities must be assigned to the various branches of the decision tree. These include for example, the expected reliability of the 3-D data as well as the subsequent interpretation. Anticipated drilling results with and without 3-D are assessed after reviewing historical data and interviewing experts to obtain 10th, 50th and 90th percentile results for various scenarios. In this way the expected value, or cumulative distribution of the expected value of the 3-D can be computed and risk can be assessed. The Minas 3-D survey is the largest (450 square kilometers) of over 25 development 3-D surveys recorded by Caltex Pacific Indonesia (CPI) over its fields in Central Sumatra. This survey was conducted after nearly 50 years of production from more than 750 wells. CPI`s 3-D experience has shown that increasing the subsurface resolution within complex high angle faulted areas results in new drilling locations in older mature fields such as Minas. Better knowledge of the oil producing reservoirs can also be used to optimize pattern waterflood locations, horizontal drilling and other tertiary recovery studies.

  17. Identifying High Potential Well Targets with 3D Seismic and Mineralogy

    Mellors, R. J. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

    2015-10-30

    Seismic reflection the primary tool used in petroleum exploration and production, but use in geothermal exploration is less standard, in part due to cost but also due to the challenges in identifying the highly-permeable zones essential for economic hydrothermal systems [e.g. Louie et al., 2011; Majer, 2003]. Newer technology, such as wireless sensors and low-cost high performance computing, has helped reduce the cost and effort needed to conduct 3D surveys. The second difficulty, identifying permeable zones, has been less tractable so far. Here we report on the use of seismic attributes from a 3D seismic survey to identify and map permeable zones in a hydrothermal area.

  18. Mature field development using 3-D seismic in Minas field, Central Sumatra, Basin, Indonesia

    Fitris, F.; Whitacre, T.P. (PT Caltex Pacific Indonesia, Sumatra (Indonesia))

    1996-01-01

    The Minas field 3-D seismic survey (450 km2) was conducted after nearly 4000 MMBO of production. Interpretation of the data has shown a significant number of remaining prospects within the field. Delineation of these locations was impossible using the 2D seismic grid. The new 3D data cube has revealed locations in simple anticlines, en-echelon folds, faulted monoclines, rollover folding, complex faulted zones, stratigraphic traps and within a fractured basement high. The new interpretation has identified more than 75 development and OSWC well locations. These wells will produce an estimated 75 MMBO. As the result of a $5MM seismic survey, the asset value will increase by nearly $1000MM.

  19. Mature field development using 3-D seismic in Minas field, Central Sumatra, Basin, Indonesia

    Fitris, F.; Whitacre, T.P. [PT Caltex Pacific Indonesia, Sumatra (Indonesia)

    1996-12-31

    The Minas field 3-D seismic survey (450 km2) was conducted after nearly 4000 MMBO of production. Interpretation of the data has shown a significant number of remaining prospects within the field. Delineation of these locations was impossible using the 2D seismic grid. The new 3D data cube has revealed locations in simple anticlines, en-echelon folds, faulted monoclines, rollover folding, complex faulted zones, stratigraphic traps and within a fractured basement high. The new interpretation has identified more than 75 development and OSWC well locations. These wells will produce an estimated 75 MMBO. As the result of a $5MM seismic survey, the asset value will increase by nearly $1000MM.

  20. 3D seismic travel time surveying - a comparison of the time-term method and tomography (an example from an archaeological site)

    Valenta, Jan; Dohnal, J.

    2007-01-01

    Roč. 63, č. 1 (2007), s. 46-58. ISSN 0926-9851 Institutional research plan: CEZ:AV0Z30460519 Keywords : shallow seismic * tomography * archaeology Subject RIV: DC - Siesmology, Volcanology, Earth Structure Impact factor: 0.938, year: 2007

  1. 3D seismic data for shallow aquifers characterisation

    Giustiniani, Michela; Accaino, Flavio; Picotti, Stefano; Tinivella, Umberta

    2009-07-01

    We present the results obtained from conventional and non-conventional analysis of 3D high-resolution seismic data acquired nearby the water spring line, which separates the upper from the lower Friuli-Venezia Giulia plain (Italy), in order to characterise an important multilayered confined aquifer. The main targets of this study were two shallow aquifers located at about 30 m and 180 m depth, respectively. The aquifer structures were reconstructed by adopting a technique consisting an iterative updating procedure, for refining and improving an initial model in depth. The method includes pre-stack depth migration, residual move-out analysis and seismic reflection tomography. In the final 3D migrated cubes, two high velocity layers were identified at about 270 m and 480 m respectively, which correspond to unknown deep aquifers, as confirmed by recent well data (stratigraphies and down-hole velocity measurements). Travel-time tomography and Amplitude Versus Offset analysis evidence that seasonal variation in the seismic response of the aquifers are not detectable. However, in this case, aquifers are well detectable by lithological changes.

  2. The use of exploration 3D seismic data for geohazards assessment in the Gulf of Mexico

    Gafford, W.T.

    1996-12-31

    3D seismic data, acquired for oil and gas exploration purposes, is now being used to supplement or even replace conventional high resolution geohazard surveys in some areas of the Gulf of Mexico. The use of 3D data has improved the identification and understanding of some types of geohazards, and the use of seismic interpretation workstations has made the data more accessible for geohazard analysis, and resulted in a more thorough interpretation of the shallow geologic section. Technologies used for exploration purposes have been adopted for use in the identification and assessment of drilling hazards in the near seafloor sedimentary section. Savings of time and money have been realized, as well as identification of areas where additional data is required. Some of the limitations and restrictions in using exploration 3D data have also been identified.

  3. Exploring the seismic expression of fault zones in 3D seismic volumes

    Iacopini, D.; Butler, R. W. H.; Purves, S.; McArdle, N.; De Freslon, N.

    2016-08-01

    Mapping and understanding distributed deformation is a major challenge for the structural interpretation of seismic data. However, volumes of seismic signal disturbance with low signal/noise ratio are systematically observed within 3D seismic datasets around fault systems. These seismic disturbance zones (SDZ) are commonly characterized by complex perturbations of the signal and occur at the sub-seismic (10 s m) to seismic scale (100 s m). They may store important information on deformation distributed around those larger scale structures that may be readily interpreted in conventional amplitude displays of seismic data. We introduce a method to detect fault-related disturbance zones and to discriminate between this and other noise sources such as those associated with the seismic acquisition (footprint noise). Two case studies from the Taranaki basin and deep-water Niger delta are presented. These resolve SDZs using tensor and semblance attributes along with conventional seismic mapping. The tensor attribute is more efficient in tracking volumes containing structural displacements while structurally-oriented semblance coherency is commonly disturbed by small waveform variations around the fault throw. We propose a workflow to map and cross-plot seismic waveform signal properties extracted from the seismic disturbance zone as a tool to investigate the seismic signature and explore seismic facies of a SDZ.

  4. DMO processing on the Ketzin 3D seismic data

    Huang, Fei; Juhlin, Christopher; Ivandic, Monika; Zhang, Fengjiao

    2014-05-01

    The Dip-moveout (DMO) correction is a process which attempts to make the finite offset data closer to zero offset data after the normal-moveout (NMO) correction. The NMO correction is then dip independent and reflections with different dips will stack coherently. DMO plays a critical role in seismic processing by enhancing the final image quality of the seismic data. In this study, we apply 3D Squeezing DMO (Hale and Artley, 1993) to seismic data from the Ketzin pilot CO2 site after NMO to study the impact of DMO on time-lapse seismic imaging and to investigate if it enhances the CO2 seismic monitoring technique. This 3D DMO method is based on an integral approach and incorporates Hale and Artley's (1993) modifications for variable velocity with time. A constant velocity algorithm is used with a gamma correction function which depends on the velocity function. An anti-alias velocity of 3000 m/s is used for the DMO. After DMO the data are stacked and F-XY deconvolution is applied. Finally, 3D finite-difference migration using the final smoothed NMO velocities is performed for each data set. We then apply a time-lapse analysis to the 3D seismic data sets and compare the results with and without DMO processing. The most important aspect of the DMO processing is determining the velocity field for the NMO step. This is done by using the initial smoothed velocity field obtained from the conventional velocity analysis before DMO as a first estimate. The data are input into the DMO process and then inverse NMO is applied. These data are then subjected to a new velocity analysis and the velocity field is updated and used as input for the NMO process. A number of iterations are generally required until the velocity field does not need further updating. In this study velocities were picked at every 20th CDP in the inline and crossline directions. Compared to the velocity spectrum without DMO processing, the velocity trend is improved and the ambiguity in the velocity picks is

  5. 3D seismic Unterhaching 2009 within hydrothermal exploration and modelling; 3D-Seismik Unterhaching 2009 im Rahmen hydrothermaler Exploration und Modellierung

    Lueschen, Ewald; Dussel, Michael; Thomas, Ruediger; Schulz, Ruediger [Leibniz-Institut fuer Angewandte Geophysik (LIAG), Hannover (Germany)

    2011-10-24

    Within the exploration of hydrothermal reservoirs, results of 3D reflexion-seismic measurements are presented. These measurements were performed in June / July 2009 according to the vibroseis method on an area of 26.3 square kilometers in the area Unterhaching (Federal Republic of Germany). The 3D seismic survey exhibits much more complex structures than previously known by 2D seismic lines. Subsequent to sinistral transtension (active in the Cretaceous to the Eocene) a short transpression impetus was performed. This is evident from graduated normal faults as well as staggered reverse fault structures and inversion structures in the Upper Jurassic. Top and base of the 600-650 m mighty Malm are well resolved. Brittle fault structures are formed linearly at the top Malm but rounded and chaotic within the Malm. This can be explained by a radical karstification / hydrothermal solution. Several circular structures are interpreted as karstified incursion structures. The seismic facies of the Malm is characterized by a shift from relatively transparent zones, layered fields, scatters and fault zones. This is an expression of smaller and larger reefs, lagoons and reef debris. Reefs are characterized by several seismic attributes. Striking low-velocity zones are oriented along the main fault zones and can be interpreted as zones that are relieved by gap porosity. Azimuth variable processing gives evidence for preferred orientations of fractures on the seismic scale. By means of the 3D seismic diverse geothermal exploration targets can be defined.

  6. Prestack reverse time migration for 3D marine reflection seismic data

    Prestack reverse time migration (RTM) is a method for imaging the subsurface using the inner product of wavefield extrapolation in shot domain and in receiver domain. It is well known that RTM is better for preserving amplitudes and phases than other prestack migrations. Since 3D seismic data is huge data volume and it needs heavy computing works, it requires parallel computing in order to have a meaningful depth image of the 3D subsurface. We implemented a parallelized version of 3D RTM for prestack depth migration. The results of numerical example for 3D SEG/EAGE salt model showed good agreement with the original geological model. We applied RTM to offshore 3D seismic reflection data. The study area is 12 × 25 km with 120 survey lines. Shot and receiver spacing is 25 m and 12.5 m. The line spacing is 100 m. Shot gathers were preprocessed to enhance signal to noise ratio and velocity model was calculated from conventional stack velocity. Both of them were used to obtain 3D image using RTM. The results show reasonable subsurface image

  7. Advancing New 3D Seismic Interpretation Methods for Exploration and Development of Fractured Tight Gas Reservoirs

    James Reeves

    2005-01-31

    In a study funded by the U.S. Department of Energy and GeoSpectrum, Inc., new P-wave 3D seismic interpretation methods to characterize fractured gas reservoirs are developed. A data driven exploratory approach is used to determine empirical relationships for reservoir properties. Fractures are predicted using seismic lineament mapping through a series of horizon and time slices in the reservoir zone. A seismic lineament is a linear feature seen in a slice through the seismic volume that has negligible vertical offset. We interpret that in regions of high seismic lineament density there is a greater likelihood of fractured reservoir. Seismic AVO attributes are developed to map brittle reservoir rock (low clay) and gas content. Brittle rocks are interpreted to be more fractured when seismic lineaments are present. The most important attribute developed in this study is the gas sensitive phase gradient (a new AVO attribute), as reservoir fractures may provide a plumbing system for both water and gas. Success is obtained when economic gas and oil discoveries are found. In a gas field previously plagued with poor drilling results, four new wells were spotted using the new methodology and recently drilled. The wells have estimated best of 12-months production indicators of 2106, 1652, 941, and 227 MCFGPD. The latter well was drilled in a region of swarming seismic lineaments but has poor gas sensitive phase gradient (AVO) and clay volume attributes. GeoSpectrum advised the unit operators that this location did not appear to have significant Lower Dakota gas before the well was drilled. The other three wells are considered good wells in this part of the basin and among the best wells in the area. These new drilling results have nearly doubled the gas production and the value of the field. The interpretation method is ready for commercialization and gas exploration and development. The new technology is adaptable to conventional lower cost 3D seismic surveys.

  8. P-Cable: New High-Resolution 3D Seismic Acquisition Technology

    Planke, Sverre; Berndt, Christian; Mienert, Jürgen; Bünz, Stefan; Eriksen, Frode N.; Eriksen, Ola K.

    2010-05-01

    We have developed a new cost-efficient technology for acquisition of high-resolution 3D seismic data: the P-Cable system. This technology is very well suited for deep water exploration, site surveys, and studies of shallow gas and fluid migration associated with gas hydrates or leaking reservoirs. It delivers unparalleled 3D seismic images of subsurface sediment architectures. The P-Cable system consists of a seismic cable towed perpendicular to a vessel's steaming direction. This configuration allows us to image an up to 150 m wide swath of the sub-surface for each sail line. Conventional 3D seismic technology relies on several very long streamers (up to 10 km long streamers are common), large sources, and costly operations. In contrast, the P-Cable system is light-weight and fast to deploy from small vessels. Only a small source is required as the system is made for relatively shallow imaging, typically above the first water-bottom multiple. The P-Cable system is particularly useful for acquisition of small 3D cubes, 10-50 km2, in focus areas, rather than extensive mapping of large regions. The rapid deployment and recovery of the system makes it possible to acquire several small cubes (10 to 30 km2) with high-resolution (50-250 Hz) seismic data in during one cruise. The first development of the P-Cable system was a cooperative project achieved by Volcanic Basin Petroleum Research (VBPR), University of Tromsø, National Oceanography Centre, Southampton, and industry partners. Field trials using a 12-streamer system were conducted on sites with active fluid-leakage systems on the Norwegian-Barents-Svalbard margin, the Gulf of Cadiz, and the Mediterranean. The second phase of the development introduced digital streamers. The new P-Cable2 system also includes integrated tow and cross cables for power and data transmission and improved doors to spread the larger cross cable. This digital system has been successfully used during six cruises by the University of Troms

  9. Seismic Wave Amplification in 3D Alluvial Basins: 3D/1D Amplification Ratios from Fast Multipole BEM Simulations

    Fajardo, Kristel C Meza; Chaillat, Stéphanie; Lenti, Luca

    2016-01-01

    In this work, we study seismic wave amplification in alluvial basins having 3D standard geometries through the Fast Multipole Boundary Element Method in the frequency domain. We investigate how much 3D amplification differs from the 1D (horizontal layering) case. Considering incident fields of plane harmonic waves, we examine the relationships between the amplification level and the most relevant physical parameters of the problem (impedance contrast, 3D aspect ratio, vertical and oblique incidence of plane waves). The FMBEM results show that the most important parameters for wave amplification are the impedance contrast and the so-called equivalent shape ratio. Using these two parameters, we derive simple rules to compute the fundamental frequency for various 3D basin shapes and the corresponding 3D/1D amplification factor for 5% damping. Effects on amplification due to 3D basin asymmetry are also studied and incorporated in the derived rules.

  10. 3-D seismic facies analysis of a reefal buildup, offshore North Sumatra

    Alexander, W.L.; Nellia, M.R. (Mobil Oil Indonesia, Jakarta (Indonesia))

    1994-07-01

    The [open quotes]A[close quotes] field is located on the Sunda shelf, offshore north Sumatra. The A-1 discovery well, drilled in 1972, found hydrocarbon gas in middle Miocene carbonate rocks of reefal origin. Six appraisal wells were subsequently drilled, the most recent in late 1990. Because of drilling problems, mainly lost circulation in the carbonate reservoir, the well data obtained from the appraisal program was generally disappointing. Prior to development of the offshore area, an extensive 3-D seismic survey was shot, a portion of which covered the [open quotes]A[close quotes] field. Interpretation of the 3-D data over the [open quotes]A[close quotes] field identified different seismic facies within the carbonate reservoir. These seismic facies have been integrated with the geological data in order to construct a depositional model for the field. The seismic facies analysis was critical for developing the model because of the inadequate geological data obtained from the wells. Three distinct facies could be identified on the 3-D seismic data and correlated with the well data: reef, near-reef and inter-reef. The main concerns this facies mapping addressed were reserve determination, areas of severe lost circulation, and the distribution of dolomite. The near-reef and inter-reef areas were found to have better reservoir properties than the reef core, thereby impacting reserve calculations. In addition, the reef facies, with zones of vuggy to near cavernous type porosity, was correlatable to wells that had experienced severe lost circulation. Finally, dolomite was found to occur only within the reef facies, enabling its distribution to be predicted.

  11. 3D seismic reverse time migration on GPGPU

    Liu, Guofeng; Liu, Yaning; Ren, Li; Meng, Xiaohong

    2013-09-01

    Reverse time migration (RTM) is a powerful seismic imaging method for the interpretation of steep-dips and subsalt regions; however, implementation of the RTM method is computationally expensive. In this paper, we present a fast and computationally inexpensive implementation of RTM using a NVIDIA general purpose graphic processing unit (GPGPU) powered with Compute Unified Device Architecture (CUDA). To accomplish this, we introduced a random velocity boundary in the source propagation kernel. By creating a random velocity layer at the left, right, and bottom boundaries, the wave fields that encounter the boundary regions are pseudo-randomized. Reflections off the random layers have minimal coherent correlation in the reverse direction. This process eliminates the need to write the wave fields to a disk, which is important when using a GPU because of the limited bandwidth of the PCI-E that is connected to the CPU and GPU. There are four GPU kernels in the code: shot, receiver, modeling, and imaging. The shot and receiver insertion kernels are simple and are computed using a GPU because the wave fields reside in GPU's memory. The modeling kernel is computed using Micikevicius's tiling method, which uses shared memory to improve bandwidth usage in 2D and 3D finite difference problems. In the imaging kernel, we also use this tiling method. A Tesla C2050 GPU with 4 GB memory and 480 stream processing units was used to test the code. The shot and receiver modeling kernel occupancy achieved 85%, and the imaging kernel occupancy was 100%. This means that the code achieved a good level of optimization. A salt model test verified the correct and effective implementation of the GPU RTM code.

  12. Fuzzy C-means Clustering for 3D Seismic Parameters Processing

    Fuqun Zhao; Liang Le

    2009-01-01

    3D seismic parameters can reflect the features of petroleum reservoir from different profiles. By analizing the3D seismic parameters, we can assess the parameters of the reservoir characterization, such as deposition,structure and growth history, fluid saturation and so on. The traditional clustering methods can’t capture thedegree of similarity between reservoir parameters very well, so we introduced in this paper the application offuzzy C-means (FCM) clustering for the processing of 3D seis...

  13. Evaluation of the 3D high resolution seismic method at the Tournemire site around the IPSN experimental station

    The IPSN experimental station of Tournemire is localized at a 200 m depth inside an abandoned railway tunnel dug in a Jurassic clayey formation. The a priori knowledge of the existing geologic structures of the clayey formations allows to test the reliability of the 3D high resolution seismic survey technique and its capability to detect these structures and discontinuities. This test study is reported in this technical note. It comprises several steps: a bibliographic synthesis and a state-of-the-art of the 3D seismic survey technique, the construction of a velocity model for the different strata of the site, a simulation of the possible seismic response of these strata with respect to the velocities chosen, the processing of the data and finally their interpretation. (J.S.)

  14. Seismic performance of a novel 3D isolation system on continuous bridges

    Ou, J. P.; Jia, J. F.

    2010-04-01

    Remarkable vertical seismic motion is one of the prominent characteristics of the near-fault earthquake motions, but the traditional and widely used base isolation system only can effectively mitigate horizontal seismic responses and structural damage. A promising three-dimensional (3D) seismic isolation bearing, consisting of laminated rubber bearing with lead core (LRB) and combined coned disc spring with vertical energy dissipation device (e.g., inner fluid viscous cylindric damper or steel damper), was proposed to mitigate horizontal and vertical structural seismic responses simultaneously and separately. Three-group seismic ground motion records were selected to validate the effectiveness of the proposed 3D seismic isolation bearing on a continuous slab bridge. The appropriate damping of the vertical damping device was presented by parametric study. The analyses results showed that the proposed 3D isolation bearing is essentially effective to mitigate vertical and horizontal structural seismic response simultaneously. Near-fault pulse-type seismic motions should be considered in seismic isolation design and evaluation. The proper damping ratio of the vertical damping device should be 20%-30% for favorable vertical isolation effectiveness. The proposed 3D seismic isolation bearing is promising to be applied to the mediate-to-short span bridge and even some building structures.

  15. Assessing a 3D smoothed seismicity model of induced earthquakes

    Zechar, Jeremy; Király, Eszter; Gischig, Valentin; Wiemer, Stefan

    2016-04-01

    As more energy exploration and extraction efforts cause earthquakes, it becomes increasingly important to control induced seismicity. Risk management schemes must be improved and should ultimately be based on near-real-time forecasting systems. With this goal in mind, we propose a test bench to evaluate models of induced seismicity based on metrics developed by the CSEP community. To illustrate the test bench, we consider a model based on the so-called seismogenic index and a rate decay; to produce three-dimensional forecasts, we smooth past earthquakes in space and time. We explore four variants of this model using the Basel 2006 and Soultz-sous-Forêts 2004 datasets to make short-term forecasts, test their consistency, and rank the model variants. Our results suggest that such a smoothed seismicity model is useful for forecasting induced seismicity within three days, and giving more weight to recent events improves forecast performance. Moreover, the location of the largest induced earthquake is forecast well by this model. Despite the good spatial performance, the model does not estimate the seismicity rate well: it frequently overestimates during stimulation and during the early post-stimulation period, and it systematically underestimates around shut-in. In this presentation, we also describe a robust estimate of information gain, a modification that can also benefit forecast experiments involving tectonic earthquakes.

  16. A 3D surface seismic pilot study at Olkiluoto, Finland. Aquisition and processing report

    Posiva carries out R and D tasks related to spent nuclear fuel disposal in Finland. Works are currently performed in the Olkiluoto Island, in western Finland. Multidisciplinary engineering and geological investigations have been carried out since 1987. Site and rock characterization continued underground after 2004, in the ONKALO underground research facility. Posiva co-operated with Andra, France (Agence Nationale des Dechets Radioactifs) to purchase the design, acquisition and processing work needed to demonstrate the capability of 3D seismic techniques to image geological and structural features in crystalline rock. A 3D seismic pilot survey was carried out at the Olkiluoto site during seven days in June 2006. Data acquisition and processing were carried out as a Finnish-Swedish joint effort of Vibrometric and Uppsala University. The seismic survey covered 600 x 650 m area to the west of ONKALO underground premises. The single-template fixed array consisted of 270 active geophones on nine NW-SE lines. Geophone line interval was 60 m and station interval 24 m (30 channels on line). Sources (469 in total) were placed on seven NE-SW lines perpendicular to the receiver lines. The shot line interval was 100 m and the shot spacing 10 m (71 shots per line). The source was a time distributed Vibsist impact source (2000-2500 J/ impact), using a hydraulic rock breaker. The receivers (28 Hz geophones) were connected to a Sercel 408 UL seismograph. The coverage and signal quality were good. The array size was small compared to imaging depth, compared with conventional 3D works. The fold (CDP bin size 12 x 12 m) was highly variable (5->100), lower at the edges of area. Reflections are already seen in the processed receiver and shot gathers. The processing comprised time stacking (impact decoding), first arrival picking, resampling, refraction and residual static corrections, spectral equalisation, band-pass filtering, NMO and DMO, stacking. Static corrections, DMO and

  17. A 3D surface seismic pilot study at Olkiluoto, Finland. Aquisition and processing report

    Juhlin, C. (Uppsala University (SE)); Cosma, C. (Vibrometric Oy, Vantaa (FI))

    2007-08-15

    Posiva carries out R and D tasks related to spent nuclear fuel disposal in Finland. Works are currently performed in the Olkiluoto Island, in western Finland. Multidisciplinary engineering and geological investigations have been carried out since 1987. Site and rock characterization continued underground after 2004, in the ONKALO underground research facility. Posiva co-operated with Andra, France (Agence Nationale des Dechets Radioactifs) to purchase the design, acquisition and processing work needed to demonstrate the capability of 3D seismic techniques to image geological and structural features in crystalline rock. A 3D seismic pilot survey was carried out at the Olkiluoto site during seven days in June 2006. Data acquisition and processing were carried out as a Finnish-Swedish joint effort of Vibrometric and Uppsala University. The seismic survey covered 600 x 650 m area to the west of ONKALO underground premises. The single-template fixed array consisted of 270 active geophones on nine NW-SE lines. Geophone line interval was 60 m and station interval 24 m (30 channels on line). Sources (469 in total) were placed on seven NE-SW lines perpendicular to the receiver lines. The shot line interval was 100 m and the shot spacing 10 m (71 shots per line). The source was a time distributed Vibsist impact source (2000-2500 J/ impact), using a hydraulic rock breaker. The receivers (28 Hz geophones) were connected to a Sercel 408 UL seismograph. The coverage and signal quality were good. The array size was small compared to imaging depth, compared with conventional 3D works. The fold (CDP bin size 12 x 12 m) was highly variable (5->100), lower at the edges of area. Reflections are already seen in the processed receiver and shot gathers. The processing comprised time stacking (impact decoding), first arrival picking, resampling, refraction and residual static corrections, spectral equalisation, band-pass filtering, NMO and DMO, stacking. Static corrections, DMO and

  18. Application Results of 3-D Seismic Exploration Technology in Coal Mines

    SUN Shenglin; WU Xizun

    2004-01-01

    This paper briefly introduces the development and present situation of China's coal seismic exploration. It focuses on analyzing the important functions of 3-D seismic exploration technology in the designing and production of coal mines, and also the results of its application.

  19. BGP Started Shooting Its Largest Oversea Seismic Survey

    2005-01-01

    @@ At 5 a.m. on November 1, 2005 of Beijing time, BGP (Bureau of Geophysical Prospecting)formally began shooting the S-51 3D seismic acquisition project in Saudi Arabia, its largest oversea seismic survey project to date. At the same time, the SRAK 2D seismic survey project was also started in the country. These two projects mark the construction of BGP's another large oversea production base.

  20. 3D seismic imaging around the 2.5 km deep COSC-1 scientific borehole, central Sweden

    Hedin, Peter; Juhlin, Christopher; Buske, Stefan

    2015-04-01

    Following the successful completion of the COSC-1 drilling campaign, a number of geophysical investigations have been performed in and around the 2.5 km deep borehole. Three different seismic experiments were conducted simultaneously in the fall of 2014 to take advantage of the same source points; 1) a Vertical Seismic Profile (VSP) in the borehole, 2) three 2D seismic profiles across the borehole, and 3) a limited 3D seismic survey (presented here). The latter is the first 3D seismic survey on land in Scandinavia to target the Caledonian Nappes and will allow mapping a small part of the Seve Nappe Complex (SNC) in 3D. Furthermore, it will allow extrapolation of results from downhole logging, core analysis and other seismic surveys to structures surrounding the borehole. A total number of 429 receivers (10 Hz single component geophones) were planted with 20 m separation along 7 lines spaced 200 m apart. The total area with receivers covered approximately 1.5 km2 and was centered on the drill site. A combination of a mechanical source (a rock breaking hydraulic hammer, near offsets) and explosive charges (0.5 kg fired at 3.5 - 5 m depth, far offsets) were used. The source points were activated along roads radiating outwards from the COSC-1 drill site in a star pattern. The nominal shot spacing was 20 m (vibrating source) or 80 m (explosives) and maximum horizontal offset was about 5.75 km. The high-grade metamorphic SNC is well known from previous 2D seismic studies to be a highly reflective unit. However, due to the complex 3D geometry and lithological variation within the unit, it has not been clearly imaged. The new 3D data provide a means to image these structures in more detail and to follow the lithological and structural interfaces observed in the core into the surrounding unit. Preliminary results from the 3D processing and correlation with borehole data will be presented.

  1. 3D, 9-C anisotropic seismic modeling and inversion

    Rusmanugroho, Herurisa

    The most complete representation of an elastic medium consists of an elastic tensor with 21 independent moduli. All 21 can be estimated from compressional and shear wave polarization and slowness vectors corresponding to wide apertures of polar and azimuth angles. In isotropic media, when seismic source and receiver components have the same orientation (such as XX and YY), the reflection amplitude contours align approximately perpendicular to the particle motions. The mixed components (such as XY and YX) have amplitude patterns that are in symmetrical pairs of either the same, or of opposite, polarity on either side of the diagonal of the 9-C response matrix. In anisotropic media, amplitude variations with azimuth show the same basic patterns and symmetries as for isotropic, but with a superimposed tendency for alignment parallel to the strike of the vertical cracks. Solutions for elastic tensor elements from synthetic slowness and polarization data calculated directly from the Christoffel equation are more sensitive to the polar angle aperture than to the azimuth aperture. Nine-component synthetic elastic vertical seismic profile data for a model with triclinic symmetry calculated by finite-differencing allows estimation of the elastic 21 tensor elements in the vicinity of a three-component borehole receiver. Wide polar angle and azimuth apertures are needed for accurately estimating the elastic tensor elements. The tensor elements become less independent as the data apertures decrease. Results obtained by extracting slowness and polarization data from the corresponding synthetic seismograms show similar results. The inversion algorithm has produced good results from field vertical seismic profile data set from the Weyburn Field in Southern Saskatchewan in Canada. Synthetic nine-component seismograms calculated from the extracted tensor are able to explain most of the significant features in the field data. The inverted stiffness elastic tensor shows orthorhombic

  2. Geo-pressure prediction from 3-D seismic data: case studies from the gulf of Mexico

    Kan, T.K.; Kilsdonk, B. [ARCO Exploration and Production Technology, Plano TX (United States)

    1999-04-01

    We have developed a hypothesis and technique to interpret overpressure in both shales and sandstones by combining 3-D seismic geo-pressure estimates with 3-D interpretation of sandstone horizons. Before we can accept them, these ideas must be tested with measured pressure data from areas with 3-D geo-pressure cubes. Additionally, Our technique depends strongly on interpreting the lateral continuity of permeable units. We hope to improve our pre-drill prediction capability by incorporating both seismic-stratigraphic and impedance based facies interpretation techniques. (author)

  3. 3-D pneumatic seismic isolation of nuclear power plants

    This paper describes the work carried at the Russian Federation Research Center of Fundamental Engineering (RCFE), in development of innovative pneumatic multicomponent low-frequency seismic isolation bearings for advanced nuclear power plants.This device incorporates both supporting spherical elements, which provide displacements in the horizontal direction, and pneumatic dampers with rubber diaphragms for displacement in the vertical direction. To decrease the relative displacements of the isolated object the system uses viscoelastic dampers. Damping devices had been specially elaborated for the reactor building seismic isolation system as a result of substantial advances in the design and operation of the HD-type hydrodampers, created at the CKTI VIBROSEISM. The procedures developed have been used for comparison of the test and computer data on model isolated steel structure (MISS) and isolated rigid mass (IRM) isolators produced by ENEA and KAERI. Most recent work has concentrated on the development of mathematical models of isolators and isolated nuclear structures. Force-deformation characteristics of the HDRB model had been calculated on the basis of a special method of non-linear elastic theory using the continual transformations method. (author)

  4. Application of Catastrophe Theory in 3D Seismic Data Interpretation of Coal Mine

    ZHAO Mu-hua; YANG Wen-qiang; CUI Hui-xia

    2005-01-01

    In order to detect fault exactly and quickly, cusp catastrophe theory is used to interpret 3D coal seismic data in this paper. By establishing a cusp model, seismic signal is transformed into standard form of cusp catastrophe and catastrophe parameters, including time-domain catastrophe potential, time-domain catastrophe time, frequency-domain catastrophe potential and frequency- domain degree, are calculated. Catastrophe theory is used in 3D seismic structural interpretation in coal mine. The results show that the position of abnormality of the catastrophe parameter profile or curve is related to the location of fault, and the cusp catastrophe theory is effective to automatically pick up geology information and improve the interpretation precision in 3D seismic data.

  5. Analysis of Paleokarst Sinkholes in the Arkoma Basin using 3-D Seismic

    Kumbalek, Michael

    Paleokarst features are important to understand, both with regards to research geologists and to the petroleum industry. In terms of geology, understanding paleokarst features can yield more information about the depositional and surface environments of past times, and how diagenetic alteration affected the environment during the formation of karst features. In the petroleum industry, paleokarst features can have positive or negative consequence resulting in a potential reservoir with enhanced porosity due to the paleokarst features, or as a geo-hazard to prepare for or avoid when drilling. Inspired by issues faced when drilling in the Ft. Worth basin, this study utilizes multiple 3-D seismic surveys and subsurface well control to map paleokarsts within the Viola Limestone in the Arkoma Basin. Calculated seismic attribute volumes used to identify paleokarst sinkholes within the Viola Group include coherency and curvature attributes. ImageJ software was used to aid in counting and measuring paleokarst sinkholes identified using seismic mapping, coherency, and curvature attribute volumes. In addition to mapping, a cumulative distribution plot was produced from the diameters of the seismically mapped paleokarst sinkholes, allowing for an estimate to be made as to what the total amount of paleokarst sinkholes are within the study area. The methods detailed in this study proved to be effective in mapping and analyzing paleokarst sinkholes within the Viola Group. The paleokarst sinkholes mapped were determined to have been formed on the outer edge of the Southern Oklahoma aulacogen, as a result of the Sylvan/Viola unconformity. In addition to this, it has been determined that these paleokarst sinkholes are linked in formation to visually similar paleokarst sinkholes located in the Ellenburger Group in the Fort Worth Basin.

  6. Interpretation of a 3D Seismic-Reflection Volume in the Basin and Range, Hawthorne, Nevada

    Louie, J. N.; Kell, A. M.; Pullammanappallil, S.; Oldow, J. S.; Sabin, A.; Lazaro, M.

    2009-12-01

    A collaborative effort by the Great Basin Center for Geothermal Energy at the University of Nevada, Reno, and Optim Inc. of Reno has interpreted a 3d seismic data set recorded by the U.S. Navy Geothermal Programs Office (GPO) at the Hawthorne Army Depot, Nevada. The 3d survey incorporated about 20 NNW-striking lines covering an area of approximately 3 by 10 km. The survey covered an alluvial area below the eastern flank of the Wassuk Range. In the reflection volume the most prominent events are interpreted to be the base of Quaternary alluvium, the Quaternary Wassuk Range-front normal fault zone, and sequences of intercalated Tertiary volcanic flows and sediments. Such a data set is rare in the Basin and Range. Our interpretation reveals structural and stratigraphic details that form a basis for rapid development of the geothermal-energy resources underlying the Depot. We interpret a map of the time-elevation of the Wassuk Range fault and its associated splays and basin-ward step faults. The range-front fault is the deepest, and its isochron map provides essentially a map of "economic basement" under the prospect area. There are three faults that are the most readily picked through vertical sections. The fault reflections show an uncertainty in the time-depth that we can interpret for them of 50 to 200 ms, due to the over-migrated appearance of the processing contractor’s prestack time-migrated data set. Proper assessment of velocities for mitigating the migration artifacts through prestack depth migration is not possible from this data set alone, as the offsets are not long enough for sufficiently deep velocity tomography. The three faults we interpreted appear as gradients in potential-field maps. In addition, the southern boundary of a major Tertiary graben may be seen within the volume as the northward termination of the strong reflections from older Tertiary volcanics. Using a transparent volume view across the survey gives a view of the volcanics in full

  7. First images and orientation of fine structure from a 3-D seismic oceanography data set

    T. M. Blacic

    2010-04-01

    Full Text Available We present 3-D images of ocean fine structure from a unique industry-collected 3-D multichannel seismic dataset from the Gulf of Mexico that includes expendable bathythermograph casts for both swaths. 2-D processing reveals strong laterally continuous reflections throughout the upper ~800 m as well as a few weaker but still distinct reflections as deep as ~1100 m. We interpret the reflections to be caused by reversible fine structure from internal wave strains. Two bright reflections are traced across the 225-m-wide swath to produce reflection surface images that illustrate the 3-D nature of ocean fine structure. We show that the orientation of linear features in a reflection can be obtained by calculating the orientations of contours of reflection relief, or more robustly, by fitting a sinusoidal surface to the reflection. Preliminary 3-D processing further illustrates the potential of 3-D seismic data in interpreting images of oceanic features such as internal wave strains. This work demonstrates the viability of imaging oceanic fine structure in 3-D and shows that, beyond simply providing a way visualize oceanic fine structure, quantitative information such as the spatial orientation of features like fronts and solitons can be obtained from 3-D seismic images. We expect complete, optimized 3-D processing to improve both the signal to noise ratio and spatial resolution of our images resulting in increased options for analysis and interpretation.

  8. First images and orientation of internal waves from a 3-D seismic oceanography data set

    T. M. Blacic

    2009-10-01

    Full Text Available We present 3-D images of ocean finestructure from a unique industry-collected 3-D multichannel seismic dataset from the Gulf of Mexico that includes expendable bathythermograpgh casts for both swaths. 2-D processing reveals strong laterally continuous reflectors throughout the upper ~800 m as well as a few weaker but still distinct reflectors as deep as ~1100 m. Two bright reflections are traced across the 225-m-wide swath to produce reflector surface images that show the 3-D structure of internal waves. We show that the orientation of internal wave crests can be obtained by calculating the orientations of contours of reflector relief. Preliminary 3-D processing further illustrates the potential of 3-D seismic data in interpreting images of oceanic features such as internal wave strains. This work demonstrates the viability of imaging oceanic finestructure in 3-D and shows that, beyond simply providing a way to see what oceanic finestructure looks like, quantitative information such as the spatial orientation of features like internal waves and solitons can be obtained from 3-D seismic images. We expect complete, optimized 3-D processing to improve both the signal to noise ratio and spatial resolution of our images resulting in increased options for analysis and interpretation.

  9. Vertical Cable Seismic Survey for Hydrothermal Deposit

    Asakawa, E.; Murakami, F.; Sekino, Y.; Okamoto, T.; Ishikawa, K.; Tsukahara, H.; Shimura, T.

    2012-04-01

    The vertical cable seismic is one of the reflection seismic methods. It uses hydrophone arrays vertically moored from the seafloor to record acoustic waves generated by surface, deep-towed or ocean bottom sources. Analyzing the reflections from the sub-seabed, we could look into the subsurface structure. This type of survey is generally called VCS (Vertical Cable Seismic). Because VCS is an efficient high-resolution 3D seismic survey method for a spatially-bounded area, we proposed the method for the hydrothermal deposit survey tool development program that the Ministry of Education, Culture, Sports, Science and Technology (MEXT) started in 2009. We are now developing a VCS system, including not only data acquisition hardware but data processing and analysis technique. Our first experiment of VCS surveys has been carried out in Lake Biwa, JAPAN in November 2009 for a feasibility study. Prestack depth migration is applied to the 3D VCS data to obtain a high quality 3D depth volume. Based on the results from the feasibility study, we have developed two autonomous recording VCS systems. After we carried out a trial experiment in the actual ocean at a water depth of about 400m and we carried out the second VCS survey at Iheya Knoll with a deep-towed source. In this survey, we could establish the procedures for the deployment/recovery of the system and could examine the locations and the fluctuations of the vertical cables at a water depth of around 1000m. The acquired VCS data clearly shows the reflections from the sub-seafloor. Through the experiment, we could confirm that our VCS system works well even in the severe circumstances around the locations of seafloor hydrothermal deposits. We have, however, also confirmed that the uncertainty in the locations of the source and of the hydrophones could lower the quality of subsurface image. It is, therefore, strongly necessary to develop a total survey system that assures a accurate positioning and a deployment techniques

  10. Astor Pass Seismic Surveys Preliminary Report

    Louie, John [UNR; Pullammanappallil, Satish [Optim; Faulds, James; Eisses, Amy; Kell, Annie; Frary, Roxanna; Kent, Graham

    2011-08-05

    In collaboration with the Pyramid Lake Paiute Tribe (PLPT), the University of Nevada, Reno (UNR) and Optim re-processed, or collected and processed, over 24 miles of 2d seismic-reflection data near the northwest corner of Pyramid Lake, Nevada. The network of 2d land surveys achieved a near-3d density at the Astor Pass geothermal prospect that the PLPT drilled during Nov. 2010 to Feb. 2011. The Bureau of Indian Affairs funded additional seismic work around the Lake, and an extensive, detailed single-channel marine survey producing more than 300 miles of section, imaging more than 120 ft below the Lake bottom. Optim’s land data collection utilized multiple heavy vibrators and recorded over 200 channels live, providing a state-of-the-art reflection-refraction data set. After advanced seismic analysis including first-arrival velocity optimization and prestack depth migration, the 2d sections show clear fault-plane reflections, in some areas as deep as 4000 ft, tying to distinct terminations of the mostly volcanic stratigraphy. Some lines achieved velocity control to 3000 ft depth; all lines show reflections and terminations to 5000 ft depth. Three separate sets of normal faults appear in an initial interpretation of fault reflections and stratigraphic terminations, after loading the data into the OpendTect 3d seismic visualization system. Each preliminary fault set includes a continuous trace more than 3000 ft long, and a swarm of short fault strands. The three preliminary normal-fault sets strike northerly with westward dip, northwesterly with northeast dip, and easterly with north dip. An intersection of all three fault systems documented in the seismic sections at the end of Phase I helped to locate the APS-2 and APS-3 slimholes. The seismic sections do not show the faults connected to the Astor Pass tufa spire, suggesting that we have imaged mostly Tertiary-aged faults. We hypothesize that the Recent, active faults that produced the tufa through hotspring

  11. Integrated approach to 3-D seismic acquisition geometry analysis: emphasizing the influence of the inhomogeneous subsurface

    van Veldhuizen, E.J.

    2006-01-01

    The seismic reflection method for imaging of the earth's interior is an essential part of the exploration and exploitation of hydrocarbon resources. A seismic survey should be designed such that the acquired data leads to a sufficiently accurate subsurface image. The survey geometry analysis method

  12. 3D seismic survey in Honjo, Akita. Problems and struggles in acquisition and processing; Akitaken Honjo koku ni okeru sanjigen jishin tansa. Genba sagyo to data shori ni okeru mondaiten

    Imahori, S.; Kotera, Y.; Nakanishi, T. [Japan Energy Development Co. Ltd., Tokyo (Japan)

    1997-05-27

    Honjo mining area where investigations are conducted is hilly and has a complicated terrain with gas pipes buried in the ground just under the access road disabling the proper positioning of shock-generating large excavators or vibrators. Auger`s shallow hole shooting method is used in this survey to execute blastings at 639 points. In this method using charge depths of 4m, different from the conventional method using deeper charge depths (20-25m), surface waves prevail in the shot records giving rise to a new problem of removing them in the stage of data processing. The 2D filter that is a powerful tool in 2D data processing is not available in a 3D-survey where the tracing intervals are irregular in the shot records. In the effort of this time, a window length as a parameter in the direction of time is specified, and the F-X dip filtering method is employed in which any event that linearly continues beyond a certain number of traces in the said window is eliminated as a linear noise. It is recommended that the weighting function be changed in the direction of space since surface wave velocities are different at different locations. 1 fig., 1 tab.

  13. 3-D seismic and reservoir modeling, ram prospect, Viosca Knoll Block 912, offshore Gulf of Mexico

    Carew, W.; Ostendorf, P.F. (Exxon Company, New Orleans, LA (United States)); Krum, G.K. (Exxon Exploration Company, Houston, TX (United States))

    1993-09-01

    Ram prospect is a large stratigraphic trap located in Viosca Knoll Block 912/956/957, offshore Alabama in 3500-ft water depth. Reservoirs are Pliocene and Miocene gas- and oil-bearing deep-water sands deposited as fan complexes in an intraslope basin. The field has been proved by a total of 12 well penetrations and is nearing the development stage. In an effort to predict reservoir performance and recovery efficiencies, we constructed three-dimensional (3-D) reservoir models Exxon's in-house 3-D modeling program (GEOSET). Reservoir simulation studies will be based upon these 3-D geological models. We used 3-D seismic data to map seismic attributes around the prospect and well control to calibrate the seismic attributes based on known reservoir characteristics, thereby deriving a facies map for the entire field. Top/base structure, gross isopach, facies polygons, porosity, and Vshale were input into GEOSET to define the overall reservoir container and fill. The paucity of well data was compensated by using the 3-D-seismic based facies as a guide to filling polygons and by creating [open quotes]pseudowells[close quotes] from the real well data. These pseudowells aided in correlating within and between polygons. The resulting 3-D models (total porosity, effective porosity, Vsand) faithfully reflect the heterogeneity inferred from both 3-D seismic data and well control and provide visualization of reservoir continuity much better than models derived from well data alone. The models serve as a framework within which one can perform reservoir simulations and run various sensitivities. Additionally, the GEOSET porosity models can provide an alternative reservoir volume calculation.

  14. Development of a 400 Level 3C Clamped Downhole Seismic Receiver Array for 3D Borehole Seismic Imaging of Gas Reservoirs

    Bjorn N. P. Paulsson

    2006-09-30

    level receiver array can be used to obtain 3D 9C data. These 9C borehole seismic data provide both compressional wave and shear wave information that can be used for quantitative prediction of rock and pore fluid types. The 400-level borehole receiver array has been deployed successfully in a number of oil and gas wells during the course of this project, and each survey has resulted in marked improvements in imaging of geologic features that are critical for oil or gas production but were previously considered to be below the limits of seismic resolution. This added level of reservoir detail has resulted in improved well placement in the oil and gas fields that have been drilled using the Massive 3D VSP{reg_sign} images. In the future, the 400-level downhole seismic receiver array is expected to continue to improve reservoir characterization and drilling success in deep and complex oil and gas reservoirs.

  15. METHODOLOGY IMPLEMENTED FOR THE 3D-SEISMIC MODELLING USING GOCAD AND NORSAR 3D SOFTWARE APPLIED TO COMPLEX AREAS IN THE LLANOS FOOTHILLS

    Piedrahita Carlos

    2007-06-01

    Full Text Available A methodology has been applied in order to use seismic modeling and generate synthetic data which are compared with results from field data. This methodology has been applied to the Sirirí area and
    subsequently it has been extended to other geographic areas. To fully develop this methodology, the GOCAD and 3D NORSAR software packages were used to build respectively the geological and 3D seismic models. This methodology involves making structural modeling by edition of the different objects (horizons and faults in the structural module of GOCAD, with the help of available seismic and well data, in order to be consistent with reality. From this result a valid seismic model is imported to NORSAR 3D for making the seismic modeling, using ray tracing to get synthetic data. Finally, this information is compared with field data that has been processed.

  16. Enhanced imaging of CO2 at the Ketzin storage site: Inversion of 3D time-lapse seismic data

    Gil, M.; Götz, J.; Ivanova, A.; Juhlin, C.; Krawczyk, C. M.; Lüth, S.; Yang, C.

    2012-04-01

    The Ketzin test site, located near Berlin, is Europe's longest-operating on-shore CO2 storage site. As of December 2011, more than 56,000 tons of food grade CO2 has been injected since June 2008 in an anticlinal structure of the Northeast German Basin. The target reservoir consists of porous, brine bearing sandstone units of the Upper Triassic Stuttgart Formation at approximately 630 to 650 m depth. In order to enhance the understanding of the structural geometry of the site and to investigate the extension of the CO2-plume, several geophysical monitoring methods are being applied at Ketzin, among these are active seismic measurements, geoelectrics and borehole measurements. Among the various seismic techniques (e.g. 2D reflection surveys, crosshole tomography, Vertical Seismic Profiling, 2D- and 3D-Moving Source Profiling) employed at this pilot site, 3D time-lapse reflection surveys are an important component. The baseline 3D survey was acquired in 2005 and the first repeat measurements were performed in 2009 after injection of about 22,000 tons of CO2. The second repeat survey is planned to be carried out in fall 2012. These measurements allow the time-lapse signature of the injected CO2 to be imaged. The time-lapse amplitude variation attributed to the injected CO2 in the reservoir matches, considering detection limits of seismic surface measurements, the expected distribution of the CO2 plume derived from reservoir simulations. Previous attempts towards a quantitative interpretation were based on integrative considerations of different types of geophysical measurements using strict assumptions and characterized by large error bars. In order to increase the resolution and reliability of the data and to improve estimation of rock properties and especially to enhance the imaging resolution of the CO2-plume, the time-lapse 3D seismic data have now been inverted for seismic impedances with different methods, which is the focus of this presentation. One difficulty

  17. Error Analysis Of 3d Polygonal Model:A Survey

    Devendra Singh Rajput

    2012-05-01

    Full Text Available Various applications of computer graphics, (like animation, scientific visualization, and virtual reality involve the manipulation of geometric models. They are generally represented by triangular meshes due to its wide acceptance to process on rendering systems. The need of realism and high visual fidelity and the latest advances on scanning devices has increased complexity and size of triangular meshes. The original 3D model gets modified because of activities like approximation, transmission, processing and storage etc. Mostly the modification occurs due to simplification approaches which primarily use geometric distance metric as their simplification criteria. But it is hard to measure a small distance error accurately whereas other geometric or appearance error (like high curvature, thin region, color, texture, normals and volumetric has greater importance. Hence it is essential to understand the applicability of various parameters to evaluate the quality of 3D model. This paper briefly surveys the various errors analysis techniques, error metrics and tools to assess the quality of 3D mesh models.

  18. Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging

    Anderson, Roger N.; Boulanger, Albert; Bagdonas, Edward P.; Xu, Liqing; He, Wei

    1996-01-01

    The invention utilizes 3-D and 4-D seismic surveys as a means of deriving information useful in petroleum exploration and reservoir management. The methods use both single seismic surveys (3-D) and multiple seismic surveys separated in time (4-D) of a region of interest to determine large scale migration pathways within sedimentary basins, and fine scale drainage structure and oil-water-gas regions within individual petroleum producing reservoirs. Such structure is identified using pattern recognition tools which define the regions of interest. The 4-D seismic data sets may be used for data completion for large scale structure where time intervals between surveys do not allow for dynamic evolution. The 4-D seismic data sets also may be used to find variations over time of small scale structure within individual reservoirs which may be used to identify petroleum drainage pathways, oil-water-gas regions and, hence, attractive drilling targets. After spatial orientation, and amplitude and frequency matching of the multiple seismic data sets, High Amplitude Event (HAE) regions consistent with the presence of petroleum are identified using seismic attribute analysis. High Amplitude Regions are grown and interconnected to establish plumbing networks on the large scale and reservoir structure on the small scale. Small scale variations over time between seismic surveys within individual reservoirs are identified and used to identify drainage patterns and bypassed petroleum to be recovered. The location of such drainage patterns and bypassed petroleum may be used to site wells.

  19. Application of 3D Seismic Data Inversion to Coal Mining Prospecting

    2005-01-01

    Seismic inversion is one of the most important methods for lithological prospecting . Seismic data with low resolution is converted into impedance data of high resolution which can reflect the geological structure by inversion The inversion technique of 3D seismic data is discussed from both methodological and theoretical aspects, and the inversion test is also carried out using actual logging data. The result is identical with the measured data obtained from roadway of coal mine. The field tests and research results indicate that this method can provide more accurate data for identifying thin coal seam and minor faults.

  20. 3D seismic experiment in difficult area in Japan; Kokunai nanchiiki ni okeru sanjigen jishin tansa jikken

    Minegishi, M.; Nakagami, K.; Tanaka, H. [Japan National Oil Corp., Tokyo (Japan). Technology Research Center

    1997-05-27

    Difficult area in this context means an exploration-difficult area supposed to store oil/gas but retarded in exploration for the lack of knowledge about the geological structure due to poor quality of available seismic survey records. Discussed in this paper is a survey conducted into an area covering the southern part of Noshiro-shi, Akita-ken, and Yamamoto-cho, Yamamoto-gun, Akita-ken. An area size suitable for data collection at a target depth of 2500m is determined using an interpretation structure compiled on the basis of available well data and 2D seismic survey data. The plan for siting shock points and receiving points is modified case by case as restrictive factors come to the surface (resulting from the complicated hilly terrain, presence of pipes for agricultural water, etc.). The peculiarities of seismic waves in the terrain are studied through the interpretation of the available well data and 2D seismic survey data for the construction of a 3D velocity model for the confirmation of the appropriateness of the plan for siting shock points and receiving points. Efforts are exerted through enhanced coordination with the contractor to acquire data so that a technologically best design may be won within the limits of the budget. The quality of the data obtained from this experiment is in general better than those obtained from previous experiments, yet many problems remain to be settled in future studies about exploration-difficult areas. 4 refs., 4 figs., 1 tab.

  1. Automatic detection of karstic sinkholes in seismic 3D images using circular Hough transform

    Heydari Parchkoohi, Mostafa; Keshavarz Farajkhah, Nasser; Salimi Delshad, Meysam

    2015-10-01

    More than 30% of hydrocarbon reservoirs are reported in carbonates that mostly include evidence of fractures and karstification. Generally, the detection of karstic sinkholes prognosticate good quality hydrocarbon reservoirs where looser sediments fill the holes penetrating hard limestone and the overburden pressure on infill sediments is mostly tolerated by their sturdier surrounding structure. They are also useful for the detection of erosional surfaces in seismic stratigraphic studies and imply possible relative sea level fall at the time of establishment. Karstic sinkholes are identified straightforwardly by using seismic geometric attributes (e.g. coherency, curvature) in which lateral variations are much more emphasized with respect to the original 3D seismic image. Then, seismic interpreters rely on their visual skills and experience in detecting roughly round objects in seismic attribute maps. In this paper, we introduce an image processing workflow to enhance selective edges in seismic attribute volumes stemming from karstic sinkholes and finally locate them in a high quality 3D seismic image by using circular Hough transform. Afterwards, we present a case study from an on-shore oilfield in southwest Iran, in which the proposed algorithm is applied and karstic sinkholes are traced.

  2. 3D Geotechnical Soil Model of Nice, France, Inferred from Seismic Noise Measurements, for Seismic Hazard Assessment.

    Bertrand, E.; Duval, A.; Castan, M.; Vidal, S.

    2007-12-01

    In seismic risk studies, the assessment of lithologic site effect is based on an accurate knowledge of mechanical properties and geometry of superficial geological formations. Therefore, we built a 3D subsurface model in the city of Nice, southeastern France, using not only geological and geotechnical data but also geophysical inputs. We used especially ambient vibration recordings to supply the lack of borehole data over the city. Nice spreads over 72 km2 and roughly 20% of the city is built upon recent alluvium deposits. Other parts of the city lie on Jurassic and Cretaceous rocks to the east and thick Pliocene conglomerates to the west. Nearly 450 boreholes located mainly in the alluvial valleys were used. Because they are essentially linked to previous planned constructions (such as road network or important building), their distribution is rather heterogeneous over the studied area. In the valleys moreover, less than 40% of the boreholes are reaching the rock basement. These boreholes have been analyzed and a representative soil column made of 9 sedimentary layers has been recognized. Shear wave velocity of these layers were obtained from Standard Penetration Test values using several empirical correlation law described in the literature. Because of its cost, an extended boring survey was not feasible to complete our data set. Traditional seismic profiling was also not intended, as it is not possible to use intensive explosive sources in town. Recent years have seen many studies using ambient vibration measurements for site effect estimation. Especially, the very simple H/V technique was proven to be suitable for microzoning studies although some limitation were pointed out when dealing with 2D or 3D structures. Nevertheless, this technique alone provides only the fundamental eigenfrequency of the site under investigation. But assuming the shear wave velocity in the sediment it can helps to constrain the depth of the bedrock thanks to the well known f0=VS/4H

  3. Zonation of North Alex Mud Volcano Highlighted by 3-D Active and Passive Seismic Data

    Bialas, J.; Lefeldt, M. R.; Klaeschen, D.; Papenberg, C. A.; Brueckmann, W.

    2010-12-01

    The West Nile Delta forms part of the source of the large turbiditic Nile Deep Sea Fan. Since the late Miocene sediments have formed an up to 10 km thick pile, which includes about 1 - 3 km of Messinian evaporates. The sediment load of the overburden implies strong overpressures and salt-related tectonic deformation. Both are favourable for fluid migration towards the seafloor guided by the fractured margin. The western deltaic system, Rosetta branch, has formed an 80 km wide continental shelf. Here at 700 m water depth the mud volcano North Alex (NA) developed his circular bathymetric feature, which proved to be an active gas and mud-expelling structure. A 3-D high-resolution multichannel seismic survey (IFM-GEOMAR P-Cable system) was completed across the mud volcano. 3-D time migration provided a 3-D data cube with a 6.25 m grid. Vertical seismic sections did reveal a large set of faults located within the main mud volcano as well as surrounding the structure. Internal faults are mainly related to episodic mud expulsion processes and continuous gas and fluid production. Deep cutting external faults surround the structure in a half circle shape. Horizontal amplitude maps (time slices) of indicate recent activity of these faults even up to the seafloor. High gas saturation of the sediments is indicated by inverted reflection events. In the centre the gas front cuts into the seafloor reflection while it dips down with increasing radius. Only with the small grid resolution inward dipping reflections become visible, which form an upward opened concave reflector plane underlying the top gas front. The interpretation assumes an oval lens shaped body (conduit) saturated with gas at the top of the mud volcano. It provides the upper termination of the mud chimney. This separation is further supported by passive seismic observations. Distant earthquakes can stimulate long-period harmonic oscillations in mud volcanoes. Such oscillations are detectable with three

  4. ) A Feasibility Study for High Resolution 3D Seismic In The Deep Offshore Nigeria

    The conventional Exploration 3D seismic in the Deep Offshore Nigeria is typically acquired with 4000m-6000m cable length at 6-8 depth and with flip-flop shooting, providing a shot point interval of 50m. the average resulting frequency content is typically between 10-60hz which is adequate for exploration interpretation. It has become common in the last few years. E.g. in Angola and the Gulf of Mexico, to re-acquire High Resolution 3D seismic, after a discovery, to improve definition of turbidite systems and accuracy of reservoir geometry for optimized delineation drilling. This feasibility study which was carried out in three different steps was due to the question on whether HR-Seismic should be acquired over TotalFinaElf AKPO discovery for optimized delineation drilling

  5. Evaluation and assessment of the seismic responses of 3-D base-isolated structures

    Hueffmann, G.K.; Sutton, W.T. [GERB Vibration Control Systems, Inc., Westmont, IL (United States)

    1995-12-01

    The 1994 Northridge earthquake offered the opportunity to evaluate and assess the seismic response of several base-isolated buildings. Assessment of 2-D base isolation is straightforward comparing separately the uncoupled horizontal and vertical building motions to the corresponding ground motions. With 3-D systems characterized by low vertical natural frequencies, the system assessment must include rocking of the structure. Neglecting this consideration leads to an erroneous conclusion that the system greatly amplifies vertical ground accelerations. The paper evaluates the seismic response of a 3-D base-isolated building as compared to the same structure on a 2-D system. The paper also shows that the vertical accelerations measured at extreme locations in the building on a 3-D base-isolation system develop mainly from rocking responses to the horizontal ground motion.

  6. Integrating 3D seismic curvature and curvature gradient attributes for fracture characterization: Methodologies and interpretational implications

    Gao, Dengliang

    2013-03-01

    In 3D seismic interpretation, curvature is a popular attribute that depicts the geometry of seismic reflectors and has been widely used to detect faults in the subsurface; however, it provides only part of the solutions to subsurface structure analysis. This study extends the curvature algorithm to a new curvature gradient algorithm, and integrates both algorithms for fracture detection using a 3D seismic test data set over Teapot Dome (Wyoming). In fractured reservoirs at Teapot Dome known to be formed by tectonic folding and faulting, curvature helps define the crestal portion of the reservoirs that is associated with strong seismic amplitude and high oil productivity. In contrast, curvature gradient helps better define the regional northwest-trending and the cross-regional northeast-trending lineaments that are associated with weak seismic amplitude and low oil productivity. In concert with previous reports from image logs, cores, and outcrops, the current study based on an integrated seismic curvature and curvature gradient analysis suggests that curvature might help define areas of enhanced potential to form tensile fractures, whereas curvature gradient might help define zones of enhanced potential to develop shear fractures. In certain fractured reservoirs such as at Teapot Dome where faulting and fault-related folding contribute dominantly to the formation and evolution of fractures, curvature and curvature gradient attributes can be potentially applied to differentiate fracture mode, to predict fracture intensity and orientation, to detect fracture volume and connectivity, and to model fracture networks.

  7. Reservoir lithofacies analysis using 3D seismic data in dissimilarity space

    Seismic data interpretation is one of the most important steps in exploration seismology. Seismic facies analysis (SFA) with emphasis on lithofacies can be used to extract more information about structures and geology, which results in seismic interpretation enhancement. Facies analysis is based on unsupervised and supervised classification using seismic attributes. In this paper, supervised classification by a support vector machine using well logs and seismic attributes is applied. Dissimilarity as a new measuring space is employed, after which classification is carried out. Often, SFA is carried out in a feature space in which each dimension stands as a seismic attribute. Different facies show lots of class overlap in the feature space; hence, high classification error values are reported. Therefore, decreasing class overlap before classification is a necessary step to be targeted. To achieve this goal, a dissimilarity space is initially created. As a result of the definition of the new space, the class overlap between objects (seismic samples) is reduced and hence the classification can be done reliably. This strategy causes an increase in the accuracy of classification, and a more trustworthy lithofacies analysis is attained. For applying this method, 3D seismic data from an oil field in Iran were selected and the results obtained by a support vector classifier (SVC) in dissimilarity space are presented, discussed and compared with the SVC applied in conventional feature space. (paper)

  8. 3D Seismic Experimentation and Advanced Processing/Inversion Development for Investigations of the Shallow Subsurface

    Levander, Alan Richard [Rice Univ., Houston, TX (United States). Earth Science Department; Zelt, Colin A. [Rice Univ., Houston, TX (United States). Earth Science Department

    2015-03-17

    The work plan for this project was to develop and apply advanced seismic reflection and wide-angle processing and inversion techniques to high resolution seismic data for the shallow subsurface to seismically characterize the shallow subsurface at hazardous waste sites as an aid to containment and cleanup activities. We proposed to continue work on seismic data that we had already acquired under a previous DoE grant, as well as to acquire additional new datasets for analysis. The project successfully developed and/or implemented the use of 3D reflection seismology algorithms, waveform tomography and finite-frequency tomography using compressional and shear waves for high resolution characterization of the shallow subsurface at two waste sites. These two sites have markedly different near-surface structures, groundwater flow patterns, and hazardous waste problems. This is documented in the list of refereed documents, conference proceedings, and Rice graduate theses, listed below.

  9. Etruscans in 3D - Surveying and 3D modeling for a better access and understanding of heritage -

    B. Jiménez Fernández-Palacios; Rizzi., A; F. Remondino

    2013-01-01

    Archaeological 3D digital documentation of monuments and historical sites should be considered a precious source of information and it can be very useful for preservation, conservation, restoration and reconstruction of Cultural Heritage. This paper reports a work dealing with 3D surveying and modeling of different Etruscan heritage sites, featuring necropolis with underground frescoed tombs dating back to VII-IV century B.C., located in the area corresponding roughly to the actual central It...

  10. Seismic Surveys in Ice-Covered Waters

    Klavenes, Jørgen Sørensen

    2013-01-01

    The master’s thesis ‘Seismic surveys in ice-covered waters’ describes several methods within the maritime seismic industry to acquire data concerning the Earth’s strata. It is the objective of this thesis to present a viable solution enabling the seismic industry to perform surveys in ice influenced waters. The hardware conducted for conventional seismic surveys are introduced and explained. Arctic area and ice issues in regards to marine seismic surveys are investigated. A suitable theor...

  11. Porosity and permeability prediction of Zechstein-2-carbonates from 3D seismic data

    Trappe, H.; Krajewski, P.; Aust, S. [Deilmann Erdoel Erdgas GmbH, Lingen (Germany)

    1995-08-01

    In this study the applicability is explored of using 3D seismic data to estimate porosity and permeability for Zechstein-2-Carbonates using advanced techniques such as Neural Networks and image Processing algorithm. The study area is located in the North German Basin next to the Dutch border. About 40 wells tested the Zechstein reservoir at depth between 3000-4000m. Production of the gas fields started in the fifties/early sixties. The wells were generally drilled on structural highs mapped by 2D seismic data. In a first step a detailed seismic modelling exercise was carried in order to find a relationship between reservoir quality and seismic attributes. Moreover the influence of overburden and hanging wall layers to the seismic signature at reservoir level was investigated. In the second phase the previously obtained results were proofed by the 3D seismic. Using amplitude and acoustic impedance porosity maps of the reservoir were calculated. While the porosity maps show a high degree of confidence unfortunately the porosity-permeability relationship within the carbonates is poor. Seismic attributes and geological parameters form the input for a selforganizing Neural Network. This Neural Network was trained to characterise facies units of variable permeabilities. After successful training the network is able to predict facies and permeability as well as porosity for a given area. The last stage of the study covers the structural framework of the study area. Seismic attributes such as two way traveltime and amplitude were processed using image processing algorithms. Edge detection operators, IHS ilumination, gradient descriptors and other procedures were applied to highlight small scale faulting which is thought to be responsible for an increase the production performance.

  12. 3D Porosity Estimation of the Nankai Trough Sediments from Core-log-seismic Integration

    Park, J. O.

    2015-12-01

    The Nankai Trough off southwest Japan is one of the best subduction-zone to study megathrust earthquake fault. Historic, great megathrust earthquakes with a recurrence interval of 100-200 yr have generated strong motion and large tsunamis along the Nankai Trough subduction zone. At the Nankai Trough margin, the Philippine Sea Plate (PSP) is being subducted beneath the Eurasian Plate to the northwest at a convergence rate ~4 cm/yr. The Shikoku Basin, the northern part of the PSP, is estimated to have opened between 25 and 15 Ma by backarc spreading of the Izu-Bonin arc. The >100-km-wide Nankai accretionary wedge, which has developed landward of the trench since the Miocene, mainly consists of offscraped and underplated materials from the trough-fill turbidites and the Shikoku Basin hemipelagic sediments. Particularly, physical properties of the incoming hemipelagic sediments may be critical for seismogenic behavior of the megathrust fault. We have carried out core-log-seismic integration (CLSI) to estimate 3D acoustic impedance and porosity for the incoming sediments in the Nankai Trough. For the CLSI, we used 3D seismic reflection data, P-wave velocity and density data obtained during IODP (Integrated Ocean Drilling Program) Expeditions 322 and 333. We computed acoustic impedance depth profiles for the IODP drilling sites from P-wave velocity and density data. We constructed seismic convolution models with the acoustic impedance profiles and a source wavelet which is extracted from the seismic data, adjusting the seismic models to observed seismic traces with inversion method. As a result, we obtained 3D acoustic impedance volume and then converted it to 3D porosity volume. In general, the 3D porosities show decrease with depth. We found a porosity anomaly zone with alteration of high and low porosities seaward of the trough axis. In this talk, we will show detailed 3D porosity of the incoming sediments, and present implications of the porosity anomaly zone for the

  13. Reconstruction of allochthonous salt emplacement from 3-D seismic reflection data, northern Gulf of Mexico

    Walters, R.D. (Univ. of Texas, Austin (United States))

    1993-05-01

    Analysis of a 110 mi[sup 2] (284 km[sup 2]) three-dimensional (3-D) seismic survey located at the northern Gulf of Mexico continental shelf edge documents the emplacement history of an allochthonous salt structure. A listric, counter-regional growth fault bounds the updip salt edge and apparently soles into a detachment surface from which salt was evacuated. The counter-regional fault system served as the feeder stock for the dome. The current salt shape is an asymmetric, elongated wedge with a diameter of 20,000 ft (6.1 km) and maximum thickness of 13,000 ft (4.0 km). A steep fault/salt contact bounded by small, oval- shaped withdrawal subbasins on the landward side contrasts with a shallow-dipping basinward salt/sediment contact. A relatively thin sediment cover onlaps the shallow basinward salt flank. Regional seismic lines support the presence of a detachment surface in the area apparently separating Upper Cretaceous to lower Tertiary( ) sediments from upper Miocene (7.0--5.4 Ma) and younger sediments. During the first stage, an allochthonous salt sheet was emplaced during the hiatus and then segmented. From the middle to late Pliocene (4.0--2.2 Ma), the salt structure evolved by downbuilding, remaining close to the ocean floor surface and covered by only a thin sediment veneer. After the beginning of the Illinoian glacial stage (0.85 Ma), the deep source salt layer was depleted and the dome was buried by up to 1,700 ft (520 m) of sediment during the third stage. Seismic amplitude anomalies associated with structural disconformities indicate a possible salt evacuation surface or salt weld along the deep feeder stock. Geohistory diagrams indicate extremely rapid local subsidence in the withdrawal basin, which contains up to 12,000 ft (3.6 km) of Pleistocene sediment. Such acceleration subsidence, in conjunction with the fault/salt geometry, contributed to the development of hydrocarbon reservoirs near the dome flanks. 69 refs., 17 figs., 1 tab.

  14. Optical 3-D-measurement techniques : a survey

    Tiziani, Hans J.

    1989-01-01

    Close range photogrammetry will be more frequently applied in industry for 3-D-sensing when real time processing can be applied. Computer vision, machine vision, robot vision are in fact synonymous with real time photogrammetry. This overview paper concentrates on optical methods for 3-D-measurements. Incoherent and coherent methods for 3-D-sensing will be presented. Particular emphasis is put on high precision 3-D-measurements. Some of the work of our laboratory will be reported.

  15. Simulation of CO2 distribution pattern in an underground CO2 injection projected calibrated by 3D seismics

    In the ongoing aquifer CO2 disposal project in the Sleipner license (North Sea), underground CO2 is being monitored by time-lapse seismic. The CO2 is being injected close to the base of a high permeable, highly porous sand unit, the Utsira Sand. In an iterative process between seismic surveys and reservoir simulations, a reservoir model featuring the major controlling heterogeneities has been developed. Well-data and seismic data prior to injection shows that the sand is divided by nearly horizontal, discontinuous shales. >From the 3-D seismic image after three years of injection, strong reflectors can be interpreted as CO2 accumulations identifying the major shale layers that control the vertical migration of CQ from the injection point to the top of the formation. By modelling this flow in reservoir simulations, it can be inferred that the CO2 is transported in distinct columns between the shales rather than as dispersed bubbles over a large area. Improvement of the geological model increases the confidence of predictions based on simulation of the long-time fate Of CO2. A possible natural aquifer flow can have a pronounced effect on the location Of CO2 accumulations due to the relatively flat topography of the trapping shales. This effect has been quantified by simulation and this phenomenon was used to adjust the localisation of the CO2 bubbles to better fir the seismic images

  16. Vertical Cable Seismic Survey for SMS exploration

    Asakawa, Eiichi; Murakami, Fumitoshi; Tsukahara, Hotoshi; Mizohata, Shigeharu

    2014-05-01

    The Vertical Cable Seismic (VCS) survey is one of the reflection seismic methods. It uses hydrophone arrays vertically moored from the seafloor to record acoustic waves generated by sea-surface, deep-towed or ocean bottom sources. Analyzing the reflections from the sub-seabed, we could look into the subsurface structure. Because the VCS is an efficient high-resolution 3D seismic survey method for a spatially-bounded area, we proposed it for the SMS survey tool development program that the Ministry of Education, Culture, Sports, Science and Technology (MEXT) started in 2009. We have been developing the VCS survey system, including not only data acquisition hardware but data processing and analysis technique. We carried out several VCS surveys combining with surface towed source, deep towed source and ocean bottom source. The water depths of these surveys are from 100m up to 2100 m. Through these experiments, our VCS data acquisition system has been also completed. But the data processing techniques are still on the way. One of the most critical issues is the positioning in the water. The uncertainty in the positions of the source and of the hydrophones in water degraded the quality of subsurface image. GPS navigation system is available on sea surface, but in case of deep-towed source or ocean bottom source, the accuracy of shot position with SSBL/USBL is not sufficient for the very high-resolution imaging. We have developed a new approach to determine the positions in water using the travel time data from the source to VCS hydrophones. In 2013, we have carried out the second VCS survey using the surface-towed high-voltage sparker and ocean bottom source in the Izena Cauldron, which is one of the most promising SMS areas around Japan. The positions of ocean bottom source estimated by this method are consistent with the VCS field records. The VCS data with the sparker have been processed with 3D PSTM. It gives the very high resolution 3D volume deeper than two

  17. INCREASING OIL RECOVERY THROUGH ADVANCED REPROCESSING OF 3D SEISMIC, GRANT CANYON AND BACON FLAT FIELDS, NYE COUNTY, NEVADA

    Makoil, Inc., of Orange, California, with the support of the U.S. Department of Energy has reprocessed and reinterpreted the 3D seismic survey of the Grant Canyon area, Railroad Valley, Nye County, Nevada. The project was supported by Dept. of Energy Grant DE-FG26-00BC15257. The Grant Canyon survey covers an area of 11 square miles, and includes Grant Canyon and Bacon Flat oil fields. These fields have produced over 20 million barrels of oil since 1981, from debris slides of Devonian rocks that are beneath 3,500 to 5,000 ft of Tertiary syntectonic deposits that fill the basin of Railroad Valley. High-angle and low-angle normal faults complicate the trap geometry of the fields, and there is great variability in the acoustic characteristics of the overlying valley fill. These factors combine to create an area that is challenging to interpret from seismic reflection data. A 3D seismic survey acquired in 1992-93 by the operator of the fields has been used to identify development and wildcat locations with mixed success. Makoil believed that improved techniques of processing seismic data and additional well control could enhance the interpretation enough to improve the chances of success in the survey area. The project involved the acquisition of hardware and software for survey interpretation, survey reprocessing, and reinterpretation of the survey. SeisX, published by Paradigm Geophysical Ltd., was chosen as the interpretation software, and it was installed on a Dell Precision 610 computer work station with the Windows NT operating system. The hardware and software were selected based on cost, possible addition of compatible modeling software in the future, and the experience of consulting geophysicists in the Billings area. Installation of the software and integration of the hardware into the local office network was difficult at times but was accomplished with some technical support from Paradigm and Hewlett Packard, manufacturer of some of the network equipment. A

  18. INCREASING OIL RECOVERY THROUGH ADVANCED REPROCESSING OF 3D SEISMIC, GRANT CANYON AND BACON FLAT FIELDS, NYE COUNTY, NEVADA

    Eric H. Johnson; Don E. French

    2001-06-01

    Makoil, Inc., of Orange, California, with the support of the U.S. Department of Energy has reprocessed and reinterpreted the 3D seismic survey of the Grant Canyon area, Railroad Valley, Nye County, Nevada. The project was supported by Dept. of Energy Grant DE-FG26-00BC15257. The Grant Canyon survey covers an area of 11 square miles, and includes Grant Canyon and Bacon Flat oil fields. These fields have produced over 20 million barrels of oil since 1981, from debris slides of Devonian rocks that are beneath 3,500 to 5,000 ft of Tertiary syntectonic deposits that fill the basin of Railroad Valley. High-angle and low-angle normal faults complicate the trap geometry of the fields, and there is great variability in the acoustic characteristics of the overlying valley fill. These factors combine to create an area that is challenging to interpret from seismic reflection data. A 3D seismic survey acquired in 1992-93 by the operator of the fields has been used to identify development and wildcat locations with mixed success. Makoil believed that improved techniques of processing seismic data and additional well control could enhance the interpretation enough to improve the chances of success in the survey area. The project involved the acquisition of hardware and software for survey interpretation, survey reprocessing, and reinterpretation of the survey. SeisX, published by Paradigm Geophysical Ltd., was chosen as the interpretation software, and it was installed on a Dell Precision 610 computer work station with the Windows NT operating system. The hardware and software were selected based on cost, possible addition of compatible modeling software in the future, and the experience of consulting geophysicists in the Billings area. Installation of the software and integration of the hardware into the local office network was difficult at times but was accomplished with some technical support from Paradigm and Hewlett Packard, manufacturer of some of the network equipment. A

  19. Construction method and application of 3D velocity model for evaluation of strong seismic motion and its cost performance

    Based on experiences of making subsurface structure models for seismic strong motion evaluation, the advantages and disadvantages in terms of convenience and cost for several methods used to make such models were reported. As for the details, gravity and micro-tremor surveys were considered to be highly valid in terms of convenience and cost. However, stratigraphy and seismic velocity structure are required to make accurate 3-D subsurface structures. To realize these, methods for directly examining subsurface ground or using controlled tremor sources (at high cost) are needed. As a result, it was summarized that in modeling subsurface structures, some sort of plan including both types of methods is desirable and that several methods must be combined to match one's intended purposes and budget. (authors)

  20. Complex Crustal Structure Beneath Western Turkey Revealed by 3D Seismic Full Waveform Inversion (FWI)

    Cubuk-Sabuncu, Yesim; Taymaz, Tuncay; Fichtner, Andreas

    2016-04-01

    We present a 3D radially anisotropic velocity model of the crust and uppermost mantle structure beneath the Sea of Marmara and surroundings based on the full waveform inversion method. The intense seismic activity and crustal deformation are observed in the Northwest Turkey due to transition tectonics between the strike-slip North Anatolian Fault (NAF) and the extensional Aegean region. We have selected and simulated complete waveforms of 62 earthquakes (Mw > 4.0) occurred during 2007-2015, and recorded at (Δ Unified Seismic Network (HUSN, Greece) and Earthquake Research Center of Turkey (AFAD-DAD). The spectral-element solver of the wave equation, SES3D algorithm, is used to simulate seismic wave propagation in 3D spherical coordinates (Fichtner, 2009). The Large Scale Seismic Inversion Framework (LASIF) workflow tool is also used to perform full seismic waveform inversion (Krischer et al., 2015). The initial 3D Earth model is implemented from the multi-scale seismic tomography study of Fichtner et al. (2013). Discrepancies between the observed and simulated synthetic waveforms are determined using the time-frequency misfits which allows a separation between phase and amplitude information (Fichtner et al., 2008). The conjugate gradient optimization method is used to iteratively update the initial Earth model when minimizing the misfit. The inversion is terminated after 19 iterations since no further advances are observed in updated models. Our analysis revealed shear wave velocity variations of the shallow and deeper crustal structure beneath western Turkey down to depths of ~35-40 km. Low shear wave velocity anomalies are observed in the upper and mid crustal depths beneath major fault zones located in the study region. Low velocity zones also tend to mark the outline of young volcanic areas. Our final 3D Earth model is tested using forward wave simulations of earthquakes (M ≥ 3.7) that were not used during the inversion process. The comparison of observed

  1. Imaging of 3-D seismic velocity structure of Southern Sumatra region using double difference tomographic method

    Southern Sumatra region has a high level of seismicity due to the influence of the subduction system, Sumatra fault, Mentawai fault and stretching zone activities. The seismic activities of Southern Sumatra region are recorded by Meteorological Climatological and Geophysical Agency (MCGA’s) Seismograph network. In this study, we used earthquake data catalog compiled by MCGA for 3013 events from 10 seismic stations around Southern Sumatra region for time periods of April 2009 – April 2014 in order to invert for the 3-D seismic velocities structure (Vp, Vs, and Vp/Vs ratio). We applied double-difference seismic tomography method (tomoDD) to determine Vp, Vs and Vp/Vs ratio with hypocenter adjustment. For the inversion procedure, we started from the initial 1-D seismic velocity model of AK135 and constant Vp/Vs of 1.73. The synthetic travel time from source to receiver was calculated using ray pseudo-bending technique, while the main tomographic inversion was applied using LSQR method. The resolution model was evaluated using checkerboard test and Derivative Weigh Sum (DWS). Our preliminary results show low Vp and Vs anomalies region along Bukit Barisan which is may be associated with weak zone of Sumatran fault and migration of partial melted material. Low velocity anomalies at 30-50 km depth in the fore arc region may indicated the hydrous material circulation because the slab dehydration. We detected low seismic seismicity in the fore arc region that may be indicated as seismic gap. It is coincides contact zone of high and low velocity anomalies. And two large earthquakes (Jambi and Mentawai) also occurred at the contact of contrast velocity

  2. Imaging of 3-D seismic velocity structure of Southern Sumatra region using double difference tomographic method

    Lestari, Titik, E-mail: t2klestari@gmail.com [Meteorological Climatological and Geophysical Agency (MCGA), Jalan Angkasa I No.2 Kemayoran, Jakarta Pusat, 10720 (Indonesia); Faculty of Earth Science and Technology, Bandung Institute of Technology, Jalan Ganesa No.10, Bandung 40132 (Indonesia); Nugraha, Andri Dian, E-mail: nugraha@gf.itb.ac.id [Global Geophysical Research Group, Faculty of Mining and Petroleum Engineering, Bandung Institute of Technology, Jalan Ganesa 10 Bandung 40132 (Indonesia)

    2015-04-24

    Southern Sumatra region has a high level of seismicity due to the influence of the subduction system, Sumatra fault, Mentawai fault and stretching zone activities. The seismic activities of Southern Sumatra region are recorded by Meteorological Climatological and Geophysical Agency (MCGA’s) Seismograph network. In this study, we used earthquake data catalog compiled by MCGA for 3013 events from 10 seismic stations around Southern Sumatra region for time periods of April 2009 – April 2014 in order to invert for the 3-D seismic velocities structure (Vp, Vs, and Vp/Vs ratio). We applied double-difference seismic tomography method (tomoDD) to determine Vp, Vs and Vp/Vs ratio with hypocenter adjustment. For the inversion procedure, we started from the initial 1-D seismic velocity model of AK135 and constant Vp/Vs of 1.73. The synthetic travel time from source to receiver was calculated using ray pseudo-bending technique, while the main tomographic inversion was applied using LSQR method. The resolution model was evaluated using checkerboard test and Derivative Weigh Sum (DWS). Our preliminary results show low Vp and Vs anomalies region along Bukit Barisan which is may be associated with weak zone of Sumatran fault and migration of partial melted material. Low velocity anomalies at 30-50 km depth in the fore arc region may indicated the hydrous material circulation because the slab dehydration. We detected low seismic seismicity in the fore arc region that may be indicated as seismic gap. It is coincides contact zone of high and low velocity anomalies. And two large earthquakes (Jambi and Mentawai) also occurred at the contact of contrast velocity.

  3. Effect of Damping and Yielding on the Seismic Response of 3D Steel Buildings with PMRF

    Alfredo Reyes-Salazar; Achintya Haldar; Ramon Eduardo Rodelo-López; Eden Bojórquez

    2014-01-01

    The effect of viscous damping and yielding, on the reduction of the seismic responses of steel buildings modeled as three-dimensional (3D) complex multidegree of freedom (MDOF) systems, is studied. The reduction produced by damping may be larger or smaller than that of yielding. This reduction can significantly vary from one structural idealization to another and is smaller for global than for local response parameters, which in turn depends on the particular local response parameter. The unc...

  4. Recovering Uniform Coverage in a 3D Survey: Case Study from Onshore Southern India

    K. Shukla

    2014-01-01

    Full Text Available We present a case study on 3D seismic acquisition under inaccessible ground conditions and heavy head-loads. The novelty in this case study is designing novel shot recovery strategies. The study area is in a densely populated region in Southern India, where a 3D survey was conducted over an area of 700 km2 to illuminate a target formation at 3000 m depth. The survey was designed in a brick pattern with six shots in a round, a line repeat interval of 2, and offset and skid grid of 300 m × 300 m. Only 40% of full fold could be obtained with conventional shot recovery methods. To increase the survey fold, two new shot recovery strategies were attempted. First, the original offset and skid grid was increased to 1100 m × 1200 m. Second, if the recovery shot did not fall in this grid, it was relocated along the swath from its original location to a distance equal to half the inline offset. The two strategies, employed together, increased the fold to 90% of full fold at the target area maintaining its uniformity. Although the illumination of the target zone was adequate for interpretation in this case, we see a need for adaptation if these strategies are to be applied in other surveys. In general, issues related to offset-limit, shot-density, and grid-spacing are a matter of continual optimization in 3D surveys.

  5. 3D Seismic Characterization of the Research Facility for Geological Storage of CO2: Hontomín (Burgos, Spain)

    Alcalde, J.; Martí, D.; Calahorrano, A.; Marzan, I.; Ayarza, P.; Carbonell, R.; Perez-Estaun, A.

    2011-12-01

    A technological research facility dedicated to the underground geological storage of CO2 is currently being developed by the Spanish research program on Carbon Capture and Storage (CCS) in Hontomin (Burgos, North of the Iberian Peninsula) This research program is being developed by the CIUDEN Foundation. CIUDEN is an initiative of 3 Spanish state departments (Science & Innovation, Environment and Industry). An extensive multidisciplinary site characterization phase has been carried out, including a multi-seismic data acquisition experiment. Within this effort a 36 km2 academic-oriented 3D seismic reflection survey was acquired in summer 2010. The aim of data acquisition effort are to provide high resolution images of the subsurface of the storage complex, constrain a baseline model for all the disciplines involved in the project. The main acquisition characteristics of this survey included: a mixed source (Vibroseis & explosive , 74% and 26% of the source points, respectively); 5000 shot points, distributed along 22 source lines (separated 250 m), 22 lines of receivers (separated 275 m); shot and receiver spacing along the source and receiver lines was of 25 m; this resulted in a nominal CDP-fold of 36 traces, with 13 m2 bins. This 3D-data was fully processed until migration. The main features within the processing sequence include static correction calculation, frequency filtering, trace amplitude equalization, rms velocity modeling, FK-domain filtering, 3D deconvolution, dip move-out corrections, residual static calculation and pre and post stack migration. The final high-resolution 3D-volume allowed to characterize the main tectonic structure of the dome complex, the fault system of the area and the feasibility of the reservoir for the storage. The target reservoir is a saline aquifer placed at 1400, approximately, within Lower Jurassic carbonates (Lias); the main seal is formed by inter-layered marls and marly limestones from Early to Middle Jurassic (Dogger

  6. Detection of coalbed fractures with P-wave azimuthal AVO in 3-D seismic exploration

    LI Guofa; PENG Suping; HE Bingshou; PENG Xiaobo; YUAN Chunfang; HU Chaoyuan

    2005-01-01

    The detection of fractures is important for production and safety in coal fields. Subsurface fractures result in azimuthal anisotropy of the seismic wave, and the amplitude of reflection wave varies with offset and azimuth.In case of weak anisotropy, the reflection coefficients of P-wave are concisely denoted as the analytic function of fracture parameters. For the purpose of predicting the coalbed fracture distribution through analyzing variation of the reflection amplitudes with offset and azimuth, 3-D seismic data with full-azimuth were acquired in a coal field in Huainan, Anhui Province. The careful analysis and process of seismic data showed that the reflection amplitude of the primary coaibed varied with azimuth in much consistent with the theoretical model. The conclusion was drawn that the coal-bed fracture in this coal field could be predicted through the method of the P-wave azimuthal AVO.

  7. Subtle traps prediction using sequence stratigraphy and 3D seismic technology: A case study from Qikou depression in Huanghua basin

    MAO Ning-bo; DAI Ta-gen; PENG Sheng-lin

    2005-01-01

    Forecasting subtle traps by sequence stratigraphy and 3D seismic data is a sensitive topic in hydrocarbon exploration. Research on subtle traps by geophysical data is the most popular and difficult. Based on the sufficiently drilling data, log data, core data and 3D seismic data, sediment sequence of Qikou depression, Huanghua basin was partitioned by using sequence stratigraphy theory. Each sediment sequence system mode was built. Sediment faces of subtle traps were pointed out. Dominating factors forming subtle traps were analyzed. Sandstone seismic rock physics and its response were studied in Tertiary System. Sandstone geophysical response and elastic modulus vary laws with pressure, temperature, porosity, depth were built. Experimental result and practice shows that it is possible using seismic information forecasting subtle traps. Integrated using geology, log, drilling data, special seismic processing technique, interpretation technique, high precision horizon calibration technique, 3D seismic visualizing interpretation, seismic coherence analysis, attribute analysis, logging-constrained inversion, time frequency analysis, subtle trapsobject is identified and interpreted. Finally, advantage object of subtle trap in this area was determined. Bottomland sand stratigraphic and lithologic reservoirs in Qinan slope zone have been founded by means of high resolution 3D seismic data field technique, high resolution 3D seismic data processing technique and seismic wave impendence inversion technique.

  8. Multi-dimensional Seismic Response Analysis of Base-Isolated Frame Structure with 3D Isolator

    Xiong Shishu; Huang Liting; Chen Jinfeng; Su Jingsu

    2005-01-01

    The three-dimensional lead-rubber dish-spring bearing (3DB) is proposed in this paper. The 3DB is composed of lead rubber bearing (LRB) and dish-spring bearing (DSB) with damper in series. The 3DB put forward in this paper is effective in the resolution of difficulties in strong vertical capacity and vertical damping of three-dimensional isolation bearings. It effectively suppresses rocking motions as well. The analytical model and motion equations of multi-dimensional seismic responses of 3D base-isolated frame structures are established. Taking a five-storey frame structure as an example, an extensive simulation analysis is carried out. The results show that the 3D base-isolated structure with the proposed 3DB is effective in 3D isolation; it can reduce seismic responses by 50 % compared to a non-isolated structure. Therefore, the 3D isolation problem in building can be solved easily and effectively with the 3DB proposed in this paper.

  9. Digital 3D Borobudur - Integration of 3D surveying and modeling techniques

    Suwardhi, D.; Menna, F.; Remondino, F.; Hanke, K.; Akmalia, R.

    2015-08-01

    The Borobudur temple (Indonesia) is one of the greatest Buddhist monuments in the world, now listed as an UNESCO World Heritage Site. The present state of the temple is the result of restorations after being exposed to natural disasters several times. Today there is still a growing rate of deterioration of the building stones whose causes need further researches. Monitoring programs, supported at institutional level, have been effectively executed to observe the problem. The paper presents the latest efforts to digitally document the Borobudur Temple and its surrounding area in 3D with photogrammetric techniques. UAV and terrestrial images were acquired to completely digitize the temple, produce DEM, orthoimages and maps at 1:100 and 1:1000 scale. The results of the project are now employed by the local government organizations to manage the heritage area and plan new policies for the conservation and preservation of the UNESCO site. In order to help data management and policy makers, a web-based information system of the heritage area was also built to visualize and easily access all the data and achieved 3D results.

  10. Seismic surveying and accelerators

    The paper deals with an investigation into the impact of earth vibrations on charged particle beams in modern colliders. It is ascertained that the displacement of accelerator magnetic elements from the perfect position results in the excitation of betatron oscillations and distortion of particle orbit position. The results of experimental investigations into seismic noises are presented for ASR, SSC, DESY and KEK. The rms orbit displacement in accelerators is estimated relying on the law of earth diffusion motion, according to which the variance of relative displacements is proportional to the distance between these points and time of observation. 6 refs., 3 figs., 2 tabs

  11. Training toward Advanced 3D Seismic Methods for CO2 Monitoring, Verification, and Accounting

    Christopher Liner

    2012-05-31

    The objective of our work is graduate and undergraduate student training related to improved 3D seismic technology that addresses key challenges related to monitoring movement and containment of CO{sub 2}, specifically better quantification and sensitivity for mapping of caprock integrity, fractures, and other potential leakage pathways. We utilize data and results developed through previous DOE-funded CO{sub 2} characterization project (DE-FG26-06NT42734) at the Dickman Field of Ness County, KS. Dickman is a type locality for the geology that will be encountered for CO{sub 2} sequestration projects from northern Oklahoma across the U.S. midcontinent to Indiana and Illinois. Since its discovery in 1962, the Dickman Field has produced about 1.7 million barrels of oil from porous Mississippian carbonates with a small structural closure at about 4400 ft drilling depth. Project data includes 3.3 square miles of 3D seismic data, 142 wells, with log, some core, and oil/water production data available. Only two wells penetrate the deep saline aquifer. In a previous DOE-funded project, geological and seismic data were integrated to create a geological property model and a flow simulation grid. We believe that sequestration of CO{sub 2} will largely occur in areas of relatively flat geology and simple near surface, similar to Dickman. The challenge is not complex geology, but development of improved, lower-cost methods for detecting natural fractures and subtle faults. Our project used numerical simulation to test methods of gathering multicomponent, full azimuth data ideal for this purpose. Our specific objectives were to apply advanced seismic methods to aide in quantifying reservoir properties and lateral continuity of CO{sub 2} sequestration targets. The purpose of the current project is graduate and undergraduate student training related to improved 3D seismic technology that addresses key challenges related to monitoring movement and containment of CO{sub 2

  12. New insights into the earliest Quaternary environments in the Central North Sea from 3D seismic

    Lamb, Rachel; Huuse, Mads; Stewart, Margaret; Brocklehurst, Simon H.

    2014-05-01

    In the past the transition between an unconformable surface in the south to a conformable horizon towards the north has made identification and mapping the base-Quaternary in the central North Sea difficult (Sejrup et al 1991; Gatliff et al 1994). However recent integration of biostratigraphy, pollen analysis, paleomagnetism and amino acid analysis in the Dutch and Danish sectors (Rasmussen et al 2005; Kuhlmann et al 2006) has allowed greater confidence in the correlation to the region 3D seismic datasets and thus has allowed the base-Quaternary to be mapped across the entire basin. The base-Quaternary has been mapped using the PGS MegaSurvey dataset from wells in the Danish Sector along the initially unconformable horizon and down the delta front into the more conformable basin giving a high degree of confidence in the horizon pick. The revised base-Quaternary surface reaches a depth of 1248 ms TWT with an elongate basin shape which is significantly deeper than the traditionally mapped surface. Using RMS amplitudes and other seismic attributes the revised base-Quaternary has been investigated along the horizon and in time slice to interpret the environments of the earliest Quaternary prior to the onset of glaciation. Combined with analysis of aligned elongate furrows over 10 km long, 100 m wide and 100 m deep suggest a deep marine environment in an almost enclosed basin with persistent strong NW-SE bottom currents in the deepest parts. Pockmarks were formed by the escape of shallow gas on the sides of a small delta in the eastern part of the basin. The progradation of large deltas from both the north and south into the basin make up the majority of the deposition of sediment into the basin. Key Words: base-Quaternary; seismic interpretation; paleoenvironments References: Gatliff, R.W, Richards, P.C, Smith, K, Graham, C.C, McCormac, M, Smith, N.J.P, Long, D, Cameron, T.D.J, Evans, D, Stevenson, A.G, Bulat, J, Ritchie, J.D, (1994) 'United Kingdom offshore regional

  13. Probabilistic Seismic Hazard Maps for Seattle, Washington, Based on 3D Ground-Motion Simulations

    Frankel, A. D.; Stephenson, W. J.; Carver, D. L.; Williams, R. A.; Odum, J. K.; Rhea, S.

    2007-12-01

    We have produced probabilistic seismic hazard maps for Seattle using over 500 3D finite-difference simulations of ground motions from earthquakes in the Seattle fault zone, Cascadia subduction zone, South Whidbey Island fault, and background shallow and deep source areas. The maps depict 1 Hz response spectral accelerations with 2, 5, and 10% probabilities of being exceeded in 50 years. The simulations were used to generate site and source dependent amplification factors that are applied to rock-site attenuation relations. The maps incorporate essentially the same fault sources and earthquake recurrence times as the 2002 national seismic hazard maps. The simulations included basin surface waves and basin-edge focusing effects from a 3D model of the Seattle basin. The 3D velocity model was validated by modeling several earthquakes in the region, including the 2001 M6.8 Nisqually earthquake, that were recorded by our Seattle Urban Seismic Network and the Pacific Northwest Seismic Network. The simulations duplicate our observation that earthquakes from the south and southwest typically produce larger amplifications in the Seattle basin than earthquakes from other azimuths, relative to rock sites outside the basin. Finite-fault simulations were run for earthquakes along the Seattle fault zone, with magnitudes ranging from 6.6 to 7.2, so that the effects of rupture directivity were included. Nonlinear amplification factors for soft-soil sites of fill and alluvium were also applied in the maps. For the Cascadia subduction zone, 3D simulations with point sources at different locations along the zone were used to determine amplification factors across Seattle expected for great subduction-zone earthquakes. These new urban seismic hazard maps are based on determinations of hazard for 7236 sites with a spacing of 280 m. The maps show that the highest hazard locations for this frequency band (around 1 Hz) are soft-soil sites (fill and alluvium) within the Seattle basin and

  14. 3D Numerical Simulation on the Sloshing Waves Excited by the Seismic Shacking

    Zhang, Lin; Wu, Tso-Ren

    2016-04-01

    In the event of 2015 Nepal earthquake, a video clip broadcasted worldwide showed a violent water spilling in a hotel swimming pool. This sloshing phenomenon indicates a potential water loss in the sensitive facilities, e.g. the spent fuel pools in nuclear power plant, has to be taken into account carefully under the consideration of seismic-induced ground acceleration. In the previous studies, the simulation of sloshing mainly focused on the pressure force on the structure by using a simplified Spring-Mass Method developed in the field of solid mechanics. However, restricted by the assumptions of plane water surface and limited wave height, significant error will be made in evaluating the amount of water loss in the tank. In this paper, the computational fluid dynamical model, Splash3D, was adopted for studying the sloshing problem accurately. Splash3D solved 3D Navier-Stokes Equation directly with Large-Eddy Simulation (LES) turbulent closure. The Volume-of-fluid (VOF) method with piecewise linear interface calculation (PLIC) was used to track the complex breaking water surface. The time series acceleration of a design seismic was loaded to excite the water. With few restrictions from the assumptions, the accuracy of the simulation results were improved dramatically. A series model validations were conducted by compared to a 2D theoretical solution, and a 3D experimental data. Good comparisons can be seen. After the validation, we performed the simulation for considering a sloshing case in a rectangular water tank with a dimension of 12 m long, 8 m wide, 8 m deep, which contained water with 7 m in depth. The seismic movement was imported by considering time-series acceleration in three dimensions, which were about 0.5 g to 1.2 g in the horizontal directions, and 0.3 g to 1 g in the vertical direction. We focused the discussions on the kinematics of the water surface, wave breaking, velocity field, pressure field, water force on the side walls, and, most

  15. 3D geological to geophysical modelling and seismic wave propagation simulation: a case study from the Lalor Lake VMS (Volcanogenic Massive Sulphides) mining camp

    Miah, Khalid; Bellefleur, Gilles

    2014-05-01

    The global demand for base metals, uranium and precious metals has been pushing mineral explorations at greater depth. Seismic techniques and surveys have become essential in finding and extracting mineral rich ore bodies, especially for deep VMS mining camps. Geophysical parameters collected from borehole logs and laboratory measurements of core samples provide preliminary information about the nature and type of subsurface lithologic units. Alteration halos formed during the hydrothermal alteration process contain ore bodies, which are of primary interests among geologists and mining industries. It is known that the alteration halos are easier to detect than the ore bodies itself. Many 3D geological models are merely projection of 2D surface geology based on outcrop inspections and geochemical analysis of a small number of core samples collected from the area. Since a large scale 3D multicomponent seismic survey can be prohibitively expensive, performance analysis of such geological models can be helpful in reducing exploration costs. In this abstract, we discussed challenges and constraints encountered in geophysical modelling of ore bodies and surrounding geologic structures from the available coarse 3D geological models of the Lalor Lake mining camp, located in northern Manitoba, Canada. Ore bodies in the Lalor lake VMS camp are rich in gold, zinc, lead and copper, and have an approximate weight of 27 Mt. For better understanding of physical parameters of these known ore bodies and potentially unknown ones at greater depth, we constructed a fine resolution 3D seismic model with dimensions: 2000 m (width), 2000 m (height), and 1500 m (vertical depth). Seismic properties (P-wave, S-wave velocities, and density) were assigned based on a previous rock properties study of the same mining camp. 3D finite-difference elastic wave propagation simulation was performed in the model using appropriate parameters. The generated synthetic 3D seismic data was then compared to

  16. Joint environmental assessment for Chevron USA, Inc. and Santa Fe Energy Resources, Inc.: Midway Valley 3D seismic project, Kern County, California

    NONE

    1996-10-01

    The proposed Midway Valley 3D Geophysical Exploration Project covers approximately 31,444 aces of private lands, 6,880 acres of Department of Energy (DOE) Lands within Naval Petroleum Reserve 2 (NPR2) and 3,840 acres of lands administered by the Bureau of Land Management (BLM), in western Kern County, California. This environmental assessment (EA) presents an overview of the affected environment within the project area using results of a literature review of biological field surveys previously conducted within or adjacent to a proposed 3D seismic project. The purpose is to provide background information to identify potential and known locations of sensitive wildlife and special status plant species within the proposed seismic project area. Biological field surveys, following agency approved survey protocols, will be conducted during October through November 1996 to acquire current resources data to provide avoidance as the project is being implemented in the field.

  17. Joint environmental assessment for Chevron USA, Inc. and Santa Fe Energy Resources, Inc.: Midway Valley 3D seismic project, Kern County, California

    The proposed Midway Valley 3D Geophysical Exploration Project covers approximately 31,444 aces of private lands, 6,880 acres of Department of Energy (DOE) Lands within Naval Petroleum Reserve 2 (NPR2) and 3,840 acres of lands administered by the Bureau of Land Management (BLM), in western Kern County, California. This environmental assessment (EA) presents an overview of the affected environment within the project area using results of a literature review of biological field surveys previously conducted within or adjacent to a proposed 3D seismic project. The purpose is to provide background information to identify potential and known locations of sensitive wildlife and special status plant species within the proposed seismic project area. Biological field surveys, following agency approved survey protocols, will be conducted during October through November 1996 to acquire current resources data to provide avoidance as the project is being implemented in the field

  18. 3D seismic imaging of the subsurface for underground construction and drilling

    3D seismic imaging of underground structure has been carried out in various parts of the world for various purposes. Examples shown below were introduced in the presentation. - CO2 storage in Ketzin, Germany; - Mine planning at the Millennium Uranium Deposit in Canada; - Planned Forsmark spent nuclear fuel repository in Sweden; - Exploring the Scandinavian Mountain Belt by Deep Drilling: the COSC drilling project in Sweden. The author explained that seismic methods provide the highest resolution images (5-10 m) of deeper (1-5 km) sub-surfaces in the sedimentary environment, but further improvement is required in crystalline rock environments, and the integration of geology, geophysics, and drilling will provide an optimal interpretation. (author)

  19. Intriguing Success in 3D Seismic Acquisition in Ecologically Critical Lawachara National Park of Bangladesh

    Bakht, Delawar; Siddique, Mohammad; Masud, Mohammad

    2010-09-15

    In-depth environmental studies were conducted in 2008 by a multi-disciplinary team of international and national specialists of SMEC International for Chevron Bangladesh for obtaining Environmental Clearance for 3D seismic acquisition in Moulvibazar Gas Field. This included Lawachara National Park which was declared as an ecologically critical area in 1996. Exclusive monitoring of potential impact mitigation mechanism identified through EIA studies resulted in to completing the project with intriguing success. This has displayed a glaring example of sharing expertise leading to successful initiative in technology transfer in the developing country like Bangladesh currently in dire quest of harnessing natural gas.

  20. High-resolution 3D seismic model of the crustal and uppermost mantle structure in Poland

    Grad, Marek; Polkowski, Marcin; Ostaficzuk, Stanisław R.

    2016-01-01

    In the area of Poland a contact between the Precambrian and Phanerozoic Europe and the Carpathians has a complicated structure and a complex P-wave velocity of the sedimentary cover, crystalline crust, Moho depth and the uppermost mantle. The geometry of the uppermost several kilometers of sediments is relatively well recognized from over 100,000 boreholes. The vertical seismic profiling (VSP) from 1188 boreholes provided detailed velocity data for regional tectonic units and for stratigraphic successions from Permian to the Tertiary and Quaternary deposits. These data, however, do not provide information about the velocity and basement depth in the central part of the Trans-European suture zone (TESZ) and in the Carpathians. So, the data set is supplemented by 2D velocity models from 32 deep seismic sounding refraction profiles which also provide information about the crust and uppermost mantle. Together with the results of other methods: vertical seismic profiling, magnetotelluric, allow for the creation of a detailed, high-resolution 3D model for the entire Earth's crust and the uppermost mantle down to a depth of 60 km. The thinnest sedimentary cover in the Mazury-Belarus anteclise is only 0.3 to 1 km thick, which increases to 7 to 8 km along the East European Craton (EEC) margin, and 9 to 12 km in the TESZ. The Variscan domain is characterized by a 1-4 km thick sedimentary cover, while the Carpathians are characterized by very thick sedimentary layers, up to about 20 km. The crystalline crust is differentiated and has a layered structure. The crust beneath the West European Platform (WEP; Variscan domain) is characterized by P-wave velocities of 5.8-6.6 km/s. The upper and middle crusts beneath the EEC are characterized by velocities of 6.1-6.6 km/s, and are underlain by a high velocity lower crust with a velocity of about 7 km/s. A general decrease in velocity is observed from the older to the younger tectonic domains. The TESZ is associated with a steep dip

  1. SEISVIZ3D: Stereoscopic system for the representation of seismic data - Interpretation and Immersion

    von Hartmann, Hartwig; Rilling, Stefan; Bogen, Manfred; Thomas, Rüdiger

    2015-04-01

    The seismic method is a valuable tool for getting 3D-images from the subsurface. Seismic data acquisition today is not only a topic for oil and gas exploration but is used also for geothermal exploration, inspections of nuclear waste sites and for scientific investigations. The system presented in this contribution may also have an impact on the visualization of 3D-data of other geophysical methods. 3D-seismic data can be displayed in different ways to give a spatial impression of the subsurface.They are a combination of individual vertical cuts, possibly linked to a cubical portion of the data volume, and the stereoscopic view of the seismic data. By these methods, the spatial perception for the structures and thus of the processes in the subsurface should be increased. Stereoscopic techniques are e. g. implemented in the CAVE and the WALL, both of which require a lot of space and high technical effort. The aim of the interpretation system shown here is stereoscopic visualization of seismic data at the workplace, i.e. at the personal workstation and monitor. The system was developed with following criteria in mind: • Fast rendering of large amounts of data so that a continuous view of the data when changing the viewing angle and the data section is possible, • defining areas in stereoscopic view to translate the spatial impression directly into an interpretation, • the development of an appropriate user interface, including head-tracking, for handling the increased degrees of freedom, • the possibility of collaboration, i.e. teamwork and idea exchange with the simultaneous viewing of a scene at remote locations. The possibilities offered by the use of a stereoscopic system do not replace a conventional interpretation workflow. Rather they have to be implemented into it as an additional step. The amplitude distribution of the seismic data is a challenge for the stereoscopic display because the opacity level and the scaling and selection of the data have to

  2. Survey of Robot 3D Path Planning Algorithms

    Liang Yang

    2016-01-01

    Full Text Available Robot 3D (three-dimension path planning targets for finding an optimal and collision-free path in a 3D workspace while taking into account kinematic constraints (including geometric, physical, and temporal constraints. The purpose of path planning, unlike motion planning which must be taken into consideration of dynamics, is to find a kinematically optimal path with the least time as well as model the environment completely. We discuss the fundamentals of these most successful robot 3D path planning algorithms which have been developed in recent years and concentrate on universally applicable algorithms which can be implemented in aerial robots, ground robots, and underwater robots. This paper classifies all the methods into five categories based on their exploring mechanisms and proposes a category, called multifusion based algorithms. For all these algorithms, they are analyzed from a time efficiency and implementable area perspective. Furthermore a comprehensive applicable analysis for each kind of method is presented after considering their merits and weaknesses.

  3. HIRE seismic reflection survey in the Olkiluoto area

    A seismic reflection survey comprising three vibroseismic lines (total length of 31.1 km) was carried out in the Olkiluoto area, western Finland, in July, 2008. The survey is a part of the project HIRE (High Resolution Reflection Seismics for Ore Exploration 2007-2010) of the Geological Survey of Finland (GTK). The Olkiluoto survey was done in co-operation with Posiva Oy. The HIRE seismic reflection survey in the Olkiluoto area revealed numerous previously unknown structures in the upper crust of the area. The most prominent structures observed are the subhorizontal strong reflectors which very probably represent Postjotnian diabase sills intruding both the Svecofennian gneisses as well as the rapakivi granites. These reflectors can be associated with the similar seismic structures recorded in marine seismic transects in the Bothnian Sea, and thus they represent a large-scale structure. The Mesoproterozoic rapakivi granites can be distinguished as homogeneous, seismically transparent domains which extend to a depth of at least 4 km. The interpreted rapakivi structures are in a good agreement with gravity modellings. On the Olkiluoto Island, reflectors could be correlated with drillhole based data on lithology and brittle fault zones. The main brittle fault zones detected in drillholes are represented as reflectors in the seismic sections, and several new structures have been interpreted. A synthetic seismogram constructed for a 1 km deep hole in the Olkiluoto Island using down-hole logs of density and P-wave velocity, suggests that the main brittle fault zones generate strong reflectors. Pegmatitic granite sometimes has a weak reflection contrast with the surrounding gneisses, and sometimes has no contrast. There seems to be a simultaneous correlation of reflectivity with fracture zones at the location of the pegmatitic granite occurrences. The HIRE 2D seismic reflection data agrees well with the earlier 3D reflection surveys in the Olkiluoto Island, and most of

  4. 3-D structure below Aevroe Island from high-resolution reflection seismic studies, southeastern Sweden

    Reflection seismology has served as a useful tool for imaging and mapping of fracture zones in crystalline rock along 2-D lines in nuclear waste disposal studies. Two 1-km-long perpendicular seismic reflection lines were acquired on Aevroe Island, southeast Sweden, in October 1996 in order to (1) test the seismic reflection method for future site investigations, (2) map known fracture zones, and (3) add to the Swedish database of reflection seismic studies of the shallow crystalline crust. An east-west line was shot with 5-m geophone and shot-point spacing, and a north-south line was shot with 10-m geophone and shotpoint spacing. An explosive source with a charge size of 100 g was used along both lines. The data clearly image three major dipping reflectors and one subhorizontal one in the upper 200 ms (600 m). The dipping reflectors (to the south, east, and northwest) intersect or project to the surface at or close to where surface-mapped fracture zones exist. The south-dipping reflector correlates with the top of a heavily fractured interval observed in a borehole (KAV01) at about 400 m. The subhorizontal zone at about 100--200 m correlates with a known fracture zone in the same borehole (KAV01). 3-D effects are apparent in the data, and only where the profiles cross can the true orientation of the reflecting events be determined. To properly orient and locate all events observed on the lines requires acquisition of 3-D data

  5. Q AS A LITHOLOGICAL/HYDROCARBON INDICATOR: FROM FULL WAVEFORM SONIC TO 3D SURFACE SEISMIC

    Jorge O. Parra; C.L. Hackert; L. Wilson; H.A. Collier; J. Todd Thomas

    2006-03-31

    The goal of this project was to develop a method to exploit viscoelastic rock and fluid properties to greatly enhance the sensitivity of surface seismic measurements to the presence of hydrocarbon saturation. To reach the objective, Southwest Research Institute scientists used well log, lithology, production, and 3D seismic data from an oil reservoir located on the Waggoner Ranch in north central Texas. The project was organized in three phases. In the first phase, we applied modeling techniques to investigate seismic- and acoustic-frequency wave attenuation and its effect on observable wave attributes. We also gathered existing data and acquired new data from the Waggoner Ranch field, so that all needed information was in place for the second phase. During the second phase, we developed methods to extract attenuation from borehole acoustic and surface seismic data. These methods were tested on synthetic data constructed from realistic models and real data. In the third and final phase of the project, we applied this technology to a full data set from the Waggoner site. The results presented in this Final Report show that geological conditions at the site did not allow us to obtain interpretable results from the Q processing algorithm for 3D seismic data. However, the Q-log processing algorithm was successfully applied to full waveform sonic data from the Waggoner site. A significant part of this project was technology transfer. We have published several papers and conducted presentations at professional conferences. In particular, we presented the Q-log algorithm and applications at the Society of Exploration Geophysicists (SEG) Development and Production Forum in Austin, Texas, in May 2005. The presentation attracted significant interest from the attendees and, at the request of the SEG delegates, it was placed on the Southwest Research Institute Internet site. The presentation can be obtained from the following link: http://www.swri.org/4org/d15/elecsys

  6. 3D elastic inversion of vertical seismic profiles in horizontally stratified media; Inversion elastique 3D de profils sismiques verticaux en milieux stratifies horizontalement

    Petit, J.L.

    1997-07-21

    This thesis is devoted to the inversion of VSP (vertical seismic profile) seismic data in order to determine the elastic properties of horizontally stratified media. The VSP records are computed using the full wave elastic modelling in isotropic and transversely isotropic media using Hankel transform, a finite difference scheme and an inverse Hankel transform algorithm, and the propagation equations are determined and numerically solved; the importance of considering a 3D wave propagation model instead of a 1 D one is emphasized. The theoretical VSP inverse problem is then considered, with the seismic waveform inversion set as a least-squares problem, consisting in recovering the distribution of physical parameters which minimize the misfit between calculated and observed VSP. The corresponding problem requires the knowledge of the source function

  7. Improving Geologic and Engineering Models of Midcontinent Fracture and Karst-Modified Reservoirs Using New 3-D Seismic Attributes

    Susan Nissen; Saibal Bhattacharya; W. Lynn Watney; John Doveton

    2009-03-31

    Our project goal was to develop innovative seismic-based workflows for the incremental recovery of oil from karst-modified reservoirs within the onshore continental United States. Specific project objectives were: (1) to calibrate new multi-trace seismic attributes (volumetric curvature, in particular) for improved imaging of karst-modified reservoirs, (2) to develop attribute-based, cost-effective workflows to better characterize karst-modified carbonate reservoirs and fracture systems, and (3) to improve accuracy and predictiveness of resulting geomodels and reservoir simulations. In order to develop our workflows and validate our techniques, we conducted integrated studies of five karst-modified reservoirs in west Texas, Colorado, and Kansas. Our studies show that 3-D seismic volumetric curvature attributes have the ability to re-veal previously unknown features or provide enhanced visibility of karst and fracture features compared with other seismic analysis methods. Using these attributes, we recognize collapse features, solution-enlarged fractures, and geomorphologies that appear to be related to mature, cockpit landscapes. In four of our reservoir studies, volumetric curvature attributes appear to delineate reservoir compartment boundaries that impact production. The presence of these compartment boundaries was corroborated by reservoir simulations in two of the study areas. Based on our study results, we conclude that volumetric curvature attributes are valuable tools for mapping compartment boundaries in fracture- and karst-modified reservoirs, and we propose a best practices workflow for incorporating these attributes into reservoir characterization. When properly calibrated with geological and production data, these attributes can be used to predict the locations and sizes of undrained reservoir compartments. Technology transfer of our project work has been accomplished through presentations at professional society meetings, peer-reviewed publications

  8. 3D seismic analysis of gravity-driven and basement influenced normal fault growth in the deepwater Otway Basin, Australia

    Robson, A. G.; King, R. C.; Holford, S. P.

    2016-08-01

    We use three-dimensional (3D) seismic reflection data to analyse the structural style and growth of a normal fault array located at the present-day shelf-edge break and into the deepwater province of the Otway Basin, southern Australia. The Otway Basin is a Late Jurassic to Cenozoic, rift-to-passive margin basin. The seismic reflection data images a NW-SE (128-308) striking, normal fault array, located within Upper Cretaceous clastic sediments and which consists of ten fault segments. The fault array contains two hard-linked fault assemblages, separated by only 2 km in the dip direction. The gravity-driven, down-dip fault assemblage is entirely contained within the 3D seismic survey, is located over a basement plateau and displays growth commencing and terminating during the Campanian-Maastrichtian, with up to 1.45 km of accumulated throw (vertical displacement). The up-dip normal fault assemblage penetrates deeper than the base of the seismic survey, but is interpreted to be partially linked along strike at depth to major basement-involved normal faults that can be observed on regional 2D seismic lines. This fault assemblage displays growth initiating in the Turonian-Santonian and has accumulated up to 1.74 km of throw. Our detailed analysis of the 3D seismic data constraints post-Cenomanian fault growth of both fault assemblages into four evolutionary stages: [1] Turonian-Santonian basement reactivation during crustal extension between Australia and Antarctica. This either caused the upward propagation of basement-involved normal faults or the nucleation of a vertically isolated normal fault array in shallow cover sediments directly above the reactivated basement-involved faults; [2] continued Campanian-Maastrichtian crustal extension and sediment loading eventually created gravitational instability on the basement plateau, nucleating a second, vertically isolated normal fault array in the cover sediments; [3] eventual hard-linkage of fault segments in both fault

  9. Subglacial Landforms and Processes: new Information From 3D Seismic Technology

    Andreassen, K.

    2007-12-01

    Three-dimensional (3D) seismic interpretation and imaging techniques provide a unique means of investigating submarine geomorphic features produced by former ice sheets. An extensive two-dimensional (2D) and 3D seismic data base is here used to image the imprints left behind by glaciers that flowed out a major cross-shelf trough (Bjornoyrenna) of the north-Norwegian continental shelf during repeated glacial episodes. Mega-scale glacial lineations characterize the seafloor geomorphology of Bjornoyrenna and smaller, contributing cross-shelf troughs, where they are inferred to represent flow-lines of former ice streams that where active during the most recent (Weichselian) glacial period. Similar features are commonly observed on buried horizons. Large- scale seafloor imprints from an early readvance after the last glacial maximum are especially well preserved. Streamlined landforms and associated lobe-shaped ridges indicate that this major cross-shelf trough hosted six separate ice stream lobes that diverged fan-like at their margins, but were not all active simultaneously. A 300 km wide grounding-zone wedge results from high sediment flux within sub-ice stream deformable beds. A 2 to 3 km thick Pleistocene record is preserved at the mouth of Bjornoyrenna, in the Bjornoya Trough Mouth Fan. The preservation of up to several hundred meters of glacigenic sediments between the buried, glacially eroded surfaces, provides here the opportunity to study the internal structure of till units. 3D seismic attribute maps reveal that megablocks and rafts commonly occur within the till units. The sediments blocks are often aligned in chains that may be up to 2 km wide and over 50 km long. The largest individual megablocks have an areal extent of over 2 km2. The sediment chains are interpreted to have been eroded, transported and deposited by grounded ice, most probably fast-flowing ice streams. This is based on the relationship between the sediment chains and the horizons revealing

  10. Imaging 3D anisotropic upper mantle shear velocity structure of Southeast Asia using seismic waveform inversion

    Chong, J.; Yuan, H.; French, S. W.; Romanowicz, B. A.; Ni, S.

    2011-12-01

    Southeast Asia as a special region in the world which is seismically active and is surrounded by active tectonic belts, such as the Himalaya collision zone, western Pacific subduction zones and the Tianshan- Baikal tectonic belt. Seismic anisotropic tomography can shade light on the complex crust and upper mantle dynamics of this region, which is the subject of much debate. In this study, we applied full waveform time domain tomography to image 3D isotropic and anisotropic upper mantle shear velocity structure of Southeast Asia. Three component waveforms of teleseismic and far regional events (15 degree ≤ Δ≤ 165 degree) with magnitude ranges from Mw6.0 to Mw7.0 are collected from 91 permanent and 438 temporary broadband seismic stations in SE Asia. Wavepackets of both fundamental and overtone modes, filtered between 60 and 400 sec, are selected automatically according to the similarity between data and synthetic waveforms (Panning & Romanowicz, 2006). Wavepackets corresponding to event-station paths that sample the region considered are weighted according to path redundancy and signal to noise ratio. Higher modes and fundamental mode wavepackets are weighted separately in order to enhance the contribution of higher modes which are more sensitive to deeper structure compared to the fundamental mode. Synthetic waveforms and broadband sensitivity kernels are computed using normal mode asymptotic coupling theory (NACT, Li & Romanowicz, 1995). As a starting model, we consider a global anisotropic upper mantle shear velocity model based on waveform inversion using the Spectral Element Method (Lekic & Romanowicz, 2011), updated for more realistic crustal thickness (French et al., 2011) as our starting model, we correct waveforms for the effects of 3D structure outside of the region, and invert them for perturbations in the 3D structure of the target region only. We start with waveform inversion down to 60sec and after several iterations, we include shorter period

  11. 9C-3D seismic interpretation of the Bakken Formation, Banner Field, North Dakota

    Comegys, Lillian R.

    The Bakken Petroleum System is a multi-reservoir play with estimated total undiscovered resources of 3.649 BBO oil and 1.85 TCF natural gas in the United States portion of the Williston Basin (Pollastro 2008). The presence of natural fractures in all three members of the Bakken Formation have been linked to high initial production (IP) and cumulative production from the Antelope Field and better reservoir potential in the Elm Coulee Field and Sanish Fields (Sturm and Gomez 2009; Honsberger 2012; Theloy 2011). Therefore, the ability of seismic data to determine the presence, orientation, and density of natural fractures is an important achievement for petroleum exploration and exploitation. The STAMPEDE 9-component seismic survey is located in Mountrail County, North Dakota, in the Banner Field, southeast of the Parshall and Sanish Fields. It is the goal of the Reservoir Characterization Project to analyze the structural influences on reservoir properties in the STAMPEDE survey area using the compressional and pure shear seismic volumes supplemented by the public well information available on the North Dakota Industrial Commission website. Fracturing induced by basement faulting and lithology changes is detectable using multicomponent seismic data in the Stampede seismic survey. Shear wave splitting analysis delineates zones of different fracture orientation and density. These areas correlate to interpreted fault intersections and the predicted area of increased fracture frequency based on facies changes in the Middle Bakken Member and its mechanical stratigraphy. Wrench fault mechanics are at work in the study area, creating isolated convergent and divergent stress regimes in the separate fault blocks. Main fault interpretations are based on shear wave isochron mapping, wireline log mapping, seismic panel observations. Fracture interpretations were made on the analysis of shear time and amplitude anisotropy maps and the correlation of a P-wave Velocity Variation

  12. 3-D seismic response of a base-isolated fast reactor

    This paper describes a 3-D response analysis methodology development and its application to a base-isolated fast breeder reactor (FBR) plant. At first, studies on application of a base-isolation system to an FBR plant were performed to identify a range of appropriate characteristics of the system. A response analysis method was developed based on mathematical models for the restoring force characteristics of several types of the systems. A series of shaking table tests using a small scale model was carried out to verify the analysis method. A good agreement was seen between the test and analysis results in terms of the horizontal and vertical responses. Parametric studies were then made to assess the effects of various factors which might be influential to the seismic response of the system. Moreover, the method was applied to evaluate three-dimensional response of the base-isolated FBR. (author)

  13. How an independent put inexpensive 3-D seismic to good use in New Mexico

    Nester, D.C. (Landmark/Concurrent Solutions, Houston, TX (US)); Emdsley, D. (Merrion Oil and Gas Corp., Farmington, NM (US))

    1992-03-23

    This paper reports that as major oil companies focus their attention offshore and overseas, independents in the U.S. continue to find and develop onshore reserves using their traditional strengths. Those strengths have been low overhead, nimble decisionmaking, and hard-won experience within a particular geologic region. Today many of these companies are emerging as even tougher competitors by applying 3-D seismic along with the multidisciplinary know-how and the interactive workstations needed to interpret it. The recent experience of Merrion Oil and Gas Corp., Farmington, N.M., offers a case in point. When oil prices collapsed in 1986, Merrion saw the chance to venture into exploration at low cost. Ever since, the company has devoted its resources to exploring for oil in the San Juan basin's Entrada sand dunes, a complex stratigraphic play found at 6,000 ft.

  14. Delineation of Gulf of Mexico Miocene reservoirs with 3-D seismic data

    Wellborn, J.E. (Conoco Inc., Houston, TX (United States)); Brown, A.R.

    1991-03-01

    Seismic amplitude anomalies associated with Miocene plays in the Mobile area, offshore Alabama, are well documented and well drilled. Three-dimensional seismic data permit a much better understanding of these sandstone reservoirs than the 2-D data available over most prospects. In this paper Digicon's Mobile Phase 1 3-D data have been studied and the extent of two hydrocarbon-bearing sand bodies recognized on both vertical and horizontal sections. The first example shows a sand body interpreted as an offshore bar. The second shows a more complex reservoir geometry and is interpreted as a deltaic fan deposit. Both reservoirs are stratigraphically controlled. Horizon slices were constructed for each reservoir by amplitude extraction using Conoco's interactive interpretation system. The resultant contours defined the areal extent of the reservoirs. After successful gas wells were drilled into both features, pay thicknesses and reserve calculations were made by integrating the amplitudes on the horizon slices and the well information in preparation for production.

  15. Refining seismic parameters in low seismicity areas by 3D trenching: The Alhama de Murcia fault, SE Iberia

    Ferrater, Marta; Ortuño, Maria; Masana, Eulàlia; Pallàs, Raimon; Perea, Hector; Baize, Stephane; García-Meléndez, Eduardo; Martínez-Díaz, José J.; Echeverria, Anna; Rockwell, Thomas K.; Sharp, Warren D.; Medialdea, Alicia; Rhodes, Edward J.

    2016-06-01

    Three-dimensional paleoseismology in strike-slip faults with slip rates less than 1 mm per year involves a great methodological challenge. We adapted 3D trenching to track buried channels offset by the Alhama de Murcia seismogenic left-lateral strike-slip fault (SE Iberia). A fault net slip of 0.9 ± 0.1 mm/yr was determined using statistical analysis of piercing lines for one buried channel, whose age is constrained between 15.2 ± 1.1 ka and 21.9-22.3 cal BP. This value is larger and more accurate than the previously published slip rates for this fault. The minimum number of five paleo-earthquakes identified since the deposition of dated layers suggests a maximum average recurrence interval of approximately 5 ka. The combination of both seismic parameters yields a maximum slip per event between 5.3 and 6.3 m. We show that accurately planned trenching strategies and data processing may be key to obtaining robust paleoseismic parameters in low seismicity areas.

  16. Modeling the Coast Mountains Batholith, British Columbia, Canada with 3D Seismic Tomography

    Quinonez, S. M.; Olaya, J. C.; Miller, K. C.; Romero, R.; Velasco, A. A.; Harder, S. H.; Cerda, I.

    2011-12-01

    The Coast Mountains Batholith on the west coast of British Columbia, Canada comprises a series of granitic to tonalitic plutons; where felsic continental crust is generated from the subduction of mafic oceanic crust by partial melting and fractionation, leaving ultra-mafic roots. In July of 2009, a large controlled-source experiment was conducted along a 400km east - west transect from Bella Bella into central British Columbia. Student volunteers from multiple universities deployed 1,800 one-component and 200 three-component geophones plus 2400 Texan data recorders with 200-m spacing intervals and shot spacing at 30-km. The 18-point sources ranged from 160 to 1,000 kg of high explosive. The geoscience component of the NSF-funded Cyber-ShARE project at UTEP focuses on fusing models developed from different data sets to develop 3-D Earth models. Created in 2007, the Cyber-ShARE Center brings together experts in computer science, computational mathematics, education, earth science, and environmental science. We leverage the Cyber-ShARE work to implement an enhanced 3-D finite difference tomography approach for P-wave delays times (Hole, 1992) with a graphical user interface and visualization framework. In particular, to account for model sensitivity to picked P-wave arrival times, we use a model fusion approach (Ochoa et al., 2010) to generate a model with the lowest RMS residual that a combination of a set of Monte Carlo sample models. In order to make the seismic tomography process more interactive at many points, visualizations of model perturbation at each iteration will help to troubleshoot when a model is not converging to highlight where the RMS residual values are the highest to pinpoint where changes need to be made to achieve model convergence. Finally, a model of the upper mantle using 3-D P-wave tomography will be made to determine the location of these ultra-mafic roots.

  17. Full Waveform 3D Synthetic Seismic Algorithm for 1D Layered Anelastic Models

    Schwaiger, H. F.; Aldridge, D. F.; Haney, M. M.

    2007-12-01

    Numerical calculation of synthetic seismograms for 1D layered earth models remains a significant aspect of amplitude-offset investigations, surface wave studies, microseismic event location approaches, and reflection interpretation or inversion processes. Compared to 3D finite-difference algorithms, memory demand and execution time are greatly reduced, enabling rapid generation of seismic data within workstation or laptop computational environments. We have developed a frequency-wavenumber forward modeling algorithm adapted to realistic 1D geologic media, for the purpose of calculating seismograms accurately and efficiently. The earth model consists of N layers bounded by two halfspaces. Each layer/halfspace is a homogeneous and isotropic anelastic (attenuative and dispersive) solid, characterized by a rectangular relaxation spectrum of absorption mechanisms. Compressional and shear phase speeds and quality factors are specified at a particular reference frequency. Solution methodology involves 3D Fourier transforming the three coupled, second- order, integro-differential equations for particle displacements to the frequency-horizontal wavenumber domain. An analytic solution of the resulting ordinary differential system is obtained. Imposition of welded interface conditions (continuity of displacement and stress) at all interfaces, as well as radiation conditions in the two halfspaces, yields a system of 6(N+1) linear algebraic equations for the coefficients in the ODE solution. An optimized inverse 2D Fourier transform to the space domain gives the seismic wavefield on a horizontal plane. Finally, three-component seismograms are obtained by accumulating frequency spectra at designated receiver positions on this plane, followed by a 1D inverse FFT from angular frequency ω to time. Stress-free conditions may be applied at the top or bottom interfaces, and seismic waves are initiated by force or moment density sources. Examples reveal that including attenuation

  18. Probabilistic Seismic Hazard Maps of Seattle, Washington, Including 3D Sedimentary Basin Effects and Rupture Directivity: Implications of 3D Random Velocity Variations (Invited)

    Frankel, A. D.; Stephenson, W. J.; Carver, D.; Odum, J.; Williams, R. A.; Rhea, S.

    2010-12-01

    We have produced probabilistic seismic hazard maps of Seattle for 1 Hz spectral acceleration, using over five hundred 3D finite-difference simulations of earthquakes on the Seattle fault, Southern Whidbey Island fault, and Cascadia subduction zone, as well as for random deep and shallow earthquakes at various locations. The 3D velocity model was validated by modeling the observed waveforms for the 2001 M6.8 Nisqually earthquake and several smaller events in the region. At these longer periods (≥ 1 sec) that are especially important to the response of buildings of ten stories or higher, seismic waves are strongly influenced by sedimentary basins and rupture directivity. We are investigating how random spatial variations in the 3D velocity model affect the simulated ground motions for M6.7 earthquakes on the Seattle fault. A fractal random variation of shear-wave velocity with a Von Karman correlation function produces spatial variations of peak ground velocity with multiple scale lengths. We find that a 3D velocity model with a 10% standard deviation in shear-wave velocity in the top 1.5 km and 5% standard deviation from 1.5-10 km depth produces variations in peak ground velocities of as much as a factor of two, relative to the case with no random variations. The model with random variations generally reduces the peak ground velocity of the forward rupture directivity pulse for sites near the fault where basin-edge focusing of S-waves occurs. It also tends to reduce the peak velocity of localized areas where basin surface waves are focused. However, the medium with random variations also causes small-scale amplification of ground motions over distances of a few kilometers. We are also evaluating alternative methods of characterizing the aleatory uncertainty in the probabilistic hazard calculations.

  19. 3D spectroscopic surveys: Exploring galaxy evolution mechanisms

    Epinat, Benoît

    2011-01-01

    I review the major surveys of high redshift galaxies observed using integral field spectroscopy techniques in the visible and in the infrared. The comparison of various samples has to be done with care since they have different properties linked to their parent samples, their selection criteria and the methods used to study them. I present the various kinematic types of galaxies that are identified within these samples (rotators, mergers, etc.) and summarize the discussions on the mass assembly processes at various redshifts deduced from these classifications: at intermediate redshift (z~0.6) merger may be the main mass assembly process whereas the role of cold gas accretion along cosmic web filaments may increase with redshift. The baryonic Tully-Fisher relation is also discussed. This relation seems to be already in place 3 Gyr after the Big-Bang and is then evolving until the present day. This evolution is interpreted as an increase of the stellar mass content of dark matter haloes of a given mass. The dis...

  20. 3D Seismic Velocity Structure in the Rupture Area of the 2010 Maule Mw=8.8 Earthquake

    Hicks, S. P.; Rietbrock, A.; Ryder, I. M.; Nippress, S.; Haberland, C. A.

    2011-12-01

    The 2010 Mw=8.8 Maule, Chile earthquake is one of the largest subduction zone earthquakes ever recorded. Up to now numerous co-seismic and some post-seismic slip models have been published based entirely on seismological, geodetic, or tsunami run-up heights, or combinations of these data. Most of these models use a simplified megathrust geometry derived mainly from global earthquake catalogues, and also simplified models of seismic parameters (e.g. shear modulus). By using arrival times for a vast number of aftershocks that have been recorded on a temporary seismic array, we present a new model for the slab geometry based on earthquake locations together with a new 3D seismic velocity model of the region, for both vp and vp/vs. We analyzed 3552 aftershocks that occurred between 18 March and 24 May 2011, recorded by the International Maule Aftershock Dataset (IMAD) seismic network. Event selection from a catalogue of automatically-determined events was based on 20 or more arrival times, from which at least 10 are S-wave observations. In total over 170,000 arrival times (~125,000 and 45,000 P and S wave arrival times respectively) are used for the tomographic reconstructions. Initially, events were relocated in a 2D velocity model based on a previously published model for the southern end of the rupture area (Haberland et al., 2009). Afterwards a staggered inversion scheme is implemented, starting with a 2D inversion followed by a coarse 3D and a subsequent fine 3D inversion. Based on our preliminary inversions we conclude that aftershock seismicity is mainly concentrated between 20 and 35 km depth along the subduction interface. A second band of seismicity between 40 and 50 km depth is also observed. Low seismic velocities and an increased vp/vs ratio characterize the marine forearc. The obtained velocity model will be discussed.

  1. Network level pavement evaluation with 1 mm 3D survey system

    Wang, Kelvin C.P.; Qiang Joshua Li; Guangwei Yang; You Zhan; Yanjun Qiu

    2015-01-01

    The latest iteration of PaveVision3D Ultra can obtain true 1 mm resolution 3D data at full-lane coverage in all 3 directions at highway speed up to 60 mph. This paper introduces the PaveVision3D Ultra technology for rapid network level pavement survey on approximately 1280 center miles of Oklahoma interstate highways. With sophisticated automated distress analyzer (ADA) software interface, the collected 1 mm 3D data provide Oklahoma Department of Transportation (ODOT) with comprehensive solut...

  2. 3D multicomponent seismic characterization of a clastic reservoir in the Middle Magdalena Valley Basin, Colombia

    Velasquez-Espejo, Antonio Jose

    The main goal of this research is to characterize the combined structural-stratigraphic trap of the Tenerife Field in the Middle Magdalena Valley Basin (MMVB), Colombia. For the first time in Colombia the structural and quantitative interpretation of modern three-dimensional multicomponent (3D-3C) seismic imaging enables a geometric description, a kinematic interpretation of the structural styles, and the facies distribution of the reservoir. A seismic petrophysics work-flow to better achieve the seismic well-tie. Edited and check-shot calibrated P-wave sonic logs were obtained and coefficients of the Gardner and Castagna equations were calibrated to match the density and shear-wave velocity depth trends for the basin. Seismic modeling was performed to evaluate the PP and PS seismic response of the reservoir interval (Mugrosa Formation). The structural interpretation methodology involves a 3D fault-correlation and horizon picking for both PP- and PS-PSTM data volumes. Geometric attributes such as coherence and curvature were used to enhance the structural discontinuities. The main unconformity of the Middle Eocene (MEU) was interpreted, and an attribute-assisted interpretation of the reservoir was conducted in detail. While P-wave data provided most of the structural interpretation, converted-wave data provide a better understanding of the faults. Traditionally, compressive thrust-propagation folds and tectonic inversion have been considered as the main mechanisms controlling the deformation in the MMVB. However, the new interpretation shown in this work provides a different structural concept that involves two major structural styles: 1. Under the MEU the Late Cretaceous and Early Paleocene deformation, dominated by east-verging thrust and partially inverted Mesozoic normal faults, is preserved. Associated folds exhibit a north-south strike, and their structural development is controlled by a long-lived structural element that dominates the area (the Infantas

  3. Comparative velocity structure of active Hawaiian volcanoes from 3-D onshore-offshore seismic tomography

    Park, J.; Morgan, J.K.; Zelt, C.A.; Okubo, P.G.; Peters, L.; Benesh, N.

    2007-01-01

    We present a 3-D P-wave velocity model of the combined subaerial and submarine portions of the southeastern part of the Island of Hawaii, based on first-arrival seismic tomography of marine airgun shots recorded by the onland seismic network. Our model shows that high-velocity materials (6.5-7.0??km/s) lie beneath Kilauea's summit, Koae fault zone, and the upper Southwest Rift Zone (SWRZ) and upper and middle East Rift Zone (ERZ), indicative of magma cumulates within the volcanic edifice. A separate high-velocity body of 6.5-6.9??km/s within Kilauea's lower ERZ and upper Puna Ridge suggests a distinct body of magma cumulates, possibly connected to the summit magma cumulates at depth. The two cumulate bodies within Kilauea's ERZ may have undergone separate ductile flow seaward, influencing the submarine morphology of Kilauea's south flank. Low velocities (5.0-6.3??km/s) seaward of Kilauea's Hilina fault zone, and along Mauna Loa's seaward facing Kao'iki fault zone, are attributed to thick piles of volcaniclastic sediments deposited on the submarine flanks. Loihi seamount shows high-velocity anomalies beneath the summit and along the rift zones, similar to the interpreted magma cumulates below Mauna Loa and Kilauea volcanoes, and a low-velocity anomaly beneath the oceanic crust, probably indicative of melt within the upper mantle. Around Kilauea's submarine flank, a high-velocity anomaly beneath the outer bench suggests the presence of an ancient seamount that may obstruct outward spreading of the flank. Mauna Loa's southeast flank is also marked by a large, anomalously high-velocity feature (7.0-7.4??km/s), interpreted to define an inactive, buried volcanic rift zone, which might provide a new explanation for the westward migration of Mauna Loa's current SWRZ and the growth of Kilauea's SWRZ. ?? 2007 Elsevier B.V. All rights reserved.

  4. Research on 3-D base isolation system applied to new power reactor 3-D seismic isolation device with rolling seal type air spring: Pt.2

    A three dimensional seismic base isolation device was developed for heavy structures and buildings such as nuclear power reactor buildings. The device realizes 3-D isolation by combining a LRB (laminated rubber bearing) for horizontal isolation with an air spring for vertical isolation in series. In this study, scale models of the 3-D base isolation device were made and were tested to examine the dynamic properties and ultimate strengths of the device. The performance of the device under earthquake excitation was examined through shaking table tests of 1/7 scale models. As the results, it was confirmed that the device worked smoothly under the horizontal and vertical excitations, and that the theoretical formulae of the orifice damping could explain the test results. The high-pressure air springs of trial production were forced to burst to find out which factor influenced ultimate strength. It was confirmed from results of the burst test that the strength of the air spring depended upon the diameter of rolling part of the bellows and the number of layers of the reinforcing fibers. Judging from the results of the shaking table test and the burst test, the developed 3-D base isolation device was applicable to a nuclear power plant building. (authors)

  5. 3D Seismic Reflection Data: Has the Geological Hubble Retained Its Focus?

    Jackson, Christopher

    2016-04-01

    In their seminal paper in 2002, Joe Cartwright and Mads Huuse referred to 3D seismic reflection data as the 'Geological Hubble', illustrating how these data had the potential to revolutionise our understanding of the genesis and evolution of sedimentary basins. 14 years on, I will here outline just some of the key recent advances made in our understanding of basin structure and stratigraphy, focusing on: (i) the intrusion and extrusion of igneous rocks; (ii) salt tectonics, with particular emphasis on intrasalt structure and the kinematics and mechanics of diapirism; (iii) the geometry and growth of normal faults; and (iv) the structure and emplacement of mass-transport complexes (MTCs). I will stress that future advances at least partly relies on hydrocarbon exploration companies and government agencies continuing to make their data freely available via easy-to-access data portals. I will issue a clarion call to academics, stressing that 'geodynamicists', sedimentologists, structural geologists and geomorphologists, amongst others, can benefit from utilising what I believe are currently an underused data type.

  6. Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M = 9.2 Sumatra earthquake

    Pollitz, F.; Banerjee, P.; Grijalva, K.; Nagarajan, B.; Burgmann, R.

    2008-01-01

    The 2004 M=9.2 Sumatra-Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ???1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress-strain relation; an aspherical perturbation in viscoelastic structure; a 'static'mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra-Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation. ?? Journal compilation ?? 2008 RAS.

  7. 3-D Seismic Images of Mud Volcano North Alex, West-Nile Delta, Egypt

    Bialas, Joerg; Klaeschen, Dirk; Papenberg, Cord; Gehrmann, Romina; Sommer, Malte

    2010-05-01

    Mud volcanoes within shelf areas are the bathymetric expression of mobilized overpressured sediments causing a feature of possible instability within the slope. Such a scene is given in the West-Nile Delta offshore Alexandria, Egypt at 700 m water depth, ,which was studied during a RWE Dea funded research project. The West Nile Delta forms part of the source of the large turbiditic Nile Deep Sea Fan. Since the late Miocene sediments have formed an up to 10 km thick pile, which includes about 1 - 3 km of Messinian evaporates. The sediment load of the overburden implies strong overpressures and salt-related tectonic deformation. Both are favourable for fluid migration towards the seafloor guided by the fractured margin. Deep-cutting channel systems like the Rosetta channel characterize the continental slope. Bathymetric expressions of slides and numerous mud volcanoes in the area are expressions of active processes, which contribute to the ongoing modification of the slope. The western deltaic system, Rosetta branch, has formed an 80 km wide continental shelf. Here at 700 m water depth the mud volcano North Alex developed his circular bathymetric feature, which proved to be an active gas and mud-expelling structure. A grid of 2-D seismic profiles did reveal a large set of faults located within the main mud volcano as well as surrounding the structure. Internal faults are mainly related to episodic mud expulsion processes and continuous gas and fluid production. Deep cutting external faults surround the structure in a half circle shape. They can be tracked up to the seafloor indicating ongoing tectonic activity of the slope area. A recently build 3-D acquisition system (funded by RWE Dea) suitable for mid-size research vessels was applied to collect an active seismic cube of the mud volcano. Based on the P-Cable design 11 parallel streamers (each 12.5 m long with 1.5 m group interval) were used to record shots of a single 210 cinch GI airgun. Based on GPS positions of

  8. Characterization of fracture reservoirs using static and dynamic data: From sonic and 3D seismic to permeability distribution. Annual report, March 1, 1996--February 28, 1997

    Parra, J.O.; Collier, H.A.; Owen, T.E. [and others

    1997-06-01

    In low porosity, low permeability zones, natural fractures are the primary source of permeability which affect both production and injection of fluids. The open fractures do not contribute much to porosity, but they provide an increased drainage network to any porosity. They also may connect the borehole to remote zones of better reservoir characteristics. An important approach to characterizing the fracture orientation and fracture permeability of reservoir formations is one based on the effects of such conditions on the propagation of acoustic and seismic waves in the rock. The project is a study directed toward the evaluation of acoustic logging and 3D-seismic measurement techniques as well as fluid flow and transport methods for mapping permeability anisotropy and other petrophysical parameters for the understanding of the reservoir fracture systems and associated fluid dynamics. The principal application of these measurement techniques and methods is to identify and investigate the propagation characteristics of acoustic and seismic waves in the Twin Creek hydrocarbon reservoir owned by Union Pacific Resources (UPR) and to characterize the fracture permeability distribution using production data. This site is located in the overthrust area of Utah and Wyoming. UPR drilled six horizontal wells, and presently UPR has two rigs running with many established drill hole locations. In addition, there are numerous vertical wells that exist in the area as well as 3D seismic surveys. Each horizontal well contains full FMS logs and MWD logs, gamma logs, etc.

  9. Characterization of the Hontomín Research Facility for Geological Storage of CO2: 3D Seismic Imaging Results

    Alcalde, J.; Martí, D.; Juhlin, C.; Malehmir, A.; Calahorrano, A.; Ayarza, P.; Pérez-Estaún, A.; Carbonell, R.

    2012-04-01

    A technological research facility dedicated to the underground geological storage of CO2 is currently being developed by the Spanish research program on Carbon Capture and Storage (CCS) in Hontomin (Burgos). This research program is being developed by the CIUDEN Foundation, an initiative launched by 3 Spanish state departments (Science & Innovation, Environment and Industry). An extensive multidisciplinary site characterization phase has been carried out, including a multiseismic data acquisition experiment. Within this effort, a 36 km2 3D seismic reflection survey was acquired in the summer of 2010. Its aim was to provide high resolution images of the subsurface of the storage complex, as well as to provide a baseline model for all the disciplines involved in the project. The target reservoir is a saline aquifer located at 1400 m, approximately, within Lower Jurassic carbonates (Lias). The main seal is formed by inter-layered marls and marly limestones of Early to Middle Jurassic age (Dogger and Lias). The main acquisition characteristics of the survey included (1) a mixed source of vibroseis and explosives with 74% and 26% of each used, respectively, (2) 5000 source points distributed along 22 source lines (separated 250 m) and (3) 22 lines of receivers (separated 275 m). Shot and receiver spacing along the source and receiver lines was 25 m, resulting in a nominal CDP-fold of 36 for 13 m2 bins. The 3D-data have been fully processed to post stack migration. The most critical processing steps included static correction calculations, time variant frequency filtering, rms velocity analysis, F-XY deconvolution, dip move-out correction, residual statics calculations and post stack migration. The final high-resolution 3D-volume shows the shape and depth of the primary reservoir-seal system, the main faults of the area and the secondary reservoir-seal sequence. It allows us to characterize the main tectonic structure of the dome complex, the fault system of the area and

  10. The internal geometry of salt structures - A first look using 3D seismic data from the Zechstein of the Netherlands

    Van Gent, Heijn; Urai, Janos L.; de Keijzer, Martin

    2011-03-01

    We present a first look at the large-scale, complexly folded and faulted internal structure of Zechstein salt bodies in NW Europe using 3D reflection seismic reflection data from two surveys on the Groningen High and the Cleaver Bank High. We focus on a relatively brittle, folded and boudinaged, claystone-carbonate-anhydrite layer (the Z3 stringer) enclosed in ductile salt. A first classification of the structures is presented and compared with observations from salt mines and analogue and numerical models. Z3 stringers not only are reservoirs for hydrocarbons but can also present a serious drilling problem in some areas. Results of this study could provide the basis for better prediction of zones of drilling problems. More generally, the techniques presented here can be used to predict the internal structure of salt bodies, to estimate the geometry of economic deposits of all kinds and locate zones suitable for storage caverns. Structures observed include an extensive network of zones with increased thickness of the stringer. These we infer to have formed by early diagenesis, karstification, gravitational sliding and associated local sedimentation. Later, this template was deformed into large-scale folds and boudins during salt tectonics. Salt flow was rarely plane strain, producing complex fold and boudin geometries. Deformation was further complicated by the stronger zones of increased thickness, which led to strongly non-cylindrical structures. We present some indications that the thicker zones also influence the locations of later suprasalt structures, suggesting a feedback between the early internal evolution of this salt giant and later salt tectonics. This study opens the possibility to study the internal structure of the Zechstein and other salt giants in 3D using this technique, exposing a previously poorly known structure which is comparable in size and complexity to the internal parts of some orogens.

  11. The role of 3-D interactive visualization in blind surveys of HI in galaxies

    Punzo, D; Roerdink, J B T M; Oosterloo, T A; Ramatsoku, M; Verheijen, M A W

    2015-01-01

    Upcoming HI surveys will deliver large datasets, and automated processing using the full 3-D information (two positional dimensions and one spectral dimension) to find and characterize HI objects is imperative. In this context, visualization is an essential tool for enabling qualitative and quantitative human control on an automated source finding and analysis pipeline. We discuss how Visual Analytics, the combination of automated data processing and human reasoning, creativity and intuition, supported by interactive visualization, enables flexible and fast interaction with the 3-D data, helping the astronomer to deal with the analysis of complex sources. 3-D visualization, coupled to modeling, provides additional capabilities helping the discovery and analysis of subtle structures in the 3-D domain. The requirements for a fully interactive visualization tool are: coupled 1-D/2-D/3-D visualization, quantitative and comparative capabilities, combined with supervised semi-automated analysis. Moreover, the sourc...

  12. Next-generation seismic experiments - II: wide-angle, multi-azimuth, 3-D, full-waveform inversion of sparse field data

    Morgan, Joanna; Warner, Michael; Arnoux, Gillean; Hooft, Emilie; Toomey, Douglas; VanderBeek, Brandon; Wilcock, William

    2016-02-01

    3-D full-waveform inversion (FWI) is an advanced seismic imaging technique that has been widely adopted by the oil and gas industry to obtain high-fidelity models of P-wave velocity that lead to improvements in migrated images of the reservoir. Most industrial applications of 3-D FWI model the acoustic wavefield, often account for the kinematic effect of anisotropy, and focus on matching the low-frequency component of the early arriving refractions that are most sensitive to P-wave velocity structure. Here, we have adopted the same approach in an application of 3-D acoustic, anisotropic FWI to an ocean-bottom-seismometer (OBS) field data set acquired across the Endeavour oceanic spreading centre in the northeastern Pacific. Starting models for P-wave velocity and anisotropy were obtained from traveltime tomography; during FWI, velocity is updated whereas anisotropy is kept fixed. We demonstrate that, for the Endeavour field data set, 3-D FWI is able to recover fine-scale velocity structure with a resolution that is 2-4 times better than conventional traveltime tomography. Quality assurance procedures have been employed to monitor each step of the workflow; these are time consuming but critical to the development of a successful inversion strategy. Finally, a suite of checkerboard tests has been performed which shows that the full potential resolution of FWI can be obtained if we acquire a 3-D survey with a slightly denser shot and receiver spacing than is usual for an academic experiment. We anticipate that this exciting development will encourage future seismic investigations of earth science targets that would benefit from the superior resolution offered by 3-D FWI.

  13. 3D SURVEYING AND MODELING OF ARCHAEOLOGICAL SITES – SOME CRITICAL ISSUES –

    Gonizzi Barsanti, S.; F. Remondino; Visintini, D

    2013-01-01

    The goal of the reported project is to test and evaluate 3D surveying and modelling methods to document the remaining ancient byzantine city walls of the archaeological site of Aquileia in Friuli Venezia Giulia, Italy. The objectives are threefold: (1) to use 3D data to create maps, façades and sections that provide information useful for archaeological purposes such as the investigation of architectural construction techniques or construction phases, (2) to evaluate and compare phot...

  14. Network level pavement evaluation with 1 mm 3D survey system

    Kelvin C.P. Wang

    2015-12-01

    Full Text Available The latest iteration of PaveVision3D Ultra can obtain true 1 mm resolution 3D data at full-lane coverage in all 3 directions at highway speed up to 60 mph. This paper introduces the PaveVision3D Ultra technology for rapid network level pavement survey on approximately 1280 center miles of Oklahoma interstate highways. With sophisticated automated distress analyzer (ADA software interface, the collected 1 mm 3D data provide Oklahoma Department of Transportation (ODOT with comprehensive solutions for automated evaluation of pavement surface including longitudinal profile for roughness, transverse profile for rutting, predicted hydroplaning speed for safety analysis, and cracking and various surface defects for distresses. The pruned exact linear time (PELT method, an optimal partitioning algorithm, is implemented to identify change points and dynamically determine homogeneous segments so as to assist ODOT effectively using the available 1 mm 3D pavement surface condition data for decision-making. The application of 1 mm 3D laser imaging technology for network survey is unprecedented. This innovative technology allows highway agencies to access its options in using the 1 mm 3D system for its design and management purposes, particularly to meet the data needs for pavement management system (PMS, pavement ME design and highway performance monitoring system (HPMS.

  15. Pseudo 3-D P wave refraction seismic monitoring of permafrost in steep unstable bedrock

    Krautblatter, Michael; Draebing, Daniel

    2014-02-01

    permafrost in steep rock walls can cause hazardous rock creep and rock slope failure. Spatial and temporal patterns of permafrost degradation that operate at the scale of instability are complex and poorly understood. For the first time, we used P wave seismic refraction tomography (SRT) to monitor the degradation of permafrost in steep rock walls. A 2.5-D survey with five 80 m long parallel transects was installed across an unstable steep NE-SW facing crestline in the Matter Valley, Switzerland. P wave velocity was calibrated in the laboratory for water-saturated low-porosity paragneiss samples between 20°C and -5°C and increases significantly along and perpendicular to the cleavage by 0.55-0.66 km/s (10-13%) and 2.4-2.7 km/s (>100%), respectively, when freezing. Seismic refraction is, thus, technically feasible to detect permafrost in low-porosity rocks that constitute steep rock walls. Ray densities up to 100 and more delimit the boundary between unfrozen and frozen bedrock and facilitate accurate active layer positioning. SRT shows monthly (August and September 2006) and annual active layer dynamics (August 2006 and 2007) and reveals a contiguous permafrost body below the NE face with annual changes of active layer depth from 2 to 10 m. Large ice-filled fractures, lateral onfreezing of glacierets, and a persistent snow cornice cause previously unreported permafrost patterns close to the surface and along the crestline which correspond to active seasonal rock displacements up to several mm/a. SRT provides a geometrically highly resolved subsurface monitoring of active layer dynamics in steep permafrost rocks at the scale of instability.

  16. Mapping the North Sea base-Quaternary: using 3D seismic to fill a gap in the geological record

    Lamb, Rachel; Huuse, Mads; Stewart, Margaret; Brocklehurst, Simon H.

    2014-05-01

    The identification and mapping of the base-Quaternary boundary in the central parts of the North Sea is problematic due to the change from an unconformable transition between Pliocene and Pleistocene deltaic deposits in the southern North Sea to a conformable one further north (Sejrup et al 1991; Gatliff et al 1994). The best estimates of the transition use seismic reflection data to identify a 'crenulated reflector' (Buckley 2012), or rely on correlating sparse biostratigraphy (Cameron et al 1987). Recent integration of biostratigraphy, pollen analysis, paleomagnetism and amino acid analysis in the Dutch and Danish sectors (Rasmussen et al 2005; Kuhlmann et al 2006) allows greater confidence in the correlation to a regional 3D seismic dataset and show that the base-Quaternary can be mapped across the entire basin. The base-Quaternary has been mapped using the PGS MegaSurvey dataset from wells in the Danish Sector along the initially unconformable horizon and down the delta front into the more conformable basin giving a high degree of confidence in the horizon pick. The mapped horizon is presented here alongside the difference between this new interpretation and the previously interpreted base-Quaternary (Buckley 2012). The revised base-Quaternary surface reaches a depth of 1248 ms TWT or approximately 1120 m (assuming average velocity of 1800 m/s) showing an elongate basin shape that follows the underlying structure of the Central Graben. The difference between the revised base-Quaternary and the traditional base-Quaternary reaches a maximum of over 600 ms TWT or approximately 540 m in the south-west with over 300 ms TWT or approximately 270 m at the Josephine well (56° 36.11'N, 2° 27.09'E) in the centre of the basin. Mapping this new base-Quaternary allows for the interpretation of the paleo-envionrment during the earliest Quaternary. Seismic attribute analysis indicates a deep water basin with sediment deposition from multiple deltas and redistribution by deep

  17. The role of 3-D interactive visualization in blind surveys of H I in galaxies

    Punzo, D.; van der Hulst, J. M.; Roerdink, J. B. T. M.; Oosterloo, T. A.; Ramatsoku, M.; Verheijen, M. A. W.

    2015-09-01

    Upcoming H I surveys will deliver large datasets, and automated processing using the full 3-D information (two positional dimensions and one spectral dimension) to find and characterize H I objects is imperative. In this context, visualization is an essential tool for enabling qualitative and quantitative human control on an automated source finding and analysis pipeline. We discuss how Visual Analytics, the combination of automated data processing and human reasoning, creativity and intuition, supported by interactive visualization, enables flexible and fast interaction with the 3-D data, helping the astronomer to deal with the analysis of complex sources. 3-D visualization, coupled to modeling, provides additional capabilities helping the discovery and analysis of subtle structures in the 3-D domain. The requirements for a fully interactive visualization tool are: coupled 1-D/2-D/3-D visualization, quantitative and comparative capabilities, combined with supervised semi-automated analysis. Moreover, the source code must have the following characteristics for enabling collaborative work: open, modular, well documented, and well maintained. We review four state of-the-art, 3-D visualization packages assessing their capabilities and feasibility for use in the case of 3-D astronomical data.

  18. High-resolution 3D seismic data characterize fluid flow systems in the SW Barents Sea

    Bünz, Stefan; Mienert, Jürgen; Rajan, Anupama

    2010-05-01

    The flow of fluids through marine sediments is one of the most dominant and pervasive processes in continental margins. These processes control the evolution of a sedimentary basin and its seafloor environment, and have implications for hydrocarbon exploration and seabed ecosystems. Many seep sites at the seafloor are associated with large but complex faunal communities that have received significant attention in recent years. However, there is a need for a better understanding of the driving mechanism of fluid flow in various geological settings, the accumulation of fluids in the subsurface and their focused flow through conduits and/or faults to the seabed. The Barents Sea is a large hydrocarbon-prone basin of the Norwegian Arctic region. A significant portion of the hydrocarbons has leaked or migrated into the shallow subsurface and is now trapped in gas-hydrate and shallow-gas reservoirs. Furthermore, there are few places in the Barents Sea, where methane gas is leaking from the seafloor into the oceanosphere. Accumulations of free gas in the shallow subsurface are considered a geohazard. They constitute a risk for safe drilling operations and they may pose a threat to global climate if the seal that is trapping them is breached. P-Cable 3D high-resolution seismic data from the Ringvassøya Fault Complex and the Polheim Sub-Platform provide new and detailed insight into fluid flow controls and accumulation mechanisms. The data shows a wide variety of fluid flow features, mostly in the form of pockmarks, bright spots, wipe-out zones or vertical zones of disturbed reflectivity. Fluids migrate by both diapiric mechanism and channelized along sedimentary layers. Glacigenic sediments generally form a strong boundary for fluid flow in the very shallow section. However, we can recognize pockmarks not only at the seafloor but also at one subsurface layer approximately 50 m below sea floor indicating a former venting period in the SW Barents Sea. At few locations high

  19. Characterization of gas hydrate distribution using conventional 3D seismic data in the Pearl River Mouth Basin, South China Sea

    Wang, Xiujuan; Qiang, Jin; Collett, Timothy S.; Shi, Hesheng; Yang, Shengxiong; Yan, Chengzhi; Li, Yuanping; Wang, Zhenzhen; Chen, Duanxin

    2016-01-01

    A new 3D seismic reflection data volume acquired in 2012 has allowed for the detailed mapping and characterization of gas hydrate distribution in the Pearl River Mouth Basin in the South China Sea. Previous studies of core and logging data showed that gas hydrate occurrence at high concentrations is controlled by the presence of relatively coarse-grained sediment and the upward migration of thermogenic gas from the deeper sediment section into the overlying gas hydrate stability zone (BGHSZ); however, the spatial distribution of the gas hydrate remains poorly defined. We used a constrained sparse spike inversion technique to generate acoustic-impedance images of the hydrate-bearing sedimentary section from the newly acquired 3D seismic data volume. High-amplitude reflections just above the bottom-simulating reflectors (BSRs) were interpreted to be associated with the accumulation of gas hydrate with elevated saturations. Enhanced seismic reflections below the BSRs were interpreted to indicate the presence of free gas. The base of the BGHSZ was established using the occurrence of BSRs. In areas absent of well-developed BSRs, the BGHSZ was calculated from a model using the inverted P-wave velocity and subsurface temperature data. Seismic attributes were also extracted along the BGHSZ that indicate variations reservoir properties and inferred hydrocarbon accumulations at each site. Gas hydrate saturations estimated from the inversion of acoustic impedance of conventional 3D seismic data, along with well-log-derived rock-physics models were also used to estimate gas hydrate saturations. Our analysis determined that the gas hydrate petroleum system varies significantly across the Pearl River Mouth Basin and that variability in sedimentary properties as a product of depositional processes and the upward migration of gas from deeper thermogenic sources control the distribution of gas hydrates in this basin.

  20. 3D geology in a 2D country : Perspectives for geological surveying in the Netherlands

    Meulen, M.J. van der; Doornenbal, J.C.; Gunnink, J.L.; Stafleu, J.; Schokker, J.; Vernes, R.W.; Geer, F.C. van; Gessel, S.F. van; Heteren, S. van; Leeuwen, R.J.W. van; Bakker, M.A.J.; Bogaard, P.J.F.; Busschers, F.S.; Griffioen, J.; Gruijters, S.H.L.L.; Kiden, P.; Schroot, B.M.; Simmelink, H.J.; Berkel, W.O. van; Krogt, R.A.A. van der; Westerhoff, W.E.; Daalen, T.M. van

    2013-01-01

    Over the last ten to twenty years, geological surveys all over the world have been entangled in a process of digitisation. Their paper archives, built over many decades, have largely been replaced by electronic databases. The systematic production of geological map sheets is being replaced by 3D sub

  1. Towards a 3d Based Platform for Cultural Heritage Site Survey and Virtual Exploration

    Seinturier, J.; Riedinger, C.; Mahiddine, A.; Peloso, D.; Boï, J.-M.; Merad, D.; Drap, P.

    2013-07-01

    This paper present a 3D platform that enables to make both cultural heritage site survey and its virtual exploration. It provides a single and easy way to use framework for merging multi scaled 3D measurements based on photogrammetry, documentation produced by experts and the knowledge of involved domains leaving the experts able to extract and choose the relevant information to produce the final survey. Taking into account the interpretation of the real world during the process of archaeological surveys is in fact the main goal of a survey. New advances in photogrammetry and the capability to produce dense 3D point clouds do not solve the problem of surveys. New opportunities for 3D representation are now available and we must to use them and find new ways to link geometry and knowledge. The new platform is able to efficiently manage and process large 3D data (points set, meshes) thanks to the implementation of space partition methods coming from the state of the art such as octrees and kd-trees and thus can interact with dense point clouds (thousands to millions of points) in real time. The semantisation of raw 3D data relies on geometric algorithms such as geodetic path computation, surface extraction from dense points cloud and geometrical primitive optimization. The platform provide an interface that enables expert to describe geometric representations of interesting objects like ashlar blocs, stratigraphic units or generic items (contour, lines, … ) directly onto the 3D representation of the site and without explicit links to underlying algorithms. The platform provide two ways for describing geometric representation. If oriented photographs are available, the expert can draw geometry on a photograph and the system computes its 3D representation by projection on the underlying mesh or the points cloud. If photographs are not available or if the expert wants to only use the 3D representation then he can simply draw objects shape on it. When 3D

  2. tomo3d: a new 3-D joint refraction and reflection travel-time tomography code for active-source seismic data

    Meléndez, A.; Korenaga, J.; Sallares, V.; Ranero, C. R.

    2012-12-01

    We present the development state of tomo3d, a code for three-dimensional refraction and reflection travel-time tomography of wide-angle seismic data based on the previous two-dimensional version of the code, tomo2d. The core of both forward and inverse problems is inherited from the 2-D version. The ray tracing is performed by a hybrid method combining the graph and bending methods. The graph method finds an ordered array of discrete model nodes, which satisfies Fermat's principle, that is, whose corresponding travel time is a global minimum within the space of discrete nodal connections. The bending method is then applied to produce a more accurate ray path by using the nodes as support points for an interpolation with beta-splines. Travel time tomography is formulated as an iterative linearized inversion, and each step is solved using an LSQR algorithm. In order to avoid the singularity of the sensitivity kernel and to reduce the instability of inversion, regularization parameters are introduced in the inversion in the form of smoothing and damping constraints. Velocity models are built as 3-D meshes, and velocity values at intermediate locations are obtained by trilinear interpolation within the corresponding pseudo-cubic cell. Meshes are sheared to account for topographic relief. A floating reflector is represented by a 2-D grid, and depths at intermediate locations are calculated by bilinear interpolation within the corresponding square cell. The trade-off between the resolution of the final model and the associated computational cost is controlled by the relation between the selected forward star for the graph method (i.e. the number of nodes that each node considers as its neighbors) and the refinement of the velocity mesh. Including reflected phases is advantageous because it provides a better coverage and allows us to define the geometry of those geological interfaces with velocity contrasts sharp enough to be observed on record sections. The code also

  3. A university-developed seismic source for shallow seismic surveys

    Yordkayhun, Sawasdee; Na Suwan, Jumras

    2012-07-01

    The main objectives of this study were to (1) design and develop a low cost seismic source for shallow seismic surveys and (2) test the performance of the developed source at a test site. The surface seismic source, referred to here as a university-developed seismic source is based upon the principle of an accelerated weight drop. A 30 kg activated mass is lifted by a mechanical rack and pinion gear and is accelerated by a mounted spring. When the mass is released from 0.5 m above the surface, it hits a 30 kg base plate and energy is transferred to the ground, generating a seismic wave. The developed source is portable, environmentally friendly, easy to operate and maintain, and is a highly repeatable impact source. To compare the developed source with a sledgehammer source, a source test was performed at a test site, a study site for mapping a major fault zone in southern Thailand. The sledgehammer and the developed sources were shot along a 300 m long seismic reflection profile with the same parameters. Data were recorded using 12 channels off-end geometry with source and receiver spacing of 5 m, resulting in CDP stacked sections with 2.5 m between traces. Source performances were evaluated based on analyses of signal penetration, frequency content and repeatability, as well as the comparison of stacked sections. The results show that both surface sources are suitable for seismic studies down to a depth of about 200 m at the site. The hammer data are characterized by relatively higher frequency signals than the developed source data, whereas the developed source generates signals with overall higher signal energy transmission and greater signal penetration. In addition, the repeatability of the developed source is considerably higher than the hammer source.

  4. 2D and 3D numerical modeling of seismic waves from explosion sources

    Over the last decade, nonlinear and linear 2D axisymmetric finite difference codes have been used in conjunction with far-field seismic Green's functions to simulate seismic waves from a variety of sources. In this paper we briefly review some of the results and conclusions that have resulted from numerical simulations and explosion modeling in support of treaty verification research at S-CUBED in the last decade. We then describe in more detail the results from two recent projects. Our goal is to provide a flavor for the kinds of problems that can be examined with numerical methods for modeling excitation of seismic waves from explosions. Two classes of problems have been addressed; nonlinear and linear near-source interactions. In both classes of problems displacements and tractions are saved on a closed surface in the linear region and the representation theorem is used to propagate the seismic waves to the far-field

  5. Seismic ground motion amplification in a 3D sedimentary basin: the effect of the vertical velocity gradient

    Ground motion amplification in sedimentary basins has been observed in some moderate or large earthquakes, such as the 1994 Northridge and 1999 Chi-Chi event. Many numerical studies with simplified 2D models have shown significant effects of the vertical velocity gradient of sediment on basin amplification. However, we need to consider a more realistic 3D model and solve wave equations with 3D numerical methods in order to improve our understanding of basin amplification. In this study, we extend a 2D pseudospectral and finite difference hybrid method to a 3D case and investigate the effects of the vertical velocity gradient for a 3D basin model. Numerical simulations were performed for four basin models with increasing vertical velocity gradients on a PC cluster using 64 processors for 67 108 864 discretized grids. The results show that the vertical velocity gradient enhances basin amplification through strong secondary surface waves and basin trapped waves. The 3D geometry of the basin causes a wave-front focusing effect that contributes significantly to a localized strong amplification with the maximum peak ground velocity in the basin. The results of this study suggest that it is important to consider the detailed properties of sedimentary basins in seismic ground motion studies. (paper)

  6. Reducing disk storage of full-3D seismic waveform tomography (F3DT) through lossy online compression

    Lindstrom, Peter; Chen, Po; Lee, En-Jui

    2016-08-01

    Full-3D seismic waveform tomography (F3DT) is the latest seismic tomography technique that can assimilate broadband, multi-component seismic waveform observations into high-resolution 3D subsurface seismic structure models. The main drawback in the current F3DT implementation, in particular the scattering-integral implementation (F3DT-SI), is the high disk storage cost and the associated I/O overhead of archiving the 4D space-time wavefields of the receiver- or source-side strain tensors. The strain tensor fields are needed for computing the data sensitivity kernels, which are used for constructing the Jacobian matrix in the Gauss-Newton optimization algorithm. In this study, we have successfully integrated a lossy compression algorithm into our F3DT-SI workflow to significantly reduce the disk space for storing the strain tensor fields. The compressor supports a user-specified tolerance for bounding the error, and can be integrated into our finite-difference wave-propagation simulation code used for computing the strain fields. The decompressor can be integrated into the kernel calculation code that reads the strain fields from the disk and compute the data sensitivity kernels. During the wave-propagation simulations, we compress the strain fields before writing them to the disk. To compute the data sensitivity kernels, we read the compressed strain fields from the disk and decompress them before using them in kernel calculations. Experiments using a realistic dataset in our California statewide F3DT project have shown that we can reduce the strain-field disk storage by at least an order of magnitude with acceptable loss, and also improve the overall I/O performance of the entire F3DT-SI workflow significantly. The integration of the lossy online compressor may potentially open up the possibilities of the wide adoption of F3DT-SI in routine seismic tomography practices in the near future.

  7. High-Resolution 3D-Seismic Investigations Indicate Focused Fluid Flow Systems in Hydrated Sediments at the Vestnesa Ridge off the W-Svalbard Margin.

    Petersen, C.; Buenz, S.; Hustoft, S.; Mienert, J.

    2007-12-01

    High-resolution seismic data were acquired using the 3D seismic P-Cable system of the University of Tromsoe to investigate how the fluid flow penetrates gas hydrate systems of the Vestnesa Ridge. The ridge represents a current-controlled sediment drift on the continental margin offshore western Svalbard. The survey area is located at the northwestern part of the Vestnesa Ridge and centered at the ridge crest that resembles an anticline in a water depth of 1250-1320 m. The seafloor morphology at the crest is characterized by an abundance of pockmarks with a diameter between 50-500 m indicating recent fluid flow activity. Since the area is within the gas hydrate stability zone (GHSZ), it is an ideal site to understand where and how fluids escape through a hydrated sediment drift. 35 reflection seismic profiles with a spacing of about 40-60 m were shot resulting in a seismic cube covering an area of approximately 22 km2. In addition, regional single channel streamer (SCS) seismic lines were acquired across the ridge perpendicular to the crest to connect the 3D area with the regional structural setting. The seismic data provide images of the subsurface to about 500 ms TWT (two-way time) below the seafloor (bsf), where gas accumulations cause acoustic attenuations that hinder deeper acoustic signal penetration. The well-stratified sediments exhibit a bottom simulating reflector (BSR) at about 200 ms TWT bsf at the base of the GHSZ. The BSR is difficult to identify due to the stratification, but it is accompanied by the onset of an ubiquitous band of strong reflectivity indicating free gas accumulation zones beneath the GHSZ. Fluid flow activity is evident from a link between gas accumulations (bright spots), gas wipeouts and disturbed reflectivity in the seismic data. These features are observed not only beneath the pockmark structures, but also in the sediment without seafloor expressions of fluid venting. The fluid source might be related to deep tectonic processes at

  8. 3D frequency modeling of elastic seismic wave propagation via a structured massively parallel direct Helmholtz solver

    Wang, S.; De Hoop, M. V.; Xia, J.; Li, X.

    2011-12-01

    We consider the modeling of elastic seismic wave propagation on a rectangular domain via the discretization and solution of the inhomogeneous coupled Helmholtz equation in 3D, by exploiting a parallel multifrontal sparse direct solver equipped with Hierarchically Semi-Separable (HSS) structure to reduce the computational complexity and storage. In particular, we are concerned with solving this equation on a large domain, for a large number of different forcing terms in the context of seismic problems in general, and modeling in particular. We resort to a parsimonious mixed grid finite differences scheme for discretizing the Helmholtz operator and Perfect Matched Layer boundaries, resulting in a non-Hermitian matrix. We make use of a nested dissection based domain decomposition, and introduce an approximate direct solver by developing a parallel HSS matrix compression, factorization, and solution approach. We cast our massive parallelization in the framework of the multifrontal method. The assembly tree is partitioned into local trees and a global tree. The local trees are eliminated independently in each processor, while the global tree is eliminated through massive communication. The solver for the inhomogeneous equation is a parallel hybrid between multifrontal and HSS structure. The computational complexity associated with the factorization is almost linear with the size of the Helmholtz matrix. Our numerical approach can be compared with the spectral element method in 3D seismic applications.

  9. High-resolution 3D seismic imaging of a pull-apart basin in the Gulf of Cadiz

    Crutchley, G.; Berndt, C.; Klaeschen, D.; Gutscher, M.

    2009-12-01

    In 2006, high-resolution 3D seismic data were acquired in the Gulf of Cadiz and the Mediterranean Sea aboard the RRS Charles Darwin as part of the HERMES (Hotspot Ecosystem Research on the Margins of European Seas) project. The P-Cable system, a cost-efficient set-up for fast acquisition of 3D seismic data on 12 single-channel streamers, was utilized to acquire seismic cubes at four different targets. Here, we present results from the second target - a WNW-ESE-oriented pull-apart basin in the southeastern Gulf of Cadiz. Initial processing has included: 1) spatial positioning of each recording channel from GPS data acquired on the outer two channels, 2) improved positioning of shot points and channels from the inversion of first arrival times, 3) application of a swell filter to improve reflection coherency, 4) CDP binning and stacking and 5) migration. The new data confirm that the southeastern Gulf of Cadiz north of the Rharb submarine valley is structurally controlled by numerous strike slip faults that were active until quite recently (within the resolution of the data). Given the location of this basin, between the extensional domain on the upper slope and the compressional toe of the accretionary wedge, we interpret the origin to be gravitational sliding on a detachment layer, possibly containing salt, but at this stage not imaged by our profiles.

  10. PLANETARY NEBULAE DETECTED IN THE SPITZER SPACE TELESCOPE GLIMPSE 3D LEGACY SURVEY

    We used the data from the Spitzer Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) to investigate the mid-infrared (MIR) properties of planetary nebulae (PNs) and PN candidates. In previous studies of GLIMPSE I and II data, we have shown that these MIR data are very useful in distinguishing PNs from other emission-line objects. In the present paper, we focus on the PNs in the field of the GLIMPSE 3D survey, which has a more extensive latitude coverage. We found a total of 90 Macquarie-AAO-Strasbourg (MASH) and MASH II PNs and 101 known PNs to have visible MIR counterparts in the GLIMPSE 3D survey area. The images and photometry of these PNs are presented. Combining the derived IRAC photometry at 3.6, 4.5, 5.8, and 8.0 μm with the existing photometric measurements from other infrared catalogs, we are able to construct spectral energy distributions (SEDs) of these PNs. Among the most notable objects in this survey is the PN M1-41, whose GLIMPSE 3D image reveals a large bipolar structure of more than 3 arcmin in extent.

  11. Record of Subducting Topography revealed in 3D Seismic Imaging of Pleistocene unconformities, offshore Southern Costa Rica

    Edwards, J. H.; Kluesner, J. W.; Silver, E. A.

    2015-12-01

    3D seismic reflection data (CRISP) collected across the southern Costa Rica forearc reveals broad, survey-wide erosional events in the upper ~1 km of slope sediments in the mid-slope to outer shelf. The upper 0-280 m of continuous, weakly deformed sediments, designated by IODP Expedition 344 as structural domain I, is bounded by a major erosional event, (CRISP-U1, dated near 1 Ma), suggesting wave-plain erosion from the present shelf break out to 25 km seaward, to a present-day water depth of 900-1300 m. The eastern toe of its surface is characterized by a large drainage system, likely including submarine channels that eroded to depths >1500 m below present-day water depth. CRISP-U1 is variably uplifted by a series of fault propagation folds and cut by an intersecting array of normal faults. Another, major erosional event, (CRISP-M1, approximately 2 Ma) extended from the outer shelf to the mid slope and removed 500-1000 m of material. Overlying CRISP-M1 is up to 1 km of sediments that are more deformed by fault propagation folds, back thrusts, and intersecting arrays of normal faults. Unconformities with smaller areal extent are variably found in these overlying sediments across the mid-slope to outer shelf, at present-day water depths >220 m. Below CRISP-M1, sediments are more densely deformed and also contain major unconformities that extend survey-wide. Both unconformities, CRISP-U1 and CRISP-M1, are encountered in well U1413 and are demarcated by major benthic foraminifera assemblage changes at 149 mbsf and ~504 mbsf (Harris et al., 2013, Proceeding of the IODP, Volume 344).CRISP-M1 is likely correlative to the major sediment facies and benthic foraminifera assemblage change found in U1379 at ~880 mbsf (Vannuchi et al., 2013). The unconformities and intersecting array of normal faults may demarcate the passing of topography on the downgoing Cocos plate, episodically lifting and then subsiding the Costa Rica margin, with amplitudes up to about 1 km.

  12. Automated fault extraction and classification using 3-D seismic data for the Ekofisk field development

    Signer, C.; Nickel, M.; Randen, T.; Saeter, T.; Soenneland, H.H.

    1998-12-31

    Mapping of fractures is important for the prediction of fluid flow in many reservoir types. The fluid flow depends mainly on the efficiency of the reservoir seals. Improved spatial mapping of the open and closed fracture systems will allow a better prediction of the fluid flow pattern. The primary objectives of this paper is to present fracture characterization at the reservoir scale combined with seismic facies mapping. The complexity of the giant Ekofisk field on the Norwegian continental shelf provides an ideal framework for testing the validity and the applicability of an automated seismic fault and fracture detection and mapping tool. The mapping of the faults can be based on seismic attribute grids, which means that attribute-responses related to faults are extracted along key horizons which were interpreted in the reservoir interval. 3 refs., 3 figs.

  13. Vertical Cable Seismic Survey for SMS Exploration in Izena Cauldron, Okinawa-Trough

    Asakawa, E.; Murakami, F.; Tsukahara, H.; Mizohata, S.; Tara, K.

    2014-12-01

    The VCS survey is one of the reflection seismic methods. It uses hydrophone arrays vertically moored from the seafloor to record acoustic waves generated by seismic sources. Because the VCS is an efficient high-resolution 3D seismic survey method for a spatially-bounded area, we proposed it for the SMS survey tool development program started by Japanese government. In 2010, we manufactured the autonomous VCS data acquisition systems. Through several experimental surveys, our VCS is successfully completed. In 2011 and 2013, we carried out the two VCS surveys using GI gun and high-voltage sparker respectively in the Izena Cauldron, Okinawa Trough, which is one of the most promising SMS areas around Japan. Because seismic survey is not proven to be effective for SMS exploration, no seismic surveys have been conducted there so far. Our strategy for SMS exploration consists of two stages. In the first stage, we carried out VCS survey with the lower frequency GI gun (but higher compared to the convebtional oil/gas exploration) and explored deeper (up to 1,500m) structure to obtain the fault system of hydrothermal flow. Next, using a high frequency (about 1 kHz higher) and high-voltage sparker, we explored very shallow (up to 200m) part to delineate the very thin SMS deposits. These two VCS dataset have been processed with 3D Prestack Depth Migration. These results are consistent with geological information from the borehole drilled nearby and give useful information to SMS exploration.

  14. 3-D seismic acquisition geometry design and analysis: Investigation of the requirements to include illumination from all multiples

    Kumar, A.

    2015-01-01

    A seismic survey should be designed such that imaging of the acquired data leads to a sufficiently accurate subsurface image. For that purpose, methods for acquisition geometry analysis and design are available. These methods are used to judge whether an acquisition geometry is suited for the specif

  15. 2D and 3D imaging of the metamorphic carbonates at Omalos plateau/polje, Crete, Greece by employing independent and joint inversion on resistivity and seismic data

    Pangratis Pangratis

    2012-07-01

    Full Text Available A geophysical survey carried out at Omalos plateau in Chania, Western Crete, Greece employed seismic as well as electrical tomography methods in order to image karstic structures and the metamorphic carbonates (Tripali unit and Plattenkalk group which are covered by post-Mesozoic deposits (terra rossa, clays, sands and gravels. The geoelectrical sections image the metamorphic carbonates which exhibit a highly irregular relief. At the central part of the plateau the thickness of post-Mesozoic deposits (terra rossa, clays, sands and gravels ranges from 40-130 m. A 3D resistivity image was generated by inverting resistivity data collected on a grid to the south west at the Omalos plateau. The 3D resistivity image delineated a karstic structure at a depth of 25 to 55 m. On the same grid the depth to the top of the karstified carbonates ranges from 25-70 m. This is also verified on the resistivity sections and seismic velocity sections along lines 5 and 7 of the above mentioned grid which are derived from the cross-gradients joint inversion.

  16. Using 3D Simulation of Elastic Wave Propagation in Laplace Domain for Electromagnetic-Seismic Inverse Modeling

    Petrov, P.; Newman, G. A.

    2010-12-01

    -Fourier domain we had developed 3D code for full-wave field simulation in the elastic media which take into account nonlinearity introduced by free-surface effects. Our approach is based on the velocity-stress formulation. In the contrast to conventional formulation we defined the material properties such as density and Lame constants not at nodal points but within cells. This second order finite differences method formulated in the cell-based grid, generate numerical solutions compatible with analytical ones within the range errors determinate by dispersion analysis. Our simulator will be embedded in an inversion scheme for joint seismic- electromagnetic imaging. It also offers possibilities for preconditioning the seismic wave propagation problems in the frequency domain. References. Shin, C. & Cha, Y. (2009), Waveform inversion in the Laplace-Fourier domain, Geophys. J. Int. 177(3), 1067- 1079. Shin, C. & Cha, Y. H. (2008), Waveform inversion in the Laplace domain, Geophys. J. Int. 173(3), 922-931. Commer, M. & Newman, G. (2008), New advances in three-dimensional controlled-source electromagnetic inversion, Geophys. J. Int. 172(2), 513-535. Newman, G. A., Commer, M. & Carazzone, J. J. (2010), Imaging CSEM data in the presence of electrical anisotropy, Geophysics, in press.

  17. Seismic response of free standing fuel rack construction to 3-D floor motion

    Seismic analysis of free standing submerged racks is complicated by the presence of water and structural non-linearities such as fuel assembly cell impact and floor interface friction. A direct time integration technique has been proposed to analyze this class of structures. Application of the time integration technique on a fourteen degree of freedom lumped mass model of the rack reveals some heretofore unpublished quirks in the structure's behavior. The method of analysis is utilized to compare the seismic response of some representative rack designs. Results show wide differences in the structural response, depending on the fabrication details of racks

  18. Assessing earthquake source models using 3-D forward modelling of long-period seismic data: application to the SCARDEC method

    Ferreira, Ana Mg; Vallee, Martin; Charlety, Jean

    2010-05-01

    Accurate earthquake point source parameters (e.g. seismic moment, depth and focal mechanism) provide key first-order information for detailed studies of the earthquake source process and for improved seismic and tsunami hazard evaluation. In order to objectively assess the quality of seismic source models, it is important to go beyond classical resolution checks. In particular, it is desirable to apply sophisticated modelling techniques to quantify inaccuracies due to simplified theoretical formulations and/or Earth structure employed to build the source models. Moreover, it is important to verify how well the models explain data not used in their construction. In this study we assess the quality of the SCARDEC method (see joint abstracts), which is a new automated technique that retrieves simultaneously the seismic moment, focal mechanism, depth and source time functions of large earthquakes. Because the SCARDEC method is based on body-wave deconvolution using ray methods in a 1D Earth model, we test how well SCARDEC source parameters explain long-period seismic data (surface waves and normal modes). We calculate theoretical seismograms using two forward modelling techniques (full ray theory and spectral element method) to simulate the long-period seismic wavefield for the 3D Earth model S20RTS combined with the crust model CRUST2.0, and for two point source models: (i) the SCARDEC model; and (ii) the Global CMT model. We compare the synthetic seismograms with real broadband data from the FDSN for the major subduction earthquakes of the last 20 years. We show that SCARDEC source parameters explain long-period surface waves as well as Global CMT solutions. This can be explained by the fact that most of the differences between SCARDEC and Global CMT solutions are linked to correlated variations of the seismic moment and dip of the earthquakes, and it is theoretically known that for shallow earthquakes it is difficult to accurately resolve these two parameters using

  19. Modeling seismic wave propagation and amplification in 1D/2D/3D linear and nonlinear unbounded media

    Semblat, Jean-François

    2011-01-01

    To analyze seismic wave propagation in geological structures, it is possible to consider various numerical approaches: the finite difference method, the spectral element method, the boundary element method, the finite element method, the finite volume method, etc. All these methods have various advantages and drawbacks. The amplification of seismic waves in surface soil layers is mainly due to the velocity contrast between these layers and, possibly, to topographic effects around crests and hills. The influence of the geometry of alluvial basins on the amplification process is also know to be large. Nevertheless, strong heterogeneities and complex geometries are not easy to take into account with all numerical methods. 2D/3D models are needed in many situations and the efficiency/accuracy of the numerical methods in such cases is in question. Furthermore, the radiation conditions at infinity are not easy to handle with finite differences or finite/spectral elements whereas it is explicitely accounted in the B...

  20. Analysis of allochthonous salt and salt welds in the northern Gulf of Mexico utilizing 3D seismic data

    Pritchett, J.A.; House, W.M. [Amoco Production Co., Houston, TX (United States)

    1996-12-31

    The widespread availability of 3D seismic data sets in the northern Gulf of Mexico has resulted in the emergence of new play concepts such as subsalt exploration. Time and depth migrations of these 3D data allow interpreters to develop a detailed understanding of the geological processes that contribute to the structural and stratigraphic framework of the Gulf. These data provide excellent imaging of structural features, and result in the correct spatial positioning of those structural elements. Analysis of the geometrical relationships between allochthonous salt, salt welds and subsalt reflectors aids in the development of salt emplacement models. These models are subsequently tied to other elements of the hydrocarbon system such as fluid migration and reservoir development. Salt sheets and horizontal salt welds often separate distinct structural domains in the supra salt and subsalt section, and complex structural deformation above salt or a salt weld may not translate into the subsalt section.

  1. Analysis of allochthonous salt and salt welds in the northern Gulf of Mexico utilizing 3D seismic data

    Pritchett, J.A.; House, W.M. (Amoco Production Co., Houston, TX (United States))

    1996-01-01

    The widespread availability of 3D seismic data sets in the northern Gulf of Mexico has resulted in the emergence of new play concepts such as subsalt exploration. Time and depth migrations of these 3D data allow interpreters to develop a detailed understanding of the geological processes that contribute to the structural and stratigraphic framework of the Gulf. These data provide excellent imaging of structural features, and result in the correct spatial positioning of those structural elements. Analysis of the geometrical relationships between allochthonous salt, salt welds and subsalt reflectors aids in the development of salt emplacement models. These models are subsequently tied to other elements of the hydrocarbon system such as fluid migration and reservoir development. Salt sheets and horizontal salt welds often separate distinct structural domains in the supra salt and subsalt section, and complex structural deformation above salt or a salt weld may not translate into the subsalt section.

  2. 3D seismic modeling and reverse‐time migration with the parallel Fourier method using non‐blocking collective communications

    Chu, Chunlei

    2009-01-01

    The major performance bottleneck of the parallel Fourier method on distributed memory systems is the network communication cost. In this study, we investigate the potential of using non‐blocking all‐to‐all communications to solve this problem by overlapping computation and communication. We present the runtime comparison of a 3D seismic modeling problem with the Fourier method using non‐blocking and blocking calls, respectively, on a Linux cluster. The data demonstrate that a performance improvement of up to 40% can be achieved by simply changing blocking all‐to‐all communication calls to non‐blocking ones to introduce the overlapping capability. A 3D reverse‐time migration result is also presented as an extension to the modeling work based on non‐blocking collective communications.

  3. The role of 3-D interactive visualization in blind surveys of HI in galaxies

    Punzo, D.; van der Hulst, J. M.; Roerdink, J. B. T. M.; Oosterloo, T. A.; Ramatsoku, M.; Verheijen, M. A. W.

    2015-01-01

    Upcoming HI surveys will deliver large datasets, and automated processing using the full 3-D information (two positional dimensions and one spectral dimension) to find and characterize HI objects is imperative. In this context, visualization is an essential tool for enabling qualitative and quantitative human control on an automated source finding and analysis pipeline. We discuss how Visual Analytics, the combination of automated data processing and human reasoning, creativity and intuition,...

  4. 3D Reconstruction of Coronary Arteries from Angiographic Images: A Survey

    Mahima Goyal; Jian Yang; Arun PrakashAgrawal

    2013-01-01

    X-Ray angiographyis considered to be the “golden standard” ofall times in the medical imaging field due to its wide range ofapplications. In this paper, a survey is performed on the basictechniques and the methodologies that have already beenproposed by various researchers in the field of 3D reconstructionof coronary arteries using one or more angiograms as inmonoplane or biplane angiographic systems or a volume of dataas in rotational angiographic systems. Also, the variousprocedures that ne...

  5. Incorporating higher order WINKLER springs with 3-D finite element model of a reactor building for seismic SSI analysis

    In order to fulfill the seismic safety requirements, in the frame of seismic requalification activities for NPP Muehleberg, Switzerland, detailed seismic analysis performed on the Reactor Building and the results are presented previously. The primary objective of the present investigation is to assess the seismic safety of the reinforced concrete structures of reactor building. To achieve this objective requires a rather detailed 3-D finite element modeling for the outer shell structures, the drywell, the reactor pools, the floor decks and finally, the basemat. This already is a complicated task, which enforces need for simplifications in modelling the reactor internals and the foundation soil. Accordingly, all internal parts are modelled by vertical sticks and the Soil Structure Interaction (SSI) effects are represented by sets of transitional and higher order rotational WINKLER springs, i.e. avoiding complicated finite element SSI analysis. As a matter of fact, the availability of the results of recent investigations carried out on the reactor building using diversive finite element SSI analysis methods allow to calibrate the WINKLER springs, ensuring that the overall SSI behaviour of the reactor building is maintained

  6. Mantle fabric of western Bohemian Massif (central Europe) constrained by 3D seismic P and S anisotropy

    Babuška, Vladislav; Plomerová, Jaroslava; Vecsey, Luděk

    2008-01-01

    Roč. 462, č. 1-4 (2008), s. 149-163. ISSN 0040-1951 R&D Projects: GA ČR GA205/04/0748; GA ČR GA205/07/1088; GA AV ČR(CZ) KJB300120605 Institutional research plan: CEZ:AV0Z30120515 Keywords : Bohemian Massif * P and S seismic anisotropy * 3D fabric of mantle lithosphere Subject RIV: DC - Siesmology, Volcanology, Earth Structure Impact factor: 1.677, year: 2008

  7. A Robust MEMS Based Multi-Component Sensor for 3D Borehole Seismic Arrays

    Paulsson Geophysical Services

    2008-03-31

    The objective of this project was to develop, prototype and test a robust multi-component sensor that combines both Fiber Optic and MEMS technology for use in a borehole seismic array. The use such FOMEMS based sensors allows a dramatic increase in the number of sensors that can be deployed simultaneously in a borehole seismic array. Therefore, denser sampling of the seismic wave field can be afforded, which in turn allows us to efficiently and adequately sample P-wave as well as S-wave for high-resolution imaging purposes. Design, packaging and integration of the multi-component sensors and deployment system will target maximum operating temperature of 350-400 F and a maximum pressure of 15000-25000 psi, thus allowing operation under conditions encountered in deep gas reservoirs. This project aimed at using existing pieces of deployment technology as well as MEMS and fiber-optic technology. A sensor design and analysis study has been carried out and a laboratory prototype of an interrogator for a robust borehole seismic array system has been assembled and validated.

  8. 3D Subsoil Model of the San Biagio `Salinelle' Mud Volcanoes (Belpasso, Sicily) derived from Geophysical Surveys

    Imposa, S.; Grassi, S.; De Guidi, G.; Battaglia, F.; Lanaia, G.; Scudero, S.

    2016-07-01

    Mud volcanoes are common in active mountain fronts. At Mt. Etna, located just between the Apennine front in Sicily and its foredeep, there are some manifestations of mud volcanism in the lower border of the volcanic edifice. The activity of these mud volcanoes is characterized by persistent emission of muddy water mixed with salts, which rises to the surface due to the gas pressure in the subsoil. The San Biagio Salinelle is one of the three mud volcano fields located around the Paternò eruptive monogenic apparatus; this old volcanic structure was one of the first subaerial volcanic manifestations that formed in the pre-Etnean phase. It is not fully clear whether and how the activity of the mud fields is connected with the volcanic activity of Mt. Etna. Noninvasive geophysical surveys were carried out in the area of the active cone of the San Biagio Salinelle, in order to identify the probable ascent path of the emitted products. Seismic ambient noise records were collected at the nodes of a specially designed grid and, subsequently, the V s values were obtained from an active seismic survey. A digital elevation model (DEM) of the area was obtained by a topographic survey, carried out with the GNSS technique (global navigation satellite system), in real-time kinematic mode. The DEM and the topographic benchmark installed will represent the reference surface for future periodic monitoring of the ongoing deformation in the area. Our results provide an accurate and detailed 3D subsurface model showing the shallower feeding system of the investigated mud volcano.

  9. Delivery mechanisms of 3D geological models - a perspective from the British Geological Survey

    Terrington, Ricky; Myers, Antony; Wood, Ben; Arora, Baneet

    2013-04-01

    The past decade has seen the British Geological Survey (BGS) construct over one hundred 3D geological models using software such as GOCAD®, GSI3D, EarthVision and Petrel across the United Kingdom and overseas. These models have been produced for different purposes and at different scales and resolutions in the shallow and deep subsurface. Alongside the construction of these models, the BGS and its collaborators have developed several options for disseminating these 3D geological models to external partners and the public. Initially, the standard formats for disseminating these 3D geological models by the BGS comprised of 2D images of cross-sections, GIS raster data and specialised visualisation software such as the LithoFrame Viewer. The LithoFrame Viewer is a thick-client software that allows the user to explore the 3D geometries of the geological units using a 3D interface, and generate synthetic cross-sections and boreholes on the fly. Despite the increased functionality of the LithoFrame Viewer over the other formats, the most popular data formats distributed remained 2D images of cross-sections, CAD based formats (e.g. DWG and DXF) and GIS raster data of surfaces and thicknesses, as these were the types of data that the external partners were most used too. Since 2009 software for delivering 3D geological models has advanced and types of data available have increased. Feature Manipulation Engine (FME) has been used to increase the number of outputs from 3D geological models. These include: • 3D PDFs (Adobe Acrobat) • KMZ/KML (GoogleEarth) • 3D shapefiles (ESRI) Alongside these later outputs, the BGS has developed other software such as GroundhogTM and Geovisionary (in collaboration with Virtalis). Groundhog is fully a web based application that allows the user to generate synthetic cross-sections, boreholes and horizontal slices from 3D geological models on the fly. Geovisionary provides some of the most advanced visualisation of 3D geological models in

  10. A seismic survey at Colima volcano (Mexico)

    Del Pezzo, E.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione OV, Napoli, Italia; La Rocca, M.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione OV, Napoli, Italia; Galluzzo, D.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione OV, Napoli, Italia; Petrosino, S.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione OV, Napoli, Italia; Cusano, P.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione OV, Napoli, Italia; Bianco, F.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione OV, Napoli, Italia; Breton, M.; Observatorio Vulcanologico de Colima – Università di Colima (Mexico); Orozco-Rojas, J.; Observatorio Vulcanologico de Colima – Università di Colima, (Mexico); Ibanez, J.; Instituto Andaluz de Geofisica - Universidad de Granada (Spain); Veneruso, M.; Centro Regionale di Competenza “Analisi e Monitoraggio del Rischio Ambientale” (AMRA)

    2008-01-01

    In the period 2-6 April 2007 a seismic survey was carried out at Solfatara Volcano, (Campi Flegrei, Southern Italy) with the aim of inferring the shallow structure and evaluating local site effects. Five circular seismic arrays equipped with 1-Hz 3-component Mark LE3Dlite sensors, were installed in the Solfatara crater. Each array consisted of 4 sensors, 3 of them evenly spaced (120°) around the circumference and the fourth placed at its center. The arrays were designed with radii of 5, 10...

  11. Development of Vertical Cable Seismic System for Hydrothermal Deposit Survey (2) - Feasibility Study

    Asakawa, E.; Murakami, F.; Sekino, Y.; Okamoto, T.; Mikada, H.; Takekawa, J.; Shimura, T.

    2010-12-01

    In 2009, Ministry of Education, Culture, Sports, Science and Technology(MEXT) started the survey system development for Hydrothermal deposit. We proposed the Vertical Cable Seismic (VCS), the reflection seismic survey with vertical cable above seabottom. VCS has the following advantages for hydrothermal deposit survey. . (1) VCS is an effective high-resolution 3D seismic survey within limited area. (2) It achieves high-resolution image because the sensors are closely located to the target. (3) It avoids the coupling problems between sensor and seabottom that cause serious damage of seismic data quality. (4) Various types of marine source are applicable with VCS such as sea-surface source (air gun, water gun etc.) , deep-towed or ocean bottom sources. (5) Autonomous recording system. Our first experiment of 2D/3D VCS surveys has been carried out in Lake Biwa, JAPAN. in November 2009. The 2D VCS data processing follows the walk-away VSP, including wave field separation and depth migration. The result gives clearer image than the conventional surface seismic. Prestack depth migration is applied to 3D data to obtain good quality 3D depth volume. Uncertainty of the source/receiver poisons in water causes the serious problem of the imaging. We used several transducer/transponder to estimate these positions. The VCS seismic records themselves can also provide sensor position using the first break of each trace and we calibrate the positions. We are currently developing the autonomous recording VCS system and planning the trial experiment in actual ocean to establish the way of deployment/recovery and the examine the position through the current flow in November, 2010. The second VCS survey will planned over the actual hydrothermal deposit with deep-towed source in February, 2011.

  12. High resolution processing of 3D seismic data for thin coal seam in Guqiao coal mine

    Li, Qiaoling; Peng, Suping; Zou, Guangui

    2015-04-01

    Accurate identification of small faults for coal seams is very important for coal-field exploration, which can greatly improve mining efficiency and safety. However, coal seams in China are mostly thin layers, ranging from 2-5 m. Moreover, the shallow coal seam with strong reflection forms a shield underneath thin coal seam which is only about 40 m deeper. This causes great difficulty in seismic processing and interpretation. The primary concern is to obtain high-resolution seismic image of underneath thin coal seam for mining safety. In this paper, field data is carefully analyzed and fit-for-purpose solutions are adopted in order to improve the quality of reprocessed data and resolution of target coal seam. Identification of small faults has been enhanced significantly.

  13. 3D FEM Numerical Simulation of Seismic Pile-supported Bridge Structure Reaction in Liquefying Ground

    Ling XianZhang, Tang Liang and Xu Pengju

    2011-04-01

    Full Text Available This study examines the establishment of liquefied ground pile-soil-bridge seismic interaction analysis of three-dimensional finite element analysis method for the simulation of liquefied ground shaking table test of pile-soil seismic interaction analysis, undertake OpenSees finite element based numerical simulation platform, for the shaking table test based on two-phase saturated porous media, Comparative numerical and experimental results, detailed test pile dynamic response of bridge structure and dynamic properties, especially liquefaction pore pressure, liquefaction of pile foundation and the dynamic response of the free field. Finite element method can reasonably predict the site of pore pressure, dynamic response; despite the conventional beam element simulation of pile, pile dynamic response can still accurately simulated.

  14. 3-D surface wave tomography of the Piton de la Fournaise volcano using seismic noise correlations

    Brenguier, Florent; M. Shapiro, Nikolai; Campillo, Michel; Nercessian, Alexandre; Ferrazzini, Valérie

    2007-01-01

    [1] We invert Rayleigh waves reconstructed from cross-correlations of 18 months of ambient seismic noise recorded by permanent seismological stations run by the Piton de la Fournaise Volcanological Observatory. By correlating noise records between 21 receivers, we reconstruct Rayleigh waves with sufficient signal-to-noise ratio for 210 inter-station paths. We use the reconstructed waveforms to measure group velocity dispersion curves at periods between 1.5 and 4.5 s. The obtained measurements...

  15. A comparison of methods for 3D target localization from seismic and acoustic signatures

    ELBRING,GREGORY J.; GARBIN,H. DOUGLAS; LADD,MARK D.

    2000-04-03

    An important application of seismic and acoustic unattended ground sensors (UGS) is the estimation of the three dimensional position of an emitting target. Seismic and acoustic data derived from UGS systems provide the taw information to determine these locations, but can be processed and analyzed in a number of ways using varying amounts of auxiliary information. Processing methods to improve arrival time picking for continuous wave sources and methods for determining and defining the seismic velocity model are the primary variables affecting the localization accuracy. Results using field data collected from an underground facility have shown that using an iterative time picking technique significantly improves the accuracy of the resulting derived target location. Other processing techniques show little advantage over simple crosscorrelation along in terms of accuracy, but may improve the ease with which time picks can be made. An average velocity model found through passive listening or a velocity model determined from a calibration source near the target source both result in similar location accuracies, although the use of station correction severely increases the location error.

  16. 3D Digital Surveying and Modelling of Cave Geometry: Application to Paleolithic Rock Art

    Diego González-Aguilera

    2009-02-01

    Full Text Available 3D digital surveying and modelling of cave geometry represents a relevant approach for research, management and preservation of our cultural and geological legacy. In this paper, a multi-sensor approach based on a terrestrial laser scanner, a high-resolution digital camera and a total station is presented. Two emblematic caves of Paleolithic human occupation and situated in northern Spain, “Las Caldas” and “Peña de Candamo”, have been chosen to put in practise this approach. As a result, an integral and multi-scalable 3D model is generated which may allow other scientists, pre-historians, geologists…, to work on two different levels, integrating different Paleolithic Art datasets: (1 a basic level based on the accurate and metric support provided by the laser scanner; and (2 a advanced level using the range and image-based modelling.

  17. Magma Migration Through the Continental Crust - 3-D Seismic and Thermo-mechanical Constraints on Sites of Crustal Contamination

    Wilson, M.; Wheeler, W.

    2002-12-01

    Current understanding of the processes and pathways by which magma travels from its mantle source, through the crust to the Earth's surface is limited by the lack of continuously exposed sections through "fossil" magmatic systems. We report results from a 50 x 30 km 3-D seismic reflection survey of part of the Voring rifted continental margin of Norway which provide the first detailed images of an entire crustal magmatic plumbing system, from a Moho-level magma chamber, through complexes of sills and dykes in the mid to upper crust, to lavas and vent fields extruded at the early Tertiary paleosurface. The Voring margin of Norway formed during a period of Late Cretaceous to early Tertiary (Eocene) continental break-up when Greenland rifted away from Eurasia, resulting in the opening the NE Atlantic Ocean. Rifting was accompanied by widespread magmatic activity, inferred to be related to the impingement of the Iceland mantle plume on the base of the continental lithosphere. Regionally, magma migration occurred in at least two pulses: 62-59 Ma (main initial phase) and 57-54 Ma (continental break-up phase). Wide-angle seismic experiments indicate the presence of a laccolith-like "high-velocity body" (HVB) in the lower crust beneath most of the outer Voring Basin with P-wave velocities (Vp 7.1-7.4 km/s) characteristic of basaltic igneous rocks, overlying typical mantle rocks with Vp of over 8 km/s. The HVB locally reaches 8 km thickness and at break-up (54 Ma) measured 300 km x 500 km - corresponding to a volume of 450,000 cubic km of basaltic magma. It is interpreted as a magmatic underplate formed over a period of several million years as rising basaltic magmas ponded at the Moho at their level of neutral buoyancy. A laterally extensive sill complex (1000 m thick) occurs at the interface between thinned crystalline basement and the overlying Mesozoic sedimentary sequence. This is interpreted as one of the main intra-crustal magma storage reservoirs and is the most

  18. 3D Finite-Difference Modeling of Acoustic Radiation from Seismic Sources

    Chael, E. P.; Aldridge, D. F.; Jensen, R. P.

    2013-12-01

    Shallow seismic events, earthquakes as well as explosions, often generate acoustic waves in the atmosphere observable at local or even regional distances. Recording both the seismic and acoustic signals can provide additional constraints on source parameters such as epicenter coordinates, depth, origin time, moment, and mechanism. Recent advances in finite-difference (FD) modeling methods enable accurate numerical treatment of wave propagation across the ground surface between the (solid) elastic and (fluid) acoustic domains. Using a fourth-order, staggered-grid, velocity-stress FD algorithm, we are investigating the effects of various source parameters on the acoustic (or infrasound) signals transmitted from the solid earth into the atmosphere. Compressional (P), shear (S), and Rayleigh waves all radiate some acoustic energy into the air at the ground surface. These acoustic wavefronts are typically conical in shape, since their phase velocities along the surface exceed the sound speed in air. Another acoustic arrival with a spherical wavefront can be generated from the vicinity of the epicenter of a shallow event, due to the strong vertical ground motions directly above the buried source. Images of acoustic wavefields just above the surface reveal the radiation patterns and relative amplitudes of the various arrivals. In addition, we compare the relative effectiveness of different seismic source mechanisms for generating acoustic energy. For point sources at a fixed depth, double-couples with almost any orientation produce stronger acoustic signals than isotropic explosions, due to higher-amplitude S and Rayleigh waves. Of course, explosions tend to be shallower than most earthquakes, which can offset the differences due to mechanism. Low-velocity material in the shallow subsurface acts to increase vertical seismic motions there, enhancing the coupling to acoustic waves in air. If either type of source breaks the surface (e.g., an earthquake with surface rupture

  19. Modeling 3-D flow in the mantle wedge with complex slab geometries: Comparisons with seismic anisotropy

    Kincaid, C. R.; MacDougall, J. G.; Druken, K. A.; Fischer, K. M.

    2010-12-01

    Understanding patterns in plate scale mantle flow in subduction zones is key to models of thermal structure, dehydration reactions, volatile distributions and magma generation and transport in convergent margins. Different patterns of flow in the mantle wedge can generate distinct signatures in seismological observables. Observed shear wave fast polarization directions in several subduction zones are inconsistent with predictions of simple 2-D wedge corner flow. Geochemical signatures in a number of subduction zones also indicate 3-D flow and entrainment patterns in the wedge. We report on a series of laboratory experiments on subduction driven flow to characterize spatial and temporal variability in 3-D patterns in flow and shear-induced finite strain. Cases focus on how rollback subduction, along-strike dip changes in subducting plates and evolving gaps or tears in subduction zones control temporal-spatial patterns in 3-D wedge flow. Models utilize a glucose working fluid with a temperature dependent viscosity to represent the upper 2000 km of the mantle. Subducting lithosphere is modeled with two rubber-reinforced continuous belts. Belts pass around trench and upper/lower mantle rollers. The deeper rollers can move laterally to allow for time varying dip angle. Each belt has independent speed control and dip adjustment, allowing for along-strike changes in convergence rate and the evolution of slab gaps. Rollback is modeled using a translation system to produce either uniform and asymmetric lateral trench motion. Neutral density finite strain markers are distributed throughout the fluid and used as proxies for tracking the evolution of anisotropy through space and time in the evolving flow fields. Particle image velocimetry methods are also used to track time varying 3-D velocity fields for directly calculating anisotropy patterns. Results show that complex plate motions (rollback, steepening) and morphologies (gaps) in convergent margins produce flows with

  20. Full 3D Microwave Tomography enhanced GPR surveys: a case study

    Catapano, Ilaria; Soldovieri, Francesco; Affinito, Antonio; Hugenschmidt, Johannes

    2014-05-01

    Ground Penetrating Radar (GPR) systems are well assessed non-invasive diagnostic tools capable of providing high resolution images of the inner structure of the probed spatial region. Owing to this capability, GPR systems are nowadays more and more considered in the frame of civil engineering surveys since they may give information on constructive details as well as on the aging and risk factors affecting the healthiness of an infrastructure. In this frame, accurate, reliable and easily interpretable images of the probed scenarios are mandatory in order to support the management of maintenance works and assure the safety of structures. Such a requirement motivates the use of different and sophisticated data processing approaches in order to compare more than one image of the same scene, thus improving the reliability and objectiveness of the GPR survey results. Among GPR data processing procedures, Microwave Tomography approaches based on the Born approximation face the imaging as the solution of a linear inverse problem, which is solved by using the Truncated Singular Value Decomposition as a regularized inversion scheme [1]. So far, an approach exploiting a 2D scalar model of the scattering phenomenon have been adopted to process GPR data gathered along a single scan. In this case, 3D images are obtained by interpolating 2D reconstructions (this is referred commonly as pseudo 3D imaging). Such an imaging approach have provided valuable results in several real cases dealing with not only surveys for civil engineering but also archeological prospection, subservice monitoring, security surveys and so on [1-4]. These encouraging results have motivated the development of a full 3D Microwave Tomography approach capable of accounting for the vectorial nature of the wave propagation. The reconstruction capabilities of this novel approach have been assessed mainly against experimental data collected in laboratory controlled conditions. The obtained results corroborate

  1. Long-range laser scanning and 3D imaging for the Gneiss quarries survey

    Schenker, Filippo Luca; Spataro, Alessio; Pozzoni, Maurizio; Ambrosi, Christian; Cannata, Massimiliano; Günther, Felix; Corboud, Federico

    2016-04-01

    In Canton Ticino (Southern Switzerland), the exploitation of natural stone, mostly gneisses, is an important activity of valley's economies. Nowadays, these economic activities are menaced by (i) the exploitation costs related to geological phenomena such as fractures, faults and heterogeneous rocks that hinder the processing of the stone product, (ii) continuously changing demand because of the evolving natural stone fashion and (iii) increasing administrative limits and rules acting to protect the environment. Therefore, the sustainable development of the sector for the next decades needs new and effective strategies to regulate and plan the quarries. A fundamental step in this process is the building of a 3D geological model of the quarries to constrain the volume of commercial natural stone and the volume of waste. In this context, we conducted Terrestrial Laser Scanning surveys of the quarries in the Maggia Valley to obtain a detailed 3D topography onto which the geological units were mapped. The topographic 3D model was obtained with a long-range laser scanning Riegl VZ4000 that can measure from up to 4 km of distance with a speed of 147,000 points per second. It operates with the new V-line technology, which defines the surface relief by sensing differentiated signals (echoes), even in the presence of obstacles such as vegetation. Depending on the esthetics of the gneisses, we defined seven types of natural stones that, together with faults and joints, were mapped onto the 3D models of the exploitation sites. According to the orientation of the geological limits and structures, we projected the different rock units and fractures into the excavation front. This way, we obtained a 3D geological model from which we can quantitatively estimate the volume of the seven different natural stones (with different commercial value) and waste (with low commercial value). To verify the 3D geological models and to quantify exploited rock and waste volumes the same

  2. Research on 3-D terrain correction methods of airborne gamma-ray spectrometry survey

    The general method of height correction is not effectual in complex terrain during the process of explaining airborne gamma-ray spectrometry data, and the 2-D terrain correction method researched in recent years is just available for correction of section measured. A new method of 3-D sector terrain correction is studied. The ground radiator is divided into many small sector radiators by the method, then the irradiation rate is calculated in certain survey distance, and the total value of all small radiate sources is regarded as the irradiation rate of the ground radiator at certain point of aero- survey, and the correction coefficients of every point are calculated which then applied to correct to airborne gamma-ray spectrometry data. The method can achieve the forward calculation, inversion calculation and terrain correction for airborne gamma-ray spectrometry survey in complex topography by dividing the ground radiator into many small sectors. Other factors are considered such as the un- saturated degree of measure scope, uneven-radiator content on ground, and so on. The results of for- ward model and an example analysis show that the 3-D terrain correction method is proper and effectual. (authors)

  3. Development of a State-Wide 3-D Seismic Tomography Velocity Model for California

    Thurber, C. H.; Lin, G.; Zhang, H.; Hauksson, E.; Shearer, P.; Waldhauser, F.; Hardebeck, J.; Brocher, T.

    2007-12-01

    We report on progress towards the development of a state-wide tomographic model of the P-wave velocity for the crust and uppermost mantle of California. The dataset combines first arrival times from earthquakes and quarry blasts recorded on regional network stations and travel times of first arrivals from explosions and airguns recorded on profile receivers and network stations. The principal active-source datasets are Geysers-San Pablo Bay, Imperial Valley, Livermore, W. Mojave, Gilroy-Coyote Lake, Shasta region, Great Valley, Morro Bay, Mono Craters-Long Valley, PACE, S. Sierras, LARSE 1 and 2, Loma Prieta, BASIX, San Francisco Peninsula and Parkfield. Our beta-version model is coarse (uniform 30 km horizontal and variable vertical gridding) but is able to image the principal features in previous separate regional models for northern and southern California, such as the high-velocity subducting Gorda Plate, upper to middle crustal velocity highs beneath the Sierra Nevada and much of the Coast Ranges, the deep low-velocity basins of the Great Valley, Ventura, and Los Angeles, and a high- velocity body in the lower crust underlying the Great Valley. The new state-wide model has improved areal coverage compared to the previous models, and extends to greater depth due to the data at large epicentral distances. We plan a series of steps to improve the model. We are enlarging and calibrating the active-source dataset as we obtain additional picks from investigators and perform quality control analyses on the existing and new picks. We will also be adding data from more quarry blasts, mainly in northern California, following an identification and calibration procedure similar to Lin et al. (2006). Composite event construction (Lin et al., in press) will be carried out for northern California for use in conventional tomography. A major contribution of the state-wide model is the identification of earthquakes yielding arrival times at both the Northern California Seismic

  4. 3D Digital Surveying and Modelling of Cave Geometry: Application to Paleolithic Rock Art

    Diego González-Aguilera; Angel Muñoz-Nieto; Javier Gómez-Lahoz; Jesus Herrero-Pascual; Gabriel Gutierrez-Alonso

    2009-01-01

    3D digital surveying and modelling of cave geometry represents a relevant approach for research, management and preservation of our cultural and geological legacy. In this paper, a multi-sensor approach based on a terrestrial laser scanner, a high-resolution digital camera and a total station is presented. Two emblematic caves of Paleolithic human occupation and situated in northern Spain, “Las Caldas” and “Peña de Candamo”, have been chosen to put in practise this approach. As a result, an i...

  5. Calibration of 3D Upper Mantle Structure in Eurasia Using Regional and Teleseismic Full Waveform Seismic Data

    Barbara Romanowicz; Mark Panning

    2005-04-23

    Adequate path calibrations are crucial for improving the accuracy of seismic event location and origin time, size, and mechanism, as required for CTBT monitoring. There is considerable information on structure in broadband seismograms that is currently not fully utilized. The limitations have been largely theoretical. the development and application to solid earth problems of powerful numerical techniques, such as the Spectral Element Method (SEM), has opened a new era, and theoretically, it should be possible to compute the complete predicted wavefield accurately without any restrictions on the strength or spatial extent of heterogeneity. This approach requires considerable computational power, which is currently not fully reachable in practice. We propose an approach which relies on a cascade of increasingly accurate theoretical approximations for the computation of the seismic wavefield to develop a model of regional structure for the area of Eurasia located between longitudes of 30 and 150 degrees E, and latitudes of -10 to 60 degrees North. The selected area is particularly suitable for the purpose of this experiment, as it is highly heterogeneous, presenting a challenge for calibration purposes, but it is well surrounded by earthquake sources and, even though they are sparsely distributed, a significant number of high quality broadband digital stations exist, for which data are readily accessible through IRIS (Incorporated Research Institutions for Seismology) and the FDSN (Federation of Digital Seismic Networks). The starting models used will be a combination of a-priori 3D models recently developed for this region, combining various geophysical and seismological data, and a major goal of this study will be to refine these models so as to fit a variety of seismic waveforms and phases.

  6. 3D Seismic Reflection Imaging of Crustal Formation Processes on the East Pacific Rise, 9°57-42'N

    Purdy, G. M.; Mutter, J. C.; Carbotte, S. M.; Canales, J. P.; Nedimovic, M. R.; Carton, H.; Newman, K. R.; Marjanovic, M.; Xu, M.; Aghaei, O.; Stowe, L. C.

    2008-12-01

    Between June 29th and August 19th 2008 the research vessel Marcus G Langseth carried out its first multi- streamer 3D seismic reflection imaging cruise, MGL08-12, by conducting a program research on the East Pacific Rise centered around 9°50'N. The primary goals were to create an accurate 3D seismic reflection image of the magmatic-hydrothermal system at this Integrated Study Site of the Ridge2000 program by imaging the structure of the axial magma chamber (AMC) lid and oceanic crust at a resolution, accuracy, and scale comparable to seafloor observations. The vessel acquired data with four, 6-kilometer solid streamers each comprising 468 active channels deployed with a total separation of 450 meters. Four gun strings with total volume of 3300 cubic inches in two groups fired alternately provide the source for a shot spacing of 37.5 meters. This configuration yields eight CMP lines for each of the sail lines that were spaced 300 m apart, and a static bin size of 6.25 m × 37.5 m in the along-track and across-track directions, respectively, providing a nominal fold of 40. The cruise accomplished the acquisition of ~3,782 km of sail line data. There are 111 across axis lines that required 10 repeated lines and 14 infills. Average feathering during the cruise was 0° ± 5° (one standard deviation), with maximum values of up to 11°. This means that 18% of the total cross axis acquisition was needed for reshoots and infilling. A 25% multiplier on planned lines for a 3D grid is probably a useful figure to use in cruise planning and is fairly standard in the seismic industry. Data quality meets or exceeds industry standards. 3D coverage was achieved in two areas. The larger comprises a set of 93 equally spaced lines forming the 3D grid between 9°57'N and 9°42'N. This grid is made up of lines from all of racetracks #1 and #2 and the northern lines of racetrack#3 and covers two principal hydrothermal vent areas in a continuous fashion. The second 3D area is comprised

  7. The Role of Faulting on the Growth of a Carbonate Platform: Evidence from 3D Seismic Analysis and Section Restoration

    Nur Fathiyah Jamaludin, Siti; Pubellier, Manuel; Prasad Ghosh, Deva; Menier, David; Pierson, Bernard

    2014-05-01

    Tectonics in addition to other environmental factors impacts the growth of carbonate platforms and plays an important role in shaping the internal architecture of the platforms. Detailed of faults and fractures development and healing in carbonate environment have not been explored sufficiently. Using 3D seismic and well data, we attempt to reconstruct the structural evolution of a Miocene carbonate platform in Central Luconia Province, offshore Malaysia. Luconia Province is located in the NW coast of Borneo and has become one of the largest carbonate factories in SE Asia. Seismic interpretations including seismic attribute analysis are applied to the carbonate platform to discern its sedimentology and structural details. Detailed seismic interpretations highlight the relationships of carbonate deposition with syn-depositional faulting. Branching conjugate faults are common in this carbonate platform and have become a template for reef growth, attesting lateral facies changes within the carbonate environments. Structural restoration was then appropriately performed on the interpreted seismic sections based on sequential restoration techniques, and provided images different from those of horizon flattening methods. This permits us to compensate faults' displacement, remove recent sediment layers and finally restore the older rock units prior to the fault motions. It allows prediction of platform evolution as a response to faulting before and after carbonate deposition and also enhances the pitfalls of interpretation. Once updated, the reconstructions allow unravelling of the un-seen geological features underneath the carbonate platform, such as paleo-structures and paleo-topography which in turn reflects the paleo-environment before deformations took place. Interestingly, sections balancing and restoration revealed the late-phase (Late Oligocene-Early Miocene) rifting of South China Sea, otherwise difficult to visualize on seismic sections. Later it is shown that

  8. Assessing the quality of earthquake source models using 3-D forward modelling of long-period seismic data

    Ferreira, A. M.; Vallée, M.; Lentas, K.

    2010-12-01

    Accurate earthquake point source parameters (e.g. seismic moment, depth and focal mechanism) provide key first-order information for detailed studies of the earthquake source process and for improved seismic and tsunami hazard evaluation. In order to objectively assess the quality of seismic source models, it is important to go beyond classical resolution/misfit checks. In particular, it is desirable to apply sophisticated modeling techniques to quantify uncertainties due to simplified theoretical formulations and/or Earth structure employed to build the source models. Moreover, it is important to verify how well the models explain data not used in their construction for a complete, quantitative assessment of the earthquake source models. In this study we compare the quality of the surface-wave Centroid Moment Tensor (CMT) method with that of the SCARDEC method, which is a new automated body-wave technique for the fast simultaneous determination of the seismic moment, focal mechanism, depth and source time functions of large earthquakes. We focus on the major shallow subduction earthquakes of the last 20 years, for which there are some systematic differences between SCARDEC and CMT source parameters, notably in fault dip angle and moment magnitude. Because the SCARDEC method is based on body-wave deconvolution using ray methods in a 1D Earth model, we test how well SCARDEC source parameters explain long-period seismic data (surface waves and normal modes) compared to the CMT method. We calculate theoretical seismograms using two forward modelling techniques (full ray theory and spectral element method) to simulate the long-period seismic wavefield for the 3D Earth model S20RTS combined with the crust model CRUST2.0, and for two point source models: (i) the SCARDEC model; and (ii) the Global CMT model. We compare the synthetic seismograms with real broadband data from the FDSN for the major subduction earthquakes of the last 20 years. We show that SCARDEC source

  9. Investigation of data acquisition parameters for Minami Noshiro 3D experiment using 3D seismic modeling; Sanjigen hado denpa modeling wo riyoshita Minami Noshiro sanjigen jishin tansa data shutoku parameter no kento

    Tanaka, H.; Nakagami, K.; Minegishi, M. [Japan National Oil Corp., Tokyo (Japan). Technology Research Center; Kano, R. [Teikoku Oil Co. Ltd., Tokyo (Japan)

    1996-10-01

    This paper describes how to grasp the characteristics of reflected waves which means how the reflected wave can be received in its magnitude depending on the complicated subsurface structure. Data acquisition parameters were also investigated, as to offset distance distribution and mute setting during the HMO correction. A velocity structure model was prepared for an area of 9km{times}8km including the given 3D seismic exploration area at Minami Noshiro, Akita Prefecture. For the geological sequence used for the velocity structure model, three formations, i.e., Katsurane Formation, Funakawa Formation, and Onagawa Formation, were inputted in an order from the shallower depth on the basis of the depth structure profile obtained from the previous data. Ray tracing was calculated by means of the two-dot dashed line tracing method. For this method, amplitude as well as travel time of waves can be calculated at the same time. This was effective for grasping the magnitude of reflected wave on simulating the traced data. For the velocity structure model used in this study, existing data inputted were old, which limited the quantity of information. However, this model would sufficiently contribute to the verification of survey design and the determination of optimal layout. 1 ref., 3 figs.

  10. Center Line Coordinates Survey for Existing Railway by 3-D Constraints Method

    Hui Wei

    2013-05-01

    Full Text Available Since the maintenance of HSR depends both on inner geometries and outer geometries information of track, it’s necessary that the survey technologies could measure the smoothness and coordinates of track in a high-density manner. But the existing technologies for track survey are either inaccurate or inefficient for HSR maintenance, the problem about track survey arises with the rapid development of high-speed rail (HSR in China. To meet the demands on track smoothness survey of existing line, a method denoted as 3-D Constraints Method is introduced in this paper. The general principle of the method is obtaining center line locus with location information by means of survey trolley and CP-Ⅲ stakes. By measuring the distances between trolley and CP-Ⅲ stakes, the boundary condition could be determined. So, the problem of center line calculation is transformed into the problem of definite integral for FOG signals. Accordingly, an algorithm of coordinate calculation is given. With the constraints form CP-Ⅲ stakes, the coordinates of center line could be obtained in a defined accuracy. Finally, the simulation result shows that the method is effective.

  11. 3D seismic expression of fluid migration and mud remobilization on the Gjallar Ridge, offshore mid-Norway

    Hansen, J.P.V.; Cartwright, J.A.; Huuse, M.;

    2005-01-01

    anomalies and possible mud volcanoes at the base Pleistocene indicating their long-term significance as vertical fluid conduits. Type C structures comprise discrete mound features that seem to jack up the Top Palaeocene (Top Brygge) horizon. These are similar to hydrothermal mounds found elsewhere on the......This paper presents a three-dimensional (3D) seismic analysis of sediment remobilization and fluid migration in a 2000-km2 area above the Gjallar Ridge located in the Vøring Basin, offshore Norway. Three distinct types of mounded structures have been identified as resulting from focused fluid....../gas migration and associated mud remobilization and intrusion. Type A structures are gently mounded, and we infer that these structures formed because of in situ remobilization of Middle Eocene to LowerMiddle Oligocene fine-grained sediments in response to fluid and minor sediment injection via deep...

  12. 3D Reconstruction of Coronary Arteries from Angiographic Images: A Survey

    Mahima Goyal

    2013-06-01

    Full Text Available X-Ray angiographyis considered to be the “golden standard” ofall times in the medical imaging field due to its wide range ofapplications. In this paper, a survey is performed on the basictechniques and the methodologies that have already beenproposed by various researchers in the field of 3D reconstructionof coronary arteries using one or more angiograms as inmonoplane or biplane angiographic systems or a volume of dataas in rotational angiographic systems. Also, the variousprocedures that need to be performed are stated and describedshowing the contribution of each step in the reconstructionprocess. These procedures include distortion correction, motioncompensation, feature extraction, background removal,projection geometry optimization and topology identification,surface visualization etc.X-Ray angiographyis considered to be the “golden standard” ofall times in the medical imaging field due to its wide range ofapplications. In this paper, a survey is performed on the basictechniques and the methodologies that have already beenproposed by various researchers in the field of 3D reconstructionof coronary arteries using one or more angiograms as inmonoplane or biplane angiographic systems or a volume of dataas in rotational angiographic systems. Also, the variousprocedures that need to be performed are stated and describedshowing the contribution of each step in the reconstructionprocess. These procedures include distortion correction, motioncompensation, feature extraction, background removal,projection geometry optimization and topology identification,surface visualization etc.

  13. Investigation of drilling failure of Well A1, E-Field, onshore Niger Delta, Nigeria, using 3-D seismic data

    Akinmosin, A.; Oladele, S.; Oriade, O. F.

    2016-02-01

    This study aimed at investigating reasons for failure of Well A1 that is surrounded by hydrocarbon discoveries in onshore Niger delta with a view to propose optimal location for a new well through interpretation of new 3-D seismic data. Sands encountered by Well A1 were delineated and tied to seismic. Structural closure was mapped and reservoirs at various depths were stacked and sectioned. Porosity, Net to Gross, and Gross Rock Volume of the reservoirs were computed. Well A1 was correlated to a nearby Well K4 and a good correlation was observed. A fault assisted multi reservoirs Prospect-E with south-westerly shift with depth was delineated on the hanging wall of structure building E-Fault whose closing contours is expected to trap hydrocarbon. Petrophysical properties of the reservoirs range from fair to good. Well A1 either perforated prospect-E beneath the oil water contact, completely missed prospect-E or punctured the wet foot wall of E-fault. Consequently, Well A1 could not impact the objective sands and only able to produce water. Well A1 failed because of wrong surface positioning which unsuccessfully targeted the shifting reservoirs. A successful exploratory well in E-Field would be a gently south-westerly deviated well whose vertical section would encounter the shallower sands and deviated section targeting the deeper sands with surface location at about 1 km southeast of Well A1.

  14. Effects of eccentricity and order of vibration modes on the inelastic seismic response of 3D steel structures

    P. Sharifi, Naser; Sakulich, Aaron R.

    2014-04-01

    In torsionally coupled buildings, the total response of the structure is the result of the translational displacement of the story's center of stiffness and the displacement due to the roof's rotation. In structures with high eccentricity, the effect of the floor's rotation in the total response is considerable. The order of vibration modes is another important parameter that changes the contribution of the different translational and rotational modes in the total response. To explore the effects of eccentricity and the order of vibration modes on the total response, a number of 3-D steel moment-resistant frames with 4, 8, and 12 stories, with different eccentricities and plans, were considered. The structures were subjected to bidirectional seismic inputs so that their peak ground accelerations were scaled to 0.4g, 0.6g, and 0.8g. Increasing the eccentricity of the structure increases the participation of rotation in the total response. Furthermore, in torsionally flexible structures, where the first or second mode of vibration is a torsional mode, the contribution of the floor's rotation can be even greater. In some cases, the displacement of exterior columns is primarily the result of the floor's rotation. This suggests that to efficiently dampen the seismic displacement of such structures, the rotational mode of the building should be controlled.

  15. Identification of kimberlite bodies in Brazil from a 3D audio-magnetotelluric survey

    De Lugao, P. P.; Eric, C. D. O.; Loureiro, F. O.; Arantes, P. R.; Pastana, A. F.

    2015-12-01

    We report on a succesfull identification of kimberlite bodies in Brazil through the use of the electromagnetic technique audio-magnetotelluric (AMT). Macnae (1979) writes that "In one large survey in South Africa, electromagnetic (EM) techniques have proven to be remarkably effective in detecting the presence of weathered clays or epiclastic kimberlite contained within the pipes." Full tensor AMT data were acquired at 65 points (stations) in a 3D configuration with frequencies ranging from 10kHz to 1Hz. The survey was located in the NW portion of the Mato Grosso state, Brazil, in na area of thick jungle coverage. During the AMT survey, few outcrops were seen because of the dense forest cover. Usually, the occurrences found were of sand deposits, indicating the occurence of Fazenda Casa Branca and Utiariti Formations and gravel from Salto das Nuvens Formation, widely used in paving trails n this region. In the area of the survey, three main targets were confirmed/identified: Kimberlite Area 1 - a classic kimberlite in the region, with the crater facies with different clasts and distinct size. We noted the occurrence of a red-brown soil and an unusual vegetation in this area. The resistivity model provided confirmed the presence of Kimberlite Area 1 and was used to identify other two areas. Area of Interest 1 - area with atypical vegetation along a trail. There is an excavation that displays soil of white color with several blocks present, there are small quartz crystal agglomerates in these blocks. The resistivity model cleary shows a conductive body here, indicative of the presence of a kimberlite. Area of Interest 2 - the presence of a kimberlite was confirmed, not exactly where the targeted Area 2 was, but the southwest of it. Close to this area, there was a very fine rock and a few blocks of pure silica, probably indicating a kimberlitic intrusion. In summary, the 3D resistivity model in depth obtained from inversion of the AMT data confirmed and identified

  16. A seismic modeling analysis of wide and narrow 3D observation systems for channel sand bodies%河道砂体宽/窄方位三维观测系统地震物理模型分析

    狄帮让; 徐秀仓; 魏建新

    2008-01-01

    The effect of the wide and narrow azimuth 3D observation systems on seismic imaging precision is becoming a hot area for studies of high precision 3D seismic acquisition methods in recent years. In this paper we utilize 3D physical seismic modeling experiments. A 3D channel sand body physical seismic model is constructed and two acquisition systems are designed with wide azimuth (16 lines) and narrow azimuth (8 lines) to model 3D seismic data acquisition and processing seismic work flows. From analysis of migrated time slice data with high quality and small size, we conclude that when the overlying layers are smooth and lateral velocities have little change, both wide and narrow azimuth observation systems in 3D acquisition can be used for obtaining high precision imaging and equivalent resolution of the channel sand body.

  17. A 3D Seismic Velocity Model Offshore Southern California from Ambient Noise Tomography of the ALBACORE OBS Array

    Kohler, M. D.; Bowden, D. C.; Tsai, V. C.; Weeraratne, D. S.

    2015-12-01

    The Pacific-North America plate boundary in Southern California extends far west of the coastline, and a 12-month ocean bottom seismometer (OBS) array spanned the western side of the plate boundary to image lithospheric seismic velocities. Velocities are modeled through stacked cross correlations of ambient noise data. Twelve months of continuous data were used from 22 OBS stations and ~30 coastal and island Southern California Seismic Network stations. Particular attention has been paid to improving signal-to-noise ratios in the noise correlations with OBS stations by removing the effects of instrument tilt and infragravity waves. Different applications of preprocessing techniques allow us to distinguish the fundamental and first higher order Rayleigh modes, especially in deep water OBS pairs where the water layer dominates crustal sensitivity of the fundamental mode. Standard time domain and frequency domain methods are used to examine surface wave dispersion curves for group and phase velocities between 5 and 50 second periods, and these are inverted for 3D velocity structure. The results define the transition in three dimensions from continental lithospheric structure in the near-shore region to oceanic structure west of the continental borderland. While the most prominent features of the model relate to thinning of the crust west of the Patton Escarpment, other notable anomalies are present north-to-south throughout the continental borderland and along the coast from the Los Angeles Basin to the Peninsular Ranges. The velocity model will help describe the region's tectonic history, as well as provide new constraints for determination of earthquake relocations and rupture styles.

  18. Newly discovered abundant fluid seep indicators off southern Costa Rica, imaged from overlapping multibeam swaths and 3D seismic data

    Kluesner, J. W.; Silver, E. A.; Gibson, J. C.; Bangs, N. L.; McIntosh, K.; von Huene, R.; Orange, D.; Ranero, C. R.

    2012-12-01

    Offshore southern Costa Rica we have identified 161 potential fluid seepage sites on the shelf and slope regions within an 11 x 55 km strip where no fluid indicators had been reported previously using conventional deep-water mutlibeam bathymetry (100 m grid cell size) and deep towed side scan sonar. Evidence includes large and small pockmarks, mounds, ridges, and slope failure features with localized anomalous high-amplitude backscatter strength. The majority of seepage indicators are associated with shallow sub-bottom reversed polarity bright spots and flat spots imaged within the CRISP 3D seismic grid. Data were collected ~50 km west of Osa Peninsula, Costa Rica onboard the R/V Marcus G. Langseth during the spring of 2011. We obtained EM122 multibeam data using fixed, closely spaced receiver beams and 9-10 times swath overlap, which greatly improved the signal-to-noise ratio and sounding density and allowed for very small grid and mosaic cell sizes (2-10 m). A gas plume in the water column, seen on a 3.5 kHz profile, is located along a fault trace and above surface and subsurface seep indicators. Fluid indicators on the outer shelf occur largely on a dense array of faults, some of which cut through the reflective basement. Seismic flat spots commonly underlie axes of large anticlines on the shelf and slope. Pockmarks are also located at the foot of mid-slope canyons, very near to the upper end of the BSR. These pockmarks appear to be associated with canyon abandonment and folded beds that channel fluids upward, causing hydrate instability. Our findings suggest that significant amounts of methane are venting into ocean and potentially into the atmosphere across the heavily deformed shelf and slope of Costa Rica.

  19. A Survey of Recent View-based 3D Model Retrieval Methods

    Liu, Qiong

    2012-01-01

    Extensive research efforts have been dedicated to 3D model retrieval in recent decades. Recently, view-based methods have attracted much research attention due to the high discriminative property of multi-views for 3D object representation. In this report, we summarize the view-based 3D model methods and provide the further research trends. This paper focuses on the scheme for matching between multiple views of 3D models and the application of bag-of-visual-words method in 3D model retrieval....

  20. 3D high resolution mineral phase distribution and seismic velocity structure of the transition zone: predicted by a full spherical-shell compressible mantle convection model

    Geenen, T.; Heister, T.; Van Den Berg, A. P.; Jacobs, M.; Bangerth, W.

    2011-12-01

    We present high resolution 3D results of the complex mineral phase distribution in the transition zone obtained by numerical modelling of mantle convection. We extend the work by [Jacobs and van den Berg, 2011] to 3D and illustrate the efficiency of adaptive mesh refinement for capturing the complex spatial distribution and sharp phase transitions as predicted by their model. The underlying thermodynamical model is based on lattice dynamics which allows to predict thermophysical properties and seismic wave speeds for the applied magnesium-endmember olivine-pyroxene mineralogical model. The use of 3D geometry allows more realistic prediction of phase distribution and seismic wave speeds resulting from 3D flow processes involving the Earth's transition zone and more significant comparisons with interpretations from seismic tomography and seismic reflectivity studies aimed at the transition zone. Model results are generated with a recently developed geodynamics modeling application based on dealII (www.dealii.org). We extended this model to incorporate both a general thermodynamic model, represented by P,T space tabulated thermophysical properties, and a solution strategy that allows for compressible flow. When modeling compressible flow in the so called truncated anelastic approximation framework we have to adapt the solver strategy that has been proven by several authors to be highly efficient for incompressible flow to incorporate an extra term in the continuity equation. We present several possible solution strategies and discuss their implication in terms of robustness and computational efficiency.

  1. Preliminary 3d depth migration of a network of 2d seismic lines for fault imaging at a Pyramid Lake, Nevada geothermal prospect

    Frary, R.; Louie, J. [UNR; Pullammanappallil, S. [Optim; Eisses, A.

    2016-08-01

    Roxanna Frary, John N. Louie, Sathish Pullammanappallil, Amy Eisses, 2011, Preliminary 3d depth migration of a network of 2d seismic lines for fault imaging at a Pyramid Lake, Nevada geothermal prospect: presented at American Geophysical Union Fall Meeting, San Francisco, Dec. 5-9, abstract T13G-07.

  2. Hypocenter relocation using a fast grid search method and a 3-D seismic velocity model for the Sumatra region

    Determination of earthquake hypocenter in Indonesia conducted by the Meteorological, Climatological, and Geophysical Agency (MCGA) has still used a 1-D seismic velocity model. In this research, we have applied a Fast Grid Search (FGM) method and a 3-D velocity model resulting from tomographic imaging to relocate earthquakes in the Sumatran region. The data were taken from the MCGA data catalog from 2009 to 2011 comprising of subduction zone and on land fault earthquakes with magnitude greater than 4 Mw. Our preliminary results show some significant changes in the depths of the relocated earthquakes which are in general deeper than the depths of hypocenters from the MCGA data catalog. The residual times resulting from the relocation process are smaller than those prior to the relocation. Encouraged by these results, we will continue to conduct hypocenter relocation for all events from the MCGA data catalog periodically in order to produce a new data catalog with good quality. We hope that the new data catalog will be useful for further studies

  3. Hypocenter relocation using a fast grid search method and a 3-D seismic velocity model for the Sumatra region

    Nugroho, Hendro [Study Program of Earth Sciences, Faculty of Earth Sciences and Technology, Institute of Technology Bandung, Jl. Ganesha No. 10, Bandung 40132, Indonesia and Meteorological, Climatological, and Geophysical Agency, Jl. Angkasa 1 No. 2, Kemayoran, Jakar (Indonesia); Widiyantoro, Sri [Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Jl. Ganesha No. 10, Bandung 40132 (Indonesia); Nugraha, Andri Dian [Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institute of Technologyc Bandung, Jl. Ganesha No. 10, Bandung 40132 (Indonesia)

    2013-09-09

    Determination of earthquake hypocenter in Indonesia conducted by the Meteorological, Climatological, and Geophysical Agency (MCGA) has still used a 1-D seismic velocity model. In this research, we have applied a Fast Grid Search (FGM) method and a 3-D velocity model resulting from tomographic imaging to relocate earthquakes in the Sumatran region. The data were taken from the MCGA data catalog from 2009 to 2011 comprising of subduction zone and on land fault earthquakes with magnitude greater than 4 Mw. Our preliminary results show some significant changes in the depths of the relocated earthquakes which are in general deeper than the depths of hypocenters from the MCGA data catalog. The residual times resulting from the relocation process are smaller than those prior to the relocation. Encouraged by these results, we will continue to conduct hypocenter relocation for all events from the MCGA data catalog periodically in order to produce a new data catalog with good quality. We hope that the new data catalog will be useful for further studies.

  4. Potential Geophysical Field Transformations and Combined 3D Modelling for Estimation the Seismic Site Effects on Example of Israel

    Eppelbaum, Lev; Meirova, Tatiana

    2015-04-01

    It is well-known that the local seismic site effects may have a significant contribution to the intensity of damage and destruction (e.g., Hough et al., 1990; Regnier et al., 2000; Bonnefoy-Claudet et al., 2006; Haase et al., 2010). The thicknesses of sediments, which play a large role in amplification, usually are derived from seismic velocities. At the same time, thickness of sediments may be determined (or defined) on the basis of 3D combined gravity-magnetic modeling joined with available geological materials, seismic data and borehole section examination. Final result of such investigation is a 3D physical-geological model (PGM) reflecting main geological peculiarities of the area under study. Such a combined study needs in application of a reliable 3D mathematical algorithm of computation together with advanced methodology of 3D modeling. For this analysis the developed GSFC software was selected. The GSFC (Geological Space Field Calculation) program was developed for solving a direct 3-D gravity and magnetic prospecting problem under complex geological conditions (Khesin et al., 1996; Eppelbaum and Khesin, 2004). This program has been designed for computing the field of Δg (Bouguer, free-air or observed value anomalies), ΔZ, ΔX, ΔY , ΔT , as well as second derivatives of the gravitational potential under conditions of rugged relief and inclined magnetization. The geological space can be approximated by (1) three-dimensional, (2) semi-infinite bodies and (3) those infinite along the strike closed, L.H. non-closed, R.H. on-closed and open). Geological bodies are approximated by horizontal polygonal prisms. The program has the following main advantages (besides abovementioned ones): (1) Simultaneous computing of gravity and magnetic fields; (2) Description of the terrain relief by irregularly placed characteristic points; (3) Computation of the effect of the earth-air boundary by the method of selection directly in the process of interpretation; (4

  5. Seismic methods for the characterisation of reservoirs in developing old natural gas fields in Germany; 3D Seismische Verfahren zur Reservoircharakterisierung bei der Entwicklung alter Erdoelfelder in Deutschland

    Krajewski, P.; Stahl, E.; Bischoff, R. [Preussag Energie GmbH, Lingen (Germany); Guderian, K.; Hasse, G.; Schmiermann, I. [BEB Erdoel und Erdgas GmbH, Hannover (Germany); Groot, P. de [De Groot-Bril Earth Sciences BV, Enschede (Netherlands)

    1998-12-31

    Two examples are chosen to describe the possiblities and limitations of using 3D seismic data for the interpretation of structures and the seismic characterisation of reservoirs. New techniques of seismic classification offer a great deal of possibilities, especially if - as in the case of Ruehme - there is a sufficiency of data from many borehole locations which enables the training of algorithms.(orig.) [Deutsch] Anhand zweier Beispiele wurden die Moeglichkeiten aber auch die Grenzen des Einsatzes 3D seismischer Daten bei der strukturellen Interpretation und der seismischen Reservoircharakterisierung aufgezeigt. Neuartige Techniken der seismischen Klassifizierung erweitern die Moeglichkeiten dabei betraechtlich, insbesondere, wenn - wie beim Beispiel Ruehme - durch die vielen Bohrlokationen ausreichend Daten zum Trainieren der Algorithmen zur Verfuegung stehen. (orig.)

  6. 3D seismic, geochemical and biostratigraphical analysis of Paleogene remobilized sand in the Norwegian-Danish Basin

    Andresen, Katrine Juul; Rønø Clausen, Ole; Heilmann-Clausen, Claus; Friis, Henrik

    2013-04-01

    This study describes remobilized Paleogene sand occurring on the hanging-wall segment north of the major D-1 normal fault in the Norwegian-Danish Basin, eastern North Sea. The remobilized sand is observed on 3D seismic data in fine-grained Eocene host-strata as cross-cutting reflections with a typical tabular, V-shaped or wing-like geometry in the seismic cross-sections and a pronounced jack-up of the overlying succession onto which onlap can be observed. In map view the remobilized sand in certain areas have a channel-like appearance. The seismic observations indicate that the sand has a remobilized origin which may be partly depositional. Particularly the observed wings and jack-up on the seismic cross-sections indicate remobilization which potentially could be generated by two different processes: a) remobilization of depositional channel sand resulting in the formation of injected wings along the sides of the channel, or b) injection of remobilized sand from the deeper Paleocene strata causing jack-up and typically V-shaped intrusions. Injection of Paleocene sand into Eocene host strata is a well-known phenomenon from the nearby Paleogene Siri Canyon located c. 15 km north of the study area. In order to acquire more information about the intrusions a geochemical study and a detailed biostratigraphical dating of cuttings and sidewall core samples from the Floki-1 well was carried out. The Floki-1 well penetrates the remobilized sand and was drilled to test an apparent 4-way closure on prospect Eocene sand which by then was interpreted to be 100 % depositional. The geochemical study of the samples from the sand identified the Floki-sandstone as a very fine grained sand and silt with a matrix of very angular silt grains. The sand does not contain clays. The matrix appears to have formed by crushing of the sand grains. Thus, heavy minerals appear to have disintegrated by crushing but still most parts of the mineral grain is found together. Glauconite grains are

  7. High-resolution seismic reflection survey near SPR surface collapse feature at Weeks Island, Louisiana

    Miller, R.D.; Xia, J. [Kansas Geological Survey, Lawrence, KS (United States); Harding, R.S. Jr. [Sandia National Labs., Albuquerque, NM (United States); Steeples, D.W. [Univ. of Kansas, Lawrence, KS (United States)

    1994-12-31

    Shallow high resolution 2-D and 3-D seismic reflection techniques are assisting in the subsurface delineation of a surface collapse feature (sinkhole) at Weeks Island, Louisiana. Seismic reflection surveys were conducted in March 1994. Data from walkaway noise tests were used to assist selection of field recording parameters. The top of the salt dome is about 180 ft below ground surface at the sinkhole. The water table is an estimated 90 ft below the ground surface. A single coherent reflection was consistently recorded across the entire area of the survey, although stacking velocity and spectral content of the event varied. On the basis of observed travel times and stacking velocities, the coherent reflection event appears to originate above the top of the salt, possibly at or near the water table. Identification of this reflector will be made form borehole investigations currently planned for the sinkhole site. A depression or time sag in this reflection event is clearly evident in both the 2-D and 3-D seismic data in the immediate vicinity of the sinkhole. The time sag appears to be related to the subsurface structure of the reflector and not to near surface topography or velocity effects. Elsewhere in the survey area, observed changes in reflection travel times and wavelet character appear to be related to subsurface geologic structure. These seismic observations may assist in predicting where future sinkholes will develop after they have been tied to borehole data collected at the site.

  8. Fault and dyke detectability in high resolution seismic surveys for coal: a view from numerical modelling*

    Zhou, Binzhong 13Hatherly, Peter

    2014-10-01

    Modern underground coal mining requires certainty about geological faults, dykes and other structural features. Faults with throws of even just a few metres can create safety issues and lead to costly delays in mine production. In this paper, we use numerical modelling in an ideal, noise-free environment with homogeneous layering to investigate the detectability of small faults by seismic reflection surveying. If the layering is horizontal, faults with throws of 1/8 of the wavelength should be detectable in a 2D survey. In a coal mining setting where the seismic velocity of the overburden ranges from 3000 m/s to 4000 m/s and the dominant seismic frequency is ~100 Hz, this corresponds to a fault with a throw of 4-5 m. However, if the layers are dipping or folded, the faults may be more difficult to detect, especially when their throws oppose the trend of the background structure. In the case of 3D seismic surveying we suggest that faults with throws as small as 1/16 of wavelength (2-2.5 m) can be detectable because of the benefits offered by computer-aided horizon identification and the improved spatial coherence in 3D seismic surveys. With dykes, we find that Berkhout's definition of the Fresnel zone is more consistent with actual experience. At a depth of 500 m, which is typically encountered in coal mining, and a 100 Hz dominant seismic frequency, dykes less than 8 m in width are undetectable, even after migration.

  9. Exposure to seismic survey alters blue whale acoustic communication.

    Di Iorio, Lucia; Clark, Christopher W

    2010-02-23

    The ability to perceive biologically important sounds is critical to marine mammals, and acoustic disturbance through human-generated noise can interfere with their natural functions. Sounds from seismic surveys are intense and have peak frequency bands overlapping those used by baleen whales, but evidence of interference with baleen whale acoustic communication is sparse. Here we investigated whether blue whales (Balaenoptera musculus) changed their vocal behaviour during a seismic survey that deployed a low-medium power technology (sparker). We found that blue whales called consistently more on seismic exploration days than on non-exploration days as well as during periods within a seismic survey day when the sparker was operating. This increase was observed for the discrete, audible calls that are emitted during social encounters and feeding. This response presumably represents a compensatory behaviour to the elevated ambient noise from seismic survey operations. PMID:19776059

  10. Dynamical analysis of a PWR internals using super-elements in an integrated 3-D model model. Part 2: dynamical tests and seismic analysis

    The results of the test analysis (frequencies) for the isolated super-elements and for the developed 3-D model of the internals core support structures of a PWR research reactor are presented. Once certified of the model effectiveness for this type of analysis the seismic spectral analysis was performed. From the results can be seen that the structures are rigid for this load, isolated or together with the other in the 3-D model, and there are no impacts among them during the earthquake (OBE). (author)

  11. Techniques for Surveying Urban Active Faults by Seismic Methods

    Xu Mingcai; Gao Jinghua; Liu Jianxun; Rong Lixin

    2005-01-01

    Using the seismic method to detect active faults directly below cities is an irreplaceable prospecting technique. The seismic method can precisely determine the fault position. Seismic method itself can hardly determine the geological age of fault. However, by considering in connection with the borehole data and the standard geological cross-section of the surveyed area, the geological age of reflected wave group can be qualitatively (or semi-quantitatively)determined from the seismic depth profile. To determine the upper terminal point of active faults directly below city, it is necessary to use the high-resolution seismic reflection technique.To effectively determine the geometric feature of deep faults, especially to determine the relation between deep and shallow fracture structures, the seismic reflection method is better than the seismic refraction method.

  12. The application of high-resolution 3D seismic data to model the distribution of mechanical and hydrogeological properties of a potential host rock for the deep storage of radioactive waste in France

    Mari, Jean-Luc

    2014-01-01

    In the context of a deep geological repository of high-level radioactive wastes, the French National Radioactive Waste Management Agency (Andra) has conducted an extensive characterization of the Callovo-Oxfordian argillaceous rock and surrounding formations in the Eastern Paris Basin. As part of this project, an accurate 3D seismic derived geological model is needed. The paper shows the procedure used for building the 3D seismic constrained geological model in depth by combining time-to-depth conversion of seismic horizons, consistent seismic velocity model and elastic impedance in time. It also shows how the 3D model is used for mechanical and hydrogeological studies. The 3D seismic field data example illustrates the potential of the proposed depth conversion procedure for estimating density and velocity distributions, which are consistent with the depth conversion of seismic horizons using the Bayesian Kriging method. The geological model shows good agreement with well log data obtained from a reference we...

  13. The Representation of Cultural Heritage from Traditional Drawing to 3d Survey: the Case Study of Casamary's Abbey

    Canciani, M.; Saccone, M.

    2016-06-01

    In 3D survey the aspects most discussed in the scientific community are those related to the acquisition of data from integrated survey (laser scanner, photogrammetric, topographic and traditional direct), rather than those relating to the interpretation of the data. Yet in the methods of traditional representation, the data interpretation, such as that of the philological reconstruction, constitutes the most important aspect. It is therefore essential in modern systems of survey and representation, filter the information acquired. In the system, based on the integrated survey that we have adopted, the 3D object, characterized by a cloud of georeferenced points, defined but their color values, defines the core of the elaboration. It allows to carry out targeted analysis, using section planes as a tool of selection and filtering data, comparable with those of traditional drawings. In the case study of the Abbey of Casamari (Veroli), one of the most important Cistercian Settlement in Italy, the survey made for an Agreement with the Ministry of Cultural Heritage and Activities and Tourism (MiBACT) and University of RomaTre, within the project "Accessment of the sismic safety of the state museum", the reference 3D model, consisting of the superposition and geo-references data from various surveys, is the tool with which yo develop representative models comparable to traditional ones. It provides the necessary spatial environment for drawing up plans and sections with a definition such as to develop thematic analysis related to phases of construction, state of deterioration and structural features.

  14. A questionnaire-based survey on 3D image-guided brachytherapy for cervical cancer in Japan. Advances and obstacles

    The purpose of this study is to survey the current patterns of practice, and barriers to implementation, of 3D image-guided brachytherapy (3D-IGBT) for cervical cancer in Japan. A 30-item questionnaire was sent to 171 Japanese facilities where high-dose-rate brachytherapy devices were available in 2012. In total, 135 responses were returned for analysis. Fifty-one facilities had acquired some sort of 3D imaging modality with applicator insertion, and computed tomography (CT) and magnetic resonance imaging (MRI) were used in 51 and 3 of the facilities, respectively. For actual treatment planning, X-ray films, CT and MRI were used in 113, 20 and 2 facilities, respectively. Among 43 facilities where X-ray films and CT or MRI were acquired with an applicator, 29 still used X-ray films for actual treatment planning, mainly because of limited time and/or staffing. In a follow-up survey 2.5 years later, respondents included 38 facilities that originally used X-ray films alone but had indicated plans to adopt 3D-IGBT. Of these, 21 had indeed adopted CT imaging with applicator insertion. In conclusion, 3D-IGBT (mainly CT) was implemented in 22 facilities (16%) and will be installed in 72 (53%) facilities in the future. Limited time and staffing were major impediments. (author)

  15. A questionnaire-based survey on 3D image-guided brachytherapy for cervical cancer in Japan: advances and obstacles.

    Ohno, Tatsuya; Toita, Takafumi; Tsujino, Kayoko; Uchida, Nobue; Hatano, Kazuo; Nishimura, Tetsuo; Ishikura, Satoshi

    2015-11-01

    The purpose of this study is to survey the current patterns of practice, and barriers to implementation, of 3D image-guided brachytherapy (3D-IGBT) for cervical cancer in Japan. A 30-item questionnaire was sent to 171 Japanese facilities where high-dose-rate brachytherapy devices were available in 2012. In total, 135 responses were returned for analysis. Fifty-one facilities had acquired some sort of 3D imaging modality with applicator insertion, and computed tomography (CT) and magnetic resonance imaging (MRI) were used in 51 and 3 of the facilities, respectively. For actual treatment planning, X-ray films, CT and MRI were used in 113, 20 and 2 facilities, respectively. Among 43 facilities where X-ray films and CT or MRI were acquired with an applicator, 29 still used X-ray films for actual treatment planning, mainly because of limited time and/or staffing. In a follow-up survey 2.5 years later, respondents included 38 facilities that originally used X-ray films alone but had indicated plans to adopt 3D-IGBT. Of these, 21 had indeed adopted CT imaging with applicator insertion. In conclusion, 3D-IGBT (mainly CT) was implemented in 22 facilities (16%) and will be installed in 72 (53%) facilities in the future. Limited time and staffing were major impediments. PMID:26265660

  16. Fluid Substitution Modeling to Determine Sensitivity of 3D Vertical Seismic Profile Data to Injected CO­2­ at an active Carbon Capture, Utilization and Storage Project, Farnsworth field, TX.

    Haar, K. K.; Balch, R. S.

    2015-12-01

    The Southwest Regional Partnership on Carbon Sequestration monitors a CO2 capture, utilization and storage project at Farnsworth field, TX. The reservoir interval is a Morrowan age fluvial sand deposited in an incised valley. The sands are between 10 to 25m thick and located about 2800m below the surface. Primary oil recovery began in 1958 and by the late 1960's secondary recovery through waterflooding was underway. In 2009, Chaparral Energy began tertiary recovery using 100% anthropogenic CO2 sourced from an ethanol and a fertilizer plant. This constitutes carbon sequestration and fulfills the DOE's initiative to determine the best approach to permanent carbon storage. One purpose of the study is to understand CO­2 migration from injection wells. CO2­ plume spatial distribution for this project is analyzed with the use of time-lapse 3D vertical seismic profiles centered on CO2 injection wells. They monitor raypaths traveling in a single direction compared to surface seismic surveys with raypaths traveling in both directions. 3D VSP surveys can image up to 1.5km away from the well of interest, exceeding regulatory requirements for maximum plume extent by a factor of two. To optimize the timing of repeat VSP acquisition, the sensitivity of the 3D VSP surveys to CO2 injection was analyzed to determine at what injection volumes a seismic response to the injected CO­2 will be observable. Static geologic models were generated for pre-CO2 and post-CO2 reservoir states through construction of fine scale seismic based geologic models, which were then history matched via flow simulations. These generated static states of the model, where CO2­ replaces oil and brine in pore spaces, allow for generation of impedance volumes which when convolved with a representative wavelet generate synthetic seismic volumes used in the sensitivity analysis. Funding for the project is provided by DOE's National Energy Technology Laboratory (NETL) under Award No. DE-FC26-05NT42591.

  17. 3D seismic response of a limited valley via BEM using 2.5D analytical Green's functions for an infinite free-rigid layer

    António, Julieta; Tadeu, António

    2002-01-01

    This paper presents analytical solutions for computing the 3D displacements in a flat solid elastic stratum bounded by a rigid base, when it is subjected to spatially sinusoidal harmonic line loads. These functions are also used as Greens functions in a boundary element method code that simulates the seismic wave propagation in a confined or semi-confined 2D valley, avoiding the discretization of the free and rigid horizontal boundaries.

  18. 3D seismic analysis of the Collyhurst Sandstone: implications for CO2 sequestration in the East Irish Sea Basin

    Gamboa, Davide; Williams, John; Kirk, Karen; Gent, Christopher; Bentham, Michelle; Fellgett, Mark; Schofield, David

    2016-04-01

    Carbon Capture and Storage (CCS) is a vital technology towards low-carbon energy resources and the mitigation of global warming trends induced by rising CO2 levels in the atmosphere. The East Irish Sea Basin (EISB) is a key area for CCS in the western UK, having high CO2 storage potentials in explored hydrocarbon fields and in saline aquifers within the Permo-Triassic Sherwood Sandstone Formation. However, the theoretical storage potential of the EISB could be poorly estimated as the reservoir-prone Lower Permian formations are not considered in detail by current estimations. This work aims to fill this gap, focusing on the characterisation of the Lower Permian Collyhurst Sandstone Formation as a viable storage unit. The potential for CO2 storage is estimated as the total volume/area of suitable closures that are isolated by structural traps, occurring at depths suitable for CO2 injection and containment (>800m). Detailed structural and stratigraphic interpretations were made using 3D seismic data to assess the storage potential of the Collyhurst Sandstone Formation in the southern EISB. The basin strata is compartmentalised by numerous N-S trending faults. A higher degree of compartmentalisation occurs within regional anticlines where elongated tilted blocks are observed, bound by predominantly west-dipping faults that induce a variable offset of the Collyhurst Sandstone strata. Contrastingly, higher lateral continuity of this formation is observed within graben basins were faults are less frequent and with minor offset, thus potentially creating larger storage closures. Fault dip orientation in the grabens is variable, with west and east dipping faults occurring as a function of large east-dipping listric faults. This study was complemented by the stress modelling of the interpreted faults in order to assess the risk of CO2 leakage. Analysis of borehole breakouts observed in four approximately vertical wells in the EISB suggest a maximum horizontal stress

  19. Exposure to seismic survey alters blue whale acoustic communication

    Di Iorio, Lucia; Clark, Christopher W.

    2009-01-01

    The ability to perceive biologically important sounds is critical to marine mammals, and acoustic disturbance through human-generated noise can interfere with their natural functions. Sounds from seismic surveys are intense and have peak frequency bands overlapping those used by baleen whales, but evidence of interference with baleen whale acoustic communication is sparse. Here we investigated whether blue whales (Balaenoptera musculus) changed their vocal behaviour during a seismic survey th...

  20. Earthquake relocation using a 3D a-priori geological velocity model from the western Alps to Corsica: Implication for seismic hazard

    Béthoux, Nicole; Theunissen, Thomas; Beslier, Marie-Odile; Font, Yvonne; Thouvenot, François; Dessa, Jean-Xavier; Simon, Soazig; Courrioux, Gabriel; Guillen, Antonio

    2016-02-01

    The region between the inner zones of the Alps and Corsica juxtaposes an overthickened crust to an oceanic domain, which makes difficult to ascertain the focal depth of seismic events using routine location codes and average 1D velocity models. The aim of this article is to show that, even with a rather lose monitoring network, accurate routine locations can be achieved by using realistic 3D modelling and advanced location techniques. Previous earthquake tomography studies cover the whole region with spatial resolutions of several tens of kilometres on land, but they fail to resolve the marine domain due to the absence of station coverage and sparse seismicity. To overcome these limitations, we first construct a 3D a-priori P and S velocity model integrating known geophysical and geological information. Significant progress has been achieved in the 3D numerical modelling of complex geological structures by the development of dedicated softwares (e.g. 3D GeoModeller), capable at once of elaborating a 3D structural model from geological and geophysical constraints and, possibly, of refining it by inversion processes (Calcagno et al., 2008). Then, we build an arrival-time catalogue of 1500 events recorded from 2000 to 2011. Hypocentres are then located in this model using a numerical code based on the maximum intersection method (Font et al., 2004), updated by Theunissen et al. (2012), as well as another 3D location technique, the NonLinLoc software (Lomax and Curtis, 2001). The reduction of arrival-time residuals and uncertainties (dh, dz) with respect to classical 1D locations demonstrates the improved accuracy allowed by our approach and confirms the coherence of the 3D geological model built and used in this study. Our results are also compared with previous works that benefitted from the installation of dense temporary networks surrounding the studied epicentre area. The resulting 3D location catalogue allows us to improve the regional seismic hazard assessment

  1. Historical Buildings Models and Their Handling via 3d Survey: from Points Clouds to User-Oriented Hbim

    Chiabrando, F.; Sammartano, G.; Spanò, A.

    2016-06-01

    This paper retraces some research activities and application of 3D survey techniques and Building Information Modelling (BIM) in the environment of Cultural Heritage. It describes the diffusion of as-built BIM approach in the last years in Heritage Assets management, the so-called Built Heritage Information Modelling/Management (BHIMM or HBIM), that is nowadays an important and sustainable perspective in documentation and administration of historic buildings and structures. The work focuses the documentation derived from 3D survey techniques that can be understood like a significant and unavoidable knowledge base for the BIM conception and modelling, in the perspective of a coherent and complete management and valorisation of CH. It deepens potentialities, offered by 3D integrated survey techniques, to acquire productively and quite easilymany 3D information, not only geometrical but also radiometric attributes, helping the recognition, interpretation and characterization of state of conservation and degradation of architectural elements. From these data, they provide more and more high descriptive models corresponding to the geometrical complexity of buildings or aggregates in the well-known 5D (3D + time and cost dimensions). Points clouds derived from 3D survey acquisition (aerial and terrestrial photogrammetry, LiDAR and their integration) are reality-based models that can be use in a semi-automatic way to manage, interpret, and moderately simplify geometrical shapes of historical buildings that are examples, as is well known, of non-regular and complex geometry, instead of modern constructions with simple and regular ones. In the paper, some of these issues are addressed and analyzed through some experiences regarding the creation and the managing of HBIMprojects on historical heritage at different scales, using different platforms and various workflow. The paper focuses on LiDAR data handling with the aim to manage and extract geometrical information; on

  2. Comment on 'The role of 3-D interactive visualization in blind surveys of HI in galaxies'

    Taylor, Rhys

    2015-01-01

    Punzo et al. (2015) recently reported on the state of the art for visualisation of H I data cubes. I here briefly describe another program, FRELLED, specifically designed for dealing with H I data. Unlike many 3D viewers, FRELLED can handle astronomical world coordinates, easily and interactively mask and label specific volumes within the data, overlay optical data from the SDSS, generate contour plots and renzograms, make basic spectral profile measurements via an interface with MIRIAD, and can switch between viewing the data in 3D and 2D. The code is open source and can potentially be extended to include any astronomical function possible with Python, displaying the result in an interactive 3D environment.

  3. Vision-Based Long-Range 3D Tracking, applied to Underground Surveying Tasks

    Mossel, Annette; Gerstweiler, Georg; Vonach, Emanuel; Kaufmann, Hannes; Chmelina, Klaus

    2014-04-01

    To address the need of highly automated positioning systems in underground construction, we present a long-range 3D tracking system based on infrared optical markers. It provides continuous 3D position estimation of static or kinematic targets with low latency over a tracking volume of 12 m x 8 m x 70 m (width x height x depth). Over the entire volume, relative 3D point accuracy with a maximal deviation ≤ 22 mm is ensured with possible target rotations of yaw, pitch = 0 - 45° and roll = 0 - 360°. No preliminary sighting of target(s) is necessary since the system automatically locks onto a target without user intervention and autonomously starts tracking as soon as a target is within the view of the system. The proposed system needs a minimal hardware setup, consisting of two machine vision cameras and a standard workstation for data processing. This allows for quick installation with minimal disturbance of construction work. The data processing pipeline ensures camera calibration and tracking during on-going underground activities. Tests in real underground scenarios prove the system's capabilities to act as 3D position measurement platform for multiple underground tasks that require long range, low latency and high accuracy. Those tasks include simultaneously tracking of personnel, machines or robots.

  4. 3D SURVEY AND AUGMENTED REALITY FOR CULTURAL HERITAGE. THE CASE STUDY OF AURELIAN WALL AT CASTRA PRAETORIA IN ROME

    M. Canciani

    2016-06-01

    Full Text Available The development of close-range photogrammetry has produced a lot of new possibility to study cultural heritage. 3D data acquired with conventional and low cost cameras can be used to document, investigate the full appearance, materials and conservation status, to help the restoration process and identify intervention priorities. At the same time, with 3D survey a lot of three-dimensional data are collected and analyzed by researchers, but there are a very few possibility of 3D output. The augmented reality is one of this possible output with a very low cost technology but a very interesting result. Using simple mobile technology (for iPad and Android Tablets and shareware software (in the case presented “Augment” it is possible to share and visualize a large number of 3D models with your own device. The case study presented is a part of an architecture graduate thesis, made in Rome at Department of Architecture of Roma Tre University. We have developed a photogrammetric survey to study the Aurelian Wall at Castra Praetoria in Rome. The surveys of 8000 square meters of surface have allowed to identify stratigraphy and construction phases of a complex portion of Aurelian Wall, specially about the Northern door of Castra. During this study, the data coming out of 3D survey (photogrammetric and topographic, are stored and used to create a reverse 3D model, or virtual reconstruction, of the Northern door of Castra. This virtual reconstruction shows the door in the Tiberian period, nowadays it's totally hidden by a curtain wall but, little and significative architectural details allow to know its original feature. The 3D model of the ancient walls has been mapped with the exact type of bricks and mortar, oriented and scaled according to the existing one to use augmented reality. Finally, two kind of application have been developed, one on site, were you can see superimposed the virtual reconstruction on the existing walls using the image

  5. Computing 3-D wavefields in mantle circulations models to test hypotheses on the origin of lower mantle heterogeneity under Africa directly against seismic observations

    Schuberth, Bernhard; Zaroli, Christophe; Nolet, Guust

    2015-04-01

    Of particular interest for the tectonic evolution of the Atlantic region is the influence of lower mantle structure under Africa on flow in the upper mantle beneath the ocean basin. Along with its Pacific counterpart, the large African anomaly in the lowermost mantle with strongly reduced seismic velocities has received considerable attention in seismological and geodynamic studies. Several seismological observations are typically taken as an indication that these two anomalies are being caused by large-scale compositional variations and that they are piles of material with higher density than normal mantle rock. This would imply negative buoyancy in the lowermost mantle under Africa, which has important implications for the flow at shallower depth and inferences on the processes that led to the formation of the Atlantic Ocean basin. However, a large number of recent studies argue for a strong thermal gradient across the core-mantle boundary that might provide an alternative explanation for the lower mantle anomaly through the resulting large lateral temperature variations. Recently, we developed a new joint forward modeling approach to test such geodynamic hypotheses directly against the seismic observations: Seismic heterogeneity is predicted by converting the temperature field of a high-resolution 3-D mantle circulation model into seismic velocities using thermodynamic models of mantle mineralogy. 3-D global wave propagation in the synthetic elastic structures is then simulated using a spectral element method. Being based on forward modelling only, this approach allows us to generate synthetic wavefields and seismograms independently of seismic observations. The statistics of observed long-period body wave traveltime variations show a markedly different behaviour for P- and S-waves: the standard deviation of P-wave delay times stays almost constant with ray turning depth, while that of the S-wave delay times increases strongly throughout the mantle. In an

  6. Evaluating integration of inland bathymetry in the U.S. Geological Survey 3D Elevation Program, 2014

    Miller-Corbett, Cynthia

    2016-01-01

    Inland bathymetry survey collections, survey data types, features, sources, availability, and the effort required to integrate inland bathymetric data into the U.S. Geological Survey 3D Elevation Program are assessed to help determine the feasibility of integrating three-dimensional water feature elevation data into The National Map. Available data from wading, acoustic, light detection and ranging, and combined technique surveys are provided by the U.S. Geological Survey, National Oceanic and Atmospheric Administration, U.S. Army Corps of Engineers, and other sources. Inland bathymetric data accessed through Web-hosted resources or contacts provide useful baseline parameters for evaluating survey types and techniques used for collection and processing, and serve as a basis for comparing survey methods and the quality of results. Historically, boat-mounted acoustic surveys have provided most inland bathymetry data. Light detection and ranging techniques that are beneficial in areas hard to reach by boat, that can collect dense data in shallow water to provide comprehensive coverage, and that can be cost effective for surveying large areas with good water clarity are becoming more common; however, optimal conditions and techniques for collecting and processing light detection and ranging inland bathymetry surveys are not yet well defined.Assessment of site condition parameters important for understanding inland bathymetry survey issues and results, and an evaluation of existing inland bathymetry survey coverage are proposed as steps to develop criteria for implementing a useful and successful inland bathymetry survey plan in the 3D Elevation Program. These survey parameters would also serve as input for an inland bathymetry survey data baseline. Integration and interpolation techniques are important factors to consider in developing a robust plan; however, available survey data are usually in a triangulated irregular network format or other format compatible with

  7. Ryukyu Subduction Zone: 3D Geodynamic Simulations of the Effects of Slab Shape and Depth on Lattice-Preferred Orientation (LPO) and Seismic Anisotropy

    Tarlow, S.; Tan, E.; Billen, M. I.

    2015-12-01

    At the Ryukyu subduction zone, seismic anisotropy observations suggest that there may be strong trench-parallel flow within the mantle wedge driven by complex 3D slab geometry. However, previous simulations have either failed to account for 3D flow or used the infinite strain axis (ISA) approximation for LPO, which is known to be inaccurate in complex flow fields. Additionally, both the slab depth and shape of the Ryukyu slab are contentious. Development of strong trench-parallel flow requires low viscosity to decouple the mantle wedge from entrainment by the sinking slab. Therefore, understanding the relationship between seismic anisotropy and the accompanying flow field will better constrain the material and dynamic properties of the mantle near subduction zones. In this study, we integrate a kinematic model for calculation of LPO (D-Rex) into a buoyancy-driven, instantaneous 3D flow simulation (ASPECT), using composite non-Newtonian rheology to investigate the dependence of LPO on slab geometry and depth at the Ryukyu Trench. To incorporate the 3D flow effects, the trench and slab extends from the southern tip of Japan to the western edge of Taiwan and the model region is approximately 1/4 of a spherical shell extending from the surface to the core-mantle boundary. In the southern-most region we vary the slab depth and shape to test for the effects of the uncertainties in the observations. We also investigate the effect of adding locally hydrated regions above the slab that affect both the mantle rheology and development of LPO through the consequent changes in mantle flow and dominate (weakest) slip system. We characterize how changes in the simulation conditions affect the LPO within the mantle wedge, subducting slab and sub-slab mantle and relate these to surface observations of seismic anisotropy.

  8. Determination of porosity and facies trends in a complex carbonate reservoir, by using 3-D seismic, borehole tools, and outcrop geology

    Zacharakis, T.G. Jr.; Comet, J.N.; Murillo, A.A. [Respol Exploracion, S.A., Madrid (Spain)] [and others

    1996-08-01

    Mesozoic carbonate reservoirs are found in the Mediterranean Sea, off the east coast of Spain. A wide variation of porosities are found in the core samples and logs: vuggy, breccia, fractures, and cavern porosity. In addition, complex Tertiary carbonate geometries include olistostromes, breccia bodies, and reef buildups, which are found on top of Mesozoic carbonates. Predicting the porosity trends within these oil productive reservoirs requires an understanding of how primary porosity was further enhanced by secondary processes, including fractures, karstification, and dolomitization in burial conditions. Through an extensive investigation of field histories, outcrop geology, and seismic data, a series of basic reservoir styles have been identified and characterized by well log signature and seismic response. The distribution pattern of the different reservoirs styles is highly heterogeneous, but by integrating subsurface data and outcrop analogs, it is possible to distinguish field-scale and local patterns of both vertical and local variations in reservoir properties. Finally, it is important to quantify these reservoir properties through the study of seismic attributes, such as amplitude variations, and log responses at the reservoir interval. By incorporating 3-D seismic data, through the use of seismic inversion, it is possible to predict porosity trends. Further, the use of geostatistics can lead to the prediction of reservoir development within the carbonate facies.

  9. Aftershock Sequences Modeled with 3-D Stress Heterogeneity and Rate-State Seismicity Equations: Implications for Crustal Stress Estimation

    Smith, Deborah Elaine; Dieterich, James H.

    2010-01-01

    In this paper, we present a model for studying aftershock sequences that integrates Coulomb static stress change analysis, seismicity equations based on rate-state friction nucleation of earthquakes, slip of geometrically complex faults, and fractal-like, spatially heterogeneous models of crustal stress. In addition to modeling instantaneous aftershock seismicity rate patterns with initial clustering on the Coulomb stress increase areas and an approximately 1/t diffusion back to the pre-mains...

  10. Development of a New Parallel Code for 3-D Joint Refraction and Reflection Travel-Time Tomography of Wide-Angle Seismic Data. Synthetic and Real Data Applications to the Study of Subduction Zones

    Meléndez, Adrià

    2014-01-01

    [eng] This dissertation is devoted to seismic tomography. I have implemented a new modelling tool for 3-D joint refraction and reflection travel-time tomography of wide-angle seismic data (TOMO3D). The reason behind this central objective is the evidence that the information based on 2-D seismic data does not allow to capture the structural complexity of many 3-D targets, and in particular that of the seismogenic zone in subduction margins. The scientific rationale for this statement, which j...