论三维等时线叠加反偏移中的有关问题

利用几何射线理论和一种考虑了积分域边缘贡献的二维稳想法，研究了三维等时线叠加反偏移中的反偏移场，反偏移加权函及反偏移孔径。结果证明，在有限孔径条件下，反偏移场的结构比已发表的结果要复杂和多。即使在最简单的条件下，反偏移场也是由反偏移信号的反偏移噪声（孔径效应）叠加而成的。而且，反偏移场的结构强烈地依赖于反偏移孔径边缘的位置。为了消除孔径效应，应采用最佳反偏移孔径和针对有限孔径的加权函数。根据定义，最佳反偏移孔径为稳相点的第一和第二Fresnel区之和。在第一Fresnel区内，有限孔径和加权函数与文献中所给出的加权函数完全一致。而在第二Fresnel区内，有限孔径加权函数为常规的加权函数与一窗函数的乘积。在反偏移孔径的边缘上，有限孔径加权函恒取零值。     

Interpretation of seismic and potential field data from western New York State and Lake Ontario

Lithoprobe and industry seismic profiles have furnished evidence of major zones of easterly dipping Grenville deformed crust extending southwest from exposed Grenville rocks north of Lake Ontario. Additional constraints on subsurface structure limited to the postulated Clarendon–Linden fault system south of Lake Ontario are provided by five east–west reflection lines recorded in 1976. Spatial correlations between seismic structure and magnetic anomalies are described from both Lake Ontario and the newly reprocessed New York lines.In the Paleozoic to Precambrian upper crust, the New York seismic sections show: (1) An easterly thickening wedge of subhorizontal Paleozoic strata unconformably overlying a Precambrian basement whose surface has an apparent regional easterly dip of 1–2°. Minor apparent normal offsets, possibly on the order of tens of meters, occur within the Paleozoic section. The generally poorly reflective unconformity may be locally characterized by topographic relief on the order of 100 m; (2) Apparent local displacement on the order of 90 m at the level of the Black River Group diminishes upward to little or no apparent offset of Queenston Shale; (3) Within the limited seismic sections, there appears to be no evidence that the complete upper crustal section is vertically or subvertically offset; (4) Dipping structure in the Paleozoic strata (15° to 35°) resembles some underlying Precambrian basement elements; (5) The surface continuity of inferred faults constituting the Clarendon–Linden system is not strongly supported by the seismic data.Beneath the Paleozoic strata, the seismic sections show both linear and arcuate reflector geometry with easterly apparent dips of 15° to 35° similar to the deep structures imaged on seismic lines from nearby Lake Ontario and on Lithoprobe lines to the north. The similarity supports an extension of easterly dipping Central Metasedimentary Belt structures of the Grenville orogen from southern Ontario to beneath western New York State.From a comparison of the magnetic and gravity fields with the New York seismic sections, we suggest: (1) The largely nonmagnetic Paleozoic strata appear to contribute negligibly to magnetic anomalies. Seismically imaged fractures in the New York Paleozoic strata appear to lie mainly west of a positive gravity anomaly. The relationship between magnetic and gravity anomalies and the changes in the geometry of interpreted Precambrian structures remains enigmatic; (2) North to northeast trending curvilinear magnetic and gravity anomalies parallel, but are not restricted to the principal trend of the postulated Clarendon–Linden fault system. Paleozoic fractures of the Clarendon–Linden system may partly overlie a southward extension of the Composite Arc Belt boundary zone.     

Numerical Analysis of Space–Time Finite Element Formulations for Miscible Displacements

A finite element formulation is proposed to approximate a nonlinear system of partial differential equations, composed by an elliptic subsystem for the pressure–velocity and a transport equation (convection–diffusion) for the concentration, which models the incompressible miscible displacement of one fluid by another in a rigid porous media. The pressure is approximated by the classical Galerkin method and the velocity is calculated by a post-processing technique. Then, the concentration is obtained by a Galerkin/least-squares space–time (GLS/ST) finite element method. A numerical analysis is developed for the concentration approximation. Then, stability, convergence and numerical results are presented confirming the a priori error estimates.     

Structure sismique du socle paléozoı̈que du bassin des Doukkala,Môle côtier,Maroc occidental. Indication en faveur de l'existence d'une phase éo-varisqueSeismic structure of the Doukkala basin,Palaeozoic basement,western Morocco: a hint for an Eovariscan fold-and-thrust belt

