Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation          下载免费PDF全文

Tongbi Tu  Ali Ercan  M. Levent Kavvas 《水文研究》2018,32(10):1406-1419

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.  相似文献   

Geochemistry of S-type granitic rocks from the reversely zoned Castelo Branco pluton (central Portugal)     

I.M.H.R. Antunes  A.M.R. Neiva  M.M.V.G. Silva  F. Corfu 《Lithos》2008,103(3-4):445-465

The zoned pluton from Castelo Branco consists of Variscan peraluminous S-type granitic rocks. A muscovite>biotite granite in the pluton's core is surrounded successively by biotite>muscovite granodiorite, porphyritic biotite>muscovite granodiorite grading to biotite=muscovite granite, and finally by muscovite>biotite granite. ID-TIMS U–Pb ages for zircon and monazite indicate that all phases of the pluton formed at 310 ± 1 Ma. Whole-rock analyses show slight variation in ⁸⁷Sr/⁸⁶Sr_310 Ma between 0.708 and 0.712, Nd_310 Ma values between − 1 and − 4 and δ¹⁸O values between 12.2 and 13.6. These geological, mineralogical, geochemical and isotopic data indicate a crustal origin of the suite, probably from partial melting of heterogeneous Early Paleozoic pelitic country rock. In detail there is evidence for derivation from different sources, but also fractional crystallization linking some of internal plutonic phases. Least-squares analysis of major elements and modelling of trace elements indicate that the porphyritic granodiorite and biotite=muscovite granite were derived from the granodiorite magma by fractional crystallization of plagioclase, quartz, biotite and ilmenite. By contrast variation diagrams of major and trace elements in biotite and muscovite, the behaviours of Ba in microcline and whole-rock δ¹⁸O, the REE patterns of rocks and isotopic data indicate that both muscovite-dominant granites were probably originated by two distinct pulses of granite magma.  相似文献   

滇西北保山核桃坪铅锌矿床成矿流体来源的同位素地球化学证据   总被引：11，自引：2，他引：9  

薛传东  韩润生  杨海林  杨志明  田世洪  刘勇强  郝百武 《矿床地质》2008,27(2):243-252

滇西北核桃坪铅锌矿是保山地块内沿澜沧江断裂带发育的一个大型铅锌多金属矿床,矿体严格受近SN向的断裂带及层间破碎带控制,赋存于上寒武统核桃坪组和沙河厂组大理岩化灰岩及矽卡岩中。文章在分析矿床地质特征的基础上,通过对V1矿体的氢、氧、碳、硫、硅和铅同位素地球化学的示踪研究,探讨了成矿溶液中水、碳、硅、硫和铅的来源以及成矿流体的演化问题。研究表明,V1矿体中矿石的石英δ18O水值变化范围为6·1‰~7·6‰,平均为6·7‰,其δDV-SMOW值变化范围为-100‰~-108‰,平均为-104‰,说明成矿流体主要来自深部岩浆分异水,并在后期成矿作用过程中有地层建造水和大气降水的加入;矿石中方解石的δ13CV-PDB和δ18O水值分别为-6·6‰~-5·9‰和5·0‰~5·2‰,反映成矿流体的碳、氧具有深部流体与围岩发生交代反应的特征。硅、硫和铅同位素具深部岩浆或地幔来源的特点,但遭受后期壳源物质的混染。据此提出该矿床是与深部岩浆热液或幔源流体贯入有关的构造热液脉型矿床,其成矿作用受控于深部过程的驱动和控制,可为指导地质找矿提供依据。  相似文献   

安徽铜陵狮子山矿田岩浆岩锆石SHRIMP定年及其成因意义   总被引：36，自引：4，他引：32  

徐晓春  陆三明  谢巧勤  柏林  储国正 《地质学报》2008,82(4):500-509

