3D Subsurface Research Articles

Use of Deep-Learning-Accelerated Gradient Approximation for Reservoir Geological Parameter Estimation

The estimation of space-varying geological parameters is often not computationally affordable for high-dimensional subsurface reservoir modeling systems. The adjoint method is generally regarded as an efficient approach for obtaining analytical gradient and, thus, proceeding with the gradient-based iteration algorithm; however, the infeasible memory requirement and computational demands strictly prohibit its generic implementation, especially for high-dimensional problems. The autoregressive neural network (aNN) model, as a nonlinear surrogate approximation, has gradually received increasing popularity due to significant reduction of computational cost, but one prominent limitation is that the generic application of aNN to large-scale reservoir models inevitably poses challenges in the training procedure, which remains unresolved. To address this issue, model-order reduction could be a promising strategy, which enables us to train the neural network in a very efficient manner. A very popular projection-based linear reduction method, i.e., propel orthogonal decomposition (POD), is adopted to achieve dimensionality reduction. This paper presents an architecture of a projection-based autoregressive neural network that efficiently derives an easy-to-use adjoint model by the use of an auto-differentiation module inside the popular deep learning frameworks. This hybrid neural network proxy, referred to as POD-aNN, is capable of speeding up derivation of reduced-order adjoint models. The performance of POD-aNN is validated through a synthetic 2D subsurface transport model. The use of POD-aNN significantly reduces the computation cost while the accuracy remains. In addition, our proposed POD-aNN can easily obtain multiple posterior realizations for uncertainty evaluation. The developed POD-aNN emulator is a data-driven approach for reduced-order modeling of nonlinear dynamic systems and, thus, should be a very efficient modeling tool to address many engineering applications related to intensive simulation-based optimization.

Quasi-3D Geoelectrical Imaging as A New Application for Landslide Investigations: A Tunnel Case Induced by Blasting Activity

Landslides are a significant hazard in mountainous regions, especially when influenced by construction activities such as tunnel excavation. In this paper, we aim to conduct a slope stability analysis as a result of tunnel blasts using quasi-3D subsurface models based on resistivity values. The study site is a construction area for the Jakarta-Bandung High-Speed Train tunnel, located in a mountainous region undergoing drill-and-blast excavation. This excavation method makes the area susceptible to landslides, which pose a threat to settlements in the Padalarang subdistrict, West Bandung Regency, Indonesia. Data was collected along four lines in 2D, and the dipole-dipole array was used to enhance resolution. Data modeling was carried out using ResIPy v3.2.3 software to create 2D and quasi-3D subsurface models based on resistivity values. The study findings indicate that the study area exhibits three resistivity ranges: low resistivity (0-30 Ωm), medium resistivity (31-49 Ωm), and high resistivity (>50 Ωm). Utilizing quasi-3D imaging, we were able to identify the dimensions and presence of slip surfaces, which can be categorized as shallow (1.5-5 m) and deep (5-20 m) criteria. This study successfully applied the quasi-3D geoelectrical approach in a susceptible environment to detect potential landslide zones.

Influence of 3D subsurface flow on slope stability for unsaturated soils

The occurrence of rainfall-induced slope failures has become more frequent due to the effect of climate change. Hence, various studies have been conducted to analyse the effect of rainfall infiltration on slope stability. Physically-based hydrological models have been commonly used with slope stability models such as the infinite slope model to develop slope susceptibility maps. However, a combination of three-dimensional (3D) water balance model with 3D limit equilibrium method (LEM) has not been commonly used. Hence, in this study, a water balance model, GEOtop was used to investigate the influence of subsurface flow in unsaturated soil under extreme rainfall conditions on regional slope stability in 3D directions. The results from the GEOtop model were used as inputs for 3D LEM slope stability analysis performed using the Scoops3D software to obtain the factor of safety (FOS) map for the region. Four slopes within the region were then selected to be modelled in the two-dimensional (2D) seepage and slope stability analyses, SEEP/W and SLOPE/W. Results from the detailed study showed that the pore-water pressures (PWPs) from the 3D water balance analyses were found to be higher than the 2D seepage analyses. Under similar PWP conditions, the FOS from the 2D slope stability analysis was observed to be lower than the 3D analysis for two out of the four slopes. However, the combined 3D water balance and slope stability analyses produced lower FOS compared to the 2D seepage and slope stability analyses due to the higher PWPs in the 3D water balance analyses. Therefore, this study highlights the importance of considering the 3D subsurface flow in unsaturated soil given that it has a significant influence on the FOS of slopes.

