Stochastic analysis of immiscible displacement of the fluids with arbitrary viscosities and its dependence on support scale of hydrological data

Stochastic analysis is commonly used to address uncertainty in the modeling of flow and transport in porous media. In the stochastic approach, the properties of porous media are treated as random functions with statistics obtained from field measurements. Several studies indicate that hydrological properties depend on the scale of measurements or support scales, but most stochastic analysis does not address the effects of support scale on stochastic predictions of subsurface processes. In this work we propose a new approach to study the scale dependence of stochastic predictions. We present a stochastic analysis of immiscible fluid–fluid displacement in randomly heterogeneous porous media. While existing solutions are applicable only to systems in which the viscosity of one phase is negligible compare with the viscosity of the other (water–air systems for example), our solutions can be applied to the immiscible displacement of fluids having arbitrarily viscosities such as NAPL–water and water–oil. Treating intrinsic permeability as a random field with statistics dependant on the permeability support scale (scale of measurements) we obtained, for one-dimensional systems, analytical solutions for the first moments characterizing unbiased predictions (estimates) of system variables, such as the pressure and fluid–fluid interface position, and we also obtained second moments, which characterize the uncertainties associated with such predictions. Next we obtained empirically scale dependent exponential correlation function of the intrinsic permeability that allowed us to study solutions of stochastic equations as a function of the support scale. We found that the first and second moments converge to asymptotic values as the support scale decreases. In our examples, the statistical moments reached asymptotic values for support scale that were approximately 1/10000 of the flow domain size. We show that analytical moment solutions compare well with the results of Monte Carlo simulations for moderately heterogeneous porous media, and that they can be used to study the effects of heterogeneity on the dynamics and stability of immiscible flow.     

内蒙古科右前旗一带梅勒图组火山岩U-Pb测年及地球化学特征

通过对科右前旗一带梅勒图组火山岩LA-ICP-MS锆石U-Pb定年和元素地球化学分析,获得梅勒图组安山质晶屑熔结凝灰岩及安山岩的喷发年龄分别为125.6±1.6Ma和120.2±2.2Ma,指示梅勒图组火山岩的形成时代介于120～126Ma之间,形成于早白垩世中期。岩石具有高的SiO_2(60.4%～65.3%)和全碱(7.0%～9.5%)含量,铝饱和指数A/CNK为0.89～0.98,属于钾玄岩及高钾钙碱性系列;稀土元素总量为124×10~(-6)～169×10~(-6),(La/Yb)N值为5.83～13.4,负Eu异常微弱(δEu=0.84～1.05);在原始地幔标准化微量元素蛛网图上,样品富集Rb、Ba、Th、U、K、La、Nd,亏损Nb、Ta、P和Ti,Sr/Y值介于19.4～29.7之间。其岩石源区为富集地幔,岩浆上升过程中受到地壳物质的混染,岩石形成于板内伸展环境,其地球动力学背景可能与古太平洋板块的俯冲有关。     

IS THE GRENVILLE PROVINCE AN ANCIENT ANALOGUE OF THE HIMALAYAN BELT?

IS THE GRENVILLE PROVINCE AN ANCIENT ANALOGUE OF THE HIMALAYAN BELT?1　All埁ｇｒｅＣＪ ,2 4others.StructureandevolutionoftheHimalayan Tibetorogenicbelt[J].Nature ,1984,30 7:17～ 2 2 . 2　BurgJP ,DavyP ,MartinodJ .Shorteningofanaloguemodelsofthecontinentallithosphere :Newhypothesisforthefor mationoftheTibetanplateau[J].Tectonics ,1994,13:475～ 483. 3　BurtmanVS ,MolnarP .GeologicalandgeophysicalevidencefordeepsubductionofcontinentalcrustbeneaththePamir[C]…     

象山港盐度分布和水体混合I.盐度分布和环流结构

利用１９８１—１９９０年的实测水文资料分析了象山港与混合密切相关的盐度分布和环流结构,并对盐度锋面出现成因进行了探讨。研究结果表明,象山港盐度受长江冲淡水活动影响较大,冬、夏口门附近盐度分别为２４—２５左右和３１以上,湾顶附近盐度季节变化不大。冬季湾内盐度空间梯度较小。夏季湾内有两个盐度锋面出现,水体垂向略有层化但并无明显的跃层出现。此外,象山港狭湾内以重力环流为主,狭湾外以水平环流为主。     

The Southland Front,New Zealand: Variability and ENSO correlations

