伪谱法地震波正演模拟的多线程并行计算   总被引：16，自引：7，他引：9       下载免费PDF全文

谢桂生  刘洪  赵连功 《地球物理学进展》2005,20(1):17-23

地震波的正演模拟,尤其是3D正演模拟,往往涉及大规模的数据存储和计算,问题的规模往往超出计算机的物理内存,或者计算时间让问题的求解者难以忍受,即使采用目前存储和计算能力很强的计算机,其计算费用仍然是十分昂贵的.本文提出一种基于多线程协同使用多CPU和计算域分割的正演模拟并行计算技术,使得问题的求解过程得以加快,大大地缩短了用户等待的时间.为了检验我们的并行算法的可行性,文中以傅利叶正演模拟技术为例,给出了声波和3D各向异性弹性波模拟的例子,并对不同版本（串行、并行）运算效率进行了比较,证实了方法的有效性.  相似文献   

三维GIS环境下的地质体可视化和特征分析   总被引：4，自引：0，他引：4  

曾新平  吴健生  杨自安  郑跃鹏  朱谷昌  张普斌 《地质与勘探》2005,41(1):72-76

地质体的三维可视化建模是三维GIS研究的一个重要方向,具有重要的应用价值.三维GIS环境下的地质体可视化的两个关键方面是地质体的表面模拟和可视化实现,文章介绍了地质体三维可视化的基本框架,以带约束Deluanay剖分法为例阐述了地质体表面模拟的原理和方法,并从取景变换、消隐、光照模型、纹理映射等方面论述了可视化的具体实现方法.最后以西部某铜矿为例,建立了三维可视化地质体模型(一个地层矿体综合模型和一个矿体三维模型),阐述了三维可视化建模的基本过程和方法,并对两个模型进行了地质特征分析.  相似文献   

坡麓相斜坡软土特性及其地质灾害防治研究   总被引：1，自引：0，他引：1  

卿三惠  黄润秋 《地质灾害与环境保护》2005,16(2):113-118

结合工程实例．提出“坡麓相斜坡软土”的概念,论述西南煤系地层地区斜坡软土的形成机理、分布规律、物理力学特征及其工程特性,研究斜坡软土对工程的危害及其地质灾害防治对策,最后列举三个典型斜坡软土工程整治实例,对指导工程实践具有重要意义。  相似文献   

Preliminary estimation of the tsunami hazards associated with the Makran subduction zone at the northwestern Indian Ocean   总被引：3，自引：3，他引：0  

Mohammad Heidarzadeh  Moharram D. Pirooz  Nasser H. Zaker  Ahmet C. Yalciner 《Natural Hazards》2009,48(2):229-243

We present a preliminary estimation of tsunami hazard associated with the Makran subduction zone (MSZ) at the northwestern Indian Ocean. Makran is one of the two main tsunamigenic zones in the Indian Ocean, which has produced some tsunamis in the past. Northwestern Indian Ocean remains one of the least studied regions in the world in terms of tsunami hazard assessment. Hence, a scenario-based method is employed to provide an estimation of tsunami hazard in this region for the first time. The numerical modeling of tsunami is verified using historical observations of the 1945 Makran tsunami. Then, a number of tsunamis each resulting from a 1945-type earthquake (M _w 8.1) and spaced evenly along the MSZ are simulated. The results indicate that by moving a 1945-type earthquake along the MSZ, the southern coasts of Iran and Pakistan will experience the largest waves with heights of between 5 and 7 m, depending on the location of the source. The tsunami will reach a height of about 5 m and 2 m in northern coast of Oman and eastern coast of the United Arab Emirates, respectively.  相似文献   

End Bearing Capacity of Drilled Shafts in Sand: A Numerical Approach     

M. M. Ahmadi  M. Khabbazian 《Geotechnical and Geological Engineering》2009,27(2):195-206

