基于并行CPU-GPU平台高温高压岩石内水输运机理的研究.pdfVIP

中国地球物理201l ·299· (8)地球内部结构及其动力学 基于并行CPU．GPU平台高温高压岩石内水输运机理研究 state inrock ofionic mechanism Study water(HydrogenOxidation)transport underultra and on CPUandGPU hightemperaturepressureparallel platform 朱伯靖石耀霖 Zhu ShiYaolin BoJing 中国科学院计算地球动力学重点实验室中国科学院研究生院地球科学学院 北京100049 地球深部结构及物质转化输运机理地球动力学领域核心问题，也是研究其它地球动力学问题理论基础 和根本出发点之一。但鉴于地球内部结构、物质构成、输运方式极端复杂性，虽然众多学者一直从地质学、 地球化学、地球物理学、地球动力学等角度做了众多理论、实验、数值等研究工作，也得到很多奠基性的 卓越成果理论；但这个问题至今还远没解决。本研究尝试在地球物理学和地球化学理论实验成果基础上， 究：利用CT扫描、数字重建，得到多尺度虚拟岩石模型，在并行CPU．GPU平台上对虚拟模型进行高温 高压数值模拟；通过对分子尺度、单晶尺度、多晶尺度、微观尺度、宏观尺度“虚拟实验”，研究高温高压 下岩石内部非自由态／自由态水输运机理，以期为地球深部板块运动过程中脱水作用提供理论支持与帮助。 (1)利用基于第一性原理的materialstudio分子结构分析软件，对橄榄石晶体结构建模，得到橄榄石 单晶体空间结构，应用分子动力学理论，模拟空间晶体结构变形与能量变化关系，及晶体结构出现缺陷和 夹杂时能量变化(缺陷：镁离子和铁离子被氢离子与氧离子替代，晶体结构出现歧变；夹杂：空间晶体结 构内部，有氢离子和氧离子，铁镁离子偏离正常位置，晶体结构出现歧变)，得到歧变能大小；把正常晶 Boltzmann 体结构、缺陷歧变晶体结构和夹杂歧变晶体结构作为同为基于第一性原理的lattice method输入 模型，以铁离子和镁离子为固体骨架，晶格结构上的氢离子和氧离子为缺陷“孔洞”，进一步非晶格节点 上的氢离子和氧气离子作为“流体”通过三种晶体结构时所需压力／温度大小，进而计算“氢一氧离子流 体”通过三种晶体结构时，晶体晶格结构改变的能量大小。 晶体尺度下的不同组分典型橄榄岩／石空间结构(橄榄石和夹杂)数值模型，应用LBM方法，把橄榄石 多晶体空间结构作为固体骨架，其它矿物夹杂作为灰色地带处理，把氢离子和氧离子作为“流体”，模拟 其通过规则多晶体晶格内部、及晶体无规则生长的晶体外部所需的压力／温度；得到“氢一氧离子流体” 通过晶体内部与晶体外部所需压力／温度大小，进而计算所需能量的大小。 日本高知岩芯中心)，利用数字重建技术，得到不同分辨率下的虚拟数字岩芯模型，利用LBM方法、模 计算水流体通过所需能量大小。 利用数字 (4)在厘米尺度上，利用典型橄榄岩50pro分辨率CT扫描数据(日本高知岩芯中心))， 构裂隙／缺陷所需的压力／温度大小，进而计算水流体通过所需能量大小。 (5)根据以上不同尺度模型，建立多尺度MS．LBM模型，模拟非自由态／自由态水通过晶体内部／ 晶体间／裂隙间隙力的大小，进而计算所需能量的大小量级。 初步研究结果表明：高温高压下岩石内部水输运时，通过晶间所需能量较通过晶体内部所需能量要小 输运方式主要是以离子态通过晶体间隙进行输运的；结合实验研究，数值模拟可为地球内部板块脱水作用 的研究提供理论支持和帮助。 A(NOYl5101KY00)支持。