Stability Zone Research Articles

Impact of gas arising from gas hydrate decomposition on the effective permeability of bottom sediments of the Arctic shelf: a laboratory simulation

Relevance. The necessity to investigate the mechanism of gas (methane) flow occurring during decomposition of gas hydrate formations within the layers of bottom sediments. This process is evident through substantial releases of methane bubbles in extensive shallow regions of the Russian Arctic shelf and has the potential to alter the equilibrium of atmospheric methane, a critical greenhouse gas. Aim. To quantitatively investigate, within the context of gas transport in bottom sediments, the impact of free and bound gas within the pores on the nonlinearity of filtration flows. Methods. The model samples with filtration properties similar to those of the bottom rocks of the East Siberian Arctic shelf. During the experiments, a specific amount of gas was injected into the saturated model sample, followed by the measurement of its effective permeability during a slow decrease in pore pressure gradient. The obtained curve was interpreted within the framework of the threshold gradient model. Results. The experiments yielded curves showing the dependency of the threshold gradient magnitude on gas saturation for several samples. It was found that the threshold gradient linearly increases with the growth in gas fraction within the pore space. Already at a gas fraction of approximately 0.02, this value in the experiments reached 0.01 MPa/m, corresponding to the hydrostatic pressure gradient in water. This suggests that even areas with relatively low gas saturation may be impermeable to convective fluid flow. This possibility should be considered when creating models for bubble gas transporting through the rock layers of bottom sediments. Furthermore, the existence of gas-saturated zones with threshold gradients can significantly impact the vertical profile of pore pressure and lead to a reassessment of the depth of gas hydrate stability zones.

Dynamics of vortex and anti-vortex solitons in a vectorial cubic-quintic complex Ginzburg-Landau equation

In this paper, we present a study of vortex and anti-vortex dynamics within a complex cubic-quintic Ginzburg-Landau vector equation (CCQGLVE). We employ a variational approach to address the analytical aspects, and the results obtained are subsequently confirmed numerically. The vortex vector (VV) and the anti-vortex vector (anti-VV) are defined with topological charges: m = 1 for VV and m = − 1 for anti-VV. Our investigation reveals that the stability zone map corresponds to the region where greater stability can be achieved for the two studied solutions. Notably, the radius of the vortex craters experiences variations either an increase or decrease depending on the competition between the coupling parameters associated with cubic and quintic cross-phase modulation (XPM). During the propagation, the interaction between a fundamental soliton and anti-VV transforms the soliton into a vortex after a short time, but both finally undergo self-confinement which probably will generates solitons. In the case of the interaction between a VV and a fundamental soliton, we observed a self-confinement and a transformation into solitons. Considering the interaction between a VV and an anti-VV, we found that both solutions are also self-confined but the anti-VV solution turns into a soliton faster than the VV solution. This confirms that the anti-VV is the better solution that can be managed with system coupling parameters than the VV one.

Dynamic accumulation of a high-grade gas hydrate system: insights from the trial production gas hydrate reservoir in the Shenhu area, northern South China Sea

The ultimate enrichment level and quantity of gas hydrate resources are influenced by the dynamic process of accumulation and preservation. High-resolution 3-D seismic data, logging while drilling (LWD), pressured coring, and in situ testing were used to characterize the dynamic accumulation and preservation of the trial production high-grade gas hydrate reservoir (HGGHR) in the Shenhu area. Through seismic variance analysis and ant-tracking, we found that newly identified mud diapir-associated faults with three development stages controlled the migration and accumulation of gas hydrate and shifted the base of the gas hydrate stability zone (BGHSZ), resulting in dynamic accumulation and dissociation of gas hydrates. The recognized double bottom simulating reflectors (BSRs) were concluded to have been formed due to the shift of the BGHSZ caused by the variational equilibrium conditions. The interval between the double BSRs was inferred to be a disequilibrium zone where gas recycling occurred, contributing to the coexistence of gas hydrates and free gas and the dynamic formation of the HGGHR. Multiple gliding faults formed within the GHSZ in the late period have altered the HGGHR and control the present thickness and distribution of the gas hydrates and free gas in the hanging wall and footwall. Under the influence of geothermal fluids and the fault system associated with the mud diapir, the HGGHR experienced dynamic accumulation with three stages, including early accumulation, medium-term adjustment, and late alteration and preservation. We conclude that four factors affected the formation, distribution, and occurrence of the HGGHR: the geothermal fluids accompanying the deep mud diapir below the reservoir, the dual supply of thermogenic gas and biogenic gas, the recycling of hydrate gas beneath the BGHSZ, and the post-gas hydrate faults developed within the GHSZ. A geological model illustrating the dynamic formation of the trial production HGGHR was proposed, providing a reference for future exploration of HGGHRs with a great production potential in deepwater settings.