铜陵狮子山矿田发育大量岩浆岩,且与矿田中的铜金多金属成矿关系密切。锆石SHRIMP同位素精确定年表明,矿田中的岩浆侵位年龄在132.4～142.9Ma之间,即晚侏罗世—早白垩世,属燕山早期晚阶段。矿田岩浆岩体是在同期岩浆活动中多次侵位形成的,岩浆侵入活动可以划分为分别起始于140Ma前后和约136Ma的早晚两次。从岩浆上升侵位到冷却结晶的时间间隔均较短,但其中白芒山辉石二长闪长岩冷却史相对较长,且经历了早期深部岩浆房中的分离结晶作用和后期构造脉动、岩浆上升侵位、减压受热、早期晶体再熔蚀及冷却结晶的过程。结合主量元素和微量元素地球化学研究认为,狮子山矿田岩浆演化的后期,即起源于上地幔或下地壳的原生岩浆在同化了壳源物质并聚集到岩浆房中以后,在滞留的过程中发生了一定程度的分离结晶作用,但尚未固结,成分上显示了一定的带状分布,在区域构造应力松弛及构造事件诱发下,随机地沿发育的构造裂隙先后上升侵位,冷凝结晶。  相似文献   

力马河镍矿Re-Os同位素研究   总被引：10，自引：0，他引：10  

陶琰 《地质学报》2008,82(9):1292-1304

四川力马河镍矿是峨眉山大火成岩省一个重要的岩浆硫化物矿床。本文通过对其主要岩、矿石类型Re、Os及其同位素组成的分析,综合探讨了成矿岩体原始岩浆性质、矿石硫化物成因、成矿机制及Re-Os同位素等时线年龄。结果表明,力马河镍矿不同类型岩矿石样品初始Os同位素组成是不均一的,富硫化物的网脉状矿石及其选纯硫化物Os同位素组成初值差异较小,其等时线年龄为265±35 Ma、与岩体锆石SHRIMP年龄263±3 Ma基本相当;硫化物含量较低的岩、矿石样品间初始Os同位素组成差异较大,其表观等时线年龄大于成矿年龄。分析认为,岩矿样品初始Os同位素组成的不均一是由含较高放射成因187Os丰度的硫化物熔体和含较低放射成因187Os丰度的硅酸盐熔体不同比例混合造成的。混合模型分析表明,硫化物含量超过30％的矿石样品初始187Os/188Os基本接近,硫化物含量低于30％的岩矿石样品初始187Os/188Os随硫化物含量上的不同差异很大,为岩浆硫化物矿床Re-Os等时线年龄可能出现多组年龄解的现象提供了一种可能的解释。成矿岩体中含放射成因187Os丰度最低的岩石样品γOs（t=260Ma）在5左右、Cu/Pd比值在7000左右,表明是基本没有受到地壳混染及硫化物熔离影响的原始岩浆结晶分异产物,估计原始岩浆Os含量在1×10－9左右,为苦橄质岩浆。矿石硫化物Re/Os比值显著高于任何赋矿橄榄岩,γOs（t=260Ma）高达110左右,综合分析揭示了力马河镍矿硫化物为二次熔离成因,模式分析认为,矿石硫化物是由原始岩浆经历R=2000左右的硫化物熔离后、其亏损岩浆再经R=200左右的硫化物熔离形成,与二次熔离相对应,成矿岩浆也经历了两次混染作用,分别为上、下地壳7％左右的混染。  相似文献   

Experimental constraints on the effects of pressure and H2O on the fractional crystallization of high-Mg island arc basalt     

Morihisa Hamada  Toshitsugu Fujii 《Contributions to Mineralogy and Petrology》2008,155(6):767-790

