分段光滑曲线边界波动方程数值模拟研究

矩形网格有限差分法在地震波传播数值模拟方面具有计算速度快的显著优势,但该方法在处理复杂边界问题上存在着效率低的严重缺陷.本文针对分段光滑曲线边界定义了尖点处的一种正则导数,给出了矩形网格情形分段光滑曲线网格边界点法向导数的一种插值计算方法.采用矩形网格有限差分法对复杂边界地球介质模型进行地震波场数值模拟,并采用波场系列快照技术揭示地震波在起伏地表和复杂介质中的传播规律.模拟结果表明：法向导数插值计算方法为矩形网格有限差分法处理复杂边界提供了有效途径,采用波场系列快照技术可以清晰地展现地震波在反射界面的反射和透射规律、在尖点的绕射规律以及在自由表面的直达波和多次反射规律.     

土层分层参数不确定时对二维工程场址地震动影响分析

工程场地地震安全性评价中计算二维复杂场地地震反应分析时,如采用一维等效线性化分析模型会带来较大的误差,而直接采用二维的非线性模型在技术上还存在一定的困难和不合理性。目前工程中多采用对一维分析结果进行二维修正的思想给出设计地震动及反应谱。然而在建立二维分析模型时,由于勘测条件的限制使模型建立出现很多不确定性。基于以往提出的二维复杂工程场址设计地震动的修正分析思想,建立了几种可能且工程认可的二维复杂场地模型,主要研究不同分层特征模型及土层剪切波速这一物理参数不确定时对设计地震动的影响,进一步考虑不同场地类别下,不同二维分层模型及土体物理参数对地震动的影响。根据分析结果提出了不同类别场地下,方便且合理建立二维复杂场地地震动分析模型的方法,为实际工程中模型的建立及参数的选取提供一些参考。     

频率域全波形反演中关于复频率的研究

Laplace-Fourier域全波形反演可以利用简单的初始模型,从缺失低频信息的地震数据中得到长波长速度模型.Laplace-Fourier域全波形反演等价于本文的复频率全波形反演,但二者的实现方式不同,因此研究复频率全波形反演,可以为二者的对比研究并发展更有效的方法奠定重要基础.本文首先比较用线性增加模型作为初始模型时几个包含不同高低频成分的频率组的反演效果,再比较结合复频率之后各个频率组的反演效果,从简单模型和复杂模型的测试中都可以看出这种复频率+频率反演的方式对反演效果有明显改善.     

山东沂水杂岩岩石化学及锆石Hf同位素研究

山东沂水杂岩主要由高角闪岩相至麻粒岩相变质的变基性岩和(紫苏)花岗岩岩体组成.本文主要研究了三个含紫苏辉石的黑云斜长角闪岩(YS06-19、YS06-41和YS06-29),三个含石榴子石的角闪二辉斜长麻粒岩(YS06-40、YS06-45、YS06-49),一个含尖晶石和石榴子石角闪二辉麻粒岩(YS06-31)和两个英灵山花岗片麻岩样品(YS06-30和YS06-48)的岩石化学和锆石Hf同位素特征.结果认为,(1)含紫苏辉石的黑云斜长角闪岩和含石榴子石的角闪二辉斜长麻粒岩对Sr、K、Rb、Ba、Ce、Th等大离子亲石元素和轻稀土元素的富集程度不同,指示了其原岩经历了不同程度的部分熔融;(2)认为英灵山花岗片麻岩是由来自于亏损地幔的基性岩部分熔融产生,这一观点同沈其韩等(2000)认识一致;(3)所研究的变基性岩的锆石Hf亏损地幔模式年龄均小于英灵山花岗片麻岩Hf亏损地幔模式年龄,指示了该变基性岩可能不是英灵山花岗片麻岩的母岩,沂水地区应该存在更古老的变基性岩石,也可能反映了这两类岩石对Hf同位素体系的保存能力不同;(4)该地区地壳生长在30亿年左右启动,大规模的地壳生长出现在2530～2740Ma.     

内蒙古锡林浩特岩群岩相学、地球化学特征及构造环境分析