Influence of the shape of a conducting chamber on the stability of rigid ballooning modes in a mirror trap

MHD stabilization of flute and ballooning modes in an axisymmetric mirror trap is studied under the assumption of strong finite Larmor radius effect that suppresses all perturbations with azimuthal numbers m⩾2 and makes the m = 1 mode ‘rigid’. The rigid mode can be effectively suppressed using perfectly conducting lateral wall without any additional means of stabilization or in combination with end MHD anchors. Numerical calculations were carried out for an anisotropic plasma produced in the course of neutral beam injection into the minimum of the magnetic field at the right angle to the trap axis. The stabilizing effect of the conducting shell made of a straightened cylinder is compared with a proportional chamber, which, on an enlarged scale, repeats the shape of the plasma column. It is confirmed that for convincing wall stabilization of the rigid modes, the plasma beta (β, the ratio of the plasma pressure to the magnetic field pressure) must exceed some critical value βcr2 . When conducting lateral wall is combined with conducting end plates imitating MHD end anchors, there are two critical betas and, respectively, two stability zones β<βcr1 and β>βcr2 that can merge, making the entire range 0<β<1 of betas allowable for stable plasma confinement. The dependence of the critical betas on the plasma anisotropy, mirror ratio, width of the vacuum gap between the plasma column and the lateral wall, radial pressure profile and the axial magnetic field profile is examined.

High-frequency harmonics suppression in high-speed railway through magnetic integrated LLCL filter.

The traction converter modulation generates switching-frequencies current harmonics. The trapped filters can eliminate these switching harmonics, reducing total inductance and filter size. Nonetheless, in comparison with the typical inductor-capacitor-inductor (LCL) filter, the trap inductor needs a larger magnetic core. Moreover, the trapped filter has not been analyzed in the traction systems. This paper proposes a magnetic integrated inductor-trap-inductor (LLCL) filter to decrease the filter's size and investigate its application in traction converters. In fact, the application range of this filter is quite broad, and it can be used in various electrical power systems, including industrial power systems, renewable energy systems, transportation systems, and building power systems. The LC-trap may be formed by connecting the equivalent trap inductor, introduced through the magnetic coupling between inverter-side and grid-side inductors, in series with the filter capacitor. Furthermore, for H-bridge unipolar pulse width modulation (PWM) traction converters, the prominent switching harmonics are concentrated at the double switching frequencies. Therefore, the stability zone is expanded by moving the resonance above the Nyquist frequency. The presented filter's features and design are thoroughly analyzed. The proposed method is finally validated by the MATLAB/Simulink simulation and hardware-in-the-loop (HIL) experimental results. Compared to the discrete windings, the integrated ones can save two magnetic cores. Furthermore, the proposed filter can meet IEEE criteria with 0.3% for all the harmonics and total harmonic distortion (THD) of 2.15% of the grid-side current.

Plastic zone characteristics and rapid stability evaluation indexes of circular tunnels subjected to non-hydrostatic stress: A novel complex variable solution incorporating 3D rock mass strength