Seismic profiles and well data from the Doukkala basin unravel the structure of the Palaeozoic basement and suggest that this coastal zone of western Morocco was affected by a compressive phase during the Frasnian. This resulted in the formation of upright, plurikilometric folds associated with reverse faults (North Doukkala), and of asymmetrical folds associated with mostly west verging ramps (South Doukkala). Folding involved all pre-Upper Frasnian formations and caused partial or total hiatus of Upper Frasnian–Strunian strata. This event can be correlated with the orogenic phase reported from more internal domains of the Morocco Hercynian belt, where it is referred to as the ‘Bretonne’ or ‘Eovariscan’ phase. To cite this article: H. Echarfaoui et al., C. R. Geoscience 334 (2002) 13–20     

震源时间函数与震源破裂过程

简要介绍了现代数字地震学的两个重要概念。一个是震源时间函数，另一个是地震震源的时空破裂过程。首先，从地震断层位移表示定量出发，介绍震源时间函数和震源破裂过程的基本概念。然后，分别介绍两种从远场地震记录中提取震源时间函数和获取有限断层面上时空破裂过程图像的反演方法。     

Modes of raising northeastern Tibet probed by explosion seismology

New wide-angle reflection and refraction seismic data provide constraints on the structure of the upper lithosphere, and test models of its evolution to raise the northeastern part of the Tibetan Plateau. Amplitudes observed for reflections from the crust-mantle boundary are sufficiently large to suggest that there is no significant partial melt in the deep crust. The data show an increase of the crustal thickness between terranes from north of the Kun Lun Fault into the Qang Tang of central Tibet, and a contrast among their intracrustal images and compositions. In the north, P and S velocities are consistent with a dominantly felsic composition and show that only the upper crust thickened. South of the Kun Lun Fault a thicker crust made of two layers could result from the superposition of the originally thin crust of the Bayan Har terrane on the lower part of the crust of the domain to the north, which upper crust it shoved and thickened. Different modes of crustal thickening, either by thickening of individual layers or superpositions and imbrication among them appear to work jointly to raise the topography.     

Most likely locations of large earthquakes in the Kanto and Tokai areas, Japan, based on the local recurrence times

We examined the hypothesis that minima in local recurrence time, T_L, or equivalently maxima in local probability, P_L, may map asperities in the Kanto and Tokai areas of Japan, where the earthquake catalog of the National Research Institute for Earth Science and Disaster Prevention (NIED) is complete at the M=1.5 (M1.5) level. We mapped T_L (P_L) based on the a- and b-values of the nearest earthquakes within 20 km of every node of a grid spaced 0.01° for M7 target events. Only earthquakes within the top 33 km were used. The b-values increase strongly with depth, in several areas. Therefore, some of the T_L (P_L) anomalies are not revealed if data from the entire crustal seismogenic zone are mixed. Thus, we mapped T_L (P_L) separately for the top 15 km and the rest of the depth range, as well as for the entire seismogenic crust. The resulting T_L- and P_L-maps show that approximately 12% of the total area shows anomalously short recurrence times. Out of six shallow target events with M≥6.5 and which occurred since 1890, five are located within the anomalous areas with T_L <450 years. We interpret this to mean that areas with anomalously short T_L map asperities, which are more likely than other areas to generate future target events. The probability that this result is due to chance is vanishingly small. The great Kanto rupture of 1923 appears to have initiated in the most significant asperity we mapped in the study area. One anomaly is located in the northeastern part of the area of the proposed future rupture of the Tokai earthquake, and another one at its southwestern corner. The absolute values of T_L calculated are uncertain because they depend on the size of the volume used for the calculation.     

2001年1月2 6日印度古杰拉特(Gujarat)Ｍem>Ssub>7.8地震时空破裂过程

