孔隙尺度多孔介质流体流动与溶质运移高性能模拟          下载免费PDF全文

周鸿翔  郑延丰  吴劳生  陈铖  曾令藻 《水科学进展》2020,31(3):422-432

深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10^-13m^2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。  相似文献   

Geology, Mineralogy, and Formation Environment of the Disseminated Gold-Silver Telluride Bulawan Deposit, Negros Occidental, Philippines     

Victor B. MAGLAMBAYAN  †  Daizo ISHIYAMA  Toshio MIZUTA  Akira IMAI  Yohei ISHIKAWA 《Resource Geology》1998,48(2):87-104

Abstract: The Bulawan deposit is located in the porphyry copper belt of southwest Negros island, Philippines. Propylitic, K–feldspar, sericitic, and carbonate alteration types can be distinguished in the deposit. Propylite alteration occurs mainly in Cretaceous-Eocene andesitic lavas and agglomerates while K–feldspar, sericite and carbonate alteration types occur mostly in the Middle Miocene dacite porphyry breccia pipes and stocks which were intruded into the andesites. K-feldspar zones occur in the inner parts of the sericitized zone. Sericite alteration overprinted the propylitized and K-feldspar alteration zones, at lower temperature than epidote and chlorite in the propylitized zone. Carbonate alteration is associated with the mineralization in the center of the breccia pipes and along faults. Mineralization consists of gold-silver telluride ores that are hosted by the carbonate– and sericite-altered dacite porphyry breccia pipes. The Bulawan ores occur mainly as disseminations, but unlike many epithermal gold deposits, lack classical epithermal colloform and crustiform quartz veins. The ore minerals are sphalerite, galena, chalcopyrite, pyrite and tetrahedite-tennantite with minor amounts of electrum, calaverite, petzite, sylvanite, hessite, tellurobismuthite, coloradoite, altaite, and rucklidgeite. Electrum and telluride minerals are associated mostly with calcite and dolomite-ankerite minerals. Fluid inclusions in quartz and calcite in clasts of propylitized andesite in the breccia pipes homogenize from about 300° to 400°C while fluid inclusions in quartz, calcite and sphalerite within the dacite porphyry breccia pipes homogenize between 300° to 310°C. The ores were formed around 300°C from hydrothermal solutions with salinity of about 6. 6 wt % NaCl equivalent. The presence of sylvanite and calaverite as intergrowths with each other, and the Ag content of calaverite are consistent with the above temperature estimate. Based on paragenesis, the Bulawan deposit formed in a pyrite-stable environment, with pH between 3. 4 and 5. 5, f_O2 between 10^-32 to 10^-30 atm, f_S2 between 10^-9.8 to 10^-7.8 atm, f_Te2 between 10^-8.9 to 10^-6.5 atm, and total sulfur content about 10^-2.8 molal. The dominant reduced sulfur species in the ore solutions may have been H₂S_(aq), and the likely aqueous tellurium species were H₂Te_(aq) and H₂TeO_3(aq). The ore minerals in the Bulawan deposit were probably formed by mixing of slightly saline and low salinity fluids.  相似文献   

A study of overflow simulations using MPAS-Ocean: Vertical grids,resolution, and viscosity     

《Ocean Modelling》2015

MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. Numerical simulations are carried out using three of the vertical coordinate types available in MPAS-Ocean, including z-star with partial bottom cells, z-star with full cells, and sigma coordinates. The results are first benchmarked against other models, including the MITgcm’s z-coordinate model and HIM’s isopycnal coordinate model, which are used to set the base case used for this work. A full parameter study is presented that looks at how sensitive overflow simulations are to vertical grid type, resolution, and viscosity. Horizontal resolutions with 50 km grid cells are under-resolved and produce poor results, regardless of other parameter settings. Vertical grids ranging in thickness from 15 m to 120 m were tested. A horizontal resolution of 10 km and a vertical resolution of 60 m are sufficient to resolve the mesoscale dynamics of the DOME configuration, which mimics real-world overflow parameters. Mixing and final buoyancy are least sensitive to horizontal viscosity, but strongly sensitive to vertical viscosity. This suggests that vertical viscosity could be adjusted in overflow water formation regions to influence mixing and product water characteristics. Lastly, the study shows that sigma coordinates produce much less mixing than z-type coordinates, resulting in heavier plumes that go further down slope. Sigma coordinates are less sensitive to changes in resolution but as sensitive to vertical viscosity compared to z-coordinates.  相似文献   