The morphological characteristics and distribution range of the plastic zone around tunnels, influenced by the true three-dimensional high geostress environment in deep engineering, constitute fundamental reference data for both tunnel design and the control of surrounding rock stability. Based on the GZZ three-dimensional rock mass strength criteria and the complex variable method, a novel semi-analytical method for solving the plastic zone of a circular tunnel subjected to non-hydrostatic stress is presented. The conformal mapping is employed to intricately transform the complex boundaries of plastic zones characterized by diverse shapes onto the unit circle boundary within the image plane. This framework proves invaluable for the precise analytical process of the complex plastic zone morphology and the plastic influence range within the surrounding rock under complex stress conditions. In contrast to traditional analytical methods, the efficacy and versatility of our proposed approach are substantiated. Notably, it demonstrates unique advantages, particularly in the context of non-hydrostatic stress fields and varied rock mass parameters, with respect to predicting both surrounding rock stability and plastic zone morphology. Through the analysis of numerous cases, the study reveals the systematic patterns of influence that in-situ stress and surrounding rock characteristic parameters have on the shape and maximum depth of the plastic zone. The results indicate that (1) there are notable disparities in the morphology and extent of the plastic zone between the GZZ and GHB criteria; (2) the transition of the plastic zone in surrounding rock from U-shaped to V-shaped is significantly influenced by the geostress-to-strength ratio Kσ; (3) there is a notable linear correlation observed between the maximum plastic zone thickness rmax and both Kσ and the stress concentration factor KSCF, importantly, this correlation is independent of far-field stress and the brittleness characteristics of the surrounding rock. These findings provide a novel discriminative method for the rapid analysis of surrounding rock stability in deep circular tunnels.

Performance Assessment of Bismuth Ferrite Nanoparticles in Enhancing Oilwell Cement Properties: Implications for Sustainable Construction

Summary Oilwell cement ensures wellbore stability and isolates zones while bearing casing load and formation pressure. Its properties, crucial in extreme downhole conditions, include compressive strength, fluid loss resistance, and durability. In the present work, bismuth ferrite nanoparticles (BFO NPs) were synthesized using the sol-gel method and used as an additive in oilwell cement. The synthesized BFO NPs were characterized using Fourier transform infrared (FTIR), X-ray diffraction (XRD), field-emission scanning electron microscope (FESEM), and dynamic light scattering (DLS) techniques to analyze the functional groups, crystalline structure, morphological features, and hydrodynamic size distribution. Tests at 70°C and 2,000 psi revealed that 1% by weight of cement (BWOC) BFO NPs increased compressive strength by ~136% and reduced fluid loss to ~64% compared with base cement. It can be conjectured that the exposed facets of BFO NPs containing oxygen act as nucleating sites that promote the ordering of the silicate tetrahedra, thereby increasing the strength and crystallinity and reducing the water loss. The experimental results confirm that the BFO NPs can improve the properties of oilwell cement slurry at high-pressure, high-temperature (HPHT) conditions. This research underscores the potential of BFO NPs as sustainable additives for optimizing oilwell cement performance under challenging HPHT conditions, paving the way for advancements in sustainable construction practices.

Contrasting Methane Seepage Dynamics in the Hola Trough Offshore Norway: Insights From Two Different Summers

AbstractThis study investigates the temporal variations in methane concentration and flare activity in the Hola trough (offshore Norway) during May 2018 and June 2022. Between these time periods, methane seep activity exhibits 3.5 times increase, as evidenced by hydroacoustic measurements. As the seep area in the Hola trough is constantly within the hydrate stability zone, the observed increase cannot be attributed to migration of its shallow boundary due to temperature increase. However, a combination of low tide conditions resulting in a lower sediment pore pressure and a bottom water temperature increase resulting in a lower methane solubility is likely to explain the increase in the number of seeps observed in June 2022. The hypothesis of tide influence is supported by data collected from a piezometer deployed and recovered during the cruise showing that the tidal effect was observed 3 m below the seafloor. Despite the numerous methane seeps detected, methane concentration and gas flow rates near the seafloor were low (&lt;19 nM and &lt;70 mL min−1, respectively) compared to other areas with methane seep activity. This is likely due to strong currents rapidly dispersing methane in the water column. Sub‐seafloor investigations identified pathways for gas migration in methane seep areas, influenced by topography. This study provides valuable insights into the temporal dynamics of methane concentrations, flare activity, and gas distribution in the Hola trough, contributing to our understanding of offshore methane dynamics in the region.

Slow Dynamics of Hydrate Systems Revealed by a Double BSR

