3D Inversion Research Articles

Shallow crustal electrical structure of the Qingchengzi orefield in Liaodong area revealed by three-dimensional inversion of Magnetotelluric data

Dense array Magnetotelluric sounding (MT) was applied in the study of the emergence of massive gold deposits relating to the tectonic activity responsible for the creation of the North China Craton. This study focuses on understanding the shallow crustal electrical structure characteristics (>1.5 km) of PbZn and AuAg deposits in the Qingchengzi orefield, through the analysis of shallow electrical resistivity distribution at depths less than 5 km, obtained by 3D inversion of full impedance at a frequency range of 320–0.107 Hz. The results showed that the electrical characteristics of Xiaotongjiapuzi and Taoyuan deposits in the 3D resistivity model are consistent with the known 2D resistivity model, and the Jianshanzi fault plays an important role in the gold mineralization, acting as a migration channel and precipitation zone for ore-forming materials. It is inferred that the fault extends at least 5 km deep with deep mineralization prospects. The strata of Liaohe Group in the north of the orefield, the intrusive granite in the south, and the NW-trending Jianshanzi and Yushanggou faults show distinguishable electrical properties, which are either high resistance, low resistance, or high and low resistance transition zone on the depth slice map. Furthermore, the Xinling rock mass has a NE trending extension, which confirms the inference that the Jianshanzi fault extends to a depth of 5 km into the crust.

Correction to: A parallel improved PSO algorithm with genetic operators for 2D inversion of resistivity data

Correction to: A parallel improved PSO algorithm with genetic operators for 2D inversion of resistivity data

3D imaging of the subsurface electrical resistivity structure in West Bohemia/Upper Palatinate covering mofettes and Quaternary volcanic structures by using Magnetotellurics

The region of West Bohemia and Upper Palatinate belongs to the West Bohemian Massif. The study area is situated at the junction of three different Variscan tectonic units and hosts the ENE-WSW trending Ohře Rift as well as many different fault systems. The entire region is characterized by ongoing magmatic processes in the intra-continental lithospheric mantle expressed by a series of phenomena, including e.g. the occurrence of repeated earthquake swarms and massive degassing of mantle derived CO2 in form of mineral springs and mofettes. Ongoing active tectonics is mainly manifested by Cenozoic volcanism represented by different Quaternary volcanic structures. All these phenomena make the Ohře Rift a unique target area for European intra-continental geo-scientific research. With magnetotelluric (MT) measurements we image the subsurface distribution of the electrical resistivity and map possible fluid pathways. Two-dimensional (2D) inversion results by Muñoz et al. (2018) reveal a conductive channel in the vicinity of the earthquake swarm region that extends from the lower crust to the surface forming a pathway for fluids into the region of the mofettes. A second conductive channel is present in the south of their model; however, their 2D inversions allow ambiguous interpretations of this feature. Therefore, we conducted a large 3D MT field experiment extending the study area towards the south. The 3D inversion result matches well with the known geology imaging different fluid/magma reservoirs at crust-mantle depth and mapping possible fluid pathways from the reservoirs to the surface feeding known mofettes and spas. A comparison of 3D and 2D inversion results suggests that the 2D inversion results are considerably characterized by 3D and off-profile structures. In this context, the new results advocate for the swarm earthquakes being located in the resistive host rock surrounding the conductive channels; a finding in line with observations e.g. at the San Andreas Fault, California.

3D Phase tensor inversion and Joint inversion using impedance and phase tensor components aiming better subsurface mapping and interpretation in Dhanjori basin, Singhbhum craton, eastern India

Dhanjori formation, a significant member of Singhbhum Craton, Eastern India, is well-known for its complex mineralogical character and huge content of Au-mineralization. This 2.1 Ga old Paleoproterozoic formation is composed of siliciclastic metasedimentary rocks interlayered with mafic/ultramafic volcanic rocks. Traces of gold has been spotted by many researchers in the sulfide-hosted mafic-ultramafic volcanic rocks and quartz pebble conglomerate units by different studies. Regional Audio Magnetotelluric (AMT) survey was conducted in 10 kHz–10 Hz frequency range over various lithological units of Dhanjori and adjacent Iron Ore Group. The key challenge faced during the analysis of magnetotelluric (MT) data is due to the presence of near-surface heterogeneities. The amplitude of the MT response gets distorted by galvanic effects which further create difficulties in modelling and inversion. If the effect is not removed in advance, it might lead to erroneous subsurface modelling and interpretations. On the contrary, the phase tensor data does not show any change in values even in the presence of galvanic distortion. Though phase tensors are the computed values from impedance tensor, still the response is free from any shallow level ambiguities unlike impedance data. We have performed 3-D inversion of phase tensor data to generate a conductivity model of Dhanjori basin. The objective was to verify whether the inclusion of phase tensor components in the inversion routine was able to provide some extra information that might add value to the analysis. Model from phase tensor inversion indicates only relative variation in the subsurface conductivity structure due to lack of amplitude information in the data. Hence, we also performed 3D joint inversion including both impedance and phase tensor components for obtaining the final conductivity model of Dhanjori to recover the absolute subsurface conductivity structure unaffected by galvanic distortion. Correlation of the conductivity model with available geoscientific information in Dhanjori and adjacent geological units helped us to relate the conductors to the sulfide mineralization prospect.

