Velocity-Interface Structure of the Southwestern Ryukyu Subduction Zone from EW9509-1 OBS/MCS Data

A wide-angle seismic survey, combining ocean-bottom seismometers (OBS) and multi-channel seismic (MCS) profiling, was implemented in the southwestern Ryukyu subduction zone during August and September 1995. In this paper, we present the data analysis of eight OBSs and the corresponding MCS line along profile EW9509-1 from this experiment. Seismic data modeling includes identification of refracted and reflected arrivals, initial model building from velocity analysis of the MCS data, and simultaneous and layer-stripping inversions of the OBS and MCS arrivals. The velocity-interface structure constructed along profile EW9509-1 shows that the northward subduction of the Philippine Sea Plate has resulted in a northward thickening of the sediments of the Ryukyu Trench and the Yaeyama accretionary wedge north of the trench. The boundary between the subducting oceanic crust and the overriding continental crust (represented by a velocity contour of 6.75 km/s) and a sudden increase of the subducting angle (from 5 degrees to 25 degrees) are well imaged below the Nanao Basin. Furthermore, velocity undulation and interface variation are found within the upper crust of the Ryukyu Arc. Therefore, the strongest compression due to subduction and a break-off of the slab may have occurred and induced the high seismicity in the forearc region. This revised version was published online in November 2006 with corrections to the Cover Date.     

井中激电地-井方式井旁球体正反演

以激发极化法中的地-井方式中的体极化球形体作为研究对象,采用解析法对地下半空间的电位场进行求解,利用提取出的井中二次场异常电位差进行反演,得出极化球体的球心离井的距离和球心的埋深。通过计算机进行正反演模拟,验证了进行正反演计算的解析式的正确性,并从物理意义上对正反演过程中出现的现象进行分析。最后给出该方法的适用条件以及注意事项。研究对激发极化法地下勘探目标的定量分析具有明确意义。     

Probabilistic modeling of seismically triggered landslides using Monte Carlo simulations

The 2004 Chuetsu earthquakes of Niigata (Japan) triggered numerous landslides, and the most widespread types of landslides were highly disrupted, relatively shallow slides and soil (debris) flows. This paper presented a method to evaluate slope instability using Newmark displacement on a pixel-by-pixel basis in a given area. The proposed method was able to integrate Newmark displacement modeling and Monte Carlo simulations within geographical information systems. In the modeling, an empirical attenuation relationship was utilized to calculate Arias intensity over this study area, and the variability of geotechnical parameters was taken into account to calculate coseismic landslide displacement. Before deriving the displacement from related inputs, the Monte Carlo simulations ran 1,500 times and generated 1,500 displacement values for each grid cell, and then means and standard deviations of displacement were calculated and probabilistic distributions can be obtained. Finally, given 10 cm as a threshold value of displacement, estimated probabilities of displacement exceeding 10 cm were shown as a map of seismic landslide hazards. The resulting hazard map was classified into four categories from very low to high level.     

顾及DEM误差自相关的坡度计算模型精度分析

基于DEM的坡度计算,其误差来源于DEM误差、DEM结构和坡度计算模型。在顾及DEM误差自相关的前提下,对四种DEM坡度计算模型进行了分析和评价。研究表明,三阶不带权差分能给出较高的坡度计算精度;在局部窗口中,格网点数量越多,坡度计算越准确;等权比不等权的坡度计算模型更准确;DEM误差自相关结构形式对坡度计算无影响。进一步的理论分析和试验分析还表明：DEM误差自相关性的存在,不仅能够改善地形分析的精度,也能改善DEM自身精度。     

地面核磁共振模型约束反演含水层参数

导电模型的地面核磁共振感应电动势是含水量的非线性函数.引入模型约束的迭代反演方法求解该非线性问题的反演.在反演过程中,根据理论公式计算迭代过程中灵敏度矩阵,并采用平坦模型和光滑模型两种约束.对均匀半空间、层状导电模型和实际数据进行了反演模拟,结果表明,模型约束迭代反演方法能从地面核磁共振感应电动势获得含水层较为合理的含水量及分布,且该结果可以从作为初始模型的均匀含水量分布反演得到.对无噪音数据,平坦模型和光滑模型约束对反演的含水量分布影响不大;但当数据存在一定的噪音时,平坦模型约束将比光滑模型约束获得更为精确的含水层参数.     