AbstractDetermining how gas hydrate distribution evolved along continental margins in the past is essential to understanding its evolution in the future. Moreover, hydrate decomposition has been linked to several catastrophic events, including some of the largest submarine landslides on Earth and the massive release of greenhouse gases into the ocean. Offshore Romania, the presence of a second bottom‐simulating reflector (BSR) provides an opportunity to gain valuable insights into hydrate dynamics since the Last Glacial Period (LGP). We conducted transient modeling of hydrate thermodynamic stability by merging in‐situ observations with indirect assessments of sea‐bottom temperature, thermal conductivity, salinity, sedimentation rate, and sea‐level variations. We reveal a strong correlation between the BSRs and the base of the Gas Hydrate Stability Zone (GHSZ) during both the present and LGP periods. The gradual evolution of the GHSZ over the past 34 ka presented here supports a conceptual model that excludes catastrophic environmental scenarios.

Geometric modeling of phase ordering for the isotropic–smectic A phase transition

BackgroundLiquid crystal (LC) mesophases have an orientational and positional order that can be found in both synthetic and biological materials. These orders are maintained until some parameter, mainly the temperature or concentration, is changed, inducing a phase transition. Among these transitions, a special sequence of mesophases has been observed, in which priority is given to the direct smectic liquid crystal transition. The description of these transitions is carried out using the Landau–de Gennes (LdG) model, which correlates the free energy of the system with the orientational and positional order.MethodologyThis work explored the direct isotropic-to-smectic A transition studying the free energy landscape constructed with the LdG model and its relation to three curve families: (I) level-set curves, steepest descent, and critical points; (II) lines of curvature (LOC) and geodesics, which are directly connected to the principal curvatures; and (III) the Casorati curvature and shape coefficient that describe the local surface geometries resemblance (sphere, cylinder, and saddle).ResultsThe experimental data on 12-cyanobiphenyl were used to study the three curve families. The presence of unstable nematic and metastable plastic crystal information was found to add information to the already developed smectic A phase diagram. The lines of curvature and geodesics were calculated and laid out on the energy landscape, which highlighted the energetic pathways connecting critical points. The Casorati curvature and shape coefficient were computed, and in addition to the previous family, they framed a geometric region that describes the phase transition zone.Conclusion and significanceA direct link between the energy landscape’s topological geometry, phase transitions, and relevant critical points was established. The shape coefficient delineates a stability zone in which the phase transition develops. The methodology significantly reduces the impact of unknown parametric data. Symmetry breaking with two order parameters (OPs) may lead to novel phase transformation kinetics and droplets with partially ordered surface structures.

Vibrational dynamics of a subjected system to external torque and excitation force

This work focuses on the dynamical response of a novel auto-parametric two-degrees-of-freedom (DOF) dynamical system when it is subjected to external torque and force. It consists of a forced oscillator with nonlinear stiffness that is connected to a linear spring. In accordance with the system’s DOF, an organizing system of equations of motion (EOM) is produced using Euler–Lagrange’s equations. To obtain the analytical approximate solutions (ASs) of the equations of this system, the multiple scales technique (MST) is employed. The provided solutions are being juxtaposed with the numerical solutions (NSs) for comparison to show how closely they agree, as well as the precision of the MST. In view of removing secular terms, the system’s solvability criteria are obtained. All arising resonance cases are classified, and two of them are then examined at once. The relevant system of modulation equations is solved numerically to illustrate the advantageous influence of diverse factors on the behavior of the modified phases and amplitudes. Based on the values of these parameters, the frequency response curves (FRCs) are drawn to explore different zones of stability and instability. A wide variety of values for the system’s parameters shows that its behavior is steady. The achieved outcomes are deemed novel, as the MST has been applied on a new dynamical model. The results have broad relevance and can be implemented in various engineering contexts, including the analysis of ship movements, in which nonlinear coupling between rolling and pitching motions can be found.

Shifts in controls and abundance of particulate and mineral-associated organic matter fractions among subfield yield stability zones