从全球数字地震台网的长周期记录中，选择了震中距小于90的27个台站的54个P波震相和44个S波震相资料．首先，用波形反演方法确定了2001年1月26日印度古杰拉特(Gujarat)ＭS7.8地震的地震矩张量、震源机制、震源时间函数和时空破裂过程等震源参数．通过矩张量反演，并根据Kutch Mainland断层的走向、地震烈度的空间分布、余震震源的空间分布和震害的空间分布，确认2001年1月26日印度古杰拉特ＭS7.8地震的发震断层的走向为92、倾角为58、滑动角为62，即一走向近东-西向、断层面向南倾斜、以逆冲为主的左旋-逆断层．这次地震所释放的地震矩为3.51020 Nm，矩震级ＭW＝7.6．然后，借助合成地震图，采用频率域求谱商的方法，得到了依赖于台站方位的27个P波震源时间函数、22个S波震源时间函数以及平均的P波震源时间函数和S波震源时间函数．对震源时间函数的分析表明，这次地震是一次连续的破裂事件，开始比较急遽，但结束比较迟缓，总持续时间约19 s．最后，以所提取的P波和S波震源时间函数为资料，采用时间域的反演技术得到了断层面上滑动的时空分布．滑动量在断层面上的静态分布表明，断层面上的最大滑动量约为7 m．断层面上的最大应力降约为30 MPa，平均应力降约为7 MPa．滑动量大于0.5 m的区域在走向方向长85 km，在断层面倾斜方向宽约60 km(相应地，在深度方向约51 km)．破裂向东扩展约50 km，向西扩展约35 km．滑动量大于0.5 m的区域的主要部分呈椭圆形，其长轴取向与断层滑动方向一致．表明此区域破裂扩展的方向即是断层错动的方向．这种现象对于走滑断层情形是多见的，但对逆冲断层情形却少见．断层面上初始破裂点以东、以上部分面积大于初始破裂点以西、以下部分的面积，这是破裂非对称性的表现，表明破裂具有自西向东、自下向上单侧破裂的特征．从滑动率随时空变化的快照可以看出，滑动率在第4 s达到最大值，此时滑动率约为0.2 m/s，滑动基本上发生在破裂起始点及其周围．从第6 s开始，起始点的破裂基本结束，破裂开始向外围扩展．破裂向西的扩展速度明显小于向东的扩展速度．在第15 s，这种环形的扩展基本结束．自16 s以后，主要是一些零星的破裂点分布在破裂区的外围．从滑动量随时空变化的快照看，破裂自起始点开始后，逐渐向四周扩展．主要的破裂(滑动量大于5 m的区域)在6～10 s，具有明显的自西向东、向上的单侧破裂特征．在第11～13 s，破裂的西端向西、向下有所扩展．整个破裂过程持续约19 s．在整个破裂过程中的平均破裂速度约为3.3 km/s．      

标定大陆科学钻探孔区地震反射体

在中国大陆科学钻探孔区 (江苏省东海县南部 )进行了系统的地球物理调查 ,包括二维地震测网和专门的地震剖面 ,大地电磁法和位场方法等 .地震调查表明 ,在超高压变质岩出露区上地壳充满了反射体 ,包括倾斜反射体与上拱的反射弧 .本文介绍大陆科学钻探先行研究中地震调查的成果 .根据大陆科学钻探预先导孔的岩芯和测井资料、井旁VSP和数值模拟结果证实 ,高波速的榴辉岩体、破碎断裂带和大型韧性剪切带都可引起倾斜的地震反射 ,而上拱的弧形反射体则是由近似直立的榴辉岩体和其中的破碎带的综合反映 .由于地壳深部广泛分布着经受变质的岩石 ,上述研究结果对标定地壳中的反射地震信号具有一定意义     

Nonlinear earthquake response analysis and energy calculation for seismic slit shear wall structures

Based on the concept of structural passive control, a new type of slit shear wall, with improved seismic performance when compared to an ordinary solid shear wall, was proposed by the authors in 1996. The idea has been verified by a series of pseudo-static and dynamic tests. In this paper a macro numerical model is developed for the wall element and the energy dissipation device. Then, nonlinear time history analysis is carried out for a 10-story slit shear wall model tested on a shaking table. Furthermore, the seismic input energy and the individual energy dissipated by the components are calculated by a method based on Newmark-β assumptions for this shear wall model, and the advantages of this shear wall are further demonstrated by the calculation results from the viewpoint of energy. Finally, according to the seismic damage criterion on the basis of plastic accumulative energy and maximum response, the optimal analysis is carried out to select design parameters for the energy dissipation device.