3D Inversion of Gravity Data Modeling Using The Chi Fact Algorithm for Revealing Subsurface Structure in Semarang City

<p class="AbstractText"><span lang="EN-AU">The interpretation of subsurface in the Semarang City with 3D inversion model using the Chi Fact algorithm has been carried out to reveal the subsurface especially related to the presence of groundwater basins. The gravity data model in this research uses 80 stations with a research area of 20 km<sup>2</sup> which is specialized in the surrounding area of Kaligarang district. The results of the 3D inversion using the Chi Fact algorithm founded the distribution of subsurface rock density values in the Semarang City area ranging from 1.6 gr/cc to 2.98 gr/cc. The distribution of subsurface density values indicates the subsurface geological structure of the Semarang City area had a normal fault leading to the southeast.</span></p>

Investigation of the Inter- and Intrascanner Reproducibility and Repeatability of Radiomics Features in T1-Weighted Brain MRI.

Radiomics is the high throughput analysis of medical images using computer algorithms, which specifically assess textural features. It has increasingly been proposed as a tool for the development of imaging biomarkers. However, an important acknowledged limitation of radiomics is the lack of reproducibility of features produced. To assess reproducibility and repeatability of radiomics variables in brain MRI through a multivisit, multicenter study. Retrospective. Fourteen individuals visiting three institutions twice, 10 males with the mean age of 36.3 years and age range 25-51. 3D T1W inversion recovery on three 1.5-T General Electric scanners. Radiomics analysis by a consultant radiologist performed on the T1W images of the whole brain on all visits. All possible radiomics features were generated. Concordance correlation coefficient (CCC) and dynamic range (DR) for all variables were calculated to assess the test-retest repeatability. Intraclass correlation coefficients (ICCs) were calculated to investigate the reproducibility of features across centers. Of 1596 features generated, 57 from center 1, 15 from center 2, and 22 from center 3 had a CCC > 0.9 and DR > 0.9. Eight variables had CCC > 0.9 and DR > 0.9 in all centers. Forty-one variables had an ICC of >0.9. No variables had CCC > 0.9, DR > 0.9, and ICC > 0.9. Repeatability and reproducibility of variables is a significant limitation of radiomics analysis in 3DT1W brain MRI. Careful selection of radiomic features is required. 4 TECHNICAL EFFICACY STAGE: 2.

Drone-based magnetic and multispectral surveys to develop a 3D model for mineral exploration at Qullissat, Disko Island, Greenland