Improved Approaches for 3D Gravity and Gradient Imaging Based on Potential Field Separation: Application to the Magma Chamber in Wudalianchi Volcanic Field, Northeastern China

The gravity and gradient anomalies contain valuable information about the underground geological structures at various depths. Deep and shallow buried source bodies are able to be identified through multi-scale field separation processes, and visual comprehensions of geological structures can be obtained via 3D density inversion techniques. In this study, we propose an improved 3D imaging strategy based on gravitational field separation using the preferential continuation filter. This strategy incorporates the relationship between spectral features and buried depths of source bodies, allowing for a one-step transformation from planar gravity and full-tensor gradient field observations to a 3D density structure in the wave-number domain. Synthetic tests validate the effectiveness and robustness of the gravity and gradient imaging approaches, highlighting their advantages in high vertical resolution and low computational requirements. Nonetheless, it should be noted that the imaging effects of horizontal gradients Γxx and Γyy are unsatisfactory due to their weak noise resistance. Thus, they are not suitable for real data applications. The other imaging approaches are further applied to recover the subsurface 3D density structure beneath the Weishan cone in Wudalianchi Volcanic Field, Northeastern China. Our results provide insights into the possible location and shape of the low-density magma chamber. Also, the potential presence of partial melts is inferred and supported from a gravity perspective. The primary advantage of these approaches is their ability to generate a reasonable geological model in scenarios with limited prior information and physical property constraints. As a result, they have significant practical value in the field of applied geophysics, including mineral exploration and volcanology studies.

Characterization of subsurface geology and hydrogeology in Kribi -Cameroon using electrical resistivity soundings and 3D-Implicit modelling: Baseline for groundwater resource management

This scientific article aims to characterize the geology and hydrogeology of Kribi, Cameroon for the purposes of water resources management, infrastructure development, and environmental protection. The specific objectives include mapping the subsurface electrical resistivity structure, identifying groundwater potential zones, and constructing a hydrogeological model for groundwater modelling. The study employed the Schlumberger method to map the electrical resistivity structure of the 3D subsurface and utilized 3D implicit modelling techniques to delineate groundwater potential zones. The study revealed that Kribi exhibits a highly complex geology, characterized by three-layered electrical resistivity models, typical of sedimentary and metamorphic formations. The major lithologies identified include laterite, clay, sand, quartz, and granite. The region encompasses areas with low water potential and brackish zones alongside high groundwater production areas. The study established the presence of aquifers within Kribi, with the unconfined aquifer found at depths ranging from 2 m to 13 m, and the confined aquifer extending from 14 m to 40 m. The aquifers are notably thick and productive, indicating substantial groundwater resources. Static groundwater volumes were estimated at 50,179,381 m3 and 317,118,946 m3 for the unconfined and confined aquifers respectively. Despite this potential, water resources management in the area is non-existent. The study recommends the implementation of Integrated Water Resources Management (IWRM) to improve water resource management, with a primary focus on groundwater modelling. The northern and southern regions of Kribi exhibit high permeability, which is advantageous for groundwater recharge and storage projects but also poses challenges for environmental remediation efforts. Consequently, safeguarding the delineated area for groundwater provision is of utmost importance.

Influence of Subsurface Critical Zone Structure on Hydrological Partitioning in Mountainous Headwater Catchments

AbstractHeadwater catchments play a vital role in regional water supply and ecohydrology, and a quantitative understanding of the hydrological partitioning in these catchments is critically needed, particularly under a changing climate. Recent studies have highlighted the importance of subsurface critical zone (CZ) structure in modulating the partitioning of precipitation in mountainous catchments; however, few existing studies have explicitly taken into account the 3D subsurface CZ structure. In this study, we designed realistic synthetic catchment models based on seismic velocity‐estimated 3D subsurface CZ structures. Integrated hydrologic modeling is then used to study the effects of the shape of the weathered bedrock and the associated storage capacity on various hydrologic fluxes and storages in mountainous headwater catchments. Numerical results show that the weathered bedrock affects not only the magnitude but also the peak time of both streamflow and subsurface dynamic storage.