Basin structures and sediment accumulation in the Baikal Rift Zone: Implications for Cenozoic intracontinental processes in the Central Asian Orogenic Belt     

《Gondwana Research》2017

In this paper we present a review of sedimentological, geomorphological, lithological, geochronological and geophysical data from major, minor and satellite basins of the Baikal Rift Zone (BRZ) and discuss various aspects of its evolution. Previously, the most detailed sedimentological data have been obtained from the basins of the central BRZ, e.g., Baikal, Tunka and Barguzin, and have been used by many scientists worldwide. We add new information about the peripheral part and make an attempt to provide a more comprehensive view on BRZ sedimentation stages and environments and their relations to local and regional tectonic events. A huge body of sedimentological data was obtained many years ago by Soviet geologists and therefore is hardly accessible for an international reader. We pay tribute to their efforts to the extent as the format of a journal paper permits. We discuss structural and facial features of BRZ sedimentary sequences for the better understanding of their sedimentation environments. In addition, we review tectono-sedimentation stages, neotectonic features and volcanism of the region. Finally, we consider the key questions of the BRZ evolution from the sedimentological point of view, in particular, correlation of Mesozoic and Cenozoic basins, bilateral growth of the Baikal rift, Miocene sedimentation environment and events at the Miocene/Pliocene boundary, Pliocene and Pleistocene tectonic deformations and sedimentation rates. The data from deep boreholes and surface occurrences of pre-Quaternary sediments, the distribution of the Pleistocene sediments, and the data from the Baikal and Hovsgol lakes sediments showed that 1) BRZ basins do not fit the Mesozoic extensional structures and therefore hardly inherited them; 2) the Miocene stage of sedimentation was characterized by low topography and weak tectonic processes; 3) the rifting mode shifted from slow to fast at ca. 7–5 Ma; 4) the late Pleistocene high sedimentation rates reflect the fast subsidence of basin bottoms.  相似文献   

Efficient nonlinear modeling of rainfall-runoff process using wavelet compression     

《Journal of Hydrology》2007,332(3-4):442-455

  相似文献   

Finite element modelling of laterally loaded lime–cement columns using a damage plasticity model     

《Computers and Geotechnics》2012

The behaviour of laterally loaded lime–cement columns in a shear box was studied. Laboratory tests are presented together with numerical analyses where the columns are simulated by a concrete damage plasticity model that considers stiffness degradation. Seven model tests were investigated where the columns were installed in a single column pattern and in rows with different column overlap in order to investigate the influence of the degree of overlapping of the columns in the rows. The results of the numerical evaluations showed good agreement with the experimental shear stress–displacement relation and a good accuracy with respect to the fractures developed.  相似文献   

The variability of the large-amplitude internal wave field on the Australian North West Shelf     

Paul Van Gastel  Gregory N. Ivey  Michael J. Meuleners  Jason P. Antenucci  Oliver Fringer 《Continental Shelf Research》2009,29(11-12):1373-1383

Observations are presented of large-amplitude internal waves (LAIWs) generated by the steepening of the internal tide on the Australian North West Shelf (NWS) over a 4-month period extending from strongly stratified summer conditions to weakly stratified winter conditions. The observations are from a site in water depth of 124 m where current and temperature measurements were made from a fixed vertical mooring and a benthic L-shaped spatial array. The observations show the LAIWs at this site to be characterized by strong seasonal variability, with energetic LAIWs of depression being dominant during summer and weaker LAIWs of elevation being dominant during the winter months as the stratification weakens, the upper mixed layer deepens, and the thermocline is close to the bottom. Waves were also seen to propagate from a range of directions towards the observation site. Modeling using the Regional Ocean Modeling System (ROMS v2.1) revealed that internal tide generation in the area occurred at water depths of between 400 and 600 m along an arc of approximately 120 km in length, some 70 km to the northwest of our experimental site. The results demonstrate both the 3D nature as well as the seasonal variation of the LAIW field.  相似文献   