In this paper, a modeling procedure is carried out to numerically analyze the end bearing capacity of drilled shafts in sand. The Mohr–Coulomb elastic plastic constitutive law with stress dependent elastic parameters is used for all numerical analyses performed in this study. The numerical results are compared with the available experimental equations. It is seen that numerical results are in good agreement with experimental equations. The variation of the end bearing capacity of drilled shafts versus embedment depth is also studied. Numerical results show that with increase in pile embedment depth, the end bearing capacity increases. However, the rate of increase becomes smaller as the pile embedment depth increases. Also, numerical analyses show that, for equal settlement, the end bearing decreases with increase in the pile diameter. Finally, a sensitivity analysis is performed to obtain the separate effect of each sand parameter on the end bearing capacity of drilled shafts, and the parameters that are most influential are identified.  相似文献   

综合应用测井及地震资料建立复杂断块地质模型     

张伟  林承焰  周明晖  贾玉梅 《物探化探计算技术》2009,31(6):553-558

针对滩海油田具有完钻井少,以地震资料为主的特点,采用井震结合的方法进行地质建模。复杂断块油藏具有许多内部断层,充分利用地震解释层位和断层成果建立合理的构造模型。采用“多步建模”的思路,选用序贯高斯同位协同模拟算法建立储层参数模型。以测井资料为硬数据,地震反演数据为软数据,先后建立自然伽玛模型、泥质含量模型和砂泥岩相模型,然后以砂泥岩相模型为约束,建立储层参数模型。经研究表明,综合使用地震解释成果和三维地震数据体,有助于建立准确、合理的构造模型;而采用“多步建模”的思路,利用序贯高斯同位协同模拟算法,并综合测井资料、地震资料建立储层参数模型,有助于从地质的角度对储层参数模拟进行约束,提高储层随机建模精度,保证储层参数模型的准确和可靠性。  相似文献   

近10年来“国际古气候模拟比较计划(PMIP)”回顾和未来古气候模拟研究热点   总被引：1，自引：1，他引：0       下载免费PDF全文

靳立亚　Bette L．Otto­  Bliesner 《第四纪研究》2009,29(6):1015-1024

国际古气候模拟比较计划(PMIP)是古气候领域里的一项长期开展的大型国际研究计划。该计划的主旨是利用气候数值模式,通过开展一系列针对古气候变化的模拟试验理解全球气候变化的机理,确定气候系统内部的各种影响气候变化的关键反馈因子; 通过比较不同气候模式对不同于现代气候(古气候)的模拟能力评价模式的优劣,进而促进气候模式的发展; 同时也为全球变化背景下利用气候模式预估未来气候变化提供借鉴。PMIP自1991年酝酿开始至今经历了PMIP1(1991～2001年)、PMIP2(2002～2007年)两个阶段,并正在开展第三阶段(PMIP3)(2008～)的研究。本文简要回顾了近10年来PMIP研究进展及概述了未来PMIP将要开展研究的重点领域和方向,主要包括4个主题： 1)基准(benchmark)时段(6kaB.P.和21kaB.P.)模拟; 2)间冰期和温暖期气候模拟; 3)快速气候变化模拟; 4)不确定性研究。  相似文献   

上新世中期中国气候的数值模拟分析   总被引：2，自引：1，他引：1       下载免费PDF全文

姜大膀 《第四纪研究》2009,29(6):1033-1043

使用中国科学院大气物理研究所的全球大气环流模式,曾就上新世中期(距今约3百万年前)气候进行过数值模拟试验,在已有针对全球气候进行分析的工作基础上,集中研究了该地质时期我国气候的变化特征。数值试验结果表明,除了青藏高原地区由于地形海拔高度变化而引起的局地冷却现象以外,上新世中期我国约100°E以东地区地表气温上升4～8℃,升温中心位于江淮流域下游地区; 该经度线以西地区地表气温升幅相对要弱一些,在1～4℃之间。上新世中期年均降水在我国东部地区显著减少,平均减幅在0.5mm/天以上,特别是在长江中游地区; 新疆北部、青海和西藏大部年均降水略有增加,而新疆中部和南部年均降水略有减少。在对流层低层,上新世中期东亚冬季风系统减弱,夏季风系统总体上略有减弱; 在对流层中层,冬季东亚大槽显著减弱,夏季中高纬度地区500hPa位势高度升高。在以上结果中,上新世中期我国气候相对于现在偏暖、东亚冬季风强度减弱得到了代用资料的支持。  相似文献   

