基于MPI的LAGFD-WAM海浪数值模式并行算法研究   总被引：5，自引：0，他引：5  

王关锁  乔方利  杨永增 《海洋科学进展》2007,25(4):401-407

在LAGFD-WAM海浪串行数值模式基础上,利用MPI信息传递机制实现其并行化。通过对模拟区域合理划分,对数据采取分块加载,实现了各个节点的负载平衡;通过对算法的改进实现了粗粒度计算,大幅度减少了通信量,从而提高了程序的执行效率。对串行计算和并行计算的效率比较表明,本文建立的方法能够得到较高的加速比。对全球海浪模式,加速比和CPU数目大致呈线性关系。对高分辨率的区域海浪模式,在128 CPU条件下加速比可以达到91.9。  相似文献   

Theoretical aspects of cap-rock and fault seals for single- and two-phase hydrocarbon columns   总被引：6，自引：0，他引：6  

N.L. Watts   《Marine and Petroleum Geology》1987,4(4)

Cap-rock seals can be divided genetically into those that fail by capillary leakage (membrane seals) and those whose capillary entry pressures are so high that seal failure preferentially occurs by fracturing and/or wedging open of faults (hydraulic seals). A given membrane seal can trap a larger oil column than gas column at shallow depths, but below a critical depth (interval), gas is more easily sealed than oil. This critical depth increases with lower API gravity, lower oil GOR and overpressured conditions (for the gas phase). These observations arise from a series of modelling studies of membrane sealing and can be conveniently represented using pressure/ depth (P/D) profiles through sealed hydrocarbon columns. P/D diagrams have been applied to the more complex situation of the membrane sealing of a gas cap underlain by an oil rim; at seal capacity, such a two-phase column will be always greater than if only oil or gas occurs below the seal.These conclusions contrast with those for hydraulic seals where the seal capacity to oil always exceeds that for gas. Moreover, a trapped two-phase column, at hydraulic seal capacity will be less than the maximum-allowed oil-only column, but more than the maximum gas-only column. Unlike membrane seals, hydraulic seal capacity should be directly related to cap-rock thickness, in addition to the magnitude of the minimum effective stress in the sealing layer and the degree of overpressure development in the sequence as a whole.Fault-related seals are effectively analogous to membrane cap-rocks which have been tilted to the angle of the fault plane. Consequently, all of the above conclusions derived for membrane cap-rocks apply to both sealing faults sensu stricto (fault plane itself seals) and juxtaposition faults (hydrocarbon trapped laterally against a juxtaposed sealing unit). The maximum-allowed two-phase column trapped by a sealing fault is greater than for equivalent oil-only and gas-only columns, but less than that predicted for a horizontal membrane cap-rock under similar conditions. Where a two-phase column is present on both sides of a sealing fault (which is at two-phase seal capacity), a deeper oil/water contact (OWC) in one fault block is associated with a deeper gas/oil contact (GOC) compared with the adjacent fault block. If the fault seal is discontinuous in the gas leg, however, the deeper OWC is accompanied by a shallower GOC, whereas a break in the fault seal in the oil leg results in a common OWC in both fault blocks, even though separate GOC's exist. Schematic P/D profiles are provided for each of the above situations from which a series of fundamental equations governing single- and two-phase cap-rock and fault seal capacities can be derived. These relationships may have significant implications for exploration prospect appraisal exercises where more meaningful estimates of differential seal capacities can be made.The membrane sealing theory developed herein assumes that all reservoirs and seals are water-wet and no hydrodynamic flow exists. The conclusions on membrane seal capacity place constraints on the migration efficiency of gas along low-permeabiligy paths at depth where fracturing, wedging open of faults and/or diffusion process may be more important. Contrary to previous assertions, it is speculated that leakage of hydrocarbons through membrane seals occurs in distinct pulses such that the seal is at or near the theoretically calculated seal capacity, once this has been initially attained.Finally, the developed seal theory and P/D profile concepts are applied to a series of development geological problems including the effects of differential depletion, and degree of aquifer support, on sealing fault leakage, and the evaluation of barriers to vertical cross-flow using RFT profiles through depleted reservoirs. It is shown that imbibition processes and dynamic effects related to active cross-flow across such barriers often preclude quantitative analysis and solution of these problems for which simulation studies are usually required.  相似文献   

