特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

中国东部上第三系硬粘土的工程地质特性及边坡开挖工程环境效应研究

The engineering geological properties of Neogene hard clays and related engineering problems are frontiers in the fields of Engineering Geology, Soil Mechanics and Rock Mechanics. Recently, it has been recognized that Neogene hard clay is the intermediate type of material between the soil and the rock. Many aspects of them, such as sampling, testing, calculating and engineering process, are special, which could not be researched by means of theories and methods of traditional Soil Mechanics of Rock Mechanics. In order to get real knowledge and instruct the engineering practice, intersect studying of multiple disciplines, including Engineering Geology, Soil Mechanics and Rock Mechanics, etc., is necessary. Neogene hard clay is one of the important study objects of regional problem rocks & soils in our country, which extensively distributed in China, especially in Eastern China. Taking the related areas along the middle line of the Project of Transferring Water from the South to the North (e.g. Nanyang basin, Fangcheng-Baofeng area and Handan-Yongnian area), South-west of Shandong, Xu-Huai area and Beijing area, etc. as main study areas, the paper divided Neogene hard clays into reduction environment dominated origin and oxidation environment dominated origin, which distributed on areas western and eastern to Mount Taihangshan respectively. Intermediate types are also existed in some areas, which mainly distribute near the edges of depositional basins; they are usually of transitions between diluvial and lacustrine deposits. As to Neogene hard clays from Eastern China, the clay particle content is high, and montmorillonite or illite/montmorillonite turbostratic mineral is the dominating clay mineral. The content of effective montmorillonite is very high in each area, which is the basis for the undesirable engineering properties of Neogene hard clays. For hard clays from the same area, the content of effective montmorillonite in gray-greenish hard clay is much higher than that in purple-brownish or brown-yellowish hard clay, which is the reason why the gray-greenish hard clay usually has outstanding expansive property. On the other hand, purple-brownish or brown-yellowish hard clay has relatively less montmorillonite, so its property is better. All of these prove that the composition (clay mineral) of Neogene hard clay is the control factor for the engineering properties. Neogene hard clays have obvious properties such as fissured, overconsolidated and expansive, which are the main reasons that many engineering problems and geological harzards usually occur in Neogene hard clays. The paper systematically elaborates the engineering properties of Neogene hard clays from Eastern China, analyses the relationships between engineering properties and basic indexes. The author introduces the ANN method into the prediction of engineering property indexes of hard clays, which provides a new way for quantitatively assessment and prediction of engineering property indexes. During investigation in the field, the author found that there exists obvious seam-sheared zone between different hard clays in Miocene Xiacaowan formation in Xu-Huai area. Similar phenomenon also exists near the borderline between Neogene hard clays and underlying coal measures in the Southwest of Shandong province, which could be observed in the cores. The discovery of seam-sheard zone has important theoretical and practical significance for engineering stability analysis and revealing the origin of fissures in Neogene hard clays. The macrostructure, medium structure and microstructure together control the engineering properties of hard clays. The author analyses and summarizes the structural effects on hard clays in detail. The complex of the strength property of hard clays is mostly related to the characteristics of fissures, which is one of the main factors that affect the choice of shear strength parameters. So structure-control theory must be inseparably combined with composition-control theory during the engineering geological and rock/soil mechanics research of hard clays. The engineering properties, such as fissured, overconsolidated and expansive, control the instability of engineering behaviors of Neogene hard clays under the condition of excavation, i.e. very sensitive to the change of existence environment. Based on test data analysis, the author elaborates the effects of engineering environment change on the engineering properties. Taking Nanyang basin as example, the author utilizes FEM to study the effects of various factors on stability of cutting canal slopes, than sets forth the characteristics, development laws and formation mechanism of the deformation and failure of hard clay canal slopes, summarizes the protection and reinforcement principles, as well as the protection and remedy steps. On the basis of comparison of engineering properties of domestic and foreign Neogene muddy deposits, in the view of whole globe and associated with the geological characteristics of China, the paper demonstrates that the intermediate type of the material between the soil and the rock, named "hard clay/soft rock", which can not be separated abruptly, really exists in China. The author has given a preliminary classification based on its geological origin and distribution law, which is very significant for promoting the mixture of Engineering Geology, Soil Mechanics and Rock Mechanics. In the course of large scales engineering construction in China, many engineering experiences and testing data are gained, summarizing these testing results and automatically managing them with computer technology are very necessary. The author develops a software named "Hard Clay-Soft Rock Engineering Geological Information Management and Analysis System (HRGIMS)", realizes the automatic and visual management of geo-engineering information, on the basis of information management, the functions of test data analysis and engineering property prediction are strengthened. This system has well merits for practice and popularization.

