Effect of pore connectivity on reflection amplitudes of an inhomogeneous wave in a composite porous solid saturated by two immiscible fluids

Present paper aims to study the phenomenon of reflection and transmission when an inhomogeneous wave strikes some discontinuity in a composite porous medium saturated by two immiscible viscous fluids. The incident wave splits into six reflected and six transmitted waves at the interface. All reflected and transmitted waves are inhomogeneous in nature with different directions of propagation vector and attenuation vector. A dimensionless parameter \(\varsigma \in [0, 1]\) is introduced to represent the extent of connection among the pores at the interface. Expression of Umov–Poynting vector is derived to obtain energy flux vector. Continuity of energy flux vector at the interface gives the required boundary conditions for the system. Connecting parameter \(\varsigma \) is also employed in boundary conditions to model the partial connection of pores at the interstices of two media. For numerical discussion we consider a porous medium composed of sandstone and ice, saturated with oil and water. The effect of parameter \(\varsigma \) and angle of incidence is determined numerically on the amplitude and the energy ratios of reflected and transmitted waves.     

Development of viscous fingering between mafic and felsic magmas: evidence from the Terra Nova Intrusive Complex (Antarctica)

Summary A wide range of types of contact morphology among mafic and felsic magmas are observed in outcrops on Vegetation Island (Terra Nova Intrusive Complex, Antarctica). Image analysis and fractal geometry techniques were applied for in-depth study of the mafic/felsic interface, with the aim of studying the origin of the varied morphologies. In particular, the length (IP_N) and fractal dimension (D_box) of interfaces were measured. Results indicate that there is a close exponential dependence of IP_N on D_box.The observed morphologies are identical to those observed during viscous fingering processes induced by the displacement of a more viscous fluid by a less viscous one. To test if viscous fingering was responsible in this case too, IP_N and D_box values were measured on viscous fingering structures obtained experimentally using various viscosity ratios (V_R) from the literature. Results indicate that, as in the natural case, there is an exponential dependence of IP_N on D_box, leading to the conclusion that the varied interface morphologies between mafic and felsic magmas are the result of viscous fingering dynamics. In addition, experimental studies clearly show that there is an exponential relationship between the viscosity ratio of fluids and the interface fractal dimension (D_box), and the ratio between the two types of magma was estimated using this relationship. It is shown that viscosity contrasts between mafic and felsic magmas varied considerably, ratios ranging from ca. 6 to 49. These results, together with outcrop evidence, provide indications regarding the evolution of the magmatic system, which generated the actual mafic/felsic associations on Vegetation Island.     

Heterogeneous modelling behaviour at an interface in porous media

The fluid flow induced by an incident wave at a discontinuity separating two porous media is governed by the hydraulic permeabilities of both media and that of the interface. In the context of Biot’s theory, we derive the time-harmonic Green’s function for the two half-space problem allowing incident fast and slow dilatational waves to assess the heterogeneous modelling behaviour for diverse hydraulic conditions. It is found that when at least one of the media is permeated with inviscid fluids, heterogeneous modelling simulates open boundary conditions. On the other hand, when the model is saturated with viscous fluids, the modelling reproduces restrained fluid flow whose values correspond to sealed pore interface conditions, in agreement with the theoretical results. Therefore the numerical technique models correctly the wave diffusion and propagation phenomena attendant at the boundary.     

辽西阜新金—多金属成矿区主要断裂带构造特征及其控矿作用

通过地质构造、显微构造和构造地球化学特征分析，论述了辽西阜新金－多金属成矿区主要断裂带的构造特征及其控矿作用，并着重指出NEE向、NNE向韧－脆性剪切带是本区主要的控矿和导矿构造，而NEE向构造的NW向次级断裂带和NNE向构造的EW向次级断裂带则为本区主要的容矿构造，现已探明的金－多金属矿体也多沿这两个方向展布。     

Fluid overpressure along an Oligocene out-of-sequence thrust in the Shimanto Belt,SW Japan