Abstract. Mineral exploration in the West Greenland flood basalt province is attractive because of its resemblance to the magmatic sulfide-rich deposit in the Russian Norilsk region, but it is challenged by rugged topography and partly poor exposure for relevant geologic formations. On northern Disko Island, previous exploration efforts have identified rare native iron occurrences and a high potential for Ni–Cu–Co–PGE–Au mineralization. However, Quaternary landslide activity has obliterated rock exposure in many places at lower elevations. To augment prospecting field work under these challenging conditions, we acquire high-resolution magnetic and multispectral remote sensing data using drones in the Qullissat area. From the data, we generate a detailed 3D model of a mineralized basalt unit, belonging to the Asuk Member of the Palaeocene Vaigat Formation. Different types of legacy data and newly acquired geo- and petrophysical as well as geochemical-mineralogical measurements form the basis of an integrated geological interpretation of the unoccupied aerial system (UAS) surveys. In this context, magnetic data aim to define the location and the shape of the buried magmatic body, and to estimate if its magnetic properties are indicative for mineralization. UAS-based multispectral orthomosaics are used to identify surficial iron staining, which serves as a proxy for outcropping sulfide mineralization. In addition, UAS-based digital surface models are created for geomorphological characterization of the landscape to accurately reveal landslide features. UAS-based magnetic data suggest that the targeted magmatic unit is characterized by a pattern of distinct positive and negative magnetic anomalies. We apply a 3D magnetization vector inversion (MVI) model to the UAS-based magnetic data to estimate the magnetic properties and shape of the magmatic body. By means of introducing constraints in the inversion, (1) UAS-based multispectral data and legacy drill cores are used to assign significant magnetic properties to areas that are associated with the mineralized Asuk Member, and (2) the Earth's magnetic and the palaeomagnetic field directions are used to evaluate the general magnetization direction in the magmatic units. Our results suggest that the geometry of the mineralized target can be estimated as a horizontal sheet of constant thickness, and that the magnetization of the unit has a strong remanent component formed during a period of Earth's magnetic field reversal. The magnetization values obtained in the MVI are in a similar range to the measured ones from a drillcore intersecting the targeted unit. Both the magnetics and topography confirm that parts of the target unit were displaced by landslides. We identified several fully detached and presumably rotated blocks in the obtained model. The model highlights magnetic anomalies that correspond to zones of mineralization and is used to identify outcrops for sampling. Our study demonstrates the potential and efficiency of using high-resolution UAS-based multi-sensor data to constrain the geometry of partially exposed geological units and assist exploration targeting in difficult or poorly exposed terrain.

Three-Dimensional Gravity Inversion in the Presence of the Sediment-Basement Interface: A Case Study in Utah, USA

We introduce a novel approach to three-dimensional gravity inversion in the presence of the sediment-basement interface with a strong density contrast. This approach makes it possible to incorporate the known information about the basement depth in the inversion. It also allows the user to determine the depth-to-basement in the initial inversion phase. One can then use this interface to constrain the final inversion phase. First, the inversion generates the depth-to-basement model based on the 3D Cauchy-type integral representation of the gravity field. Then, in the second phase, full 3D voxel-type inversion applies the depth-to-basement model determined in the first phase as an a priori constraint. We use this approach to the 3D inversion of the Bouguer gravity anomaly data observed in Utah, USA. The results of inversion generated a 3D density model of the top layers of the earth’s crust, including unconsolidated sediments and the top of the crystalline basement.

A parallel improved PSO algorithm with genetic operators for 2D inversion of resistivity data

A parallel improved PSO algorithm with genetic operators for 2D inversion of resistivity data

Accelerating Non-LTE Synthesis and Inversions with Graph Networks

Abstract The computational cost of fast non-LTE synthesis is one of the challenges that limits the development of 2D and 3D inversion codes. It also makes the interpretation of observations of lines formed in the chromosphere and transition region a slow and computationally costly process, which limits the inference of the physical properties on rather small fields of view. Having access to a fast way of computing the deviation from the LTE regime through the departure coefficients could largely alleviate this problem. We propose to build and train a graph network that quickly predicts the atomic level populations without solving the non-LTE problem. We find an optimal architecture for the graph network for predicting the departure coefficients of the levels of an atom from the physical conditions of a model atmosphere. A suitable data set with a representative sample of potential model atmospheres is used for training. This data set has been computed using existing non-LTE synthesis codes. The graph network has been integrated into existing synthesis and inversion codes for the particular case of Ca ii. We demonstrate orders-of-magnitude gain in computing speed. We analyze the generalization capabilities of the graph network and demonstrate that it produces good predicted departure coefficients for unseen models. We implement this approach in Hazel2 and show how the inversions nicely compare with those obtained with standard non-LTE inversion codes. Our approximate method opens up the possibility of extracting physical information from the chromosphere on large fields of view with time evolution.

Resistivity imaging from magnetotelluric soundings of the Rio Preto Belt (NE Brazil): The role of lateral variations in the Precambrian basement on the evolution of phanerozoic basins