锡林浩特岩群是内蒙古中东部锡林郭勒微地块内原锡林郭勒杂岩的重要组成部分,由一套片麻岩夹层状斜长角闪岩、磁铁石英岩和变粒岩等变质岩组成,形成于中元古代［1］。岩群中片麻岩样品显示轻稀土元素富集的中等分异特征((La/Yb)N=4.71~7.97),斜长角闪岩具有平坦型((La/Yb)N=0.74~0.95)和轻稀土元素微弱富集型((La/Yb)N=1.06~3.63)两种稀土元素配分模式。在微量元素N MORB标准化图解中,斜长角闪岩显示出洋中脊玄武岩和岛弧玄武岩的过渡类型特征。同时结合岩相学和岩石地球化学图解分析,推断片麻岩原岩为富含泥质的亚杂砂岩,而斜长角闪岩原岩为一套具有拉斑质到钙碱性过渡型特征的玄武岩;变质沉积岩形成于活动大陆边缘大陆岛弧环境;变质基性火山岩形成于岛弧弧后盆地。二者指示锡林浩特岩群形成于活动大陆边缘体系。     

Silicate‐oxide mineral chemistry of mafic–ultramafic rocks as an indicator of the roots of an island arc: The Chilas Complex,Kohistan (Pakistan)

The Chilas Complex is a major lower crustal component of the Cretaceous Kohistan island arc and one of the largest exposed slices of arc magma chamber in the world. Covering more than 8000 km², it reaches a current tectonic width of around 40 km. It was emplaced at 85 Ma during rifting of the arc soon after the collision of the arc with the Karakoram plate. Over 85% of the Complex comprises homogeneous, olivine‐free gabbronorite and subordinate orthopyroxene–quartz diorite association (MGNA), which contains bodies of up to 30 km² of ultramafic–mafic–anorthositic association (UMAA) rocks. Primary cumulate textures, igneous layering, and sedimentary structures are well preserved in layered parts of the UMAA in spite of pervasive granulite facies metamorphism. Mineral analyses show that the UMAA is characterized by more magnesian and more aluminous pyroxene and more calcic plagioclase than those in the MGNA. High modal abundances of orthopyroxene, magnetite and ilmenite (in MGNA), general Mg–Fe–Al spatial variations, and an MFA plot of whole‐rock analyses suggest a calc‐alkaline origin for the Complex. Projection of the pyroxene compositions on the Wo–En–Fs face is akin to those of pyroxenes from island arcs gabbros. The presence of highly calcic plagioclase and hornblende in UMAA is indicative of hydrous parental arc magma. The complex may be a product of two‐stage partial melting of a rising mantle diaper. The MGNA rocks represent the earlier phase melting, whereas the UMAA magma resulted from the melting of the same source depleted by the extraction of the earlier melt phase. Some of the massive peridotites in the UMAA may either be cumulates or represent metasomatized and remobilized upper mantle. The Chilas Complex shows similarities with many other (supra)subduction‐related mafic–ultramafic complexes worldwide.     

Ophiolite complexes as collision zone markers: making the case for the Antique Ophiolite Complex, Panay Island, central Philippines

The collision between the North Palawan Block (NPB) and Philippine Mobile Belt (PMB) has been the subject of studies considering its significance in help-ing define the tectonic evolution of the Philippine is-land arc system. The geology of the western Panay island reveals the presence of a continent-related block (Buruanga Peninsula) juxtaposed to an oceanic frag-ment (Antique Ophiolite Complex). Our recent work in the Buruanga Peninsula helped us define the terrane boundary between the Peninsula and the Antique Ophiolite Complex. However, considering available published data, the Antique Ophiolite Complex has never been considered to be a part of the NPB and to mark the collision zone between Palawan and the PMB.     

Guandishan Granitoids of the Paleoproterozoic Lüliang Metamorphic Complex in the Trans-North China Orogen: SHRIMP Zircon Ages, Petrogenesis and Tectonic Implications