Abstract. Spatiotemporal yield heterogeneity presents a significant challenge to agricultural sustainability efforts and can strain the economic viability of farming operations. Increasing soil organic matter (SOM) has been associated with increased crop productivity, as well as the mitigation of yield variability across time and space. Observations at the regional scale have indicated decreases in yield variability with increasing SOM. However, the mechanisms by which this variability is reduced remain poorly understood, especially at the farm scale. To better understand the relationship between SOM and yield heterogeneity, we examined its distribution between particulate organic matter (POM) and mineral-associated organic matter (MAOM) at the subfield scale within nine farms located in the central United States. We expected that the highest SOM concentrations would be found in stable, high-yielding zones and that the SOM pool in these areas would have a higher proportion of POM relative to other areas in the field. In contrast to our predictions, we found that unstable yield areas had significantly higher SOM than stable yield areas and that there was no significant difference in the relative contribution of POM to total SOM across different yield stability zones. Our results further indicate that MAOM abundance was primarily explained by interactions between crop productivity and edaphic properties such as texture, which varied amongst stability zones. However, we were unable to link POM abundance to soil properties or cropping system characteristics. Instead, we posit that POM dynamics in these systems may be controlled by differences in decomposition patterns between stable and unstable yield zones. Our results show that, at the subfield scale, increasing SOM may not directly confer increased yield stability. Instead, in fields with high spatiotemporal yield heterogeneity, SOM stocks may be determined by interactive effects of topography, weather, and soil characteristics on crop productivity and SOM decomposition. These findings suggest that POM has the potential to be a useful indicator of yield stability, with higher POM stocks in unstable zones, and highlights the need to consider these factors during soil sampling campaigns, especially when attempting to quantify farm-scale soil C stocks.

Spatial Variability of the Methane Hydrate Stability Zone’s Upper Boundary Parameters in the Water Column of the Sea of Okhotsk

Spatial Variability of the Methane Hydrate Stability Zone’s Upper Boundary Parameters in the Water Column of the Sea of Okhotsk

Experimental Research on the Floating Amount of Shield Tunnel Based on the Innovative Cumulative Floating Amount Calculation Method

The study of shield tunnel segment flotation is crucial for controlling the precision of underground excavation projects. Based on Winkler’s beam foundation theory, the load structure method, and the equivalent continuous beam model, and by considering the mechanical and spatial conditions that cause segment flotation, a novel theoretical calculation method for cumulative flotation is proposed using a simplified equivalent stiffness model of the tunnel. Additionally, a new concept of “equivalent flotation force” is introduced. The rationality and applicability of this theoretical calculation method are verified by comparing it with on-site construction data from the Yuanjiang River Crossing Tunnel Project in Changde, Hunan Province, China. The experimental results demonstrate that the theoretical calculation closely approximates the surface deformation monitoring data of the tunnel alignment in the eastern section of the project, and their deformation patterns are similar. Near the starting shaft, there is significant settlement influenced by stratum loss due to smaller tunnel flotation, with greater settlement occurring in the upper part. However, at approximately 45 m into both sections, they enter a deformation stability zone showing significant correlation in longitudinal deformation. Through comparison and verification of on-site experiments and theoretical model analysis, we preliminarily elucidate the feasibility of this innovative cumulative flotation theoretical calculation method which provides an important theoretical basis for assessing segment flotation issues in subsequent tunnel shield construction evaluations.

Bottom simulating reflections across the northern Gulf of Mexico slope

Bottom simulating reflections (BSRs) are widely used indicators for natural gas hydrate in reflection seismic data. BSRs, however, are typically studied on a local scale at a single site and less frequently on a large, regional scale. Herein, we map and describe BSRs over a 68,000 km2 area across the northern Gulf of Mexico slope, using publicly released industry-acquired 3D seismic data. We map 74 BSR zones that cover a total area of 2060 km2. To understand the characteristics of these BSR zones, we classify them based on their BSR characteristics in the seismic data (continuous, discontinuous, and clustered) and the geological settings in which they appear (structural, stratigraphic, and venting). The most common type of BSRs are discontinuous BSRs; such BSRs were mapped over ∼60% of the total BSR area. Discontinuous BSRs occur most commonly in structural or stratigraphic settings. Discontinuous BSRs indicate that gas and hydrate are likely concentrated in discrete layers crossing the base of the hydrate stability zone (HSZ). Continuous BSRs occur in 27% of the total BSR area and they primarily occur in venting and stratigraphic settings. At vent sites, hydrate formation likely reduces sediment permeability and focuses fluid flow laterally underneath the base of the HSZ, likely resulting in a continuous BSR along the area of lateral gas flow. Clustered BSRs, chaotic sequences of high amplitude reflections, occur in only 13% of the total BSR area and are only found within structural settings. In addition, the majority of clustered BSR zones are close to paleochannels, suggesting that clustered BSRs may be more common where coarser-grained sediments exist.