Integrating 3D subsurface imaging, seismic attributes, and wireline logging analyses: Implications for a high resolution detection of deep-rooted gas escape features, eastern offshore Nile Delta, Egypt

The eastern part of the offshore Nile Delta is one of Egypt's most productive gas provinces. This study utilizes well logs dataset and 2D seismic from Port Fouad marine (PFM) Field to investigate reservoir lithofacies, rock-types as well as the associated gas chimneys, pockmarks, bright anomalies of escaping fluids in the seafloor, and high-faulted zones caused by the chimney's developments. First, structural smoothing and median filters are used to improve the continuity of the stratigraphic horizons while preserving various geological structures, in order to pick all formation tops and structural faults. Afterwards, variance and energy attributes were computed and extracted for amplitude enhancement of various geobodies such as faults, gas chimneys, and bright reflectors. After that, petrophysical analyses revealed that the Upper Miocene Wakar Formation hosts scattered, thin gas pay zones (<10m) with low hydrocarbon saturation values (average Shc = 27.9%), high water saturation (average Sw > 70%), and high shale volume (average Vsh > 30%). These poor reservoir characteristics prevented accumulation of thick gas zones, thereby several gas escape features (e. vertical chimneys) are observed in the region. Finally, the estimated petrophysical parameter logs and the structure maps derived from seismic data are crucial components of the built 3D structural model. This model reveals the presence of four major normal faults (F1, F2, F3, and F4) oriented NW-SE and WNW-ESE, creating grabens and half grabens. These major faults (F1, F2, F3, and F4) are directly related to the inverted tectonics during the Cretaceous-Paleogene compression (Syrian Arc phase). A link between these structural patterns and gas escape zones is proposed, and the absence of good reservoir facies likely promoted an active phase of gas seepage in the study region. In addition, such gas escape chimneys can create catastrophic events on the seafloor and are considered direct evidence for the presence of gas generation and failure in accumulation.

Estimation of Geotechnical Parameters for Coal Exploration from Quasi-3D Electrical Resistivity Measurements

Geotechnical parameters are crucial for mine planning and operation at different stages of development. However, estimating these parameters requires a large number of boreholes and subsequent detailed analysis of the samples, making it a cumbersome exercise. Moreover, even after conducting these studies, it is not possible to cover the entire operational area. To address this issue, this study presents an indirect method of estimating geotechnical parameters through mathematical relations using resistivity data. The present study incorporated 2D and 3D subsurface imaging techniques for exploring coal reserves and analyzing geotechnical parameters that define subsurface soil properties. Electrical resistivity tomography (ERT) was utilized for data acquisition, employing a Dipole–dipole array with a multielectrode ABEM Terrameter LS instrument. Six parallel profiles were conducted, each 400 m in length, with an inter-electrode spacing of 10 m and a spacing of 50 m between profiles. These profiles were combined into a 3D dataset referred to as quasi-3D ERT. The inversion process for both 2D and 3D data was performed using the Res2dinv and Res3dinv programs, respectively. This study overcame the challenges of 2D resistivity sections by evaluating horizontal depth slices in the x-z plane from layers 1 to 10, reaching a depth of 81.2 m. The geotechnical parameters, including cohesion, friction angle, moisture content, and plastic index, were derived from the resistivity data. The ERT method proved to be cost-effective and efficient in determining soil properties over a large area compared with traditional laboratory analysis of borehole samples. Additionally, the variation of geotechnical parameters with resistivity values exhibited unique characteristics. The results from both the 2D and quasi-3D ERT were well correlated with the borehole data. Such studies are valuable for resource exploration and mine planning purposes.

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

Immersive visualization of 3D subsurface ground model developed from sparse boreholes using virtual reality (VR)