中-新元古代海洋的硫和碳同位素的演化

Data on seawater carbon isotope in the Mesoproterozoic and Neoproterozoic is abundant. However, the sulfur isotopic age curve of seawater sulfates determined through the analysis of sulfur isotopic composition of marine evaporite is uncertain in the Mesoproterozoic and Neoproterozoic since evaporites are generally rare in Precambrian. The Mesoproterozoic and Neoproterozoic Carbonate Formations preserve not only the carbon isotopic records, but also the sulfur isotopic records of coeval seawater in the Huabei Platform and the Yangtze Platform, China. Sulfur isotopic composition can be determined by the extraction of trace sulfate from carbonate samples. Successive measurements of sulfur and carbon isotopic compositions of carbonate samples from the Mesoproterozoic and Neoproterozoic strata in the Huabei Platform and the Yangtze Platform was accomplished through the extracting of trace sulfate from carbonates. Sulfur and carbon isotopic compositions of coeval seawater were obtained from analytical results of sulfur and carbon isotopes of the same sample without diagenetic alteration. The high-resolution age curve of sulfur isotope given in this paper may reflect the trend of variations in sulfur isotope composition of seawater sulfates during the Mesoproterozoic and Neoproterozoic. It can be correlated with the characteristics of variation in age curve of carbon isotope of coeval seawater carbonates. The δ34S values of seawater varied from +10.3-37.0‰ during the Mesoproterozoic, which took on oscillated variation on the whole. The δ34S values took on high values in the Mesoproterozoic Chuanlinggou stage, Tuanshanzi stage Tieling stage and in Neoproterozoic Jing'eryu stage. The average of those was about +30‰. The sulfates have low δ34S values in the Mesoproterozoic Yangzhuang stage and Hongshuizhuang stage, The average of those was all lower than +20‰. There occured large-amplitude changs in δ34S values of seawater during the Mesoproterozoic. Large-amplitude oscillate of 534S values occured in the intervals of 1600~1400Ma and 1300~1200Ma. The δ13C values of seawater are mostly negative in Changcheng stage of late Paleoproterozoic, -0 ± 1‰ range in Jixian stage of Mesoproterozoic , and the positive 2±2‰ commonly in early Neoproterozoic Jing'eryu stage. From 1000 Ma to 900 Ma, about 108 years interval of oceanic 513C record is shortage. At the end of Paleoproterozoic (1700 - 1600 Ma), the oceanic 813C values change from -3‰ to 0‰, but strongly oscillate near 1600 Ma. Two larger variations of seawater 513C values occur in the Mesoproterozoic: one is a cycle of about 4%o happens at ca. 1400 Ma; another is rise from >2‰ to>5‰ at ca. 1250 Ma and then become stable at the near 1000 Ma. There appears a large positive excursion over +20‰ in 534S value of ancient seawater sulfates in the early Doushantuo stage. Simultaneously, 8 C values of ancient seawater occur a positive excursion reaching 10‰. These allow δ4S values and 513C values to reach high values of+51.7‰ and +6.9‰, respectively. The range of variation in 834S values of seawater is relatively narrow and 513C values are quite high in the middle Doushantuo stage. Then, δ34S values of seawater become oscillating, the same happens in δ13C values. Negative excursions in 834S values and 813C values occur simultaneously at the end of the Doushantuo stage, and the minimum of δ34S values and δ13C values dropped to -11.3‰ and -5.7‰, respectively. The ancient seawater in the Dengying stage has high δS values and δ13C values. Most of the δ34S values of the trace sulfate samples varied between +23.6‰ and +37.9‰ except two boundaries of the Dengying Formation, and the S13C values of the carbonate samples of the Dengying Formation varied between +0.5‰ and +5.0‰. There appeared large negative excursion in 834S values and δ13C values of ancient seawater at the bounder of Precambrian-Cambrian. The isotopic characteristics of sulfur and carbon implicated that the organic productivity and isotopic fractionation caused by biology were low and the palaeoceanic environment was quite unstable during the Mesoproterozoic. The increase and subsequent oscillation of seawater δ13C value occurred from 1700 to 1600 Ma and near 1300 Ma may be responsible to the two global tectonic events happened at coeval time. The characteristics of variation in sulfur and carbon isotopes of ancient seawater imply strong changes in oceanic environment, which became beneficial to inhabitation and propagation of organism. The organic production and the burial rate of organic carbon once reached a quite high level during the Doushantuo stage. However, the state of environment became unstable that means the global climate and the environment possibly were fluctuating and reiterating after the global glaciation. The negative excursions of S34S values and δ13C values occurring at the end of the Doushantuo stage represent a global event, which might be relative to the oxidation of deep seawater. The isotopic characteristics of sulfur and carbon implicated that there were a high organic productivity and a high burial rate of organic carbon in the Dengying stage. It is obvious that the palaeoceanic environment in Dengying stage was stable corresponding and beneficial for biology to inhabit and propagate except for the two boundaries. The tendency of sulfur and carbon isotopic variations maybe resulted from the gradual oxygenation of ocean environment during the Dengying stage. It has been reported that the secular variations of the sulfur isotopic compositions in seawater was negative correlated with that of carbon isotopic compositions. However, our results show that it is not the case. They were negatively correlated in some intervals and positively in some other intervals of the Mesoproterozoic and Neoproterozoic. The difference in correlation may be associated with the changes in conditions of redox in oceanic environment, e.g. sharp change of the oxidation-reduction interface. The strong changes in global environment may induce the abnormality to occur in the biogeo chemical S and C cycles in the ocean and accordingly sharp Variations in isotopic composition of seawater sulfur and carbon during the Mesoproterozoic and Neoproterozoic. Simultaneously, the global tectonism caused large changes of 87Sr/86Sr ratios. The leading factor that causes the variation in isotopic composition is different in the different intervals of the Mesoproterozoic and Neoproterozoic. Thus, there may exist different models of the biogeochemical S and C cycles in the ocean during the Mesoproterozoic and Neoproterozoic.

