Conductivity Contrast Research Articles

Integral equation based wellbore preconditioner for 3D electromagnetic response modeling

Modeling subsurface controlled-source electromagnetic (CSEM) responses using the finite-element (FE) method is challenging in the presence of highly conductive wellbore casing. The very large conductivity contrast between the casing and the host formation leads to increased computation time and potentially unstable solutions. We address this difficulty by preconditioning an FE solver with an integral equation (IE) primary solution that captures the CSEM response of a realistic-sized steel wellbore casing. Our hybrid IE-FE approach determines the primary field using 2D integral-equation forward modeling and then interpolates the IE-computed solution onto the nodes. Then using an existing FE simulator, we solve for the secondary electric and magnetic fields. This approach removes the need for an ultra-fine FE mesh around the wellbore, thereby improving FE solution stability while greatly reducing FE computation time. Our method is illustrated by modeling the CSEM responses of idealized fluid-bearing zones.

Contrast of Hydraulic Conductivity Induces Transport of Combined Pollutants in High- and Low-permeability Systems

The transport process of pollutants in the environment can be influenced by heterogeneous geologic architecture and pollutant interactions. However, there has been a lack of research on co-transport behaviors of combined pollutants in heterogeneous aquifers. In this study, a series of two-dimensional tank experiments were carried out to study the transport behavior of toluene and naphthalene in both homogeneous and heterogeneous aquifers. The results revealed that the coexisting solutes facilitated the transport of toluene and naphthalene in the homogeneous aquifers, potentially due to competitive adsorption between these compounds. In the high- and low-permeability systems, the transport rates for both toluene and naphthalene decreased while exhibiting characteristics such as early arrival, long tails, and multiple peaks. The spatial analysis of pollutant distribution indicated that hydraulic conductivity contrast played a critical role in inducing back diffusion phenomenon. Furthermore, toluene exhibited more pronounced matrix diffusion compared to naphthalene in heterogeneous aquifers, characterized by higher concentrations, wider diffusion range in low-permeability zones. And the β value for toluene is smaller than naphthalene in CTRW model, indicating that the former is more sensitive to the hydraulic conductivity contrast. This study provides novel insights into understanding the co-transport behavior of combined pollutants in heterogenous aquifers.

Numerical simulation on transient electromagnetic response of separation layer water in coal seam roof

Mining stress induces deformation and fracture of the overlaying rock, which will result in water filling the separation layer if the aquifer finds access to abscission space along the fracture channels. Accurate detection is crucial to prevent water hazards induced by water-bearing fractures. The 3-D time-domain finite-difference method with Yee’s grid was adopted to calculate full-space transient electromagnetic response; meanwhile, a typical geologic and geophysical model with a water-bearing block in an separation layer was built according to regional tectonics and stratigraphic developments. By using numerical simulation, the induced voltage and apparent resistivity for both vertical and horizontal components were acquired, and then an approximate inversion was carried out based on the “smoke ring” theory. The results indicate that the diffusion velocity of induced voltage is significantly affected by the water-bearing body in the fracture, and the horizontal velocity of induced voltage is lower than the vertical one. The induced voltage curves indicate that the horizontal response to an anomaly body is stronger than the vertical one, leading to a high apparent resistivity resolution of conductivity contrast and separation layer boundary in the horizontal direction. The results of 3-D simulation making use of a measured data model also demonstrate that the horizontal component of apparent resistivity can reflect the electrical structure in a better way; however, its ability to recognize the concealed and fine conductor is rather weak. Accordingly, the observation method or numerical interpolation method needs to be further improved for data processing and interpretation.

WELL-TO-WELL (W2W) ELECTROMAGNETIC TOMOGRAPHY MODELING ADVANCEMENT: IMPROVING PRECISION AND EFFECTIVENESS WITH REGULARIZATION