A detailed analysis of the Southland Front, a shelf-break system off the southeast coast of South Island, New Zealand is presented. The position, temperature, temperature range and width of the front are determined using a new statistical front detection algorithm and 21 years worth of Advanced Very High Resolution Radiometer satellite sea surface temperature (SST) data. Overall, the front is strongest (highest SST gradients) in the summer and winter, and the across front gradient decreases northward in all seasons, consistent with an equatorward decrease in stability and divergence of isobaths. The surface expression of the front moves further offshore during the winter months and is found closest inshore in the summer. Seasonality of the front is strongly controlled by the annual cycle of subtropical and subantarctic water mass temperatures. Both the temperature and strength of the front are interannually variable, and correlated with the El Niño-Southern Oscillation (ENSO); they both decrease during El Niño, and increase during La Niña events. ENSO indices lead changes in the fronts temperature by up to 6 months. Conversely, the gradient may change up to 6 months in advance of peak ENSO indices. The strength and sign of correlations is seasonally dependent.     

40Ar/39Ar traverse — Grenville Front Tectonic Zone to Britt Domain,Grenville Province,Ontario, Canada*

Abstract ⁴⁰Ar/³⁹Ar data (on hornblende, muscovite and K-feldspar) are presented for samples from the western Grenville Province taken along a 140-km traverse from the Grenville Front into the Britt domain. Our interpretation is based on 28 new analyses, synthesized with 20 previously reported from the traverse area. In regions where comparisons are possible, muscovite and (large domain) K-feldspar apparent ages appear similar (at c. 920–930 Ma), but throughout the traverse, these are c. 60–70 Myr younger than the hornblende ages. The inferred cooling rate over the c. 350–500°C temperature range, c.2°C Myr^-1, is appropriate for exhumation controlled by post-orogenic erosional unroofing. At the Grenville Front Tectonic Zone (GFTZ) — Britt domain boundary there is a c. 25-Myr offset in both hornblende and muscovite/K-feldspar ages. We interpret the lower ages in the Britt domain to reflect variations in crustal thickness and geothermal gradient between the flank and interior of a thick orogen. The argon data from the GFTZ are interpreted in the context of an asymmetric crustal-scale antiformal structure developed during a late episode of convergence. Hornblende from rocks on either side of the core of the antiform has an apparent age of c. 990 Ma, our estimate of the age of the compressional event. In the west, we infer that these date the short-lived thermal event associated with the development of the crustal-scale antiform previously postulated. In the east, the ages reflect the cooling of material brought toward the surface in the flank of the antiform. Hornblendes from the antiform core appear to contain excess radiogenic argon. We suggest that this was the ambient argon in rocks transported from depth that was subsequently trapped when the rocks cooled rapidly.     

A stratigraphic network across the Subtropical Front in the central South Atlantic: Multi-parameter correlation of magnetic susceptibility, density, X-ray fluorescence and δO records

Computer aided multi-parameter signal correlation is used to develop a common high-precision age model for eight gravity cores from the subtropical and subantarctic South Atlantic. Since correlations between all pairs of multi-parameter sequences are used, and correlation errors between core pairs (A, B) and (B, C) are controlled by comparison with (A, C), the resulting age model is called a stratigraphic network. Precise inter-core correlation is achieved using high-resolution records of magnetic susceptibility κ, wet bulk density ρ and X-ray fluorescence scans of elemental composition. Additional δ¹⁸O records are available for two cores. The data indicate nearly undisturbed sediment series and the absence of significant hiatuses or turbidites. After establishing a high-precision common depth scale by synchronously correlating four densely measured parameters (Fe, Ca, κ, ρ), the final age model is obtained by simultaneously fitting the aligned δ¹⁸O and κ records of the stratigraphic network to orbitally tuned oxygen isotope [J. Imbrie, J. D. Hays, D. G. Martinson, A. McIntyre, A. C. Mix, J. J. Morley, N. G. Pisias, W. L. Prell, N. J. Shackleton, The orbital theory of Pleistocene climate: support from a revised chronology of the marine δ¹⁸O record, in: A. Berger, J. Imbrie, J. Hays, G. Kukla, B. Saltzman (Eds.), Milankovitch and Climate: Understanding the Response to Orbital Forcing, Reidel Publishing, Dordrecht, 1984, pp. 269-305; D. Martinson, N. Pisias, J. Hays, J. Imbrie, T. C. Moore Jr., N. Shackleton, Age dating and the orbital theory of the Ice Ages: development of a high-resolution 0 to 300.000-Year chronostratigraphy, Quat. Res. 27 (1987) 1-29.] or susceptibility stacks [T. von Dobeneck, F.Schmieder, Using rock magnetic proxy records for orbital tuning and extended time series analyses into the super-and sub-Milankovitch Bands, in: G. Fischer, G. Wefer (Eds.), Use of proxies in paleoceanography: Examples from the South Atlantic, Springer-Verlag, Berlin (1999), pp. 601-633.]. Besides the detection and elimination of errors in single records, the stratigraphic network approach allows to check the intrinsic consistency of the final result by comparing it to the outcome of more restricted alignment procedures. The final South Atlantic stratigraphic network covers the last 400 kyr south and the last 1200 kyr north of the Subtropical Front (STF) and provides a highly precise age model across the STF representing extremely different sedimentary regimes. This allows to detect temporal shifts of the STF by mapping δMn / Fe. It turns out that the apparent STF movements by about 200 km are not directly related to marine oxygen isotope stages.     