Data errors and hydrological modelling: The role of model structure to propagate observation uncertainty     

《Advances in water resources》2013

  相似文献   

Significance of indigenous Eocene larger foraminifera Discocyclina dispansa in Western Foothills,Central Taiwan: A Paleogene marine rift basin in Chinese continental margin     

《Journal of Asian Earth Sciences》2013

The Western Foothills of Taiwan was known to be composed of Late Oligocene to Pleistocene shallow marine strata continuously deposited on the stable passive Chinese continental margin without significant stratigraphic break. Here we present multiple micropaleontological evidences, including occurrence of larger foraminifera Discocyclina dispansa ex. interc. sella-dispansa and calcareous nannoplanktons, to show that there are Middle Eocene marine strata (first named as the Chungliao Formation) exposed in the Tsukeng anticline of the Western Foothills, central Taiwan. Occurrences of intact tests with thin delicate outer rims and well-preserved embryonic chambers suggest that the Discocyclina dispansa ex. interc. sella-dispansa (Lutetian to Bartonian in the Tethys region) are buried indigenously on shallow inner shelf during an episodic transgression in the Early Middle Eocene. The conclusion is consistent with a biostratigraphy study of calcareous nannoplanktons (Zones NP14–15) in the shale/sandstone alternations overlying the Discocyclina-bearing bed of the Chungliao Formation and calcareous nannofossils of Zone NP16 integrated with an age dating of 38.8 ± 1 Ma (Late Middle Eocene) on zircon grains of the overlying Pinglin Tuff. The Middle Eocene syn-rift sequences (Chungliao Formation and Pinglin Tuff) exposed along the Tsukeng anticline are unconformably covered by the latest Oligocene–Miocene post-rift sequence, a scenario similar to what have been drilled in the East China Sea-Taiwan Strait-South China Sea. This rift basin (named as the Nantou Basin) is sitting on the Peikang Basement High margin which further extends southwestward to the Central Uplift of the Pearl River Mouth Basin in the northern slope of the South China Sea. The present work documents a hitherto unknown occurrence of the exposed early Tertiary marine rift basin sequence in the Western Foothills of Taiwan. The study extends our knowledge of the Western Foothills geohistory from the Late Oligocene downward to the Early Middle Eocene. The occurrence of the Paleogene Nantou rift basin in the Western Foothills may also suggest that there could have similar Paleogene rift sequences exposed in other parts of the Taiwan mountain belt like the Hsüehshan Range and the Central Range east of the Western Foothills.  相似文献   

The importance of including conductivity and dielectric permittivity information when processing low-frequency GPR and high-frequency EMI data sets     

《Journal of Applied Geophysics》2013

The full solution for the wavenumber equation in electromagnetic (EM) theory is available, but not routinely used when processing ground penetrating radar (GPR), and electromagnetic induction (EMI) data. The wavenumber approach is important as it is used for the development of concepts such as skin depth and phase velocity, as well as being the basis for more complete interpretation of EM data sets. Approximations that make the solution simpler are common, and sufficiently accurate, provided that the underlying assumptions are not grossly violated. With the advent of lower-frequency GPR systems (25 MHz and below) and higher-frequency EMI systems (greater than 100 kHz) such approximations need to be re-examined. This paper reviews the full wavenumber expression and then compares phase velocity and skin depth equations based on approximations with the equations for the same parameter based on the full solution. This comparison allows the conditions under which the assumptions are valid to be refined. In this paper it is shown that for GPR surveys conducted under transition band conditions, the error in phase velocity estimates based on low-loss assumptions may be 40%. Similarly, for EMI surveys the skin depth estimation errors may be more than 30% when the equation based on quasi-static assumptions is used instead of the full solution.  相似文献