Study on the Mechanism of Natural Gas Hydrate Decomposition and Seabed Seepage Triggered by Mass Transport Deposits

Previous studies indicate that mass transport deposits are related to the dynamic accumulation of natural gas hydrates and gas leakage. This research aims to elucidate the causal mechanism of seabed seepage in the western region of the southeastern Qiongdongnan Basin through the application of seismic interpretation and attribute fusion techniques. The mass transport deposits, bottom simulating reflector, submarine mounds, and other phenomena were identified through seismic interpretation techniques. Faults and fractures were identified by utilizing variance attribute analysis. Gas chimneys were identified using instantaneous frequency attribute analysis. Free gas and paleo-seepage points were identified using sweetness attributes, enabling the analysis of fluid seepage pathways and the establishment of a seepage evolution model. Research has shown that in areas where the mass transport deposits develop thicker layers, there is a greater uplift of the bottom boundary of the gas hydrate stability zone, which can significantly alter the seafloor topography. Conversely, the opposite is true. The research indicates that the upward migration of the gas hydrate stability zone, induced by the mass transport deposits in the study area, can result in the rapid decomposition of gas hydrates. The gas generated from the decomposition of gas hydrates is identified as the principal factor responsible for inducing seabed seepage. Moderate- and low-speed natural gas seepage can create spiny seamounts and domed seamounts, respectively.

Methane seeps on the U.S. Atlantic margin: An updated inventory and interpretative framework

Since the discovery of >570 methane flares on the northern U.S. Atlantic margin between Cape Hatteras and Georges Bank in the last decade, the acquisition of thousands of kilometers of additional water column imaging data has provided greater coverage at water depths between the outer continental shelf and the lower continental slope. The additional high-resolution data reveal >1400 gas flares, but the removal of probable duplicates from the combined database of new flares and those recognized in 2014 yields ∼1139 unique sites. Most of these sites occur in clusters of 5 or more seeps, leaving about 275 unique locations (including 47 clusters) for seepage along the margin. As a function of depth, seep distribution is heavily skewed toward the upper continental slope at water depths shallower than 400 m on the southern New England margin and ∼ 550 m in the Mid-Atlantic Bight, with additional seeps clustered at ∼1100 m and just deeper than ∼1400 m in both sectors. Despite little ongoing tectonic deformation or active faulting on this passive margin, a variety of processes driven from below the seafloor (e.g., migration of fluids along faults or through permeable strata, seepage above diapirs or other pre-existing structures) and from above (e.g., erosion, sapping, unroofing) contribute to the development of seeps in different settings along the margin. In addition, the prevalence of seeps on promontories overlooking shelf-breaking canyons may be directly related to the three-dimensional nature of the hydrate stability zone in these locations. As a function of depth, the parts of the slope at the contemporary landward limit of gas hydrate stability are devoid of seeps, and the upper slope zones with the most concentrated seepage were not within the gas hydrate stability zone even during the Last Glacial Maximum. Thus, if the large number of upper slope seeps is at least partially sourced in gas hydrate degradation, the gas emitted at these seeps must have migrated there from greater depths on the continental slope.

Using the LabVIEW Simulation Program to Design and Determine the Characteristics of Amplifiers

Because of the large number of parameters that interact in amplifier functions, determining dynamic regime parameters as well as the mode of function of amplifier stages is an extremely complex problem. This paper describes a LabVIEW application for studying the functioning of an amplifier in various connections. The user selects the generator’s parameters, the type of connection and its parameters, as well as the load circuit characteristics. The application can determine both the stage characteristics and the Bode characteristics. The amplifier’s stability zone, as well as its gain and phase, are determined based on these characteristics. An important advantage of this application is that the design of the amplifier stage can be created starting from some parameters that the amplifier can establish, from which the values of components can be determined. In order to validate the simulation results from the LabVIEW application, the specialized program Multisim was used, as well as experimental measurements using the Electronics Explorer Board. Both Multisim and Electronics Explorer Board can determine Bode characteristics. In both simulations and experimental amplifiers, the same schemes with the same transistor were used. The application can be used for educational purposes as well as to design an amplifier’s stage to achieve specific parameters.

