Taxila LBM: a parallel,modular lattice Boltzmann framework for simulating pore-scale flow in porous media

The lattice Boltzmann method is a popular tool for pore-scale simulation of flow. This is likely due to the ease of including complex geometries such as porous media and representing multiphase and multifluid flows. Many advancements, including multiple relaxation times, increased isotropy, and others have improved the accuracy and physical fidelity of the method. Additionally, the lattice Bolzmann method is computationally very efficient, thanks to the explicit nature of the algorithm and relatively large amount of local work. The combination of many algorithmic options and efficiency means that a software framework enabling the usage and comparison of these advancements on computers from laptops to large clusters has much to offer. In this paper, we introduce Taxila LBM, an open-source software framework for lattice Boltzmann simulations. We discuss the design of the framework and lay out the features available, including both methods in the literature and a few new enhancements which generalize methods to complex geometries. We discuss the trade-off of accuracy and performance in various methods, noting how the Taxila LBM makes it easy to perform these comparisons for real problems. And finally, we demonstrate a few common applications in pore-scale simulation, including the characterization of permeability of a Berea sandstone and analysis of multifluid flow in heterogenous micromodels.     

A new theoretical method for analyzing confined dry granular flows

A granular body may deform in a continuous fashion such that the solid particles remain in close contact. Previous research works have always used the frictional Coulomb-like continuum treatment for analyzing granular bodies. However, this approach is only applicable for quasi-static conditions and cannot capture the complicated granular contact behavior of solid particles inside a failing granular body. This paper applies a revised Savage–Hutter equation to model granular flows moving down a confined, sloping channel. The Coulomb contact friction law is modified to consider the effect of the shear rate inside a granular body. This new method also considers the confinement effect of a sloping channel on granular flow mobility. The derived depth-averaged equations of motion bear a resemblance to nonlinear shallow-water wave equations. Results computed using the derived equations are compared with measurements from flume model tests, and consistency is found between the two.     

中国西部地区方沸石胶结相与碎屑岩次生优质储集层形成机制

含油气盆地碎屑岩次生优质储集层的发育,除了与广泛发育的方解石胶结相有关外,在中国西部还有一类与沸石族自生矿物胶结的成岩相有关,其中以溶蚀方沸石胶结相形成的次生优质储集层与油气聚集成藏又有十分密切的关联。经实例分析:与自生方沸石胶结有关的成岩相主要受控于扇(辫状河)三角洲前缘相带、火山岩物源和成岩早期的碱性环境;次生优质储集层的形成在时间上主要处于成岩中期A阶段,在空间上受控于邻近的烃源岩、圈闭和输导体的沟通。方沸石矿物在物理性质上具有抗压易溶(有机酸)作用,较其他常见的自生矿物(胶结物)相对优越,就方沸石易于溶蚀的特点而言是较为理想的自生矿物(胶结物)。其溶蚀对孔隙度的贡献率在50％～75％之间,发育几率占绝对优势。所以,在方沸石发育的相带中寻找以粒间溶蚀孔为主的优质储集层,具有可靠的操作性和现实性。     

沉积岩平行层面与垂直层面方向热导率与孔隙连通性之间的关系

根据对渤海湾盆地沾化凹陷泥质岩、粉砂岩、砂岩和白云岩平行层面与垂直层面方向热导率差异及孔隙连通性差异特征的研究,探讨了沉积岩平行层面与垂直层面方向热导率与孔隙连通性之间的关系。热导率实验结果表明:沾化凹陷泥质岩、粉砂岩和砂岩的平行层面方向热导率大于垂直层面方向热导率;白云岩的平行层面方向热导率小于垂直层面方向热导率。沉积岩平行层面与垂直层面方向热导率差异和孔隙连通性差异之间具有正相关性。孔隙特征分析结果表明:孔隙连通性较好方向,热导率也较高;孔隙连通性较差的的方向,热导率也较低。     

陕北煤矿区地质环境保护与恢复治理措施——以神府矿区郭家湾煤矿为例

郭家湾煤矿地处神府矿区新民采区西北部,为生态环境脆弱区。井田七层煤三个煤组开采后引发地面塌陷、含水层疏干、土地结构变化等地质环境问题,严重制约了地方经济的发展。为了煤炭资源可持续发展,研究陕北煤矿区地质环境保护与恢复治理措施,以1∶10000野外详细调查为依据,首先进行现状评估,然后针对矿山工程建设可能遭受、加剧、引发的矿山地质环境问题进行预测评估,最后提出了以移民搬迁、回填裂缝、生态恢复、建立矿山地质环境监测体系等为主的防治措施,让煤矿开采与矿山地质环境问题恢复治理并行。旨为陕北地区乃至全国煤矿,在保护地质环境中开发利用矿产资源提供依据和示范。     