Out-of-sequence thrusts (OSTs) exposed in ancient accretionary prisms are considered as fossil analogs of present-day megasplay faults in subduction margins and can provide direct information about the conditions of deformation during thrust activity. In modern as well as in ancient accretionary prisms, first-order megasplay faults or OSTs truncate or merge with faults of lesser importance called second-order OSTs. Structural analysis of the Makinokuchi fault, a branch of an Oligocene to lower Miocene second-order OST in the Tertiary Shimanto Belt of central Kyushu, SW Japan, brings information about the conditions of deformation at the time of thrusting. The studied exposure shows that the fault footwall and, to a much lesser extent, the fault hanging-wall, consist of quartz-cemented syntectonic dilatant hydraulic breccias testifying to pore fluid pressures larger than the least principal stress component. The footwall sandstones are crossed by several centimeters thick quartz veins that merge with the footwall breccias. The continuity between the veins and the breccias suggest that the veins acted as conduits which likely collected fluids from the footwall side sandstones upward and toward the fault. Fluid inclusions indicate that the quartz cementing the breccias and that filling the feeder veins crystallized from similar fluids and under similar pressure and temperature conditions (245–285 °C and 5–8 km depth). These similarities suggest that the fluids responsible for syn-tectonic hydraulic brecciation were collected from the footwall through the conduits. The fluid inclusion trapping temperatures are close to the temperatures expected to be reached along the seismogenic zone. Our analysis shows that fluid overpressures can play a key role in the growth and activity of second-order OSTs in accretionary prisms and suggests that fluids collected along second-order OSTs or splay faults may flow upward along first-order OSTs or megasplay faults.     

Viscosity ratio estimates from natural Rayleigh-Taylor instabilities

The geometry of natural gravity-driven overturns potentially contains information on viscosity ratios of the material involved. Irregularities at the interface between two superposed horizontal fluid layers tend to amplify and grow with time with a regular wavelength if the top iayer is more dense, Published formulas (Woidt), for such Rayleigh-Taylor (R-T) instabilities have been used to construct graphs which can constrain the effective viscosity ratio at the onset of the overturn. Only the initial thicknesses of the two layers and the final wavelength of the structures need be known. The theory is based on an infinite two-layer system of Newtonian fluids between rigid boundaries with no slip or free slip at contacts.
Geological examples included are; soft-sediment dewatering on the scale of decimetres in northern Norway, an unusual instability between basalt and sand in Australia of hundreds of metres, salt diapirism of kilometres in Germany and the Persian Gulf, and granitic diapirism from northern Sweden of tens of kilometres. Viscosity ratios in the order of 1–10 are found for soft sediments with the top layer as the more viscous. Ratios inferred from patterns of salt diapirs are very large, with the overburden approximately 10³-⁻10⁵ times more viscous than the salt. The viscosities of granites are indicated as lower than their overburden and with ratios of 10³.     

The study of reflection/transmission phenomena in a corrugated interface between two magnetoelastic transversely isotropic media

In this paper, an integrated approach has been carried out with an intent to study the reflection and transmission phenomena of plane SH-type wave on a corrugated interface separating two magnetoelastic transversely isotropic half-space. Closed form expressions for reflection and transmission coefficients have been derived for both plane and corrugated surface. Rayleigh’s method of approximation have been incorporated to deduce equations for the first- and second-order approximation of corrugation. Analytical solutions for both the half-spaces have been worked out. All possible cases have been apprehended specially for anisotropic and an isotropic medium. The effects of magnetoelastic coupling parameter, angle at which the wave crosses the magnetic field, corrugation amplitude, frequency factor and wave number have been explained by collaborating with graphical analysis.     

粉质黏土-混凝土衬砌接触面冻结强度影响因素显著性分析

冻土与结构接触面的力学性能对寒区工程研究具有重要意义, 容易受到外界因素干扰。为研究不同因素对粉质黏土与混凝土衬砌接触面冻结强度的影响, 设计考虑含水率、 冻结温度、 冻结时间的三水平正交直剪试验。试验结果表明： 在冻结初期, 提高土体含水率, 降低冻结温度, 延长冻结时间有利于接触面早期冻结强度的提升, 且这种作用效果是线性的。基于显著性分析理论, 发现含水率、 冻结温度、 冻结时间对接触面早期冻结强度的显著性等级均为极显著, 三者之中, 冻结温度影响最大, 含水率影响最小。采用多元线性回归, 忽略因素之间的二阶交互作用, 建立了考虑含水率、 冻结温度、 冻结时间作用的粉质黏土与混凝土接触面的冻结强度预测模型, 为类似工程提供科学参考。     

Compaction-driven fluid flow in viscoelastic rock

Abstract