An experimental study of the flow structure in arrays of cold bidirectional swirling jets

Jets arrays have become a promising technology for industrial applications, including heat transfer and combustion processes. This paper presents the first experimental results of arrays of bidirectional swirling jets, including their flow structures and specific turbulent interactions. Particle image velocimetry of non-reacting linear arrays showed the formation of recirculation zones within each jet of the array. These regions are the most stable, with fluctuations less than 0.002 of the bulk inlet velocity Vin. In addition, jets merging locations also show a significant decrease in velocity fluctuations, with values V′/Vin &amp;lt; 0.005, which is most clearly seen at a smaller pitch distance. However, the highest turbulent fluctuations occur in shear layers, reaching values of V'/Vin ≈ 0.02 for linear arrays and V'/Vin ≈ 0.03 for planar configurations. The flow structure of planar arrays is more complex, leading to the formation of secondary vortex structures at locations of jet interaction. This results in a significant decrease in the flow swirl number, from an initial value of Sin = 2.1 in the vortex chamber, to S1 = 0.36–0.49 at z/dout2 = 1 and S2 = 0.19–0.27 at z/dout2 = 2 beyond the outlet nozzle. Strong cross-sectional motions are observed in planar arrangements, leading to the development of additional regions of negative axial velocity between the outlet nozzles. The boundaries of these regions have nearly zero axial velocity values and can provide reliable ignition and flame stabilization in the case of reactive bidirectional jets. Therefore, the flow structure in planar arrays forms a large stabilization zone, and the obtained results can be used to develop a new multipoint combustion technique.

INFLUENCE OF PLEISTOCENE GLACIATION ON PETROLEUM SYSTEMS AND GAS HYDRATE STABILITY IN THE OLGA BASIN REGION, BARENTS SEA

This study presents the results of a 2D numerical basin and petroleum systems model of the Olga Basin in the NW Barents Sea offshore northern Norway, a frontier exploration area in which there are abundant seafloor oil and gas seepages. The effects of Pleistocene ice sheet advances on rock properties and subsurface fluid migration in this area, and on seafloor hydrocarbon seepage, are not well understood. The 2D numerical model takes account of recurrent ice advances and retreats, together with related erosional and temperature effects, and investigates the influence of these parameters on fluid migration. Model results show that Pleistocene glaciations reduced the temperature in the sedimentary succession in the Olga Basin by up to 20 °C, for example in the uppermost Cretaceous and Jurassic sediments which underlie the seafloor down to a depth of 0.5 to 1 km. The decrease in temperature was in general predominantly related to the intensity of glacial erosion, which was set in this study to a depth of 600 m based on previous studies. Hydrocarbon fluids expelled from potential thermogenic source rocks of Carboniferous to Triassic ages on the SW margin of the Olga Basin migrated to the seafloor through permeable carrier beds. However, fluid migration to the surface in the NE of the study area took place along fault conduits. In a closed fault model scenario, only 0.3 Mt of hydrocarbons are modelled to have migrated along the 0.5 km wide model section; in a second scenario with partially open faults, about 22 Mt of hydrocarbons, representing about 11% of the total hydrocarbons generated by potential thermogenic source rocks in the study area, were lost to the surface during the Pleistocene. The potential for microbial methane generation in the Olga Basin was limited both during the Pleistocene and at the present day due to the significant reduction in temperature during glacial episodes, and due to the intense glacial‐related erosion of the Mesozoic to Cenozoic stratigraphy. During glacial stages, the gas hydrate stability zone beneath the ice sheet is modelled to have extended to a depth of up to 900 m for a pure methane composition, and to a depth of up to 1100 m for a possible thermogenic‐sourced mixed gas composition of 90% methane, 7% propane and 3% ethane. Gas hydrates with this mixed composition are modelled to have been stable in the Olga Basin during the last three glacial advances and into the present. These modelling results provide an insight into the key factors controlling the migration and surface leakage of hydrocarbon fluids in the Olga Basin region, and into the effects of glaciations on rock properties in a glaciated basin.