H₂O-undersaturated melting experiments of synthesized basalt (SiO₂ = 50.7 wt.%, MgO = 8.3 wt.%, Mg# = 60) were conducted at fO₂ corresponding to NNO+1 and NNO−1 to clarify the effects of pressure (2–7 kbar) and H₂O on fractional crystallization in island arcs. H₂O content was ranged from nominally anhydrous to 4.4 wt.%. Differentiation trends, namely the liquid lines of descent, change sensitively according to pressure-H₂O relations. Tholeiitic differentiation trends are reproduced with H₂O ≤ ∼2 wt.% in primary magma. With such quantities of H₂O, fractional crystallization is controlled by olivine + plagioclase at 2 kbar. Increasing the pressure from 2 to ≥4 kbar induces early crystallization of orthopyroxene instead of olivine and therefore SiO₂ enrichment in the residual melts is suppressed. Increasing H₂O (≥ ∼2 wt.% in primary magma) stabilizes clinopyroxene relative to orthopyroxene and/or magnetite. Although the phase relations and proportions strongly depend on fO₂ and H₂O content, differentiation trends are always calc-alkaline.  相似文献   

Depths of Partial Crystallization of H2O-bearing MORB: Phase Equilibria Simulations of Basalts at the MAR near Ascension Island (7 11{degrees}S)     

Almeev  Renat; Holtz  Francois; Koepke  Jurgen; Haase  Karsten; Devey  Colin 《Journal of Petrology》2008,49(1):25-45

Phase equilibria simulations were performed on naturally quenchedbasaltic glasses to determine crystallization conditions priorto eruption of magmas at the Mid-Atlantic Ridge (MAR) east ofAscension Island (7–11°S). The results indicate thatmid-ocean ridge basalt (MORB) magmas beneath different segmentsof the MAR have crystallized over a wide range of pressures(100–900 MPa). However, each segment seems to have a specificcrystallization history. Nearly isobaric crystallization conditions(100–300 MPa) were obtained for the geochemically enrichedMORB magmas of the central segments, whereas normal (N)-MORBmagmas of the bounding segments are characterized by polybariccrystallization conditions (200–900 MPa). In addition,our results demonstrate close to anhydrous crystallization conditionsof N-MORBs, whereas geochemically enriched MORBs were successfullymodeled in the presence of 0·4–1 wt% H₂O in theparental melts. These estimates are in agreement with direct(Fourier transform IR) measurements of H₂O abundances in basalticglasses and melt inclusions for selected samples. Water contentsdetermined in the parental melts are in the range 0·04–0·09and 0·30–0·55 wt% H₂O for depleted and enrichedMORBs, respectively. Our results are in general agreement (within±200 MPa) with previous approaches used to evaluate pressureestimates in MORB. However, the determination of pre-eruptiveconditions of MORBs, including temperature and water contentin addition to pressure, requires the improvement of magma crystallizationmodels to simulate liquid lines of descent in the presence ofsmall amounts of water. KEY WORDS: MORB; Mid-Atlantic Ridge; depth of crystallization; water abundances; phase equilibria calculations; cotectic crystallization; pressure estimates; polybaric fractionation  相似文献   

Pressures of Crystallization of Icelandic Magmas   总被引：1，自引：0，他引：1  

Kelley  Daniel F.; Barton  Michael 《Journal of Petrology》2008,49(3):465-492

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass  相似文献   

南海东北部下地壳高速层的成因探讨   总被引：2，自引：0，他引：2  

刘安  武国忠  吴世敏 《地质论评》2008,54(5):609-616

通过对南海北部大陆边缘地壳结构分析,指出南海东北部存在下地壳高速层,大致分布在112°E~120°E,19°N~22°N的陆坡和拉张程度大的陆架地区,呈NEE向延伸,在海底地震仪剖面上最大的厚度有8km,向南海海盆方向减薄。通过对比综合分析认为,高速层物质组成是底侵作用形成的熔岩垫,由于伸展作用,南海海底扩张（30Ma）前后底侵作用形成了熔岩垫,并促使南海北部大陆边缘地壳抬升,导致区域性抬升剥蚀。  相似文献   

内蒙阿右旗铁板井镍矿矿床特征及成因探讨     

谢成连  刘蕾  冷莹莹  文中田  闫海忠 《地质与勘探》2008,44(6):27-30

铁板井绦矿床赋存于辉橄岩体中,辉橄岩体具有全岩矿化的特点,岩体及镍矿体的形成经历了三种不同的成矿过程:首先是岩浆深源液态重力分异作用,而后发生岩浆深源熔离-贯入(成岩、成矿)作用,形成岩体和镍矿体,最后在岩浆期后热液作用的叠加下完成整个成矿过程.该矿床是典型的主要由岩浆熔离作用形成的岩浆矿床.  相似文献