内蒙古岱海DH99a岩芯沉积物碳含量变化与湖区全新世气候变迁

Daihai Lake, a graben-type closed lake, lies ca. 10 km east of Liangcheng County, Inner Mongolia, north-central China. For its location at the transition of semi-humid and semi-arid areas, and in the north edge of the East Asian monsoon, the lake is sensitive to changes in climate and environment. Based on analyses of total inorganic carbon (TIC), total organic carbon (TOC) and the ratio of total organic carbon to total nitrogen (C/N ratio) of DH99a core sediments recovered in the central part of Daihai Lake, the data suggest Holocene climatic history of the lake region is reconstructed. In this paper, the TIC and TOC contents of the lake sediments are closely related to climate changes. 1) Changes in TIC content of the lake sediments is closely related to climate, which directly reflect changes of temperature in the lake region, i.e., higher TIC content is linked with warmer temperature and stronger evaporation; 2) Changes in TOC and C/N ratio reflect the regional precipitation, i.e., relatively higher TOC content and C/N ratio indicate higher rainfall which results in stronger river flow and more organic matter entering into the lake. Data of the TIC content, TOC content and C/N ratio of DH99a core sediments suggest that climatic history of the Daihai Lake region is characterized by 4 stages. During the interval of ca. 11200-7500 a BP, higher TIC content, relatively lower TOC content and C/N ratio value indicate a warmer and slight dry condition over the lake area. From ca. 7500 to 4500 a BP, high values of TIC content, with an increase in TOC content and C/N ratio suggest the climate was warm and humid. Changes of TIC content, TOC content show that both temperature and precipitation displayed obvious fluctuations during the period, i.e., slightly cool and humid ca. 7500-6700 a BP, warm and moist ca. 6700-5300a BP, mild and comparatively humid ca. 5300-4500a BP. Between ca. 4500 and 2900 a BP, TIC content and TOC content decreased gradually while fluctuating, C/N ratio displayed a decreasing trend. These data imply that the climate generally became cooler and drier than the preceding period. The lowest values of TIC content, TOC content and C/N ratio during the interval of ca. 2900-0 a BP, demonstrate that the climate was severe, and became cool and dry. However the relative higher values of TIC content, TOC content and C/N ratio between ca. 1700 and 1300 a BP may denote an increase both in temperature and in precipitation. Data of TIC content, TOC content and C/N ratio in Daihai DH99a core sediments indicate that the warm period was asynchronous with the humid time, the warm interval began in ca. 11200 a BP, and ended in ca. 2900 a BP. The humid period was ca. 7500-2900 a BP. During ca. 7500^500 a BP, the climate was warm and humid, which was the climatic optimum of the Holocene Epoch in the Daihai Lake region. Data of TIC content, TOC content and C/N ratio in Daihai DH99a core sediments imply that the Holocene climate was unstable, the fluctuating events happened occasionally. Such as the cold and dry climate in ca. 4400-4200a BP, the warm and humid climatic condition in the period between ca.1700 and 1300a BP.