淮北矿区煤层气资源勘查潜力评价

淮北矿区地处华东腹地,煤层气资源丰富,具备大规模开发的前景。选用层次分析法与灰色聚类法评价方法,通过构建层次结构模型、计算指标权重、制定参数分类评价标准、运用灰色聚类法判别灰类、计算综合优度值,对矿区16个评价单元的山西组和下石盒子组煤层气资源潜力进行了评价,结果显示,淮北矿区山西组煤层气勘查潜力中等,下石盒子组勘查潜力较好。     

淮南煤田潘三煤矿4-2煤中矿物对岩浆侵入的响应

以淮南煤田潘三煤矿4-2煤层中的岩浆岩、天然焦和煤样品为研究对象,利用X射线衍射仪和X射线荧光光谱仪测试其矿物组成和主量元素含量,研究地下岩浆侵入对煤中矿物成分的影响。结果表明:煤变质成天然焦后,其总的矿物含量增多,并伴随着后生矿物如微斜长石和镁绿泥石的出现。天然焦中后生矿物的组成元素可能来源于岩浆或热液中的主量元素如Si、Al、Na、K、Fe和Mg。     

聚合物钻井液在河北广宗勘探区深孔施工中的应用

在河北广宗勘探区深孔施工中,由于地表钻探用水水质影响及原钻井液体系不合适,造成钻井液性能指标达不到要求,性能难以控制,尤其失水量严重超标,接连发生孔内事故。针对存在的问题,经过分析和实验,采用化学方法处理钻探用水,选用聚合物钻井液,解决了深部钻孔施工存在的难点,收到了较好的效果,单孔深度达到1937.35m,并刷新了本单位及局属单位的钻探记录。     

External seismic stability of vertical geosynthetic-reinforced soil walls using pseudo-static method

Analysis of external stability of vertical geosynthetic-reinforced soil (GRS) walls is very important in the seismic prone zone. The scope of this paper is to obtain required minimum reinforcement length, L _min, for external seismic stability of vertical GRS walls by pseudo-static limit equilibrium method. Then, L _min can be calculated to resist sliding, eccentricity, and bearing capacity failure modes. The parameters considered include both horizontal and vertical seismic coefficients (k _h and k _v ), surcharge load (q), wall height (H) and the properties of retained backfill, GRS, and foundation soil. Results show that L _min against sliding failure mode, L _min,S , increases more quickly than that against the other two failure modes with the increase in k _h , q, or unit weight of retained backfill, γ _b , while L _min,S decreases more quickly than that against the other two failure modes with increase in friction angle of retained backfill, ? _b , or unit weight of GRS, γ _r . For the different failure modes, the effect of k _v on L _min is not identical with the change of k _h , and in addition, L _min/H will tend to remain unchanged with the increase in H. In general, L _min against bearing capacity failure mode, L _min,BC, is larger than L _min against the other two failure modes. However, L _min,BC will be less than L _min against eccentricity failure mode, L _min,E , for k _h exceeding 0.35, or friction angle of foundation soil, ? _f , exceeding 37°, and L _min,BC will also be less than L _min,S for friction angle of GRS, ? _r , being no more than 26°.     

Availability of coastal groundwater discharge as an alternative water resource in a large-scale reclaimed land,Korea

This study was conducted to identify the availability of coastal groundwater discharge (CGD), subsurface fluids flowing from inland through the coastal area to sea, as an alternative water resource for a large-scale reclaimed land. The behaviors of stable isotopes indicated that groundwater originated from inland precipitation and traveled as CGD along the coast line. Most of the groundwater samples collected from domestic wells installed along the old coast line were considered to be relatively fresh from the correlation analysis among chemical constituents. The average electrical conductivity (EC) values of the samples were identified as averaging 1,125–1,297 μS cm^?1, corresponding to appropriate crop growth. A weathered-rock layer in a small catchment within the reclaimed land was proved to be a main CGD pathway, with electrical resistivity anomalies ranging from 7 to 14 Ω m. Five monitoring wells were placed in this catchment to delineate the occurrence of CGD. Long-term vertical EC profiling results for the monitoring wells indicated that CGD occurs within a depth of 30 m below the ground surface. Annual monitoring data for groundwater level and EC demonstrated that the water quality of CGD was improved by introducing fresh terrestrial groundwater. A remarkable improvement in water quality (EC decrease of 900–1,600 μS cm^?1) of CGD was observed during the saline water pumping test that explains how CGD could be an alternative water resource for the reclaimed land.