We present a resistivity model for the boundary of the Paleozoic Parnaíba Basin and the Cretaceous Sanfranciscana Basin (NE Brazil). The boundary underlies the Neoproterozoic Rio Preto belt, whose basement consists of reworked Archean and Paleoproterozoic rocks. In order to characterize the subsurface geoelectrically, 17 broadband magnetotelluric (MT) soundings were acquired and distributed along a ∼90 km profile across Rio Preto belt. The transfer functions were obtained with a robust code, thus obtaining reliable results for all four components of the impedance tensor. Some apparent resistivity and phase curves were shown to be heavily contaminated by local electromagnetic noise at low frequencies, and others showed minimal radio signals associated with the deadband frequency range. In this study, the MT data set was interpreted using 2D and 3D inversion procedures. Comparing the two MT models, the crustal-scale structure is well resolved in the 2D model, but the 3D model explains significantly better the measured data related to deep/shallow structures and proved to be more sensitive to conductive anomalies. The final resistivity image derived from the 3D inversion denotes a resistive basement with vertical conductive anomalies in the SE portion and a conductive domain in the NW portion. The homogeneous conductive block corresponds to the Archean domain of the Rio Preto belt, and the heterogeneous resistive block has the morphology and physical properties of the Paleoproterozoic block. Lateral geological correlation and gravity data support that the basement of the Parnaíba Basin in the study area is formed by the Archean block of the Rio Preto Belt. The boundary between the resistive blocks is a reverse fault that is also the tectonic front of a sub-horizontal detachment fault. In the Paleoproterozoic block, the conductors are associated with graphitic and volcano-sedimentary sequences, the dome and keel structures explain the geometry of the conductors. Additionally, graben-likes features were revealed in this block, interpreted as Cambrian basins filled by Cretaceous sediments of the Urucuia and Areado formations. These structures have a good correlation with those defined in the eastern edge of the Parnaíba Basin in previous studies. Towards the NW end of the section, the resistive blocks abruptly turn into a high conductivity domain (HCD) covered by a shallow conductive layer of Parnaíba Basin sandstones that modifies the electrical properties of the basement by interaction of basal fluids.

3D resistivity survey over mapped caves in eogenetic karst terrane, west-central Florida, USA

This study assesses the capability and practical applications of quasi-3 Dimensional (3D) Electrical Resistivity Tomography (ERT) for mapping air-filled conduits in eogenetic karst. A high-resolution quasi-3D ERT survey, consisting of multiple parallel 110m-long 2D profiles, was conducted over two mapped cave systems on the Brooksville Ridge, Florida. The irregularly shaped caves have diameters ranging up to 4 m and span depths from 3 m to 11 m below ground surface. Dipole-dipole array geometries with L2 (least squares) rather than L1 (robust) inversion produced the best fits of resistivity highs with the mapped cave locations. As expected, 3D inversions of sets of parallel lines produced better results than 2D inversions of individual transects. Better imaging was obtained of a cave over which cave-parallel profiles were run in addition to cross-cave profiles. However, even with the best acquisition and processing steps, there are significant misfits in the apparent size of the large cave sections, and narrower conduits are not imaged. Resolution decreased significantly with depth, as expected given the method limitations and the site constraints on profile lengths. 3D visualization techniques are explored, and found helpful in examining the data and comparing mapped caves and 3D resistivity datasets; however, when applied to eogenetic karst terrain, ERT has limited capacity to detect smaller cavities, which may require the additional use of other geophysical or subsurface investigative methods. With sinkholes continuing to be of concern to residential and urban development in west-central Florida, the results of this research present additional insight on the potential of quasi-3D ERT methods to map and characterize the potential hazards posed by karst terranes.

3D constrained gravity inversion to model Moho geometry and stagnant slabs of the Northwestern Pacific plate at the Japan Islands

The motivation for this work was to model the crustal thickness, and the geometry of the Mohorovičić (Moho) discontinuity, the crust-upper mantle boundary around the Japanese archipelago. The western Pacific Ocean is one of the most eminent subduction zones in the world. The subduction zone complexity is due to alteration within the geometry of the stagnant slab. We compiled a three-dimensional density distribution model for subduction plates around the Japan Islands, and estimated the depth to Moho using the constrained non-linear inversion of gravity data derived from the latest Gravity field and Ocean Circulation Explorer, or GOCE, model known as GOCO06s. The study had three steps, starting with the data corrections followed by the inversion, and ending with the modeling and comparing of the estimated results with previous models and studies. In this study, we used available receiver function results, seismic tomography models, global crustal thickness, and geological and tectonic information as constraints. The model estimated the crustal thickness, surface geometry of the Moho, and geometry of the stagnant slab in the Kuril Trench, Japan Trench, Nankai Trough, and Izu-Bonin Trough. In gravity inversion, we used 34 km and 600 kg m−3 as the Moho depth reference and density contrast, respectively. Results showed that the crustal thickness ranged from 14 to 43 km with a − 0.78 mean misfit for Moho depths derived from seismic P-wave analysis. The calculated gravity data had a 0.16 mGal misfit with the observed data. The normalized analytical signal helped to sharply delineate the different plates and their edges. Three cross-sections were conducted to model the stagnant slab along the subduction zone of the Eastern Eurasian plate, western Pacific plate, Okhotsk plate, and Philippine Sea plate. The eastern area of the entire region (Pacific plate) had a low crust thickness and high density compared to the western area (Eurasian plate), which had a high crust thickness and low density.