P-T Paths Derived from Garnet Growth Zoning in Danba Domal Metamorphic Terrain, Sichuan Province, West China     

Cheng Suhua  Lai Xingyun  You Zhendong 《中国地质大学学报(英文版)》2009,20(2)

Danba (丹巴) domal metamorphic terrain belongs to Songpan (松潘)-Ganze (甘孜) orogenic belt, where typical Barrovian and Buchan metamorphic zones are preserved. The former included chlorite, biotite, garnet, staurolite, kyanite and sillimanite zones, while the latter only developed silimanite+muscovite and sillimanite+K-feldspar zones. Integrated study has been carried on metamorphic reactions of garnet production and consumption, P-T paths and P-T-X-M phase relation and thermal tectonic model for Danba metamorphic zones. Petrological textures in thin sections show that garnet production and consumption in kyanite-sillimanite zone is mainly attributed to ChI+Ms+PI+Q=Grt+Bt+H2O and kyanite=sillimanite respectively. Based on mineral compositions, the geothermobarometry gives an average P, T condition of (4.9±0.3)×108 Pa, 543±30℃ for the first growth stage of the garnet and (5.8±0.3)±108 Pa, 534±29 ℃ for the second stage of garnet growth respectively. Anti-counter clockwise P-T paths were drawn using Gibbs method by NCMnKFMASH system for sample G98686 in the kyanite zone. The P-T-X-M modeling for the first mineral assemblages shows that the prediction is similar to the measured values in gossular, almandine and spessartine but mole fraction of pyrope and Fe/(Fe+Mg) deviated far from the contours; while that for the second mineral assemblages exhibits that the prediction is consistent with the measured value of pyrope, grossular content and Fe/(Fe+Mg) of garnet. A thermal tectonic model that there are at least three structure levels across the thrnst-decollement zones is presented according to the P-T paths, metamorphic grades and deformation styles for the staurolite-kyanite zone of the Barroviau type metamorphism, which will provide some constraints for the evolution of the nappe complex.  相似文献   

High-Frequency Rayleigh-Wave Method     

Jianghai Xia  Richard D Millerg  Xu Yixian  Luo Yinhe  Chen Chao  Liu Jiangping  Julian Ivanov  Chong Zeng 《中国地质大学学报(英文版)》2009,20(3)

High-frequency (≥2 Hz) Rayleigh-wave data acquired with a multichannei recording sys-tem have been utilized to determine shear (S)-wave velocities in near-surface geophysics since the early 1980s. This overview article discusses the main research results of high-frequency surface-wave tech-niques achieved by research groups at the Kansas Geological Survey and China University of Geosciences in the last 15 years. The multichannel analysis of surface wave (MASW) method is a nou-iuvasive acoustic approach to estimate near-surface S-wave velocity. The differences between MASW results and direct borehole measurements are approximately 15% or less and random. Studies show that simultaneous inversion with higher modes and the fundamental mode can increase model resolution and an investigation depth. The other important seismic property, quality factor (Q), can also be estimated with the MASW method by inverting attenuation coefficients of Rayleigh waves. An inverted model (S-wave velocity or Q) obtained using a damped least-squares method can be assessed by an optimal damping vector in a vicinity of the inverted model determined by an objective function, which is the trace of a weighted sum of model-resolution and model-covariance matrices. Current developments include modeling high-frequency Rayleigh-waves in near-surface media, which builds a foundation for shallow seismic or Rayleigh-wave inversion in the time-offset domain; imaging dispersive energy with high resolution in the frequency-velocity domain and possibly with data in an arbitrary acquisition geometry, which opens a door for 3D surface-wave techniques; and successfully separating surface-wave modes, which provides a valuable tool to perform S-wave velocity profiling with high-horizontal resolution.  相似文献