Analytics and visualization of multi-dimensional and complex geo-data, such as three-dimensional (3D) subsurface ground models, is critical for development of underground space and design and construction of underground structures (e.g., tunnels, dams, and slopes) in engineering practices. Although complicated 3D subsurface ground models now can be developed from site investigation data (e.g., boreholes) which is often sparse in practice, it remains a great challenge to visualize a 3D subsurface ground model with sophisticated stratigraphic variations by conventional two-dimensional (2D) geological cross-sections. Virtual reality (VR) technology, which has an attractive capability of constructing a virtual environment that links to the physical world, has been rapidly developed and applied to visualization in various disciplines recently. Leveraging on the rapid development of VR, this study proposes a framework for immersive visualization of 3D subsurface ground models in geo-applications using VR technology. The 3D subsurface model is first developed from limited borehole data in a data-driven manner. Then, a VR system is developed using related software and hardware devices currently available in the markets for immersive visualization and interaction with the developed 3D subsurface ground model. The results demonstrate that VR visualization of the 3D subsurface ground model in an immersive environment has great potential in revolutionizing the geo-practices from 2D cross-sections to a 3D immersive virtual environment in digital era, particularly for the emerging digital twins.

Evaluation of Mangroves Effectiveness in Strengthening Coastal Bund Using Geophysical Method (Non-Destructive Testing) at Tanjong Laboh, Batu Pahat, Johor

In recent years, mangrove forests have rapidly become an endangered ecosystem in the world. On the coastal shorelines, the saltwater intrusion (SI) phenomenon is a main concern that can affect the quality of water. The saltwater intrusion might occur due to sea level rise, and extreme natural disasters along with changes in mangrove ecosystem distribution. The purpose of this research is to evaluate the effectiveness of mangroves for strengthening coastal bunds in Tanjong Laboh, Batu Pahat, Johor. By using the non-destructive geophysical method, the presence of saltwater intrusion and subsurface profile of the survey line in the coastal area with and without mangroves was determined. The coastal bund conditions with and without mangroves were also compared in order to evaluate how the mangrove forest acts as coastal protection. The result shows that the resistivity values of 0.1 to 1.0 ohm are more dominant in the mangrove bund area than in the engineered coastal bund and rock outcrop areas. The saltwater intrusion integrity map was produced from the 2D subsurface image profiles for both survey lines. The comparison of the coastal bund conditions for two different scenarios, which with and without mangroves, has shown that the coastal bund without the mangroves experienced cracking and erosion compared to the coastal bund with the presence of mangroves. This indicates the importance of mangroves in protecting the coastal regions.

3D Visualization of Geothermal System Structure Based on Inversion Model of Gravity Data. Case Study: Mt. Salak Region, West Java

The geothermal power plant located in Gunung Salak plays a crucial role in increasing the electricity supply transmitted to the Java-Bali region, as the energy demand continues to rise. The objective of this research is to determine the 3D subsurface structure of Gunung Salak, specifically the distribution of the reservoir as the target for geothermal energy using the Gravity method. Gravity data, including gravity disturbance (gd), geoid, and Digital Elevation Model (DEM), were obtained from the ICGEM website with a total of 48740 data each. Based on the results of the residual anomaly map, the low anomalies beneath Gunung Salak have values from -5.15 to -1.88 mGal, which are suspected to be associated with the magma chamber. The high anomalies beneath the manifestations have values from 0.92 to 5.01 mGal, indicating andesitic basalt intrusive rocks believed to be the reservoir rock. Through the 3D inversion modeling of the subsurface structure of the Gunung Salak geothermal system, a clay cap with a density from 2.47 to 2.5 g/cc at depths of 0 to 700 m and andesitic basalt rock as the reservoir with a density from 2.74 to 2.91 g/cc at depths of 700 to 3000 m have been identified.

Subsurface geological and geophysical data from the Po Plain and the northern Adriatic Sea (north Italy)