岱海岩芯沉积的粒度组成与湖区全新世降水变化历史

Sediment cores DH99a and DH99b recovered in the central part of Daihai Lake in north-central China were analysed at 2- to 4-crn intervals for grain-size distribution. Grain-size distributions of the lake sediments are inferred to be a proxy for past changes in East Asian monsoon precipitation, such that greater silt-size percentage and higher median grain size reflect increased monsoonal precipitation rates. The grain-size record of Daihai Lake sediments spanning the last ca 11,000 yr indicates that the monsoonal precipitation in the lake region can be divided into three stages: the Early, Middle and Late Holocene. During the Early Holocene before ca 7900 cal yr BP, the median grain size (Md) and the silt-fraction content were relatively low and constant, suggesting relatively low precipitation over the lake region. The Middle Holocene between ca 7900 and 3100 cal yr BP was marked by intensified and highly variable monsoonal precipitation, as indicated by high and variable Md values and silt contents of the lake sediments. During this period, average precipitation rate gradually increased from ca 7900 to 6900 cal yr BP, displayed intense oscillations between ca 6900 and 4400 cal yr BP, and exhibited a decreasing trend while fluctuating from ca 4400 to 3100 cal yr BP Although generally high during the Middle Holocene, both the Md and the silt content assumed distinctly low values at the short intervals of ca 6500-6400, 6000-5900, 5700-5600, 4400-4200 cal yr BP, implying that monsoonal precipitation might have been significantly reduced during these intervals. During the Late Holocene since ca 3100 cal yr BP, grain-size values suggest that precipitation decreased. However, during the Late Holocene, relatively higher Md values and silt contents occurring between ca 1700 to 1000 cal yr BP may denote an intensification of hydrological cycles in the lake area. Changes in the East Asian monsoonal precipitation were not only directly linked with the changing seasonality of solar insolation resulting from progressive changes in the Earth's orbital parameters, but also may have been closely related to variations in the temperature and size of the Western Pacific Warm Pool, in the intensity of the El Nino-Southern Oscillation, and in the path and strength of the North Equatorial Current in the western Pacific.