This research describes the direct and inverse problems of cross-well electromagnetic tomography. The model geometry has azimuthal symmetry, which simplifies the forward modeling and the inversion procedure. In the direct problem, the finite element method is used in the numerical solution of the Helmholtz equation. In the inverse problem, the study discusses the use of three stabilizer functionals: Global Smoothness (GS), Total Variation (TV), and Absolute Equality (AE). The first uses a smoothing function in the L2 norm, while the latter uses smoothing in the L1 norm, for it accepts abrupt changes between adjacent parameters. The results show that the TV method generated good estimates of both geometry and conductivity of the bodies, both for small and large conductivity contrast between the targets and the surrounding environment. Through the results, one can also observe that the regularization of the Total Variation presented a better estimate of the parameters than the Global Smoothness. In most of the synthetic models used in this work, the best estimates of the proposed model occurred when Absolute Equality constraints were used on the cells at the edges of the inversion grid, in addition to the stabilizer functional.

Room temperature thermal rectification in suspended asymmetric graphene ribbon

Thermal rectifiers are essential in optimizing heat dissipation in solid-state devices to enhance energy efficiency, reliability, and overall performance. In this study, we experimentally investigate the thermal rectification phenomenon in suspended asymmetric graphene ribbons (GRs). The asymmetry within the graphene is introduced by incorporating periodic parallel nanoribbons on one side of the GR while maintaining the other side in a pristine form. Our findings reveal a substantial thermal rectification effect in these asymmetric graphene devices, reaching up to 45% at room temperature and increasing further at lower environmental temperatures. This effect is attributed to a significant thermal conductivity contrast between pristine graphene and nanoribbon graphene within the asymmetric structure. We observe that the incorporation of nanoribbons leads to a notable reduction in thermal conductivity, primarily due to phonon scattering and bottleneck effects near the nanoribbon edges. These findings suggest that graphene structures exhibiting asymmetry, facilitated by parallel nanoribbons, hold promise for effective heat management at the nanoscale level and the development of practical phononic devices.

Synergistic effects of boron and saponin in mitigating salinity stress to enhance sweet potato growth

Salinity stress significantly hinders plant growth by disrupting osmotic balance and inhibiting nutrient uptake, leading to reduced biomass and stunted development. Using saponin (SAP) and boron (B) can effectively overcome this issue. Boron decreases salinity stress by stabilizing cell walls and membranes, regulating ion balance, activating antioxidant enzymes, and enhancing water uptake. SAP are bioactive compounds that have the potential to alleviate salinity stress by improving nutrient uptake, modulating plant hormone levels, promoting root growth, and stimulating antioxidant activity. That’s why the current study was planned to use a combination of SAP and boron as amendments to mitigate salinity stress in sweet potatoes. Four levels of SAP (0%, 0.1%, 0.15%, and 0.20%) and B (control, 5, 10, and 20 mg/L B) were applied in 4 replications following a completely randomized design. Results illustrated that 0.15% SAP with 20 mg/L B caused significant enhancement in sweet potato vine length (13.12%), vine weight (12.86%), root weight (8.31%), over control under salinity stress. A significant improvement in sweet potato chlorophyll a (9.84%), chlorophyll b (20.20%), total chlorophyll (13.94%), photosynthetic rate (17.69%), transpiration rate (16.03%), and stomatal conductance (17.59%) contrast to control under salinity stress prove the effectiveness of 0.15% SAP + 20 mg/L B treatment. In conclusion, 0.15% SAP + 20 mg/L B is recommended to mitigate salinity stress in sweet potatoes.

Tidal Composition Analysis of Global Sq Current System

AbstractIn this study, we examine the spatiotemporal variability of the global equivalent ionospheric current system estimated from geomagnetic Solar‐quiet (Sq) variations at 124 mid‐ to low‐latitude ground‐based magnetometer stations during 1–31 May 2020. In this period, geomagnetic activity was particularly low, that is, the hourly geomagnetic activity index Hp60 did not exceed 4o. The spherical harmonic analysis is performed on the Sq variations to estimate the equivalent current function at each universal time hour. Hourly maps of the global Sq current system are used to evaluate, for the first time, the tidal composition of mid‐latitude Sq currents and its temporal evolution. Although tides are known to be an important driver of Sq currents, the tidal composition analysis was difficult in previous studies that focused on Sq currents at particular longitude sectors. Our results show that the migrating tidal components (DW1, SW2, and TW3) are predominant in both hemispheres. This is as expected from the strong day‐night contrast in ionospheric conductivities. The day‐to‐day variations of the migrating tidal components are, however, not strongly correlated between the two hemispheres, suggesting that these variations arise not only from the global effect of solar radiation on ionospheric conductivities but also from the local effect of neutral winds. Variations associated with non‐migrating tides are also found. Especially, eastward‐propagating diurnal and semidiurnal tides with zonal wavenumber 1 (DE1 and SE1) are the largest non‐migrating components. Their production mechanisms remain to be understood.