Compaction driven fluid flow is inherently unstable such that an obstruction to upward fluid flow (i.e. a shock) may induce fluid-filled waves of porosity, propagated by dilational deformation due to an effective pressure gradient within the wave. Viscous porosity waves have attracted attention as a mechanism for melt transport, but are also a mechanism for both the transport and trapping of fluids released by diagenetic and metamorphic reactions. We introduce a mathematical formulation applicable to compaction driven flow for the entire range of rheological behaviors realized in the lithosphere. We then examine three first-order factors that influence the character of fluid flow: (1) thermally activated creep, (2) dependence of bulk viscosity on porosity, and (3) fluid flow in the limit of zero initial connected porosity. For normal geothermal gradients, thermally activated creep stabilizes horizontal waves, a geometry that was thought to be unstable on the basis of constant viscosity models. Implications of this stabilization are that: (1) the vertical length scale for compaction driven flow is generally constrained by the activation energy for viscous deformation rather than the viscous compaction length, and (2) lateral fluid flow in viscous regimes may occur on greater length scales than anticipated from earlier estimates of compaction length scales. In viscous rock, inverted geothermal gradients stabilize vertically elongated waves or vertical channels. Decreasing temperature toward the earth’s surface can induce an abrupt transition from viscous to elastic deformation-propagated fluid flow. Below the transition, fluid flow is accomplished by short wavelength, large amplitude waves; above the transition flow is by high velocity, low amplitude surges. The resulting transient flow patterns vary strongly in space and time. Solitary porosity waves may nucleate in viscous, viscoplastic, and viscoelastic rheologies. The amplitude of these waves is effectively unlimited for physically realistic models with dependence of bulk viscosity on porosity. In the limit of zero initial connected porosity, arguably the only model relevant for melt extraction, travelling waves are only possible in a viscoelastic matrix. Such waves are truly self-propagating in that the fluid and the wave phase velocities are identical; thus, if no chemical processes occur during propagation, the waves have the capacity to transmit geochemical signatures indefinitely. In addition to solitary waves, we find that periodic solutions to the compaction equations are common though previously unrecognized. The transition between the solutions depends on the pore volume carried by the wave and the Darcyian velocity of the background fluid flux. Periodic solutions are possible for all velocities, whereas solitary solutions require large volumes and low velocities. © Elsevier, Paris     

Insights into the collection and emplacement of granitic magma based on 3D seismic images of normal fault-related salt structures

ABSTRACT
3D seismic data from salt structures in the Southern North Sea show how viscous fluids tend to migrate towards and flow upwards along the footwalls of active normal faults where zones of low stress are inherent. This process provides one potential solution for the question of how enough space is created to accommodate large granitic bodies in the crust. During episodes of regional extension, granitic melt may collect in the lower crust below normal detachments in domed areas from which magma may be sourced. Subsequently, in the upper crust, emplacement can occur in fault-bounded laccoliths at rates commensurate with other geological processes.     

东准噶尔蒙西斑岩铜钼矿床脉体特征及其形成机制

蒙西斑岩铜钼矿床位于东准噶尔伊吾县琼河坝花岗岩北侧,以发育细脉、网脉状矿化为特征。脉体类型多样,包括石英脉、石英硫化物脉和硫化物脉等。根据脉体力学成因机制,其又可划分为水压破裂充填脉和构造破裂充填脉。构造破裂充填脉体按破裂形成的位错特征有正断与逆断两种。地表石英脉产状陡立,明显受断裂控制。脉体体积分数统计结果显示流体富集区呈向北缓倾的带状分布于深100～400m范围内,并与矿化富集带有较好的对应关系。脉体富集带内,脉体以共轭形式存在,一组为倾角较小的逆断破裂充填脉,另一组为倾角较大的正断破裂充填脉,他们可能为缓倾逆冲剪切带的次级破裂与充填脉,即富矿带内脉体是沿矿区低角度逆冲断裂次级破裂面充填的。矿区地表及深部(400m以下)脉体以陡立为主,矿化较弱。矿区流体的运移具先沿水压直立破裂往上运移,进入剪切带后沿剪切带次级破裂侧向和向上运移,并在剪切带中富集成矿。低温矿物组合脉体穿切高温矿物组合脉体的特征说明脉体形成过程矿区处于隆升构造环境,这对斑岩铜矿成矿有利。     

筋-土界面力学特性的水平循环剪切试验研究

加筋材料与土界面力学特性是加筋土结构设计和稳定性分析的重要依据。在大型多功能界面剪切仪上,进行了土工格栅和标准砂界面的水平循环剪切试验,分析循环荷载作用下筋-土界面力学特性,研究循环荷载特征对界面剪切力和法向位移的影响,探讨水平循环荷载作用下筋-土相互作用规律。结果表明,随循环次数增加,剪切力-位移关系曲线逐渐重合,剪切力峰值变化不明显,界面剪切模量增大并趋于稳定;循环剪切位移幅值和剪切速率对剪切强度的影响不明显,但对试样法向位移有较大影响;水平循环剪切过程中试样发生明显的法向位移,加筋可有效限制试样的竖向变形。     

