叠前正演模拟与偏移的网络并行计算

在叠前深度偏移和非零炮检距声波方程正演计算过程中包含了大量的可并行计算的成分。作者在本文中提出叠前正演模拟与偏移的网络并行计算算法，并基于TCP／IP协议，将该算法设计成网络并行处理程序，极大地提高了计算效率。实际运算结果证明，本文提出的并行算法和技术路线是切实可行的。     

Hydrogeologic studies for nuclear-waste disposal in Sweden

Safety Report 1997 (SR 97) of the Swedish Nuclear Fuel and Waste Management Company is a comprehensive performance assessment of three hypothetical radioactive waste repositories in Sweden. It includes hydrogeologic studies (data compilation, parameter synthesis, and groundwater flow modeling) to determine groundwater flow and the associated uncertainties for the three sites. This report reviews and compares the hydrogeologic characterization programs, the site characteristics, and the groundwater flow models used in the SR 97 performance assessment. Although differences in site-characterization programs tend to mask differences in site characteristics, the sites do have notable differences that affect the results of the performance assessment. The effects of model uncertainties evaluated by the variant cases appear to be smaller than the variability of results for the base case (best estimate of site conditions) of each site. Electronic Publication     

Two-dimensional geostatistical modeling and prediction of the fracture system in the Albala Granitic Pluton, SW Iberian Massif, Spain

We use a fracture index distribution method of geostatistical modelling and prediction to characterize quantitatively the fracture system in two-dimensions (2D) in the Mina Ratones area, located in the Albalá Granitic Pluton (SW Iberian Massif). The fracture index (FI) is a quantitative estimate of the fracture density in discrete domains. To validate the results of geostatistical modeling a detailed structural map of the area was also made on a scale of 1:1000. The resulting grids, expressed as pixel-maps, describe the continuous value of the FI in 2D for the whole Mina Ratones area. Based on the modelled distribution of the FI and their correlation with mapped faults, we distinguish two structural domains in the studied area: elongated bands of fracture zones with high FI values and romboidal blocks located between them with low FI values. The separation between both domains is gradual. Though a threshold value of the FI that separate both structural domains is not clearly defined, the fracture zones generally present FI>1 and the individualized blocks FI<0.50. As a consequence, the obtained grid of the FI permits the quantitative structural classification of the granitic massif in 2D and understanding fault zone architecture in the Mina Ratones area.     

Groundwater resources evaluation case study via discrete fracture flow modeling

The Arizona Department of Transportation (ADOT) is preparing to upgrade State Route 260 between Payson and Heber. It is estimated that a total of about one million cubic meters of water will be required for embankment construction during a period of about 84 months to upgrade the first 33.8 km of the highway. ADOT is investigating various sources of construction water for use in the highway improvement project, including groundwater resources along the highway corridor. A region known as the RV site, underlain by fractured granite, is located 12.9 km east of Payson. The site includes three springs, a creek and several wells. Several boreholes and observation wells were made to a maximum depth of 157 m to obtain fracture data and to conduct pumping tests with monitoring. Fracture data recorded by acoustic televiewer logs were used to build a fracture network model for the rock mass. Results of a 24-hour and a 7-day pumping tests were used to calibrate hydraulic parameters of a finite element discrete fracture fluid flow model considering the region as a heterogeneous, anisotropic, fractured medium. A 38-day multi-well pumping test was used to validate the calibrated numerical model. The calibrated model showed the capability to provide reasonably accurate predictions for new pumping tests conducted in the same well field. The validated model was then used to simulate pumping exceeding a 7 year period under different scenarios incorporating different sets of boundary conditions and different pumping rates at multi-wells, with and without recharge, to evaluate the yield of the aquifer and to assess the effect of long-term pumping on the environment. The results indicated that (a) the combined yield of the wells in the RV site is sufficient to meet the water demand for the ADOT highway project and (b) the water levels in the well field would decline between 3.0 and 7.6 m after one year of pumping and by 12.2 to over 30.5 m during the life of the project.     

Modeling,ray tracing,and block nonlinear travel-time inversion in 3D

We describe an integrated forward and inverse three-dimensional modeling system that can deal with complex geological structures. The system has been designed to handle large-scale problems by using a distributed approach. It uses seismic ray tracing for forward simulation, time-to-depth mapping, and nonlinear travel-time inversion.A novel decomposition method is our tool to attack large-scale problems in a parallel approach. The system is fully implemented and we demonstrate its performance with synthetic examples.This research was partially supported by the National Science Foundation under SBIR Grant III-9300992.     

Thermal evolution of the Ryoke metamorphic belt, southwestern Japan: Tectonic and numerical modeling

