Forward modeling and inversion of tensor CSAMT in 3D anisotropic media

Tensor controlled-source audio-frequency magnetotellurics （CSAMT） can yield information about electric and magnetic fields owing to its multi-transmitter configuration compared with the common scalar CSAMT. The most current theories, numerical simulations, and inversion of tensor CSAMT are based on far-field measurements and the assumption that underground media have isotropic resistivity. We adopt a three-dimensional （3D） staggered-grid finite difference numerical simulation method to analyze the resistivity in axial anisotropic and isotropic media. We further adopt the limited-memory Broyden- Fletcher-Goldfarb-Shanno （LBFGS） method to perform 3D tensor CSAMT axial anisotropic inversion. The inversion results suggest that when the underground structure is anisotropic, the isotropic inversion will introduce errors to the interpretation.

2D joint inversion of CSAMT and magnetic data based on cross-gradient theory

A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm.

Compensation of body shake errors in terahertz beam scanning single frequency holography for standoff personnel screening

In the terahertz（THz） band, the inherent shake of the human body may strongly impair the image quality of a beam scanning single frequency holography system for personnel screening. To realize accurate shake compensation in imaging processing, it is quite necessary to develop a high-precision measure system. However, in many cases, different parts of a human body may shake to different extents, resulting in greatly increasing the difficulty in conducting a reasonable measurement of body shake errors for image reconstruction. In this paper, a body shake error compensation algorithm based on the raw data is proposed. To analyze the effect of the body shake on the raw data, a model of echoed signal is rebuilt with considering both the beam scanning mode and the body shake. According to the rebuilt signal model, we derive the body shake error estimated method to compensate for the phase error. Simulation on the reconstruction of point targets with shake errors and proof-of-principle experiments on the human body in the 0.2-THz band are both performed to confirm the effectiveness of the body shake compensation algorithm proposed.

RESEARCH ON TECHNOLOGY OF CONTROLLED SOURCE RADIO MAGNETOTELLURIC SYSTEM'S TRANSMITTER

The basic theory and principle of the artificial Controlled Source Radio Magneto Telluric(CSRMT) method is studied and a novel CSRMT transmitter in kH z frequency band is designed. The specific circuit and measured results of some key modules in transmitter are presented, and some outdoor experimental tests have been carried out, which shows that the completed prototype of transmitter can generate a continuous sine current with frequency up to 35.33 k Hz, peak-to-peak amplitude up to 40 A, and realize a reliable transmitting mode with multi-frequency and high-current. The transmitter has a wide operating band, large magnetic moment and high waveform fidelity, and can meet the requirements of shallow geological exploration with high-resolution.

The Realization of Terahertz Image Reconstruction with High Resolution Based on the Amplitude of the Echoed Wave by using the Phase Retrieval Algorithm

An indirect imager working at terahertz band is presented and implemented,which is suitable for high-resolution planar object detection.The proposed imager employs a simple quasi-optics design to transmit and to receive terahertz waves,and adopts incoherent detection technology to extract the intensity of echoed signal,which results in a relatively low complexity and cost.Moreover,the Fienup Fourier phase-retrieval algorithm is successfully modified and is applied to retrieve the phase of the echoed signal and reconstruct the target image.Imaging experiments on typical planar objects are performed with the imager working at 0.2 THz,and the experimental results demonstrate the good performance of the proposed imager and validate the effectiveness of the reconstruction algorithm.

Study of a millimeter-wave squint indirect holographic algorithm suitable for imaging with large field-of-view

In this paper a millimeter-wave （MMW） squint indirect holographic method is presented, which is suitable for imaging with a large field-of-view. The proposed system employs the squint operation mode to remove the background and twin- image interferences, which achieves a similar effect to off-axis holography but leaves out the large-aperture quasi-optical component. The translational scanning manner enables a large field of view and ensures the image uniformity, which is difficult to realize in off-axis holography. In addition, a corresponding imaging algorithm for the presented scheme is developed to reconstruct the image from the recorded hologram. Some imaging results on typical objects, obtained with electromagnetic simulation, demonstrate good performance of the imaging scheme and validate the effectiveness of the image reconstruction algorithm.