Three-Dimensional Inversion of Semi-Airborne Transient Electromagnetic Data Based on a Particle Swarm Optimization-Gradient Descent Algorithm

Semi-airborne transient electromagnetics (SATEM) is a geophysical survey tool known for its ability to perform three-dimensional (3D) observations and collect high-density data in large volumes. However, SATEM data processing is presently restricted to 3D model-driven inversion, which is not conducive to detailed surveys. This paper presents a new 3D model- and data-driven inversion algorithm using the particle swarm optimization (PSO) and gradient descent (GD) algorithms. PSO is used to suppress the multiplicity of solutions associated with inverse problems, and the GD algorithm is employed to accelerate the convergence of the inversion process. For the PSO-GD algorithm, a new model-updating equation is established and a cosine probability function is introduced as a weighting term for PSO and GD algorithms to ensure a smooth transition between the two algorithms in the iterative process. The α-trimmed filter function is used as a regularization constraint to smooth the model. The stability and reliability of the PSO-GD algorithm are verified through numerical simulations. Finally, the new algorithm is applied to the processing of SATEM measurements of the Qinshui coal mine in Jincheng, Shanxi Province, China.

3D elastic multisource full waveform inversion on distributed GPU clusters

Full-waveform inversion (FWI) serves as a useful tool to quantitatively investigate the properties of the media. One of the greatest challenges is its enormous computational cost, especially for the three-dimensional (3D) case. The multisource encoding strategy is a crucial method to speed up FWI. However, current 3D elastic FWI algorithms with the simultaneous sources, such as random time series and random polarity methods, may be subjected to the crosstalk noise, which severely degrades inversion quality. We present a novel 3D elastic crosstalk-free multisource full-waveform inversion method (MFWI). The encoded sources with a series of single-frequency components are excited to perform the seismic wavefield simulation. Extracting the wavefield corresponding to each single-frequency harmonic source by phase-sensitive detection (PSD) algorithm, the multisource wavefield can be decomposed. Therefore, summing the gradients with respect to the model parameters of all decomposed single frequencies can effectively avoid the crosstalk in our MFWI algorithm. Furthermore, to increase the practicability and universality of the algorithm, we adopt the graphic processing unit (GPU) based on the domain decomposition to improve the MFWI algorithm. The multi-scale strategy with a flexible encoding number in each scale helps to enhance the robustness of our algorithm. The numerical examples demonstrate that our method can achieve domain decomposition simulation and obtain 3D crosstalk-free multisource elastic wave inversion results.

An efficient multigrid solver based on a four-color cell-block Gauss-Seidel smoother for 3D magnetotelluric forward modeling

Practical application of 3D magnetotelluric inversion requires efficient forward modeling of electromagnetic (EM) fields in the earth. To resolve realistic 3D structures, large computational domains and extremely large linear systems of equations are required. The iterative solvers, which are almost exclusively used to solve these systems, can be inefficient due to the abundant null space of the curl-curl operator. Multigrid (MG) solvers are considered a potentially efficient technique for solving such problems. However, due to the abundant null solution space and existence of the air layer, MG solvers can still converge slowly or even diverge. We have developed an efficient MG solver for finite-difference frequency-domain EM solution. In this algorithm, the excellent smoothing property of an efficient four-color cell-block Gauss-Seidel (GS) is exploited to remove the short-range errors effectively, and the interpolation and prolongation operators are used to handle the long-range errors. They work as a whole to speed the convergence of our algorithm remarkably. Because all of the nodes for the four-color cell-block GS are grouped into four colors and the edge components attached to different nodes in each color are completely decoupled, this can be used to develop a highly vectorized or parallelized algorithm. Another important property is that our algorithm is locally current divergence free, effectively eliminating spurious solutions in the null space of the curl-curl operator. The accuracy and efficiency of the algorithm are verified by comparing the numerical solutions obtained with our MG solver to those from the biconjugate gradient stabilized solver with different preconditioners based on synthetic models and a model from 3D inversion. Comparisons, in terms of iteration number and computational time, indicate that our algorithm is extremely stable and efficient relative to the other solvers. Our MG algorithm will be suitable for massively parallel computing as well.