Abstract The Ryoke metamorphic belt of southwestern Japan is composed of Cretaceous Ryoke granitoids and associated metamorphic rocks of low-pressure facies series. The Ryoke granitoids are divided into sheet-like bodies (e.g. Gamano granodiorite) and stock-like bodies. The Gamano granodiorite intruded concordantly into the high-grade metamorphic rocks without development of a contact metamorphic aureole, and the intrusion ages of the granodiorite are similar to the ages of thermal peak of the low pressure (low-P) metamorphism. It is suggested that the low-P Ryoke metamorphism resulted from the intrusion of the Gamano granodiorite. In this study, a simple 1-D numerical model of conductive heat transfer was used to evaluate the thermal effects of emplacement of the Gamano granodiorite. Calculated temperature-time ( T-t ) paths are characterized by a rapid increase of metamorphic temperature and a relatively short-lived period of high temperature. For example, the T-t path at the 15-km depth is characterized by a rapid average increase in temperature of 1.4 × 10-³°C/year and high temperatures for < ca 0.5 Ma. The calculated peak temperature for each depth is nearly equal to the petrologically estimated value for each correlated metamorphic zone. The results suggest that the magma-intrusion model is one possible thermal model for low-pressure facies series metamorphism.     

用边界单元法模拟二维地形对大地电磁场的影响

本文用边界单元法模拟二维地形对大地电磁场的影响.首先用格林公式将二维大地电磁的边值问题转变成积分方程,然后用边界单元法解积分方程,得到地形上的大地电磁场和它的法向导数,由此可计算电阻率.与有限单元法相比,本方法剖分后的地形与实际地形的拟合程度高,向计算机输入原始数据的工作简单,可以在微机上计算地形引起的视电阻率.本方法适合于在野外生产现场进行大地电磁法的地形改正.计算表明,地形对H_x型波的影响比对E_x型波严重得多;随着周期的增长,地形对H_x型波的影响变得稳定.     

FORECASTING IMPACTS OF HYDRAULIC STRUCTURESON AN ALLUVIAL RIVER

In the last decades the construction of large hydraulic projects on alluvial rivers in China has spurred the development of both physical and one-dimensional mathematical modeling to the point that deposition and erosion predicted by modeling in either way conform to prototype observations over periods of 3 to 19 years. The prototype data are those of Gezhouba on the Yangtze River and Danjiangkou on the Han River. The accuracy attained is in the neighborhood of 20% which is considered good enough in sedimentation engineering. Thus these techniques of modeling may be applied to forecast the impacts of hydraulic structures on an alluvial river. Mathematical modeling in 2 space variables has also been developed, but still awaits further verification.     

The reservoir-associated earthquakes of April 1983 in western Thailand: Source modeling and implications for induced seismicity

On 22 April 1983, a very large area of Thailand and part of Burma were strongly shaken by a rare earthquake (m _b=5.8,M _s=5.9). The epicenter was located at the Srinagarind reservoir about 190 km northwest of Bangkok, a relatively stable continental region that experienced little previous seismicity. The main shock was preceded by some foreshocks and followed by numerous aftershocks. The largest foreshock ofm _b=5.2 occurred 1 week before the main shock, and the largest aftershock ofm _b=5.3 took place about 3 hours after the main shock. Focal mechanisms of the three largest events in this earthquake sequence have been studied by other seismologists using firts-motion data. However, the solutions for the main shock and the largest aftershock showed significant inconsistency with known surface geology and regional tectonics. We reexamined the mechanisms of these three events by using teleseismicP-andS-waveforms and through careful readings ofP-wave first motions. The directions of theP axes in our study range from NNW-SSE to NNE-SSW, and nodal planes strike in the NW-SE to about E-W in agreement with regional tectonics and surface geology. The main shock mechanism strikes 255°, dips 48°, and slips 63.5°. The fault motions during the main shock and the foreshock are mainly thrust, while the largest aftershock has a large strike-slip component. The seismic moment and the stress drop of the mainshock are determined to be 3.86×10²⁴ dyne-cm and 180 bars, respectively. The occurrence of these thrust events appears to correlate with the unloading of the Srinagarind reservoir. The focal depths of the largest foreshock, the main shock, and the largest aftershock are determined to be 5.4 km, 8 km, and 22.7 km, respectively, from waveform modeling and relative location showing a downward migration of hypocenters of the three largest events during the earthquake sequence. Other characteristics of this reservoir-induced earthquake sequence are also discussed.     

A NEW RECEPTOR MODEL:A DIRECT TRILINEAR DECOMPOSITION FOLLOWED BY A MATRIX RECONSTRUCTION

In many cases,monitoring data for ambient airborne particles can be organized in the form of a three-way data table with one way for chemical species,one for sampling periods and one for sites.A directtrilinear decomposition followed by a matrix reconstruction(DTDMR)is developed to analyze such adata table as a whole.The three-way data set is composed into three two-way matrices by a direct trilineardecomposition(DTD).The column vectors of each of the matrices are called‘source profiles’,‘emissionpatterns’and‘site coefficients’respectively.Particulate sources are identified by examining both theirsource profiles and emission patterns.After the sources have been identified,emission patterns and sitecoefficients are used to produce a three-way matrix that gives estimates of mass contributions of sourcesto the samples collected at every site in every period.By simulation study,not only has the method beenverified,but a good indicator has been found that shows the number of factors(i.e.sources)in thesystem.Unlike other receptor models,DTDMR does not require source profile data and does not involvetrial-and-error procedures.Since DTDMR identifies sources based on variations in two dimensions,it hasa higher potential to distinguish two sources that have similar chemical compositions.The DTDMRmodel has provided excellent results with simulated data and has been applied in a real world three-waydata set.