四川省拉拉铜矿床岩石学及地球化学研究

拉拉铜矿床位于扬子地块西缘的康滇隆起中段，赋存于经历变质的古元古代河口群海相火山-沉积岩系中，是我国重要的受变质热液叠加或改造的火山成因块状硫化物矿床（VHMS）之一，也是四川省最大的铜矿生产基地。本论文以该矿床为研究对象，在系统整理前人资料的基础上，进行了深入细致的野外地质考察和系统采样，选取代表性样品，运用电子探针分析技术、元素地球化学、同位素地球化学（Rb-Sr、Sm-Nd、S）、锆石U-Pb定年等方法手段，对矿区新元古代大规模区域变质作用的变质条件、变质原岩、广泛分布于矿区多个赋矿层位的辉长辉绿岩脉、岩石及矿石性质、硫源、成矿时代等主要岩石学及矿床学问题进行了系统研究，并结合前人研究成果，探讨了矿床的成因机制。论文主要取得以下几点认识： （1）矿区变质火山岩的成岩年龄为1695±20Ma，岩浆在上升侵位的过程中捕获了古元古代早期地壳物质，后期变质改造的时代约为700-800Ma。 （2）基于石榴石黑云母片岩中石榴石黑云母矿物对的电子探针分析数据，利用石榴石-黑云母地质温度计和多硅白云母地质压力计确定了矿区新元古代大规模区域变质作用的变质条件：变质温度为530-580℃，变质压力上限为0.66-0.76GPa，对应于高绿片岩相。 （3）矿区出露的岩石类型以云母片岩类和钠长岩类岩石为主，二者共同构成了矿区的围岩及赋矿岩石，通过对代表性样品的元素地球化学、同位素地球化学研究，采用DF函数判别法、(Al/3-K)-(Al/3-Na)图解法、La/Yb-TR图解法，并结合岩矿鉴定及结构构造特征，恢复了矿区的变质原岩。其中，矿区的云母片岩类样品的变质原岩为沉积岩，以页岩为主，主要来自大陆岛弧及上地壳物质的风化；而钠长岩类样品的变质原岩为火山岩，源自富集地幔，包括“右倾型”钠长岩和“平坦型”钠长岩，前者为分离结晶作用早期并经历一定程度地壳混染的产物，后者主要形成于分离结晶作用晚期，Nb、Ta、Zr、Hf、P、Ti等达到饱和并以副矿物晶出。 （4）综合矿区变质沉积岩及变质火山岩的构造环境判别结果，确定了二者形成于大陆岛弧的弧后盆地环境。 （5）通过对矿区变质火山岩性质、形成机制及硫同位素研究，讨论了矿床中成矿元素可能的富集机制，即地幔中大量的亲铜元素及硫元素进入母岩浆并随之上升形成区域分布的火山岩，为矿床的形成提供了必要条件。 （6）侵入河口群地层并广泛出露于矿区多个赋矿层位的辉长辉绿岩脉的化学成分与板内碱性玄武岩类似且形成于大陆裂谷环境下的类似OIB源区，是高温地幔柱部分熔融的产物，幔源岩浆在上升侵位过程中受到了地壳物质尤其是上地壳物质较小程度的混染，表明新元古代Rodinia超大陆裂解对拉拉地区的成矿作用有重要影响。 （7）拉拉铜矿床主要存在两个成矿期次，即火山喷气-沉积成矿期（1700 Ma左右）和变质热液成矿期（700-800Ma）。矿区主要存在两种类型矿石，即条纹、条带状矿石和块状、浸染状矿石，前者主要赋存于变质火山岩中，可能形成于火山喷气-沉积成矿期，后者主要赋存于变质沉积岩中，可能为变质热液成矿期的产物。