Single-layer dual-band terahertz filter with weak coupling between two neighboring cross slots

A dual-band terahertz（THz） filter consisting of two different cross slots is designed and fabricated in a single molybdenum layer. Experimental verification by THz time-domain spectroscopy indicates good agreement with the simulation results. Owing to the weak coupling between the two neighboring cross slots in the unit cell, good selectivity performance can be easily achieved, both in the lower and higher bands, by tuning the dimensions of the two crosses. The physical mechanisms of the dual-band resonant are clarified by using three differently configured filters and electric field distribution diagrams. Owing to the rotational symmetry of the cross-shaped filter, the radiation at normal incidence is insensitive to polarization. Compared with the THz dual-band filters that were reported earlier, these filters also have the advantages of easy fabrication and low cost, which would find applications in dual-band sensors, THz communication systems, and emerging THz technologies.

Experimental demonstration of an invisible cloak with irregular shape by using tensor transmission line metamaterials

We present the design and the experimental demonstration of an invisible cloak with irregular shape by using tensor transmission line（TL） metamaterials. The fabricated cloak consists of tensor TL unit cells exhibiting anisotropic effective material parameters, while the background medium consists of isotropic TL unit cells. The simulated and the measured field patterns around the cloak show a fairly good agreement, both demonstrate that the fabricated cloak can shield the cloaked interior area from electromagnetic fields without perturbing the external fields. The scattering of the cloaked perfect electric conductor（PEC） is minimized. Furthermore, the nonresonant property of the TL structure results in a relatively broad bandwidth of the realized cloak, which is clearly observed in our experiment.

AN ULTRA-WIDEBAND ANTENNA FOR THE REINFORCED CONCRETE DETECTION

In this paper,an Ultra-Wideband(UWB) planar antenna is proposed for the reinforced concrete detection,which consists of a pair of planar waterdrop arms,a microstrip to coplanar parallel-strips transition and a shallow rectangular cavity.In order to overcome the disadvantages of the shallow cavity,some absorbing material is loaded to weaken the narrow-band effect of the cavity and the crosstalk interference.The simulated and measured results show that the proposed antenna has a large bandwidth from 0.48 GHz to 3.6 GHz with Voltage Standing Wave Ratio(VSWR) below 2 and a fractional bandwidth about 200% under the center frequency of 1.6 GHz,directional radiation characteristics and small late-time ringing in the time domain,which can be suitable for nondestructive detection of the reinforced concrete.

EFFICIENT COMPUTATIONAL METHOD FOR TRANSIENT EM RESPONSE OF LARGE-LOOP SOURCE OVER LAYERED EARTH

Calculations of the ElectroMagnetic(EM) response produced by a large horizontal loop placed over layered medium are rather complex because its integral expression contains the product of two Bessel functions and has a divergent term. In this paper, an improved fast Hankel and Gaver-Stehfest transforms are introduced to solve the strong-oscillation and slow-decay properties of the integrand, where one Bessel function in the product is substituted by a carefully chosen polynomial of high accuracy and the other used as the digital filter coefficients in the convolution integral. Comparisons prove the validity and the efficiency of the proposed method.

The Absorbing Properties of Two-Dimensional Plasma Photonic Crystals

Plasma photonic crystals composed of periodic plasma and dielectric materials have attracted considerable attention because of their tunable photonic band gaps,but only their band structures or negative refractive index properties have been addressed in previous works.In this paper,through studying the transmission and reflection characteristics of two types of twodimensional plasma photonic crystals,it is found that plasma photonic crystals play an important role in absorbing waves,and they show broader band and higher amplitude absorption characteristics than bulk plasmas.Also,the absorption of plasma photonic crystals can be tuned via plasma parameters;varying the collision frequency can make the bandwidth and amplitude tunable,but cannot change the central frequency,whereas varying the plasma frequency would control both the location and the amplitude of the absorbers.These features of plasma photonic crystals have potential for terahertz tunable absorber applications.