Theoretical study on contact current density in the boundary of random rough surfaces

We theoretically calculate here how random roughness affects the contact current density and equivalently contact resistance by utilizing a perturbative expansion method. Then, By taking into consideration self-affine roughness on the surface of the solids, we find the contact current density. The obtained results illustrate the total contact current exhibits significant dependence on the properties of the self-affine rough surface such as root-mean-square of roughness and Hurst exponent. Assuming a lower value of root-mean-square roughness for a solid in a higher potential compared to another which is connected to a lower potential, results in a decrement of the total contact current compared to the current between solids without roughness. By considering similar rough parameters for both solids, the contribution of the perturbed contact current becomes zero. Moreover, repercussions by considering the conductivity contrast between components are explored.

Modeling and Analysis of a Radiative Thermal Memristor

This study presents a theoretical framework for a radiative thermal memristor (RTM), utilizing Tungsten-doped vanadium dioxide (WVO) as the phase-change material (PCM) and silicon carbide (SiC) in the far-field regime. The behavior of the RTM is depicted through a Lissajous curve, illustrating the relationship between net flux (Q) and a periodically modulated temperature difference ΔT(t). It is established that temperature variations in the memristance (M) of the RTM form a closed loop, governed by PCM hysteresis. The analysis explores the impact of thermal conductivity contrast (r) and periodic thermal input amplitude (θ) on the Q–ΔT curve and the M–ΔT curve and negative differential thermal resistance (NDTR), revealing notable effects on the curve shapes and the emergence of NDTR. An increasing r leads to changes in the Lissajous curve’s shape and enhances the NDTR influence, while variations in both r and (θ) significantly affect the Q values and Lissajous curve amplitudes. In the M–ΔT curve, the height is linked to thermal conductivity contrast (r), with increasing r resulting in higher curve heights.

Exploration of Thermal Bridging Through Shrub Branches in Alpine Snow

AbstractIn the high Arctic, thermal bridging through frozen shrub branches has been demonstrated to cool the ground by up to 4°C during cold spells, affecting snow metamorphism and soil carbon and nutrients. In alpine conditions, the thermal conductivity contrast between shrub branches and snow is much less than in the Arctic, so that the importance of thermal bridging is uncertain. We explore this effect by monitoring ground temperature and liquid water content under green alders and under nearby alpine tundra in the Alps. During a January 2022 cold spell, the ground temperature at 5 cm depth under alders is 1.3°C colder than under alpine tundra. Ground water freezing under alders is complete, while water remains liquid under tundra. Finite element simulations reproduce the observed temperature difference between both sites, showing that thermal bridging does affect ground temperature also under Alpine conditions.

Thermal conductivity switching for a Y–Mg alloy hydride thin film due to hydrogenation/dehydrogenation reactions using dilute H2 gas

Thermal switching requires a significant contrast in thermal conductivity between the on and off states. We focus on thermal conductivity switching performance and mechanism for switchable mirror materials, which changes reversible metallic and semiconductor states due to hydrogenation and dehydrogenation. A thin film of yttrium–magnesium (Y–Mg) alloy hydride covered with a Pd catalyst layer was fabricated on quartz glass substrates by dc magnetron sputtering using a 60 at. % Y and 40 at. % Mg alloy target and a mixture of 50% Ar and 50% H2 gases. The crystal structure, electrical conductivity, and thermal conductivity in each state were measured using in situ x-ray diffraction analysis, Hall effect measurement, and thermoreflectance apparatus, respectively. The Y–Mg alloy hydride film was hydrogenated and dehydrogenated on exposure to a mixture of 3% H2 in N2 gas and air, respectively. The structural change in Y hydrides due to hydrogenation and dehydrogenation was clarified, whereas Mg or Mg hydride in the film showed no apparent crystallization. The thermal conductivity of the on-state was 4.5 times larger than that of the off-state. The thermal conductivity change from hydrogenated to dehydrogenated state was ∼5.4 W m−1 K−1, and ∼2.5 W m−1 K−1 of thermal conductivity change could be attributed to electron contribution based on the estimation using Wiedemann–Franz law. The thermal conductivity changes of Y–Mg alloy hydrides due to hydrogenation/dehydrogenation resulted from both electrons and phonons.