Abstract. The Po Plain (Italy) is one of the most densely populated and productive regions of Europe, characterized by a flourishing economy (also linked to strategic subsurface resources) and several world cultural and natural heritage sites. The coupling of socio-economic interests with geological hazards (i.e. seismic, subsidence, and flooding hazards) in this area requires accurate knowledge of the subsurface geology, the active geological processes, and the impact of human activities on natural environments to mitigate the potential natural and anthropic risks. Most data unveiling the subsurface geology of this region were produced by the hydrocarbon exploration industry. Indeed, the Po Plain hosts many hydrocarbon fields that have been discovered since the early 1950s, giving rise to the subsurface exploration through extensive seismic reflection surveys and drilling of numerous deep wells. In this work, geological and geophysical data from 160 deep wells drilled for hydrocarbon exploration and/or exploitation purposes in the Po Plain and in the facing northern Adriatic Sea have been collected and digitized along with several published geological cross-sections and maps. These data have been used to reconstruct the overall subsurface 3D architecture and to extract the physical properties of the subsurface geological units. The digitized data are suitable to be imported into geo-software environments so as to derive the geophysical and mechanical properties of the geological units for a wealth of applied and scientific studies such as geomechanical, geophysical, and seismological studies. The integrated dataset may represent a useful tool in defining regional first-order strategies to ensure the safety of the urbanized areas and human activities and to reduce natural and anthropic risks that may affect this crucial region of Europe. In particular, the data collected would be useful to highlight sensible areas where data collection and more detailed studies are needed. Nowadays, such issues are particularly relevant for the underground industry development related to the increasing interest in possible CO2 and hydrogen underground storage, which can play a fundamental role in the energy transition process towards decarbonization goals. The full dataset is available at the following link: https://doi.org/10.5281/zenodo.8126519 (Livani et al., 2023).

Data-driven multi-stage sampling strategy for a three-dimensional geological domain using weighted centroidal voronoi tessellation and IC-XGBoost3D

Conventional site investigation schemes typically use an empirical sampling strategy, which involves equal sampling spacing with regular grid patterns. However, this approach assigns equal importance to the entire site and does not account for subsurface stratigraphic variations and site constraints. In this study, a data-driven multi-stage sampling strategy is proposed to adaptively optimize the borehole locations for a three-dimensional (3D) geotechnical site, considering the subsurface stratigraphic uncertainties and irregular site geometries. The initial sampling plan is determined based on weighted centroidal Voronoi tessellation, which assigns various sampling densities to zones depending on their importance. Measurements obtained in the previous sampling stage are combined with prior geological knowledge to establish and update a 3D geological domain using a physics-informed geological modelling method. The next optimal location for sampling is adaptively determined with the objective of maximizing the reduction in the stratigraphic uncertainty. The proposed method represents the first data-driven 3D site planning method that explicitly considers 3D subsurface stratigraphic variations. The performance of the proposed multi-stage sampling strategy is illustrated using a simulation study. The results indicate that the proposed method efficiently identifies the optimal sampling locations while comprehensively considering the 3D subsurface geological uncertainties and irregular site geometries.

Application of electrical resistivity tomography (ERT) for vertical high-density polyethylene (HDPE) membrane detection

At landfill sites, uncontrolled leachate is a serious cause of groundwater and soil contamination at landfills and can occur extensively when there are looploopholesin the vertical HDPE (High-density polyethylene) membrane. For this reason, detecting HDPE membranes is necessary. The conventional electrode running poles approach is limited, it reflects the distribution of resistivity in 2D subsurface profiles, which is difficult to deploy electrodes above vertical HDPE membrane. According to the structural characteristics and field conditions of vertical HDPE membranes, the two-row current and potential (CP) electrode array is proposed to be applied to the testing of vertical HDPE membranes. The forward geoelectric model of the vulnerability was built using pyGimli, and data on 1519 apparent resistivity were acquired. The previous knowledge without the vulnerability was taken into consideration as the constraint condition, and the regional control was added to conduct the inversion, all using the least square approach. The model was used to analyze data on the detection of vertical HDPE membranes in a landfill In Yancheng, and potential flaws were discovered. Results show that the CP electrical measuring device can be efficient, quick, and accurate in the positioning of vertical HDPE membrane hole detection, and it is a cost-effective method. To verify the results, in the vertical HDPE membrane, through the Wenner device, the inversion of the 2D profile, and low resistivity anomalies are detected. This is thought to be through the leachate of leakage of HDPE membrane resistance rate is about 1.8 ohm. The findings demonstrate that the vertical HDPE membrane susceptibility may be located and detected using a CP electrical measuring instrument, which is a practical and efficient approach.