Workshop Review: Digital rock physics: From pore to reservoir scale

With continuous optimization of reservoir performance in the oil and gas industry globally, there are outstanding technical challenges that require development and implementation of technologies and methodologies. One of the main challenges is how to build a static reservoir model that integrates rock-property heterogeneity to derive dynamic reservoir models and improve hydrocarbon recovery. Rock properties are measured at different scales by using advanced measurements ranging from 3D surface seismic inversion to log data and core data. The advancement of rock digital imaging has been proven to enhance interpretation results. Digital rock physics plays a vital role by utilizing new technologies such as machine learning to get the best outcome.

Probing the Southern African Lithosphere With Magnetotellurics—Part I: Model Construction

AbstractThe Southern African Magnetotelluric Experiment (SAMTEX) involved the collection of data at over 700 sites in Archean to Proterozoic southern Africa, spanning features including the Kalahari Craton, Bushveld Complex, and voluminous kimberlites. Here, we present the first 3D inversions of the full SAMTEX data set. In this paper, we focus on assessing the robustness of the 3D models by comparing two different inversion codes, jif3D and ModEM, and two different subsets of the data, one containing all acceptable data and the other containing a smaller selection of undistorted, high‐quality data. Results show that the main conductive and resistive features are imaged by all inversions, including deep resistive features in the central Kaapvaal Craton and southern Congo Craton and a lithospheric‐scale conductor beneath the Bushveld Complex. Despite this, differences exist between the jif3D and ModEM inverse models that derive mainly from the differences in regularization between the models, with jif3D producing models that are very smooth laterally and with depth, while ModEM produces models with more discrete conductive and resistive features. Analysis of the differences between these two inversions can provide a good indication of the model resolution. More minor differences are apparent between models run with different subsets of data, with the models containing all acceptable data featuring higher wavelength conductivity variations than those run with fewer stations but also demonstrating poorer data fit.

Lithologic Hydrocarbon Accumulations in the Upper Jurassic Kalazha Formation in the Foreland Slope, the South Marginal Junggar Basin, Northwest China

The foreland slope in the south marginal Junggar Basin (SMJB) is geologically favorable for mass lithologic oil and gas reservoir accumulation. Based on outcrops, thin sections, individual-well facies, and cross-well stratigraphic correlation, the stratigraphic distribution and features of the Upper Jurassic Kalazha Formation were investigated in the SMJB. There are two sedimentary sources in the southwest and southeast Kalazha Formation, which thin out from south to north and vanish in the northern slope. Fan deltas occur in the southwest, and braided river deltas occur in the southeast. Kalazha reservoir rocks and sedimentary facies were predicted using joint 2D and 3D seismic impedance inversion. According to the overall analysis, the Kalazha Formation in the foreland slope is favorable for hydrocarbon accumulation, considering structural highs as the destination of long-term hydrocarbon migration and the updip wedge-out zone in the north with potential lithologic traps. The accumulation model was established with source rocks, reservoir rocks, caprocks, and migration systems for the lithologic hydrocarbon reservoirs in the Kalazha Formation in the SMJB. As per the model and regional stratigraphic and structural features, Kalazha lithologic hydrocarbon reservoirs mainly formed in the annular updip wedge-out zone in the foreland slope. Our deliverables may provide useful information for mass lithologic hydrocarbon accumulation in the foreland slope and deep to ultra-deep hydrocarbon exploration in the SMJB.

Strongly Confining Light with Air-Mode Cavities in Inverse Rod-Connected Diamond Photonic Crystals

Three-dimensional dielectric optical crystals with a high index show a complete photonic bandgap (PBG), blocking light propagation in all directions. We show that this bandgap can be used to trap light in low-index defect cavities, leading to strongly enhanced local fields. We compute the band structure and optimize the bandgap of an inverse 3D rod-connected diamond (RCD) structure, using the plane-wave expansion (PWE) method. Selecting a structure with wide bandgap parameters, we then add air defects at the center of one of the high-index rods of the crystal and study the resulting cavity modes by exciting them with a broadband dipole source, using the finite-difference time-domain (FDTD) method. Various defect shapes were studied and showed extremely small normalized mode volumes (Veff) with long cavity storage times (quality factor Q). For an air-filled spherical cavity of radius 0.1 unit-cell, a record small-cavity mode volume of Veff~2.2 × 10−3 cubic wavelengths was obtained with Q~3.5 × 106.