Distributed Field Rotator Composed of Isolated Components

Transformation optics offers remarkable control over electromagnetic fields and has recently opened an exciting gateway to design ＇field rotator devices＇. We propose a distributed field rotator with open windows based on composite transformation optics, which consists of a central circular region and several isolated components. The number, position and size of the components can be controlled freely by the design purpose. Full-wave simulations are performed to demonstrate its function, which is equivalent to a classic field rotator. However, such a distributed rotator makes it much easier to access and make use of the rotated field in the central region, compared to the closed classic field rotator, especially in the case of 3D situations.

Laser-Chirp Controlled Terahertz Wave Generation from Air Plasma

We report the laser-chirp controlled terahertz(THz) wave generation from two-color-laser-induced air plasma.Our experimental results reveal that the THz wave is affected by both the laser energy and chirp,leading to radiation minima that can be quantitatively reconstructed using the linear-dipole-array model.The phase difference between the two colors,determined by the chirp and intensity of the laser,can account for the radiation minima.Furthermore,we observe an asynchronous variation in the generated THz spectrum,which suggests a THz frequency-dependent phase matching between the laser pulse and THz wave.These results highlight the importance of laser chirp during the THz wave generation and demonstrate the possibility of modulating the THz yields and spectrum through chirping the incident laser pulse.This work can provide valuable insights into the mechanism of plasma-based THz wave generation and offer a unique means to control THz emissions.

A new image processing method for discriminating internal layers from radio echo sounding data of ice sheets via a combined robust principal component analysis and total variation approach

Discriminating internal layers by radio echo sounding is important in analyzing the thickness and ice deposits in the Antarctic ice sheet.The signal processing method of synthesis aperture radar(SAR)has been widely used for improving the signal to noise ratio(SNR)and discriminating internal layers by radio echo sounding data of ice sheets.This method is not efficient when we use edge detection operators to obtain accurate information of the layers,especially the ice-bed interface.This paper presents a new image processing method via a combined robust principal component analysis-total variation(RPCA-TV)approach for discriminating internal layers of ice sheets by radio echo sounding data.The RPCA-based method is adopted to project the high-dimensional observations to low-dimensional subspace structure to accelerate the operation of the TV-based method,which is used to discriminate the internal layers.The efficiency of the presented method has been tested on simulation data and the dataset of the Institute of Electronics,Chinese Academy of Sciences,collected during CHINARE 28.The results show that the new method is more efficient than the previous method in discriminating internal layers of ice sheets by radio echo sounding data.

NON-DESTRUCTIVE PAVEMENT LAYER THICKNESS MEASUREMENT USING EMPIRICAL MODE DECOMPOSITION WITH GPR

Ground Penetrating Radar(GPR) is an effective Non-Destructive Testing(NDT) technique for highway pavement surveys, which is able to acquire continuous pavement data compared with traditional core drilling method. In this study, we proposed an accurate and efficient method to estimate the thickness of each pavement layer using an air-coupled GPR system. For this work, the main difficulties are estimating each pavement layer's time delay and dielectric constant. We first give the basic signal model for pavement evaluation, and then present an Intrinsic Mode Functions(IMFs) product detector to determine each pavement layer's time delay. This method is based on Empirical Mode Decomposition(EMD), which is an adaptive signal decomposition procedure and proved to be suitable for suppressing noises in GPR signal. The dielectric constant was determined by metal reflection measurement. The laboratory and highway experiments illustrate that the proposed thickness estimation method yields reasonable result, thus meets the requirements of practical highway pavement survey with massive GPR data.

Numerical Analysis of Electromagnetic Fields in Multiscale Model