Surface buoyant plumes from melting icebergs in the Labrador Sea

Canada׳s Department of Fisheries and Oceans (DFO) conducts annual surveys in the Labrador Sea along the repeat hydrography line AR7W. The occupation of the AR7W line in May 2013 was followed by the experiment aimed at resolving the imprint of melting drifting icebergs on the upper layer thermohaline characteristics in the Labrador Sea. We present high-resolution observations around two icebergs conducted with the towed undulating platform Moving Vessel Profiler (MVP). The first iceberg drifted in relatively warm water of Atlantic origin (~2.5–3.1 °C) off Greenland, while the second iceberg was on the Labrador shelf in cold water below 0 °C. Both icebergs had a lengthscale of O(100 m). In both cases surface buoyant plumes fed by melt water and attached to the iceberg were observed. The plumes were evident in the anomalous thermohaline characteristics of the seawater. Their density anomalies were sufficiently strong to produce visible frontal structures, which imply a development of the intrinsic dynamics associated with a plume. The first plume formed over a time interval of ~10 h, while the second plume formed over several days and extended for more than 1 km (tenfold the iceberg׳s size). Strong vertical displacements of the pycnocline were observed near the second iceberg. They are interpreted as the internal wave wake. This interpretation is based on the temporal scale of these oscillations (local buoyancy frequency), as well as on the spatial orientation of these waves with respect to the iceberg drift relative to the pycnocline. The observed internal waves partially overlapped with the plume and affected its structure. The saline seawater splashing by swell contributed to the surface melting of the icebergs. Scaling analysis of the second plume suggests that it could be in the “rotational” dynamic regime with recirculating anticyclonic flow.     

Geomorphic inferences from regolith thickness, chemical denudation and CRN erosion rates near the glacial limit, Boulder Creek catchment and vicinity, Colorado

Granitic regolith, developed in the Boulder Creek catchment and adjacent areas, records a history of deep weathering, some of which may predate Quaternary time. Field and well-log measurements of weathering, chemical denudation and rates of erosion derived from ¹⁰Be cosmogenic radionuclide (CRN) data help to quantify rates of landscape change in the post-orogenic Rocky Mountains. The density of oxidized, fractured bedrock ranges from 2.7 to about 2.2 g cm^− 3, saprolite and grus have densities between 2.0 and 1.8 g cm^− 3, and 30 soil samples averaged 1.6 ± 0.2 g cm^− 3. Highly weathered regolith in 540 wells averages 3.3 m thick, mean depth to bedrock in 1661 wells is 7 m, and the weathered thickness exceeds 10 m in relatively large local areas east of the late Pleistocene glacial limit. Thickness of regolith shows no simple relationship to rock type or structure, local slope, or distance from channels. Catchments in the vicinity of the Boulder Creek have an average CRN erosion rate of 2.2 ± 0.7 cm kyr^− 1 for the past 10,000 to 40,000 yr. Annual losses of cations and SiO₂ vary from about 2 to 5 g m^− 2 over a runoff range of 10 to nearly 160 cm.Using measured rates in simple box models shows that if a substantial fraction of void space is created by volume expansion in the weathering rock materials, 7 m of weathered rock materials could form in as little as 230 kyr. If density loss results mainly from chemical denudation and some volume expansion, however, the same weathering profile would take > 1340 kyr to form. Rates of erosion measured by CRN could be balanced by the rate of soil formation from saprolite if the annual solute loss from soil is 2.0 g m^− 2 and 70% of the density decrease from saprolite to grus and soil results from strain. Saprolite, however, forms from oxidized bedrock at a far slower rate and rates of saprolite formation cannot balance soil and grus losses to erosion. The zone of thick weathered regolith is likely an eroding relict landscape. The undulating surface marked by relatively low relief and tors is not literally a topographic surface of Eocene, Oligocene or Miocene age unless it was covered with deposits that were removed in Pliocene or Quaternary time.     

Interaction of Orographic Disturbance with Front

1.IntroductionTheimportanceoforographiceffectsonfrontwasrecognizedintheearly20thcentury.Butforthecomplexityofthisproblem,theinvestigationoforographiceffectsonfrontfromdynamicalviewpointisnottakenuntilthe1980s.Bannon(1983)derivedanalyticalsolutionsforthequasi--geostrophicfrontforcedbyahorizontalwinddeformationfieldthatmovesoveratwo--dimensionalmountainridge.Thesolutionsshowthatasacoldfrontapproachestheridge,itweakens,relativetotheflat--bottomsolution,andthefrontstrengthensasitmovesdowntothelees…