西藏札达盆地控盆断裂有限元数值模拟

札达盆地是喜马拉雅构造带中的一个山间断陷盆地,其演化过程与盆地两侧的控盆构造密切相关。对控盆断裂的构造应力场进行模拟计算,将有助于进一步深化对本区构造控盆的认识。因此,在对盆地构造地质进行详细调查的基础上,结合本区的深部地质与地球物理资料,对札达盆地控盆断裂的构造应力场进行了模拟。计算结果表明,札达盆地的演化明显受盆地两侧边界断裂的控制,札达盆地是在整体南北向挤压应力的作用下,不同块体差异隆升作用的结果。其南侧的控盆断裂为北倾的正断层,北侧的控盆断裂为南倾的逆断层,二者共同形成了南降北升的翘板式断陷盆地运动过程,是喜马拉雅地块在陆内汇聚挤压构造环境中构造应力场调整的一种方式。     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

On the physical model of earthquake precursor fields and the mechanism of precursors＇timespace distribution──origin and evidences of the strong body earthquake──generating model

Ｏｎｔｈｅｐｈｙｓｉｃａｌｍｏｄｅｌｏｆｅａｒｔｈｑｕａｋｅｐｒｅｃｕｒｓｏｒｆｉｅｌｄｓａｎｄｔｈｅｍｅｃｈａｎｉｓｍｏｆｐｒｅｃｕｒｓｏｒｓ＇ｔｉｍｅｓｐａｃｅｄｉｓｔｒｉｂｕｔｉｏｎ──ｏｒｉｇｉｎａｎｄｅｖｉｄｅｎｃｅｓｏｆ　ｔｈｅｓｔｒｏｎ...     

中国三个大地震的震源参数及讨论

大地震的震源机制及地震矩、应力降等参数,不(又对认识地震的破裂过程,而且对预报强地面运动都是非常重要的。一些文章对板内地震应力降的特点、地震矩与高频源谱的关系及区域特征的讨论,使我们可直接由长周期体波或面波得出的地震矩M₀估计出不同地区中大地震的高频源谱,以供强地面运动预报参考。 由于对发生在中国板内的1976年11月6日盐源地震（M_s=6.5）和1976年11月15日宁河地震（M_s=6.3）的震源参数还没有详尽的报道;另外,对1973年7月14日     

Large eddy simulation of hydrocyclone—prediction of air-core diameter and shape

Hydrocyclones are widely used in the mining and chemical industries. An attempt has been made in this study, to develop a CFD (computational fluid dynamics) model, which is capable of predicting the flow patterns inside the hydrocyclone, including accurate prediction of flow split as well as the size of the air-core. The flow velocities and air-core diameters are predicted by DRSM (differential Reynolds stress model) and LES (large eddy simulations) models were compared to experimental results. The predicted water splits and air-core diameter with LES and RSM turbulence models along with VOF (volume of fluid) model for the air phase, through the outlets for various inlet pressures were also analyzed. The LES turbulence model led to an improved turbulence field prediction and thereby to more accurate prediction of pressure and velocity fields. This improvement was distinctive for the axial profile of pressure, indicating that air-core development is principally a transport effect rather than a pressure effect.     

Increasing critical sensitivity of the Load/Unload Response Ratio before large earthquakes with identified stress accumulation pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.     

Variability in Early Amazonian Tharsis stress state based on wrinkle ridges and strike-slip faulting

We present evidence for a decrease in the magnitude of Tharsis-circumferential compressive stress during the Late Hesperian to the Middle Amazonian based on chronologic changes in the predominant style of faulting in southern Amazonis Planitia. Using high-resolution MOLA topography, we identify a population of strike-slip faults that exhibit Middle Amazonian-aged displacements of regional chrono-stratigraphic units. These strike-slip faults are adjacent to an older population of previously documented Late Hesperian-aged thrust faults (wrinkle ridges). Along-strike orientations of these thrust and strike-slip faults reveal the Tharsis-radial stress to be the area's most compressive remote principal stress and that this stress orientation and magnitude persisted throughout the Late Hesperian to the Middle Amazonian. We show that the change in the predominant style of faulting from thrust faulting to strike-slip faulting during this time requires a decrease of the Tharsis-circumferential compressive stress to a magnitude less than lithostatic load, with negligible change in stress orientation.     

Seismic activity associated with the subducting motion of the Philippine Sea plate beneath the Kanto district, Japan

The Pacific plate and the Philippine Sea plate overlap and subduct underneath the Kanto region, central Japan, causing complex seismic activities in the upper mantle. In this research, we used a map selection tool with a graphic display to create a data set for earthquakes caused by the subducting motion of the Philippine Sea plate that are easily determined. As a result, we determined that there are at least four earthquake groups present in the upper mantle above the Pacific plate. Major seismic activity (Group 1) has been observed throughout the Kanto region and is considered to originate in the uppermost part of mantle in the subducted Philippine Sea plate, judging from the formation of the focal region and comparison with the 3D structure of seismic velocity. The focal mechanism of these earthquakes is characterized by the down-dip compression. A second earthquake layer characterized by down-dip extension (Group 2), below the earthquakes in this group, is also noted. The focal region for those earthquakes is considered to be located at the lower part of the slab mantle, and the Pacific plate located directly below is considered to influence the activity. Earthquakes located at the shallowest part (Group 3) form a few clusters distributed directly above the Group 1 focal region. Judging from the characteristics of later phases in these earthquakes and comparing against the 3D structure of seismic velocity, the focal regions for the earthquakes are considered to be located near the upper surface of the slab. Another earthquake group (Group 4) originates further below Group 2; it is difficult to consider these earthquakes within a single slab. The seismic activities representing the upper area of the Philippine Sea plate are Group 3. This paper proposes a slab geometry model that is substantially different from conventional models by strictly differentiating the groups.     

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

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

Effect of a submerged bay-mouth breakwater on tsunami behavior analyzed by 2D/3D hybrid model simulation

The hybrid numerical model had been developed to simulate a complicated 3D flow around structures generated by tsunami. In the model, the conventional 2D model is adopted for the wide region far from structures and the 3D non-hydrostatic pressure model is used in the limited region adjacent to structures. The applicability of the model is shown by comparisons of the numerical results with the experimental results of the laboratory model tests and the numerical analysis results of the conventional whole 2D simulation. In addition, the effect of a submerged structure at the opening of a breakwater is discussed from the numerical simulations by the hybrid model. The submerged structure improves the stability of the rubble mound and reduces the tsunami inflow into the bay, while it increases the water surface velocity around the opening of the breakwater. The increase of surface velocity causes the increases of impulsive forces by collision with drafts and so on.