湄洲湾潮流数值计算     

孙长青  赵可胜  郭耀同 《海洋科学》2004,28(12):35-38

采用不规则三角形网格的分步杂交法,建立了湄洲湾海域二维变动边界潮流数值模型,通过计算得到最大潮流分布、同潮时线与等振幅线、潮流椭圆、潮致欧拉余流分布及不同时刻潮流场分布。  相似文献   

多点液压式波浪能海水淡化系统建模与仿真     

胡缘  杨绍辉  何宏舟  陈沪  郑松根  张迪 《海洋工程》2019,37(1):134-141

为缓解淡水资源短缺及化石能源过度使用问题,提出多点液压式波浪能海水淡化系统,该系统主要由采能装置、液压传递系统与反渗透膜海水淡化设备组成。系统的采能装置采用振荡浮子式,可将波浪能转换为浮子振荡从而被液压系统吸收达到采集波浪能的目的。为了提高液压式波浪能海水淡化系统的采能效率及淡水率,利用AMEsim软件对液压传递系统进行建模与仿真,分析了蓄能器、浮子个数及波高对液压传递系统输出响应的影响。结果表明:蓄能器能够使液压马达的输出响应更加稳定;当浮子的数量增加时,液压系统达到稳定的运行状态所需的时间更短,从而有利于提高系统的效率;波高在2 m左右时,本系统的产水量达到最大。  相似文献   

基于并行计算的海洋平台谱分析疲劳筛选技术研究     

李洛东  刘成名  梁园华 《海洋工程》2018,36(5):56-64

为了得到海洋平台疲劳失效风险最大的焊接区域,提出了以有限元网格单元应力为基础的谱分析筛选方法,并开发了程序。程序采用并行计算架构,单次完成数万个单元的累积损伤度/疲劳寿命计算,并通过有限元软件以云图方式显示疲劳寿命。采用该程序计算双浮筒半潜式平台的疲劳寿命。结果表明,该型平台疲劳失效风险较高的区域为立柱与上壳体下浮体连接对角线转角区域、浮筒中纵舱壁与立柱连接区域。并行计算方法的引入,大大缩减疲劳损伤度计算的分析时间,采用22核心的CPU,时间缩短为原来的1/18.5,17小时完成全平台外壳单元谱分析计算。  相似文献   

地下水与地表水水量交换识别及交换量计算——以新汴河宿州段为例          下载免费PDF全文

束龙仓  宫荣  栾佳文  龚建师  陶小虎  赵贵章  陶月赞 《水科学进展》2022,33(1):57-67

为提高地下水与地表水交换量计算结果的准确性,本文利用水力联系、水头差、水温、氡-222、氢氧稳定同位素构建综合识别方法(HHTRO),对新汴河宿州段地下水与地表水水量交换进行识别,并计算交换量。计算结果表明：研究河段单位河长地表水补给地下水的水量变化范围为8.69～366.82 m³/(d·m),地下水补给地表水的水量变化范围为0.72～120.90 m³/(d·m);研究河段左岸为地下水补给地表水,单位河长净补给量为45.26 m³/(d·m);河段右岸为地表水补给地下水,单位河长净补给量为214.33 m³/(d·m);研究河段地下水与地表水水量交换以地表水补给地下水为主,地表水补给地下水的比例为55.14%。本研究可推动地下水与地表水交换量计算方法的发展,为流域或区域水资源评价提供必要的理论方法。  相似文献   

基于粒径分维值检验的松散介质渗透系数空间估值——以河北平原为例     

罗晓娟  徐恒力  叶含春 《地质科技情报》2007,26(1):87-90

为了克服松散含水介质非均质性各向异性对介质渗透系数的求取所带来的困难,以水动力条件为线索,利用互相关法寻找同一时代地层及同一时代地表水的主、次流方向.在主流方向上,建立了石家庄-束鹿-衡水段介质颗粒粒径随距离递变的数学模型,并以此为基础推求无钻孔地段的粒径.同时,运用分形理论,计算了中、上更新统地层主流方向平均粒径空间分布的分维值,并将其作为检验粒径估值正确与否的判据.通过不同岩性的渗透系数-深度模型,实现了用平均粒径预测渗透系数的全过程.  相似文献   