小波边缘分析建模波阻抗反演方法在煤层解释中的应用

受限于地震波主波长的1／4,常规地震处理难以分辨较薄煤层。基于模型的波阻抗反演虽然克服了厚、薄煤层计算问题．但因难以得到准确的高频部分,很多情况下反演结果的精度和分辨率往往不能满足地质预测的要求。小波边缘分析建模波阻抗反演方法是从地震数据中提取地震特征参数,用于建模并参与迭代反演。该方法即可以弥补反演过程中高频成分的误差,减少对钻井数据的依赖程度,得到较为合理的初始模型,又可提高反演的精度和分辨率．使反演结果更好的反映实际地下地质情况。利用该方法,对大井-将军庙及新疆红沙泉2个勘探区的侏罗系中统西山窑组含煤地层进行波阻抗反演,反演结果与钻井资料进行比较,其中大井一将军庙0802及0809号孔单煤层厚度误差分别为3％、8％．新疆红沙泉多煤层厚度及煤层之间的夹层厚度预测结果除B6煤层误差较大外,其它误差均小于7％．可见该反演方法的计算结果比较准确。     

努力提高三维地震勘探精度，为西部煤炭工业做出新贡献

以国家重大产业技术开发专项“西部煤炭资源高精度三维地震勘探技术”项目的由来、意义和总体研究目标为引,概括的介绍了项目依托工程中各个专项技术研究完成情况,并对非均匀介质成像技术、高精度三维地震静校正技术、高密度采集技术、特观技术、岩性反演技术、属性体解释技术等六项重大关键技术取得的突破性进展进行了重点说明。指出随着我国煤炭生产重点的逐步西移,应加强诸如叠前、叠后深度偏移技术的研究,以解决复杂山区三维地震面元内地震反射波散射问题,提高其三维地震勘探精度,为西部煤炭工业做出新贡献!     

SOME NEW UNDERSTANDINGS OF LAND-SURFACE PROCESSES IN ARID AREA FROM THE HEIFE*

The survey for the HEIFE(Atmosphere-Land Surface Processes Experiment at Heihe River Basin,Western China) is given in the paper.The following basic subjects for land-surface processes in arid areas are studied:(1) the general characteristics of the energy budget on ground surface in arid areas;(2) the parameterization of the land surface processes;(3) the interaction between oasis and its desert circumstances,a special phenomenon in arid areas.The analysis shows that the sensible heat flux in the surface energy budget is in the majority,and the latent heat flux may be neglected.The influence of atmospheric stratification stability on the turbulent transfer of energy and substance must be considered in parameterization of land surface processes in arid areas.The "cold island effect" phenomenon in oasis and the "humidity inversion" phenomenon in desert near oasis are the result of the interaction between them.The results would improve the understanding of land surface processes in arid areas.     

Information structure analysis for quantitative assessment of mineral resources and the discovery of a silver deposit

We propose an information-structure-analysis (ISA) method to quantify the correlations between quantitative and qualitative variables as well as within each type of variable. This method is applied to the evaluation of mineral resources in the western Zheijiang Province of China. The district contains a number of silver-bearing Fe–Cu–Pb–Zn mineral deposits near igneous complexes and FeCuPbZn zones away from the complexes. Silver anomalies occur not only in the known Fe–Cu–Zn–Pb deposits, but also in the country rock, suggesting the possible existence of silver deposits far from the igneous complexes.The tonnage distribution of silver is modeled by Monte Carlo simulation. This simulation is conducted on the basis of the correlations between silver (Ag) and lead (Pb), since no known data on silver is available. The known tonnage distribution of lead in 11 control cells was used to approximate the tonnage distribution of silver in the Monte Carlo simulation. With ISA and Monte Carlo methods, the total amount of potential polymetallic resources in 49 cells in the western Zhejiang Provice is predicted. Significantly, a deposit with about 24 tonnes of silver has been found within our exploration target area.     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.