Electrical resistivity tomography through reinforced concrete floor

AbstractThe electrical resistivity tomography (ERT) method is often challenged by the presence of reinforced concrete (RC) in urban and industrial environments, because the embedded metallic wire mesh can severely distort the distribution of subsurface currents. We investigate one typical scenario in real applications, in which an RC floor overlays the natural topsoil or rock. Our synthetic forward simulations show that the embedded wire mesh behaves like a local good conductor in data of small source‐receiver separations and acts like an equal‐potential object that keeps the potential from decaying at large source‐receiver separations. Routine ERT inversions that ignore the RC cannot work properly because the thin and highly conductive wire mesh may be manifested as large uninterpretable low‐resistivity anomalies in the imaging results. Two remedies are adopted to improve the ERT resolution in such cases. First, we find a top layer with high conductivity in our model to adequately represent the wire mesh; then, we initiate the inversion with the top‐layer model as the starting and reference model. This warm‐start approach overcomes the difficulty of recovering the large conductivity contrast between metallic objects and regular earth materials. Second, underground electrodes are added to the survey array, so more information from depth can be obtained to fight against the dominance of current channelling in the wire mesh. Finally, our strategies are used to invert a real ERT dataset from an indoor manufacturing plant, where RC covers the entire floor of the building and electrodes are in contact with the soil through open holes in the floor. Our simulation and field data inversion verify our findings and demonstrate the effectiveness of our solutions in improving the resolution of ERT when the survey is carried out over RC floor in urban and industrial environments.

Limits of three‐dimensional target detectability of logging while drilling deep‐sensing electromagnetic measurements from numerical modelling

AbstractSubsurface energy resources are often found in three‐dimensional and non‐spatially continuous rock formations that exhibit electrical anisotropy. Deep‐sensing tri‐axial borehole electromagnetic measurements are currently being used to detect three‐dimensional fluid‐bearing subsurface formations, but borehole environmental, geometrical and instrument‐design factors, together with measurement noise, constrain their practical range of detection and spatial resolution. By understanding the interplay of the above factors on the detectability and sensitivity of borehole deep electromagnetic measurements, one can potentially quantify the uncertainty of both target spatial location (relative to the well trajectory) and electrical resistivity contrast, thereby improving the certainty of three‐dimensional well navigation in real time. We implement a finite‐volume method to numerically solve Maxwell's equations for three‐dimensional electrically anisotropic heterogeneous rock formations in the calculation of magnetic fields measured with deep‐sensing tri‐axial borehole electromagnetic instruments. The calculated magnetic fields at measurement locations are described as the per cent difference between measurements acquired in rock formations with and without conductive or resistive three‐dimensional targets. We quantify the maximum radial distance of detection from the well trajectory and the spatial sensitivity of a commercially available deep‐sensing electromagnetic instrument with respect to environmental factors and measurement acquisition parameters by assuming that the borehole electromagnetic instrument can reliably detect offset three‐dimensional targets if the corresponding per cent measurement difference exceeds the threshold for measurement noise. Results indicate that commercially available tri‐axial deep‐sensing borehole electromagnetic instruments can achieve maximum detection distances between a quarter of and a full transmitter–receiver spacing. In addition, we show that radial detection distances vary from 0.3 to 2 skin depths depending on the geological environment, that is, target conductivity contrast with respect to the embedding background, electrical anisotropy of the background formation and targets and measurement acquisition parameters, that is, frequency and transmitter–receiver spacings. The above findings are not only important for instrument design and measurement‐acquisition planning but also for the effective implementation of real‐time inversion‐based measurement interpretation procedures.

Flexible germanium monotelluride phase change films with ultra-high bending stability for wearable piezoresistive sensors