西秦岭八方山—二里河Pb-Zn矿区硅质岩的微区成分特征及演化

硅质岩以SiO2为主,并广泛发育于造山带内。激光Raman、SEM-EDS分析结果显示,研究区硅质岩微组构特征记录了造山作用、成矿作用和重结晶作用等的改造。SEM分析显示,研究区硅质岩中石英颗粒较小且呈紧密堆积,这些特征符合热水沉积快速结晶的特点。激光Raman分析结果中,石英颗粒自边缘向中心核部的拉曼特征峰(463cm-1)经高斯拟合(Gaussian Fitting)后尖锐程度逐渐递增,半高宽(FWHM)递减,反映了石英自身的重结晶作用。在石英颗粒边界与碳酸盐脉的接触部位,结晶程度和有序度自内向外呈递增趋势,反映了石英颗粒遭受了流体的影响,SEM分析结果中石英颗粒的絮状反应残余边界进一步支持了该结论。在SEM下,金属硫化物在裂隙和空隙中发生沉淀并呈线状和零星浸染状分布,与应力作用及杂质矿物有着密切的关系。另外,据硅质岩中的脉体穿插关系,指示西秦岭存在碳酸盐(铁白云石)热液→金属硫化物热液→硅质-碳酸盐(方解石)热液的多阶段流体演化,与东秦岭流体作用具有较好的相似性。     

Bitumens associated with lead,zinc and flnorite ore minerals in North Derbyshire,England

A bitumen deposit in north Derbyshire, England, is described and studied using the techniques of gas chromatography, infra-red spectrophotometry and elemental analysis. The bitumens are associated with lead-zinc-fluorite ore minerals concentrated along an unconformable contact between the Carboniferous Limestone and the overlying shales. Three varieties of bitumen are distinguished and are compared with the dispersed bitumens in both the Carboniferous Limestone and the Edale Shales, with crude oils believed to be derived from these shales, and with inclusions in the hydrothermal mineral fiuorite, associated with the limestone-shale contact. One of the varieties recognized, a brittle brown solid, contains aliphatic hydrocarbons very similar to those of the shales and to those of the crude oils likely derived from these shales. The other two are a viscous oil and a brittle, black solid, which geologic evidence suggests were originally derived from the limestone. They contain a very complex mixture of aliphatic hydrocarbons quite different from those in the limestone, however. Aliphatic hydrocarbons isolated from fluorite, which is hydrothermally formed, in nearby massive deposits are almost identical to those of the limestone and shale. This observation indicates that alteration of the aliphatic hydrocarbons in the sediment by hydrothermal fluids is an unlikely explanation for the origin of the complex mixtures of hydrocarbons in the viscous oil and brittle, black solid. On the basis of indirect evidence of high nitrogen content, presence of unsaturated hydrocarbons, and suitability of the environment for bacterial growth, it is suggested that selective bacterial alteration of the aliphatic hydrocarbons is the main process responsible for their present composition. A geologic history for the deposit is postulated that involves a two-phase introduction of bitumens. Each phase is suggested to have begun as a pulse of warm saline fluids migrating along the shale-limestone unconformity passed through the topographic high at Windy Knoll. Microbial oxidation of the bitumens may have taken place during the deposition process or, more likely, as a recent secondary oxidation process.     

Tectonic inheritance and kinematic strain localization as trigger for the formation of the Helvetic nappes,Switzerland

The Helvetic nappes in Switzerland consist of sediments, which have been sheared off and thrust over the crystalline basement of the European passive continental margin during Alpine orogeny. Their basal shear zones usually root above the external crystalline massifs. However, the mechanisms that initiated the shear zones and the associated nappe formation are still debated. We perform two-dimensional numerical simulations of the shearing of linear viscous fluids above a linear viscous fluid with considerably higher viscosity (quasi-undeformable). The boundary between the fluid, mimicking the sediments, and the quasi-undeformable fluid, mimicking the basement, exhibits geometrical perturbations, mimicking half-grabens. These geometrical perturbations can trigger significant strain localization and the formation of shear zones within the linear viscous fluid although no rheological softening mechanism is active. This kinematic, ductile strain localization is caused by the half-grabens and the viscosity ratio between basement and sediments. The viscosity ratio has a strong control on the kinematics of strain localization, whereas the depth of the half-grabens has a weak control. For sediment viscosities in the order of 10²¹ Pas and typical half-graben geometries of 5 km depth and 25 km width the localization generates (a) low-angle shear zones at the basement-sediment interface, but also entirely within the sediments, (b) horizontal transport >10 km associated with the shear zones, (c) shear zones with thickness in the order of 100 m, (d) an ordered stacking of model nappes and (e) shear zones that root above the basement. The results suggest that tectonic inheritance in the form of half-grabens and associated kinematic strain localization could have been the triggering mechanism for Helvetic nappe formation, and not rheological softening mechanisms, which might, however, have subsequently further intensified shear localization significantly.     