Abstract: The Paleoproterozoic Lüliang Metamorphic Complex (PLMC) is situated in the middle segment of the western margin of the Trans-North China Orogen (TNCO), North China Craton (NCC). As the most important lithological assemblages in the southern part of the PLMC, Guandishan granitoids consist of early gneissic tonalities, granodiorites and gneissic monzogranites, and younger gneissic to massive monzogranites. Petrochemical features reveal that the early gneissic tonalities and granodiorites belong to the medium-K calc-alkaline series; the early gneissic monzogranites are transitional from high-K calc-alkaline to the shoshonite series; the younger gneissic to massive monzogranites belong to the high-k calc-alkaline series, and all rocks are characterized by right-declined REE patterns and negative Nb, Ta, Sr, P, and Ti anomalies in the primitive mantle normalized spidergrams. SHRIMP zircon U–Pb isotopic dating reveals that the early gneissic tonalities and granodiorites formed at ~2.17 Ga, the early gneissic monzogranites at ~2.06 Ga, and the younger gneissic to massive monzogranites at ~1.84 Ga. Sm–Nd isotopic data show that the early gneissic tonalities and granodiorites have εNd(t) values of +0.48 to ?3.19 with Nd-depleted mantle model ages (TDM) of 2.76–2.47 Ga, and early gneissic monzogranites have εNd(t) values of ?0.53 to ?2.51 with TDM of 2.61–2.43 Ga, and the younger gneissic monzogranites have εNd(t) values of ?6.41 to ?2.78 with a TDM of 2.69–2.52 Ga.These geochemical and isotopic data indicate that the early gneissic tonalities, granodiorites, and monzogranites were derived from the partial melting of metamorphosed basaltic and pelitic rocks, respectively, in a continental arc setting. The younger gneissic to massive monzogranites were derived by partial melting of metamorphosed greywackes within the continental crust. Combined with previously regional data, we suggest that the Paleoproterozoic granitoid magmatism in the Guandishan granitoids of the PLMC may provide the best geological signature for the complete spectrum of Paleoproterozoic geodynamic processes in the Trans-North China Orogen from oceanic subduction, through collisional orogenesis, to post-orogenic extension and uplift.     

A comparative analysis of DEM‐based models to estimate the solar radiation in mountainous terrain

Daily solar radiation estimates of four up‐to‐date solar radiation models (Solar Analyst, r.sun, SRAD and Solei‐32), based on a digital elevation model (DEM), have been evaluated and compared in a Mediterranean environment characterized by a complex topography. The models' estimates were evaluated against 40 days of radiometric data collected in 14 stations. Analyzed sky conditions ranged from completely overcast conditions to clear skies. Additionally, the role of the spatial resolution of the DEM has been evaluated through the use of two different resolutions: 20 and 100 m. Results showed that, under clear‐sky conditions, the daily solar radiation variability in the study area may be reasonably estimated with mean bias errors under 10% and root mean square error values of around 15%. On the other hand, results proved that the reliability of the estimates substantially decreases under overcast conditions for some of the solar radiation models. Regarding the role of the DEM spatial resolution, results suggested that the reliability of the estimates for complex topography areas under clear‐sky conditions improves using a higher spatial resolution.     

库车坳陷复杂高陡构造地震成像研究

复杂构造地震成像主要取决于叠前地震数据品质、偏移速度可靠性和偏移算子成像精度. 库车坳陷异常复杂的近地表条件导致极低信噪比的地震采集数据. 该区逆冲推覆高陡构造刺穿盐体大面积分布, 盐层厚度变化大、顶底面形态复杂, 盐下断裂带破碎、小断块发育, 形成异常复杂的地震成像问题. 本文重点研究三个关键环节:(1)精细的叠前地震预处理研究: 根据该区地震地质复杂性和地震资料特征, 采用一些新的方法技术和技术组合从振幅与时移的大、中、小尺度变化三个层次来解决资料信噪比问题, 重建深部反射信号; (2)三级偏移速度分析研究:利用库车坳陷盐刺穿逆冲推覆构造建模理论及变速成图配套技术解决叠前时间偏移速度场时深转换问题,利用井约束低频速度地震迭代反演技术解决连井层速度场与偏移速度场的融合问题,实现从DMO速度分析、叠前时间偏移速度分析到叠前深度偏移速度分析的有机衔接,建立拓扑结构相对保持的叠前深度偏移速度模型;(3)基于退化Fourier偏移算子的半解析波动方程叠前时间和深度偏移研究, 极大地改善了地震偏移过程中高波数波的成像问题. 通过对库车坳陷大北、博孜、却勒、西秋4和西秋10等复杂高陡构造的叠前时间和深度偏移地震成像处理,取得了较好的应用效果.