Although non-stoichiometric chalcogenide phase change materials (e.g., Ge2Sb2Te5, Sb2Te3) have been used in wearable electronics, their non-stoichiometric vacancies will inevitably degrade their bending stability and electrical contrast. In this study, we demonstrated that stoichiometric germanium monotelluride films deposited on flexible polyimide substrates exhibited excellent phase transition characteristics and ultra-high bending stability under different bending states. The GeTe films showed improved bending stability under compressive strain than under tensile strain, and the GeTe films still exhibited high conductivity contrast (6 orders of magnitude) after 2000 bending cycles, with potential use in flexible memory devices. As a proof-of-concept, the bending stability of the stoichiometric GeTe films was shown to be significantly better than the non-stoichiometric Ge2Sb2Te5 films. Notably, we fabricated a flexible piezoresistive sensor based on the GeTe film, which exhibited good sensitivity (gauge factor is 38), low detection limit (0.1%), short response time (83.3 ms), good reliability and repeatability (more than 1200 cycles) and high heat resistance (200 ℃), and was used to track human pulse, as well as finger and face movements. The results showed that the piezoresistive sensor had good sensitivity and a stable time-resolved electrical response, and could be used to detect the movement of different human body parts and the small signals of different pressure points. These findings indicated that GeTe-based wearable piezoresistive sensors hold significant promise for medical health monitoring and biosensing applications.

Correlation of interfacial and dielectric characteristics in atomic layer deposited Al2O3/TiO2 nanolaminates grown with different precursor purge times

Considering the potential applications of Al2O3/TiO2 nanolaminates (ATA NLs) in storage capacitors, device-grade ATA NLs are fabricated using an ALD system, wherein the effect of precursor purging time on interfacial, and dielectric properties is thoroughly investigated. With an increase in half-cycle purging time from 2 to 4 s, the observed improvement in interface quality and sublayer density of these NLs is ascribed to the efficient removal of reaction by-products and impurities. Moreover, with an increase in purge time from 2 to 4 s, the increase in dielectric constant and concurrent decrease in dielectric loss from ∼132 to 154 and from ∼0.29 to 0.2, respectively, are primarily assigned to the improvement in sublayer conductivity contrast assisted Maxwell–Wagner interfacial polarization across Al2O3/TiO2 interfaces. The NL based devices fabricated at 4 s purging time, exhibited a capacitance density of ∼18.94 fF/μm2, low equivalent oxide thickness of ∼1.82 nm, and reduced leakage current density of ∼3.04 × 10−5 A/cm2 at 2 V applied bias, which demonstrates its suitability as high-k materials for energy storage applications. Furthermore, this study not only gives an insight of the purging time induced growth chemistry of ATA NLs but also explores the possibility of improving its dielectric performance essential for multifaceted applications.

MilliKelvin microwave impedance microscopy in a dry dilution refrigerator.

Microwave impedance microscopy (MIM) is a near-field imaging technique that has been used to visualize the local conductivity of materials with nanoscale resolution across the GHz regime. In recent years, MIM has shown great promise for the investigation of topological states of matter, correlated electronic states, and emergent phenomena in quantum materials. To explore these low-energy phenomena, many of which are only detectable in the milliKelvin regime, we have developed a novel low-temperature MIM incorporated into a dilution refrigerator. This setup, which consists of a tuning-fork-based atomic force microscope with microwave reflectometry capabilities, is capable of reaching temperatures down to 70 mK during imaging and magnetic fields up to 9T. To test the performance of this microscope, we demonstrate microwave imaging of the conductivity contrast between graphite and silicon dioxide at cryogenic temperatures and discuss the resolution and noise observed in these results. We extend this methodology to visualize edge conduction in Dirac semi-metal cadmium arsenide in the quantum Hall regime.

Application of audio-frequency magnetotelluric data to cover characterisation – validation against borehole petrophysics in the East Tennant region, Northern Australia