Modeling technique for electromagnetic fields excited by antennas is an important topic in computational electromagnetics, which is concerned with the numerical solution of Maxwell's equations. In this paper, a novel hybrid technique that combines method of moments(MoM) with finite-difference time-domain(FDTD) method is presented to handle the problem. This approach employed Huygen's principle to realize the hybridization of the two classical numerical algorithms. For wideband electromagnetic data, the interpolation scheme is used in the MoM based on the dyadic Green's function. On the other hand, with the help of equivalence principle, the scattered electric and magnetic fields on the Huygen's surface calculated by MoM are taken as the sources for FDTD. Therefore, the electromagnetic fields in the environment can be obtained by employing finite-difference time-domain method. Finally, numerical results show the validity of the proposed technique by analyzing two canonical samples.

PULSE SHRINKING TIME-TO-DIGITAL CONVERTER FOR UWB APPLICATION

A kind of architecture of Time-to-Digital Converter(TDC) for Ultra-WideBand(UWB) application is presented. The proposed TDC is based on pulse shrinking, and implemented in a Field Programmable Gate Array(FPGA) device. The pulse shrinking is realized in a loop containing two Programmable Delay Lines(PDLs) or a two-channel PDL. One line(channel) delays the rising edge and the other line(channel) delays the falling edge of a circulating pulse. Delay resolution of PDL is converted into a digital output code under known conditions of pulse width. This delay resolution measurement mechanism is different from the conventional time interval measurement mechanism based on pulse shrinking of conversion of unknown pulse width into a digital output code. This mechanism automatically avoids the influence of unwanted pulse shrinking by any circuit element apart from the lines. The achieved relative errors for four PDLs are within 0.80%–1.60%.

4-ELEMENT PLANAR ARRAY ANTENNA FOR UWB APPLICATION

A novel low-cost 4-element planar array antenna directly fed by a coaxial cable for Ultra-WideBand(UWB) application is presented. The proposed antenna consists of 2×2 bowtie elements and a simple 1:4 power divider feeding network. Compared to the basic bowtie element, the impedance bandwidth of the array antenna has a significant improvement that the low cut-off frequency is extended from 6 GHz of the bowtie element to 2 GHz. The measured results show that the proposed antenna has a large bandwidth of 2 GHz to 11 GHz for Voltage Standing Wave Ratio(VSWR)<2, and exhibits a bidirectional radiation pattern and a modest gain across the operating band and a peak gain of about 9 dBi at 11 GHz.

UWB TEM HORN ANTENNA FOR THE ASPHALT PAVEMENT INVESTIGATION

A simple Ultra-WideBand(UWB)exponentially-tapered Transverse ElectroMagnetic(TEM)horn antenna is presented for the asphalt detection based on Ground Penetrating Radar(GPR).In order to reduce the reflections from the antenna aperture,some absorbing material is loaded on the outer surface of the conductor.Comparing with the traditional TEM horn antenna,the proposed antenna has a small size and a large impedance bandwidth.Simulated and measured results show that the proposed TEM horn antenna has a low Voltage Standing Wave Ratio(VSWR)below 2 over the whole band from 0.35 GHz to 12 GHz,good radiation characteristics,and small late-time ringing,which can perfectly meet the requirements of the GPR application.

Numerical Simulation of Antennas with Improved Integral Equation Method

Simulating antennas around a conducting object is a challenge task in computational electromagnetism,which is concerned with the behaviour of electromagnetic fields. To analyze this model efficiently, an improved integral equation-fast Fourier transform(IE-FFT) algorithm is presented in this paper. The proposed scheme employs two Cartesian grids with different size and location to enclose the antenna and the other object, respectively. On the one hand, IE-FFT technique is used to store matrix in a sparse form and accelerate the matrix-vector multiplication for each sub-domain independently. On the other hand, the mutual interaction between sub-domains is taken as the additional exciting voltage in each matrix equation. By updating integral equations several times, the whole electromagnetic system can achieve a stable status. Finally, the validity of the presented method is verified through the analysis of typical antennas in the presence of a conducting object.