Finite Element-Analysis of an Overfold Using the Viscous Fluid Model

Symmetrical folding of a single layer and folding process of an overfold are analyzed by using the finite element method. In the analysis, the constituent substances are assumed to be isotropic and incompressible and to behave as Newtonian fluids. The penalty function method is used to deal with incompressibility, The results of the analysis are compared with the data measured in the field and with those obtained by the elastopla(?)e finite element analysis, It shows that the overfold geometry may be well simulated by using the viscous fluid model. Thus, a possible mechanism of the formation of overfold is proposed.     

辽宁白云金矿床成矿构造与成矿结构面特征浅析

白云金矿是辽宁省青城子矿集区大型金矿床之一,文章根据矿山资料及野外详细调查,主要对矿区的成矿构造与成矿结构面进行了梳理。成矿构造主要为成矿期活动的断裂和褶皱,利用侵入其中的含矿岩脉（石英斑岩、花岗斑岩、闪长玢岩和煌斑岩）确定成矿期断裂主要为近东西向、北西向和北东向断裂,成矿期的褶皱主要为姚家岭—天桥岭—李家堡子倒转向斜,该倒转向斜是青城子一带印支期岩浆活动形成的推覆构造形成;成矿结构面主要分为四种,即侵入岩与围岩接触面、层间断裂面或接触面、围岩断裂或裂隙面和水压致裂结构面,并从力学角度分析了这四种结构面的产生机制。最后,讨论了在强烈的由南向北挤压应力作用下,含矿热液沿不同的张性裂隙侵位形成不同矿体样式。     

中国云南金顶铅锌矿碎屑灌入体和水力压裂构造的观察及流体动力学分析

金顶铅锌矿曾被认为是形成于相对被动地受构造与沉积相控制的盆地流体,类似于许多其他沉积岩容矿的金属矿床。然而,一些近期研究显示在金顶矿区露天采场存在砂岩灌入体、侵位角砾岩和水力压裂构造,说明矿床可能形成于超压流体的爆发性释放。文中报道在地下坑道（跑马坪和峰子山）中新观察到的与超压流体有关的构造,这些构造显示比露天采场有更...     

Montecarlo DLA-type simulations of wetting effects in fluid displacement in porous media

In this work, we report diffusion-limited aggregation (DLA)-type Montecarlo computations of a stochastic model of displacement of a viscous fluid by another that preferentially wets a porous medium, for the case when both fluids are immiscible in the absence of buoyancy forces. The model has the aim to simulate cooperative invasion processes found in experiments of immiscible wetting displacement. The model considers the nonlocal effects of the Laplacian pressure field and the capillary forces via hydrodynamic equations in the Darcy regime with a boundary condition for the pressure at the interface. The boundary condition contains two different types of disorder: the capillary term, which constitutes an additive random disorder, and a term containing an effective random surface tension, which couples to a curvature (it constitutes a multiplicative random term that carries nonlocal information of the whole pressure). We generate different displacement patterns for different setting of the parameters of the model. We analyze these patterns by studying the scaling properties of the interface that separate the two fluids and calculating the fractal dimension of the interface. The results show the existence of three distinct regimes of scaling. One regime at the smallest-length scales is due to the multiplicative random disorder together with the nonlocal coupling; it reveals itself in a roughness exponent α ≈ 0.80. Additionally, we find a DLA-type scaling regime with a roughness exponent α ≈ 0.60 at the largest scales and intermediate scaling regime with α ≈ 0.70 corresponding to invasion percolation with trapping. Each regime has definite scaling ranges that depend on the capillary number and the relative wetting tendency of the fluids. The behavior of the fractal dimensions of the interfaces of the aggregates constitutes a further confirmation of the existence of three scaling regimes and the multi-self-affinity of the perimeter of the interface boundaries.