The characterisation of the thickness and geology of cover sequences significantly improves targeting for mineral exploration in buried terrains. Audio-frequency Magnetotelluric (AMT) data is applicable to characterise cover sequences, where their conductivity (inverse resistivity) can be differentiated. We present a regional study from the under-cover East Tennant region in the Northern Territory (Australia) where we have applied deterministic and probabilistic inversion methods to derive 2D and 1D resistivity models. We integrated these models with information of co-located basement penetrating boreholes (lithological and geophysical logs) to ground-truth and validate the models and to improve geophysical interpretations. In the East Tennant region, borehole lithology and wireline logging demonstrate that the modelled AMT response is largely controlled by the mineralogy of the cover and basement rocks. The bulk conductivity is due primarily to bulk mineralogy and the success of using the AMT models to predict cover thickness is shown to be dependent on whether the bulk mineralogy of cover and basement rocks are sufficiently different to provide a detectable conductivity contrast. Our investigation of a range of geological scenarios that differ in thickness, complexity and geology of the cover and basement rocks suggests that in areas where there is sufficient difference in bulk mineralogy and where the stratigraphy is relatively simple, AMT models predict the cover thickness with high certainty. In more complex scenarios interpretation of AMT models may be more ambiguous and requires integration with other data (e.g. drilling, wireline logging, potential field modelling). Overall, we conclude that the application of the method has been validated and the results compare favourably with borehole stratigraphy logs once geological (i.e. bulk mineralogical) complexity is understood. This demonstrates that the method is capable of identifying major litho-stratigraphic units with resistivity contrasts. Our results have assisted with the planning of regional drilling programs and have helped to reduce the uncertainty and risk associated with intersecting targeted stratigraphic units in covered terrains.

Split ring resonator topology based microwave induced thermoacoustic imaging (SRR-MTAI).

Microwave-induced thermoacoustic imaging (MTAI) using low-energy and long-wavelength microwave photons has great potential in detecting deep-seated diseases due to its unique ability of visualizing intrinsic electric properties of tissue in high resolution. However, the low contrast in conductivity between a target (e.g., a tumor) and the surroundings sets a fundamental limit for achieving a high imaging sensitivity, which significantly hinders its biomedical applications. To overcome this limit, we develop a split ring resonator (SRR) topology based MTAI (SRR-MTAI) approach to achieve highly sensitive detection by precise manipulation and efficient delivery of microwave energy. The in vitro experiments show that SRR-MTAI demonstrates an ultrahigh sensitivity of distinguishing a 0.4% difference in saline concentrations and a 2.5-fold enhancement of detecting a tissue target which mimicks a tumor embedded at a depth of 2 cm. The in vivo animal experiments conducted indicate that the imaging sensitivity between a tumor and the surrounding tissue is increased by 3.3-fold using SRR-MTAI. The dramatic enhancement in imaging sensitivity suggests that SRR-MTAI has the potential to open new avenues for MTAI to tackle a variety of biomedical problems that were impossible previously.

The emergent behavior of three-dimensional and three-phase thermal network model for the heterogeneous materials

The study of the heat transfer properties of heterogeneous materials has been a world-wide research focus, for example in sectors related to thermal management in aerospace, architecture, and geology. In this paper, the emergent behavior and thermal conduction characteristics of three-dimensional and three-phase heterogeneous materials thermal networks are studied using the finite element method. The results show that the existence of percolation paths of each phase has a strong impact on the effective properties of the network when the contrast in thermal conductivities of each phase is high, and percolation also affects the effective thermal conductivity of the whole thermal network system. However, when the contrast in thermal conductivity between the two phases is low, the thermal networks exhibit a more consistent and “emergent” behavior, and the effective thermal conductivity of thermal networks at the same volume fraction changes to a lower extent from network to network. This paper also demonstrates that a logarithmic mixing rule can predict the effective thermal conductivity in the low contrast emergent region in three-dimensional networks, and the modeling method provides new approaches for the design of multi-phase composites and prediction of their thermal conduction properties.

Surface entrapment of micromotors by a background temperature field

The fabrication of self-propelling micromotors and the study of their propulsion strategies have gained attention due to their wide range of applications in the medical, engineering, and environmental fields. The role of a background temperature field in the precise navigation of a self-thermophoretic micromotor near an insulated wall has been investigated by employing exact solutions to the energy equation and creeping flow. We report bound states for half-coated micromotors appearing as steady-state sliding, damped, and periodic oscillations when the dimensionless external temperature gradient (S) is in the range of 0.15≤S<0.26. The sliding height is lower with S but remains insensitive to the thermal conductivity contrast. Moreover, the stationary states for the self-propelled, asymmetrically coated micromotors transform into scattering trajectories. We highlight the combinations of S and coating coverage needed for guided swimming up or against the field along with a broad spectrum of counter-intuitive temporal variations of its navigating locations. These unique observations have been ascribed to a confinement-mediated dynamic coupling between the passive and active propulsion mechanisms.