岩体裂隙等效水力隙宽的统计确定方法   总被引：2，自引：1，他引：1  

钟嘉高  梁杏  任志刚 《地质科技情报》2007,26(4):103-106

针对天然裂隙岩体等效水力隙宽确定的复杂性,利用平行板水流的立方定律及流量等效原理,推导出光滑平行板单裂隙的等效水力隙宽公式,同时对隙宽分布的频率特性进行了探讨,认为非均匀隙宽单裂隙的等效水力隙宽可近似取40%频率的隙宽,或取隙宽的几何平均值与调和平均值的算术平均值.  相似文献   

In situ stress field of eastern Australia     

R. R. Hillis  J. R. Enever  S. D. Reynolds 《Australian Journal of Earth Sciences》2013,60(5):813-825

New in situ data based on hydraulic fracturing and overcoring have been compiled for eastern Australia, increasing from 23 to 110 the number of in situ stress analyses available for the area between and including the Bowen and Sydney Basins. The Bowen Basin displays a consistent north‐northeast maximum horizontal stress (σ_H) orientation over some 500 km. Stress orientations in the Sydney Basin are more variable than in the Bowen Basin, with areas of the Sydney Basin exhibiting north‐northeast, northeast, east‐west and bimodal σ_H orientations. Most new data indicate that the overburden stress (σ_V) is the minimum principal stress in both the Bowen and Sydney Basins. The Sydney Basin is relatively seismically active, whereas the Bowen Basin is relatively aseismic. Despite the fact that in situ stress measurements sample the stress field at shallower depth than the seismogenic zone, there is a correlation between the stress measurements and seismicity in the two areas. Mohr‐Coulomb analysis of the propensity for failure in the Sydney Basin suggests 41% of the new in situ stress data are indicative of failure, as opposed to 13% in the Bowen Basin. The multiple pre‐existing structural grains in the Sydney Basin further emphasise the difference between propensity for failure in the two areas. Previous modelling of intraplate stresses due to plate boundary forces has been less successful at predicting stress orientations in eastern than in western and central Australia. Nonetheless, stress orientation in the Bowen Basin is consistent with that predicted by modelling of stresses due to plate boundary forces. Variable stress orientations in the Sydney Basin suggest that more local sources of stress, such as those associated with the continental margin and with local structure, significantly influence stress orientation. The effect of local sources of stress may be relatively pronounced because stresses due to plate boundary forces result in low horizontal stress anisotropy in the Sydney Basin.  相似文献   

Modeling three‐dimensional hydraulic fracture propagation using virtual multidimensional internal bonds     

K. Huang  Z. Zhang  Ahmad Ghassemi 《国际地质力学数值与分析法杂志》2013,37(13):2021-2038

Propagation of fractures, especially those emanating from wellbores and closed natural fractures, often involves Mode I and Mode II, and at times Mode III, posing significant challenges to its numerical simulation. When an embedded inclined fracture is subjected to compression, the fracture edge is constrained by the surrounding materials so that its true propagation pattern cannot be simulated by 2D models. In this article, a virtual multidimensional internal bond (VMIB) model is presented to simulate three‐dimensional (3D) fracture propagation. The VMIB model bridges the processes of macro fracture and micro bond rupture. The macro 3D constitutive relation in VMIB is derived from the 1D bond in the micro scale and is implemented in a 3D finite element method. To represent the contact and friction between fracture surfaces, a 3D element partition method is employed. The model is applied to simulate fracture propagation and coalescence in typical laboratory experiments and is used to analyze the propagation of an embedded fracture. Simulation results for single and multiple fractures illustrate 3D features of the tensile and compressive fracture propagation, especially the propagation of a Mode III fracture. The results match well with the experimental observation, suggesting that the presented method can capture the main features of 3D fracture propagation and coalescence. Moreover, by developing an algorithm for applying pressure on the fracture surfaces, propagation of a natural fracture is also simulated. The result illustrates an interesting and important phenomenon of Mode III fracture propagation, namely the fracture front segmentation. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献