A machine learning approach using legacy geophysical datasets to model Quaternary marine paleotopography

High-resolution subsurface marine mapping tools, including chirp and 3D seismic, enable the reconstruction of ancient landscapes that have been buried and subsequently submerged by marine transgression. However, the established methods for paleotopographic reconstruction require time consuming field and data interpretation efforts. Here we present a novel methodology using machine learning to estimate Marine Isotope Stage 2 (MIS2) paleotopography over a large (22 000 km2) area of the Northern Gulf of Mexico with meter-scale accuracy (2.7 m mean prediction error, 4.3 m 1-σ mean uncertainty). A relatively small area (3300 km2) of high-resolution (30 × 30 m) interpreted paleotopography is used as training and validation data, while modern bathymetry and MIS2 paleovalley location (binary deep/shallow paleotopography) are used as predictors. This approach merges the high-resolution of modern mapping techniques and the broad coverage of low-resolution legacy geophysical data. Machine learning-modeled paleotopography is not a substitute for precise high-resolution paleotopography reconstruction techniques, but it can be used to reasonably approximate paleotopography over large areas with greatly reduced expense and expertise.

Building Information Modelling for Application in Geotechnical Engineering

BIM (Building Information Modelling) is used to create and manage data during design, construction, and operation. It helps to effectively manage resources and optimize processes in the construction industry. Geotechnical engineering is one of the complex disciplines that may require BIM integration. Various data types must be provided in a timely manner and require real-time feedback, fast processing, and construction guidance. The first problem presented in the paper is the use of the traditional 2D-based method used by engineers for a particular task. It seems to be impractical when some adjustments are included. Another issue is the lack of communication between the workers. It poses the problem of information exchange and misunderstanding during the interpretation of technical data. This paper aims to find different integration techniques and steps for integrating geotechnical data into the BIM process. Methods used to examine the topic are qualitative research, literature review, and case studies. These methods were useful for studying the problems and introducing the soil information into the BIM application. Firstly, a case study with I-BIM was considered, and the BIM–FEM–BIM interaction was applied to introduce geotechnical information with Plaxis 3D. The results have shown that further development of BIM in infrastructure is needed. Another case study explored the present state of the geotechnical design in BIM and potential solutions. The new frameworks were recreated: many boreholes were imported to the BIM, and a 3D geometric model of the entire hill was created for the hill fortification structure with soil clogging. The last two studies in Malaysia modeled a 3D subsurface and used two geotechnical formats, AGS and CVS. The first includes more information than the second; however, the second can be used for a more generalized model. Overall, BIM–FEM interaction can be used as a geometric model for data transfer. However, the open data format of the Industry Foundation Class (IFC) or geotechnical data format of the AGS and CVS were suggested to be used for greater flexibility. It was also found that excessive information makes the model loaded and complex. Therefore, it was recommended that big data be summarized properly with minimal loss of necessary data. Further research is needed to understand data transmission schemes of geotechnical information better. Moreover, it is recommended to put all the strategies directly into practice to create a geotechnical design.

Ambient Noise Tomography for Coral Islands

As valuable land in the ocean, coral islands are not only important bases for making use of marine resources and protecting marine rights and interests, but also important for breakthrough research in many fields of earth science. Hence, the economical and efficient determination of the underground structure of coral islands has become significant in coral island engineering geology, but remains challenging for traditional marine geophysical prospecting and drilling methods. While ambient noise tomography with dense arrays has been widely used in continental regions, its applicability to coral islands remains undetermined. In this study, based on the data recorded by a dense array on an isolated coral island in the South China Sea, we analyzed the ambient noise characteristics and obtained a 3D subsurface structure of the coral island using ambient noise tomography. We made the following findings: ① The ambient noise frequencies can be roughly categorized into three levels: < 1, 1–5, and > 5 Hz. The spectral characteristics of the noise below 5 Hz were consistent at different stations, but there were significant differences in the characteristics of the noise above 5 Hz. ② For ambient noise frequencies below 5 Hz, cross-correlation functions with high quality could be obtained with only 24 h of waveform data. However, it was difficult to extract meaningful cross-correlation functions for ambient noise frequencies above 5 Hz. ③ The S-wave velocity in the coral island was higher toward the sea and lower toward the lagoon, which was consistent with the high degree of cementation of the outer reef flat stratum on the seaward side. ④ There were two low-velocity horizons at 25–75 and 200–300 m, which were in good agreement with the high-porosity horizons that were revealed by drilling core samples, reflecting the weathering history of the reef. Our research demonstrates that ambient noise tomography is a potentially economical, efficient, and environmentally friendly method for the geological prospecting of coral reefs.

