1
|
Nicolson E, Lines D, Mohseni E, MacLeod CN. Single-Bit Reception With Coded Excitation for Lightweight Advanced Ultrasonic Imaging Systems. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2024; 71:1120-1131. [PMID: 38748526 DOI: 10.1109/tuffc.2024.3399743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/05/2024]
Abstract
The demand for an efficient and reliable ultrasonic phased array imaging system is not unique to a single industry. Today's imaging systems can be enhanced in a number of areas including; improving scanning and processing times, reducing data storage requirements, simplifying hardware, and prolonging probe lifespan. In this work, it is shown that by combining the use of coded excitation with single-bit data capture, a number of these areas can be improved. Despite using single-bit receive data, resolution can be recovered through the coded excitation pulse compression process, and shown to produce high signal-to-noise ratio (SNR) images of phase coherence imaging (PCI) and total focusing method (TFM) of tip diffraction in a carbon steel sample. Comparison with conventional single-cycle transmission pulses has shown that little imaging performance degradation is seen despite a significant reduction in data resolution and size. This has also been shown to be effective at low excitation voltages with gain compensation due to the obsolescence of signal saturation concerns when considering single-bit receive data. The ability to compute high-resolution ultrasonic images from low-resolution input data at low transmission voltages has important implications for data compression, acquisition and imaging performance, operator safety, and hardware simplification for ultrasonic imaging systems across industrial and medical fields.
Collapse
|
2
|
Yang Y, Wang P, Jia Y, Jing L, Shi Y, Sheng H, Jiang Y, Liu R, Xu Y, Li X. Rail fracture monitoring based on ultrasonic-guided wave technology with multivariate coded excitation. ULTRASONICS 2024; 136:107164. [PMID: 37748363 DOI: 10.1016/j.ultras.2023.107164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2023] [Revised: 08/21/2023] [Accepted: 09/12/2023] [Indexed: 09/27/2023]
Abstract
Steel lay the basis of railroad train traffic. Rail fracture is the most serious injury to the rail, which should be monitored in time. The conventional mono-pulse exciting ultrasonic guide wave (UGW) has low energy, the conventional Barker code has limited coding sequence length and the orthogonal complementary Golay code has the problem of low monitoring efficiency. This study proposes multivariate coded excitation (MCE) to excite UGW to monitor rail fracture, and the feasibility of coding and decoding the method is theoretically derived and verified. Ideally, the MCE generated based on the 3-bit Barker code (B3: 1,1, -1) and 4-bit orthogonal complementary Golay code (GA4: 1,1,1, -1; GB4: 1,1, -1,1) is calculated to have a main-lobe power level (MPL) gain of 8.1020 dB, which is significantly higher than the MPL of Barker and Golay codes. For the proposed method, finite element modeling simulation and experimental study are carried out respectively. Analyze and process the data, and calculate the gain of the difference in amplitude (DIA) between the amplitude of echo caused by rail fracture and the amplitude in the healthy rail. The gain of the DIA of the echoes in the MCE is above 50 dB (experimental data, the value of simulation data is 5 dB) under different degrees of rail fracture, while the gain of the DIA of the echoes caused by the other three excitation methods is below 40 dB (experimental data, the value of simulation data is 1 dB). Simulation and experimental results show that the MCE makes up for the shortcomings of the conventional Barker code and Golay code, improves the excitation energy of the monitoring system, and the high gain of the DIA of the echoes is more conducive to the identification of rail fracture damage.
Collapse
Affiliation(s)
- Yuan Yang
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| | - Ping Wang
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| | - Yinliang Jia
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| | - Lixuan Jing
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| | - Yu Shi
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| | - Hongwei Sheng
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| | - Yi Jiang
- College of Electrical Engineering, Nanjing Vocational University of Technology, Nanjing 210023, China.
| | - Renbao Liu
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| | - Yihang Xu
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| | - Xin Li
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| |
Collapse
|
3
|
Xu C, Gao G, Deng M. Lamb wave based damage imaging under nonlinear chirp excitation. ULTRASONICS 2023; 135:107108. [PMID: 37487295 DOI: 10.1016/j.ultras.2023.107108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Revised: 06/10/2023] [Accepted: 07/14/2023] [Indexed: 07/26/2023]
Abstract
Considering a trade-off between temporal-spatial resolution and multi-mode nature of Lamb waves, tone bursts with short durations are usually used as excitations in Lamb wave based damage detection. A short-duration excitation usually requires a large amplitude to carry sufficient energy so as to obtain response signals with enough signal-to-noise ratio and cover a large inspection area. In this paper, an alternative Lamb wave damage imaging method using nonlinear chirp (nonlinear frequency modulation, NLFM) excitation with a long duration and a small amplitude is proposed. The signal processing techniques of pulse compression and dispersion compensation are adopted to compress the long-duration wave packets of response signals into short ones. Compared with conventional tone burst excitations with short durations and small amplitudes, due to the long duration of the nonlinear chirp excitation and the use of pulse compression, sufficient energy can be applied to transducers under small amplitude excitations so the image contrast in imaging will not degrade. Furthermore, as large amplitude excitations are no longer required, high voltage amplifiers are not necessary so the hardware of the Lamb wave testing system is simplified. Experiments on a carbon steel plate with an artificial crack are carried out and Lamb wave signals are collected using a linear array consisting of nine PZTs. Experimental results under the NLFM signal and conventional tone bursts are provided. Experimental results show that under the condition of the same excitation amplitude, the proposed method under the NLFM excitation can achieve better imaging quality compared with methods under conventional tone bursts.
Collapse
Affiliation(s)
- Caibin Xu
- College of Aerospace Engineering, Chongqing University, Chongqing 400044, China
| | - Guangjian Gao
- Department of Physics, Army Logistics Academy, Chongqing 401331, China
| | - Mingxi Deng
- College of Aerospace Engineering, Chongqing University, Chongqing 400044, China.
| |
Collapse
|
4
|
Jiang C, Li Z, Zhang Z, Wang S. A New Design to Rayleigh Wave EMAT Based on Spatial Pulse Compression. SENSORS (BASEL, SWITZERLAND) 2023; 23:3943. [PMID: 37112283 PMCID: PMC10146980 DOI: 10.3390/s23083943] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 04/10/2023] [Accepted: 04/11/2023] [Indexed: 06/19/2023]
Abstract
The main disadvantage of the electromagnetic acoustic transducer (EMAT) is low energy-conversion efficiency and low signal-to-noise ratio (SNR). This problem can be improved by pulse compression technology in the time domain. In this paper, a new coil structure with unequal spacing was proposed for a Rayleigh wave EMAT (RW-EMAT) to replace the conventional meander line coil with equal spacing, which allows the signal to be compressed in the spatial domain. Linear and nonlinear wavelength modulations were analyzed to design the unequal spacing coil. Based on this, the performance of the new coil structure was analyzed by the autocorrelation function. Finite element simulation and experiments proved the feasibility of the spatial pulse compression coil. The experimental results show that the received signal amplitude is increased by 2.3~2.6 times, the signal with a width of 20 μs could be compressed into a δ-like pulse of less than 0.25 μs and the SNR is increased by 7.1-10.1 dB. These indicate that the proposed new RW-EMAT can effectively enhance the strength, time resolution and SNR of the received signal.
Collapse
|
5
|
Li X, Yang J, Zhang G. A Crack Size Quantification Method Using High-Resolution Lamb Waves. SENSORS 2021; 21:s21206941. [PMID: 34696154 PMCID: PMC8541407 DOI: 10.3390/s21206941] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Revised: 10/12/2021] [Accepted: 10/18/2021] [Indexed: 11/16/2022]
Abstract
Traditional tone burst excitation cannot attain a high output resolution, due to the time duration. The received signal is much longer than that of excitation during the propagation, which can increase the difficulty of signal processing, and reduce the resolution. Therefore, it is of significant interest to develop a general methodology for crack quantification through the optimal design of the excitation waveform and signal-processing methods. This paper presents a new crack size quantification method based on high-resolution Lamb waves. The linear chirp (L-Chirp) signal and Golay complementary code (GCC) signal are used as Lamb wave excitation signals. After dispersion removal, these excitation waveforms, based on pulse compression, can effectively improve the inspection resolution in plate-like structures. A series of simulations of both healthy plates and plates with different crack sizes are performed by Abaqus CAE, using different excitation waveforms. The first wave package of the S0 mode after pulse compression is chosen to extract the damage features. A multivariate regression model is proposed to correlate the damage features to the crack size. The effectiveness of the proposed crack size quantification method is verified by a comparison with tone burst excitation, and the accuracy of the crack size quantification method is verified by validation experiments.
Collapse
|
6
|
Gao F, Hua J, Lin J, Zeng L. Improved spectrum method for impact damage characterization in the composite beam using Lamb waves. ULTRASONICS 2021; 111:106336. [PMID: 33341687 DOI: 10.1016/j.ultras.2020.106336] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/09/2020] [Revised: 10/11/2020] [Accepted: 12/09/2020] [Indexed: 06/12/2023]
Abstract
Lamb wave spectral methods are one candidate to characterize invisible damages in composite structures. Unfortunately, multiple reflections resulting from geometric boundaries could distort the Lamb wave spectrum, which may cover the important signatures concerning structural integrity. To eliminate spectral interference, a cepstrum based filtering method is proposed to separate various reflection features. In particular, the fundamental spectrum contributing to structural integrity can be extracted smoothly by removing harmonic fluctuations regarding wave reflections with an optimized filter. Subsequently, to establish the quantitative relationships between fundamental spectrum features and the impact energies, damage indices involving spectrum energy mean and median frequency shift are introduced to quantify the change of spectrum intensity and distribution respectively, which shows good performance on impact damage characterization. Finally, the experiment was carried out on a T300 composite beam, in which strong reflections come after direct waves. Meanwhile, the severity of impact damage was simulated by the free droppings of a steel ball with different heights. The results show the effectiveness of the proposed method for characterizing impact damages with improved sensitivity.
Collapse
Affiliation(s)
- Fei Gao
- School of Reliability and Systems Engineering, Beihang University, Xueyuan Road No.37, Haidian District, Beijing, China; Advanced Manufacturing Center, Ningbo Institute of Technology, Beihang University, Ningbo 315100, China.
| | - Jiadong Hua
- School of Reliability and Systems Engineering, Beihang University, Xueyuan Road No.37, Haidian District, Beijing, China; Advanced Manufacturing Center, Ningbo Institute of Technology, Beihang University, Ningbo 315100, China.
| | - Jing Lin
- School of Reliability and Systems Engineering, Beihang University, Xueyuan Road No.37, Haidian District, Beijing, China; State Key Laboratory of Manufacturing System Engineering, Xi'an Jiaotong University, 28 Xianning Road, Xi'an, China.
| | - Liang Zeng
- State Key Laboratory of Manufacturing System Engineering, Xi'an Jiaotong University, 28 Xianning Road, Xi'an, China.
| |
Collapse
|
7
|
Ultrasonic Imaging of Thick Carbon Fiber Reinforced Polymers through Pulse-Compression-Based Phased Array. APPLIED SCIENCES-BASEL 2021. [DOI: 10.3390/app11041508] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The use of pulse-compression in ultrasonic non-destructive testing has assured, in various applications, a significant improvement in the signal-to-noise ratio. In this work, the technique is combined with linear phased array to improve the sensitivity and resolution in the ultrasonic imaging of highly attenuating and scattering materials. A series of tests were conducted on a 60 mm thick carbon fiber reinforced polymer benchmark sample with known defects using a custom-made pulse-compression-based phased array system. Sector scan and total focusing method images of the sample were obtained with the developed system and were compared with those reconstructed by using a commercial pulse-echo phased array system. While an almost identical sensitivity was found in the near field, the pulse-compression-based system surpassed the standard one in the far-field producing a more accurate imaging of the deepest defects and of the backwall of the sample.
Collapse
|
8
|
Gao F, Zeng L, Lin J, Shao Y. Damage assessment in composite laminates via broadband Lamb wave. ULTRASONICS 2018; 86:49-58. [PMID: 29361331 DOI: 10.1016/j.ultras.2018.01.005] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/03/2017] [Revised: 11/08/2017] [Accepted: 01/10/2018] [Indexed: 06/07/2023]
Abstract
Time of flight (ToF) based method for damage detection using Lamb waves is widely used. However, due to the energy dissipation of Lamb waves and the non-ignorable size of damage in composite structure, the performance of damage detection is restricted. The objective of this research is to establish an improved method to locate and assess damages in composite structure. To choose appropriate excitation parameters, the propagation characters of Lamb waves in quasi-isotropic composite laminates are firstly studied and the broadband excitation is designed. Subsequently, the pulse compression technique is adopted for energy concentration and high-accuracy distance estimation. On this basis, the gravity center of intersections of path loci is employed for damage localization and the convex envelop of identified damage edge points is taken for damage contour estimation. As a result, both damage location and size can be evaluated, thereby providing the information for quantitative damage detection. The experiment consisting of five different sizes of damage is carried for method verification and the identified results show the efficiency of the proposed method.
Collapse
Affiliation(s)
- Fei Gao
- Shaanxi Key Laboratory of Mechanical Product Quality Assurance and Diagnostics, School of Mechanical Engineering, Xi'an Jiaotong University, Xi'an, Shaanxi Province 710049, China.
| | - Liang Zeng
- Shaanxi Key Laboratory of Mechanical Product Quality Assurance and Diagnostics, School of Mechanical Engineering, Xi'an Jiaotong University, Xi'an, Shaanxi Province 710049, China.
| | - Jing Lin
- State Key Laboratory for Manufacturing Systems Engineering, Xi'an Jiaotong University, Shaanxi Province, China.
| | - Yongsheng Shao
- Shaanxi Key Laboratory of Mechanical Product Quality Assurance and Diagnostics, School of Mechanical Engineering, Xi'an Jiaotong University, Xi'an, Shaanxi Province 710049, China.
| |
Collapse
|
9
|
Malo S, Fateri S, Livadas M, Mares C, Gan TH. Wave Mode Discrimination of Coded Ultrasonic Guided Waves Using Two-Dimensional Compressed Pulse Analysis. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2017; 64:1092-1101. [PMID: 28504936 DOI: 10.1109/tuffc.2017.2693319] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Ultrasonic guided waves testing is a technique successfully used in many industrial scenarios worldwide. For many complex applications, the dispersive nature and multimode behavior of the technique still poses a challenge for correct defect detection capabilities. In order to improve the performance of the guided waves, a 2-D compressed pulse analysis is presented in this paper. This novel technique combines the use of pulse compression and dispersion compensation in order to improve the signal-to-noise ratio (SNR) and temporal-spatial resolution of the signals. The ability of the technique to discriminate different wave modes is also highlighted. In addition, an iterative algorithm is developed to identify the wave modes of interest using adaptive peak detection to enable automatic wave mode discrimination. The employed algorithm is developed in order to pave the way for further in situ applications. The performance of Barker-coded and chirp waveforms is studied in a multimodal scenario where longitudinal and flexural wave packets are superposed. The technique is tested in both synthetic and experimental conditions. The enhancements in SNR and temporal resolution are quantified as well as their ability to accurately calculate the propagation distance for different wave modes.
Collapse
|
10
|
Isla J, Celga F. The Use of Binary Quantization for the Acquisition of Low SNR Ultrasonic Signals: A Study of the Input Dynamic Range. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2016; 63:1474-1482. [PMID: 27244732 DOI: 10.1109/tuffc.2016.2571843] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Low-power excitation and/or low sensitivity transducers, such as electromagnetic acoustic transducers, piezoelectric paints, air-coupled transducers, and small elements of dense arrays, may produce signals below the noise threshold at the receiver. The information from those noisy signals can be recovered after averaging or pulse compression using binary (1-b) quantization only without experiencing significant losses. Hence, no analog-to-digital converter is required, which reduces the data throughput and makes the electronics faster, more compact, and energy efficient. All these are especially attractive for applications that require arrays with many channels and high sampling rates, where the sampling rate can be as high as the system clock. In this paper, the theory of binary quantization is reviewed, mainly from previous work on wireless sensor networks, and the signal-to-noise ratio (SNR) of the input signals under which binary quantization is of practical interest for ultrasound applications is investigated. The main findings are that in most practical cases binary quantization can be used with small errors when the input SNR is on the order of 8 dB or less. Moreover, the maximum SNR after binary quantization and averaging can be estimated as 10log10N-2 dB , where N is the number of averages.
Collapse
|
11
|
Xu K, Ta D, Cassereau D, Hu B, Wang W, Laugier P, Minonzio JG. Multichannel processing for dispersion curves extraction of ultrasonic axial-transmission signals: Comparisons and case studies. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2016; 140:1758. [PMID: 27914382 DOI: 10.1121/1.4962491] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Some pioneering studies have shown the clinical feasibility of long bones evaluation using ultrasonic guided waves. Such a strategy is typically designed to determine the dispersion information of the guided modes to infer the elastic and structural characteristics of cortical bone. However, there are still some challenges to extract multimode dispersion curves due to many practical limitations, e.g., high spectral density of modes, limited spectral resolution and poor signal-to-noise ratio. Recently, two representative signal processing methods have been proposed to improve the dispersion curves extraction. The first method is based on singular value decomposition (SVD) with advantages of multi-emitter and multi-receiver configuration for enhanced mode extraction; the second one uses linear Radon transform (LRT) with high-resolution imaging of the dispersion curves. To clarify the pros and cons, a face to face comparison was performed between the two methods. The results suggest that the LRT method is suitable to separate the guided modes at low frequency-thickness-product ( fh) range; for multimode signals in broadband fh range, the SVD-based method shows more robust performances for weak mode enhancement and noise filtering. Different methods are valuable to cover the entire fh range for processing ultrasonic axial transmission signals measured in long cortical bones.
Collapse
Affiliation(s)
- Kailiang Xu
- Department of Electronic Engineering, Fudan University, Handan Road No 220, 200433, Shanghai, China
| | - Dean Ta
- Department of Electronic Engineering, Fudan University, Handan Road No 220, 200433, Shanghai, China
| | - Didier Cassereau
- Sorbonne Universités, UPMC Univ Paris 06, CNRS, INSERM, Laboratoire d'Imagerie Biomédicale (LIB), 15 rue de l'école de médecine, 75006, Paris, France
| | - Bo Hu
- Department of Electronic Engineering, Fudan University, Handan Road No 220, 200433, Shanghai, China
| | - Weiqi Wang
- Department of Electronic Engineering, Fudan University, Handan Road No 220, 200433, Shanghai, China
| | - Pascal Laugier
- Sorbonne Universités, UPMC Univ Paris 06, CNRS, INSERM, Laboratoire d'Imagerie Biomédicale (LIB), 15 rue de l'école de médecine, 75006, Paris, France
| | - Jean-Gabriel Minonzio
- Sorbonne Universités, UPMC Univ Paris 06, CNRS, INSERM, Laboratoire d'Imagerie Biomédicale (LIB), 15 rue de l'école de médecine, 75006, Paris, France
| |
Collapse
|
12
|
Hua J, Lin J, Zeng L, Luo Z. Minimum variance imaging based on correlation analysis of Lamb wave signals. ULTRASONICS 2016; 70:107-122. [PMID: 27155349 DOI: 10.1016/j.ultras.2016.04.020] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/11/2015] [Revised: 03/12/2016] [Accepted: 04/24/2016] [Indexed: 06/05/2023]
Abstract
In Lamb wave imaging, MVDR (minimum variance distortionless response) is a promising approach for the detection and monitoring of large areas with sparse transducer network. Previous studies in MVDR use signal amplitude as the input damage feature, and the imaging performance is closely related to the evaluation accuracy of the scattering characteristic. However, scattering characteristic is highly dependent on damage parameters (e.g. type, orientation and size), which are unknown beforehand. The evaluation error can degrade imaging performance severely. In this study, a more reliable damage feature, LSCC (local signal correlation coefficient), is established to replace signal amplitude. In comparison with signal amplitude, one attractive feature of LSCC is its independence of damage parameters. Therefore, LSCC model in the transducer network could be accurately evaluated, the imaging performance is improved subsequently. Both theoretical analysis and experimental investigation are given to validate the effectiveness of the LSCC-based MVDR algorithm in improving imaging performance.
Collapse
Affiliation(s)
- Jiadong Hua
- State Key Laboratory of Manufacturing System Engineering, Xi'an Jiaotong University, Xi'an, Shannxi Province 710049, PR China
| | - Jing Lin
- State Key Laboratory of Manufacturing System Engineering, Xi'an Jiaotong University, Xi'an, Shannxi Province 710049, PR China; Shaanxi Key Laboratory of Mechanical Product Quality Assurance and Diagnostics, Xi'an Jiaotong University, Xi'an, Shannxi Province 710049, PR China.
| | - Liang Zeng
- State Key Laboratory of Manufacturing System Engineering, Xi'an Jiaotong University, Xi'an, Shannxi Province 710049, PR China
| | - Zhi Luo
- State Key Laboratory of Manufacturing System Engineering, Xi'an Jiaotong University, Xi'an, Shannxi Province 710049, PR China
| |
Collapse
|