弹性波的三维有限元模拟

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.     

Unified approach to amplitude attenuation and coda excitation in the randomly inhomogeneous lithosphere

A unified model is proposed for explaining the frequency dependent amplitude attenuation and the coda wave excitation on the basis of the single scattering process in the randomly inhomogeneous lithosphere. Adopting Birch's law and a direct proportion between density and wave velocity, we statistically describe the inhomogeneous medium by one random function characterized by the von Karman autocorrelation function. We calculate the amplitude attenuation from the solid angle integral of scattered wave energy on the basis of the Born approxiimation after subtracting the travel-time fluctuation effect caused by slowly varying velocity inhomogeneities. This subtraction is equivalent to neglect energy loss by scattering within a cone around the forward direction. The random inhomogeneity of the von Karman autocorrelation function of order 0.35 with the mean square fractional fluctuation of 7.2×10^–3 1.3×10^–2 and the correlation distance of 2.15.1 km well explains observed backward scattering coefficientg and the ratioQ _P ^–1 /Q _S ^–1 , and observed and partially conjecturedQ _S ^–1 for frequencies between 0.5 Hz and 30 Hz.     

Deep structure of the Nojima Fault, southwest Japan, estimated from borehole observations of fault-zone trapped waves

To estimate the deep structure of the southern part of the Nojima Fault, southwest Japan without the influence of near-surface structures, we analyzed the Love-wave-type fault-zone trapped waves (LTWs) recorded by a borehole seismometer at 1800 m depth. We examined the polarization, dispersion, and dominant frequency of the wavetrain following the direct S-wave in each seismogram to identify the LTW. We selected eight candidates for typical LTWs from 462 records. Because the duration of the LTW increases with hypocentral distance, we infer that the low velocity fault-zone of the Nojima Fault continues towards the seismogenic depth. In addition, since the duration of the LTW increases nonlinearly with hypocentral distance, we infer that the S-wave velocity of the fault-zone increases with depth. The location of events showing the LTW indicates that the fault-zone dips to the southeast at 75° and continues to a depth of approximately 10 km. We assumed a uniform low-velocity waveguide to estimate the average structure of the fault-zone. We estimated the average width, S-wave velocity, and Q_s of the fault-zone by comparing an analytical solution of the LTW with measured data. The average width, S-wave velocity, and Q_s of the fault-zone are 150 to 290 m, 2.5 to 3.2 km/s, and 40 to 90, respectively. Hence the fault-zone structure with a larger width and smaller velocity reduction than the fault-zone model estimated by previous surface observation is more suitable to represent the average fault-zone structure of the Nojima fault. The present study also indicated that the shallow layers and/or a shallow fault-zone structure drastically changes the characteristics of the LTW recorded at the surface, and therefore cause a discrepancy in the fault-zone model between the borehole observation and surface observation.     

冲击凿岩钎杆的瞬态动力学分析

以ANSYS为平台，按照真实工作状态对凿岩钎杆进行了瞬态动力学分析，并对不同活塞形状下分析结果数据进行了比较，为冲击凿岩系统设计提供了理论依据。     

Attenuation of aqueous metal transport at two natural acid rock drainage occurrences in the Yukon Territory, Canada

Naturally acidic drainage associated with pyritic black shale has been observed in many locations in the Yukon Territory. While not necessarily linked to known mineral deposits, most of these natural acid rock drainage occurrences show elevated dissolved concentrations of trace elements, especially zinc, nickel, copper, cadmium and arsenic. Based on field observations, microbial investigation, chemical analyses and geochemical modeling, the fate and transport of potentially deleterious elements at two natural acid drainage occurrences with slightly different settings are examined. The Macintosh Creek is a small, acidic stream （pH 2.98-3.40）, 2 km long, located in the Macmillan Pass area of east-central Yukon amidst known sedimentary exhalative massive sulfide mineralization but remains undisturbed by exploration activities. Its trace metal content is apparently derived from groundwater discharges, which gave as much as 5.0, 2.5, 0.7, 0.13 and 0.03 mg/L ofZn, Ni, Cu and As, respectively. Interaction and sorption reactions with algal mats, biofilms and iron oxyhydroxides appear to be the dominant mechanisms attenuating aqueous contaminant transport along the stream. Cryogenic precipitation further consolidates the ferricrete formation and reduces the mobility of the sorbed metals. The tributaries of the Engineering Creek along the Dempster Highway in northern Yukon drain through a series of dolomite, phyllite, argillite, limestone, black shale, sandstone and conglomerate with no known concentration of mineralization. In this area, the water chemistry fully reflects the local geology with acidic streams invariably associated with black shale occurrences. Groundwater seeps in the headwaters area of the km-180 Creek completely enclosed in black shale gave pH 3.0 and as much as 148, 39, 2.9 and 9.1 mg/L of Zn, Ni, Cu and As, respectively. Sorption with iron oxyhydroxide and organic matter appear to dominate the attenuation of contaminant transport along the stream. However, once entered into carbonate-dominated terrains, secondary carbonate minerals exercise additional geochemical control on the local water chemistry as a result of neutralization.     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

Elastic properties of granulite facies rocks of Mahabalipuram, Tamil Nadu, India

Compressional and shear wave velocities and attenuation measurements have been carried out in some of the borehole samples of acidic, basic and intermediate granulites of Mahabalipuram, Tamil Nadu, India. The results have been obtained at ambient conditions using ‘time-of-flight’ pulse transmission technique at 1.0 MHz frequency. The results show linear relationships between velocity and density, and velocity and attenuation properties of the rocks. The acidic granulites show lower velocities and higher attenuation than the intermediate and basic granulites. The average values of the Poisson’s ratio of acidic, intermediate and basic granulites have been found to be 0.210, 0.241 and 0.279 respectively. The variations in velocities and attenuation in these low porosity crystalline rocks are found to be strongly influenced by their mineral composition. The laboratory velocity data (extrapolated to high pressure) of the present study and the published field velocity data from deep seismic sounding studies indicate that these granulite facies rocks may belong to mid-crustal depths only.