Hybrid achromatic microlenses with high numerical apertures and focusing efficiencies across the visible

Compact visible wavelength achromats are essential for miniaturized and lightweight optics. However, fabrication of such achromats has proved to be exceptionally challenging. Here, using subsurface 3D printing inside mesoporous hosts we densely integrate aligned refractive and diffractive elements, forming thin high performance hybrid achromatic imaging micro-optics. Focusing efficiencies of 51–70% are achieved for 15μm thick, 90μm diameter, 0.3 numerical aperture microlenses. Chromatic focal length errors of less than 3% allow these microlenses to form high-quality images under broadband illumination (400–700 nm). Numerical apertures upwards of 0.47 are also achieved at the cost of some focusing efficiency, demonstrating the flexibility of this approach. Furthermore, larger area images are reconstructed from an array of hybrid achromatic microlenses, laying the groundwork for achromatic light-field imagers and displays. The presented approach precisely combines optical components within 3D space to achieve thin lens systems with high focusing efficiencies, high numerical apertures, and low chromatic focusing errors, providing a pathway towards achromatic micro-optical systems.

PO-02-212 TWISTED FIBROTICALLY-INSULATED MYOBUNDLES CREATE THE COMMON TRACK AND PIVOTING POINTS FOR REENTRANT DRIVERS MAINTAINING PERSISTENT AF IN-VIVO

Localized reentrant drivers have been suggested to maintain persistent atrial fibrillation (perAF). However, their existence, electrophysiological mechanisms, and structural substrates in-vivo are controversial. Define AF driving mechanism in-vivo and ex-vivo in a perAF canine model. A clinically relevant canine model (n=5) of self-sustained (4 months) perAF was developed to study AF maintaining mechanisms both in-vivo and ex-vivo by integrated 3D imaging, including high-resolution near-infrared optical mapping (NIOM, 0.3-0.5mm resolution), contrast-enhanced magnetic resonance imaging (CE-MRI, 0.15mm resolution). Transmural fibrosis was validated with histology. CE-MRI was used to quantitatively analyze 3D atrial structure, including transmural fiber orientation and fibrosis content in the sub-epi, intramural, and sub-endocardial layers. In all canine models (n=5), sequential in-vivo to ex-vivo sub-surface NIOM and 3D CE-MRI revealed that the self-sustained perAF was primarily maintained by localized reentrant drivers with temporal stability 53±17% (n=10), which recurrent on the same arrhythmogenic hubs. The 3D structure of the arrhythmogenic driver hubs includes a common intramural reentry track and multiple return paths formed by fibrotically-insulated bundles. The in-vivo spatio-temporal activation patterns of reentrant drivers were preserved ex-vivo (Figure A), including driver AFCL (91±7 ms vs 96±12 ms, p>0.05) and temporal driver stability (51±8% vs 52±6%, n=4, p>0.05). Targeted ablation of the driver reentrant tracks terminated or converted AF to slower tachycardia. Quantitative CE-MRI analysis of 3D myofiber orientations revealed that the transmural fiber twist along the common reentry tracks and pivoting points is greater than in non-driver regions (Figure B-C). The driver reentrant tracks had increased interstitial fibrotic strands and significantly higher intramural fibrosis than the neighboring non-driver regions (Figure D-E). Our results provide in-vivo proof-of-concept that spatially stable reentrant drivers are a mechanism of perAF maintenance. Intramural fibrotic strands and fibrotically-insulated myobundles in regions with the highest transmural myofiber twist which may form the common tracks and pivoting points for return paths of the spatially-stable reentrant AF driver.