Abstract
Although the quartz accelerometer has made great advances, the performance, in some specific applications such as tilt sensing, needs to be well compensated in high temperature environment. Based on the high temperature testing of low cost quartz accelerometers, we found that the normalized positive and negative parts are asymmetrical at high temperature and the temperature curve of zero sensor output is related to the roll angle of the sensor. Traditional temperature compensation method only considers the temperature factor and ignores the roll sensitivity, which leads to deteriorated accuracy. To solve this problem, this paper proposes a novel and simple mathematical model to obtain a more accurate expression of zero sensor output, which makes the sensor output more robust at high temperature. Experimental results on two low cost quartz accelerometers demonstrate that the proposed model is feasible and effective, which could reduce the temperature drift error of the sensor output typically from 0.01 g to 0.001 g. Furthermore, we introduce the compensated sensors in the three-axis inclinometer system for tilt sensing, and the evaluation results show that the temperature drift error of the inclination in the range (, ) is reduced typically from to compared to the traditional method.
1. Introduction
Quartz accelerometers have become popular in many applications like motion control [1] and inertial measurement units [2] due to their simple construction and high accuracy. However, along with the increasing requirement for excellent performance of sensors, especially in harsh environment, for example, high temperature, high shock, and high vibration, some recalibration or recompensation methods should be adopted to obtain a satisfying accuracy depending on the specific usage [3, 4]. Generally, the accuracy of the accelerometer can be influenced by manufacturing technique, internal structure, and surrounding environment in which the environmental temperature makes the performance much worse. A directed compensation method is to add a temperature control system [5], which could make the sensor work under the controlled temperature. However, the main drawbacks of the temperature control system include long stabilization time, large power, and increasing size, which cannot meet the actual demand in many applications. A commonly used alternative solution is to construct a mathematical model [6], which first models how the performance changes with temperature and then, according to the measured temperature of the sensor, corrects the sensor output. Apparently, such mathematical models are more suitable for most applications due to their simplicity and practicality.
Estimation methods need to be settled for the accelerometer in the thermal calibration in which iterative methods are mostly utilized to solve the calibration parameters (e.g., scale factors, misalignments, and biases of the accelerometer triad) and to achieve high estimation accuracy. Qian et al. [7] used least squares method to estimate the parameters of the proposed linear model. Yang et al. [8] presented an improved iterative nonlinear calibration method by using least squares method and sequence quadric program method. Ang et al. [9] proposed a nonlinear regression model to reduce the deterministic errors associated with scale factor, bias, and misalignment of the dual axis accelerometer, and the performance was verified by tilt and motion sensing. Won and Golnaraghi [10] proposed a mathematical model of six calibration parameters and use an iterative method to estimate the parameters of the nonlinear model. The objective of the proposed iterative method is to make the calculated gain factors and biases of each axis match the corresponding true values. Because the need of an initial rough estimate makes this algorithm inconvenient, Zhang et al. [11] proposed an improved multiposition calibration for solving the unknown parameters without any initial guess. However, the accelerometer temperature drift is still left as an unresolved problem. As a consequence, many trials have been made to establish thermal models of calibration parameters to attack this problem. Aggarwal et al. [12] explored the effects of thermal variations on biases and scale factors at different temperature values through the thermal chamber and then proposed three-order polynomial thermal models for ADI microelectromechanical system sensors. To investigate the thermal property in varying temperature conditions, Aggarwal et al. [13] considered the thermal ramp experiment from which a simple polynomial temperature model is developed for the inertial sensor biases and scale factors. After compensating the thermal errors, the inertial navigation solution was significantly improved. These two trials have to utilize a turntable with a temperature-controlled incubator that requires the precise orientation information. To improve the robustness to the turntable error, in [14], an indirect calibration technique for estimating the body-frame drift induced by the variation of the accelerometer sensitivity axes due to temperature changes is proposed. However, this estimation method depends on and is limited to the accuracy of horizontal accelerometer measurements. Moreover, Zhang et al. [15] utilized a nonvertical rotation axis observation method to attack the problem of the relationship between the gyroscope triad and the accelerometer triad in constant thermal conditions.
In this paper, we mainly focus on the usage of the quartz accelerometer in tilt sensing systems [16], such as Measurement While Drilling (MWD) and Logging While Drilling (LWD), which are used to monitor and guide the down borehole directional drilling in oil and gas exploration. Two main involved angles, inclination and roll angle, are illustrated in Figure 1. Suppose that an accelerometer is mounted and its sensitive axes are arranged along the -axis in the coordinate system, inclination is the angle that -axis makes with the down direction . It is 0° when -axis is down and 90° when -axis is horizontal. Roll angle is defined as the angle of counterclockwise rotation about the -axis (looking in the positive -axis direction) in the gravitational field. Particularly, the inclination plays a very important role in tilt sensing systems, and the precision of inclination directly determines the performance of the whole system, so the accuracy of the sensor output is the key for accurate measurement at high temperature. In this paper, we pay our attention to the temperature compensation of low cost quartz accelerometers, and the main contributions include proposing a mathematical temperature compensation model for zero sensor output, especially, when the temperature curve changes with the roll angle, presenting a feasible and effective solution for the whole compensation of quartz accelerometer, and applying the proposed temperature compensation method in the tilt sensing system.
This paper is organized as follows: Section 2 contains the problem statement derived from our collected sensor data; Section 3 gives the description of the proposed temperature compensation method; Section 4 illustrates the experimental results and introduces the proposed model into the application of tilt sensing. Discussion and conclusions are drawn in Section 5.
2. Problem Statement
Although the performance of the quartz accelerometer is quite good to satisfy many applications at room temperature or under normal temperature, the accuracy may become poor and even unacceptable for some specific usage in high temperature environment. Generally speaking, temperature compensation is needed in those cases. However, we found that some characteristics of the sensor (e.g., asymmetry and roll sensitivity) may make the commonly used temperature compensation methods less effective.
To research how the performance changes with temperature, we have analyzed some sensor data from several quartz accelerometers. We use the monoaxial quartz accelerometer JSD-II/B-MJ, provided by the manufacturer (http://langfang037637.11467.com/), to constitute the hardware platform. Since the datasheet of this kind of accelerometer is not available online now, some of the most significant parameters are reported in Table 1 (although the provided measurement range is , in our experiments, we only care about the static range of since the derived calibration model is only valid for this range, so the designed circuit only transfers the corresponding signals in the range ). To compensate the sensor’s temperature performance, the repeatability is the most important point. Fortunately, this basic assumption is still tenable in our observations. Meanwhile, we found that the data from one sensor are identical and repeatable, and the changing trend of the data is also consistent for different sensors. In our experiments, the maximum and minimum of the sensor output both change consistently with temperature when rotating the sensor to different roll angles, but the value of the zero sensor output changes differently with temperature at different roll angles (from Figure 1, when , the sensor output is close to zero; when , the sensor output is close to maximum; when , the sensor output is close to minimum. So we call them zero, maximum, and minimum sensor output, resp. To distinguish the term “bias,” we use “zero sensor output” to express the approximate value when ). Note that, in the thermal chamber, only a mechnical fixture is used for fixing the cylinder structure in which the sensor is mounted, so it is hard to make of the accelerometer accurately . We assume that the trend of the temperature curve at near inclination is similar to that at precise inclination. For simplicity, when the inclination of the sensor is close to , we also call it zero sensor output.
Some representative illustrations are shown in Figure 2. For instance, when the roll angle is and the inclination is near , the sensor output of sensor #1 changes with the rising temperature while keeping the sensor static, and the data points are plotted with red squares in Figure 2(a). However, when we turn over the sensor to make the roll angle and the inclination is still near , the change of the sensor output of sensor #1 with the rising temperature becomes different, which is shown by green circles in Figure 2(a). The same situation occurs when the roll angles are 0° and 180° for sensor #2, and the corresponding data are depicted in Figure 2(b). In the experiments, roll angles are roughly read out by a simple scaled circle, which is attached to the cylinder structure, since we found that the accuracy of roll angles has less influence on the accuracy of our model. To convert floating-point arithmetic to faster integer arithmetic in numerical calculation, the analog output voltage is first converted to digital signal by an A/D (analog to digital) converter and then transformed by the formula , where is the digital output voltage (unit: V). In all our experiments, we use the digitalized as the sensor output. The phenomenon implies that the roll angle is relevant to the change of zero sensor output with rising temperature. All external environmental settings are identical except the change of the roll angle. It is worth noting that the phenomenon in Figure 2 should not be caused by random noises, for the differences of the zero sensor output at different roll angles could not be ignored and the phenomenon is regular and repeatable. In this paper, we only discuss the observed situation and find out how to compensate the data drift caused by temperature. Therefore, we adopt a data-driven method, by assuming that all possible errors, including the misalignment, temperature drift, and some uncertain errors caused by the operating environment, are compensated together. From Figure 2, since zero sensor output is approximately linear with temperature at any roll angle, for simplicity, we can fit it using a linear function. As a result, we give an assumption as follows.
(a) Sensor #1
(b) Sensor #2
Assumption 1. The zero sensor output is approximately linear with temperature at any roll angle, but the slope may change with roll angles.
Moreover, we further investigate the relationship between the sensor output and three possible factors: roll angle, temperature, and inclination. For an ideal situation, the sensor output should satisfy . But it is impossible because our purpose is to obtain an accurate inclination by using the sensor output. Since temperature curve of the sensor output is related to the roll angle, we try to utilize the roll angle to compensate the temperature drift of the sensor output. Therefore, when the inclination is fixed in the thermal chamber, we record the sensor output of the tested accelerometer at every roll angle with 10° interval in the range at three stable temperature points (to avoid the temperature hysteresis, all experimental data are recorded after staying at least 20 minutes at an expected temperature point). As a result, as shown in Figure 3(a), at a fixed inclination and a fixed temperature, the relation between the sensor output and roll angle is an approximate sine function because of the misalignment of the rotation axis and the sensitivity axis of the accelerometer. Hence the plots in Figure 3(a) are the projection of the rotation into the and plane and is hence describing a sine curve with an offset. When the temperature changes, for example, from 15°C to 52°C and then 92°C, the sine curves become different. As a result, as illustrated in Figure 3(b), the difference of the sensor output at different temperature can be also approximated by a sine function of roll angle. The magenta squares depict the differences of sensor output between temperature points 52°C and 15°C, and the black crosses stand for that between 92°C and 52°C. Furthermore, we found that the amplitude of the sine curve in Figure 3(b) is related to the inclination, and the function could be also an approximate sine curve. Due to the limited experiment condition, we cannot give out a lot of points at different inclination angles to fit the sine curve. We tested some specific positions to give a rough estimation. For instance, the maximum amplitude of the sine curve at inclination is 37, which can be expressed as , while the maximum amplitude of the sine curve at inclination is about 27, which can be represented as . In fact, as it will be mentioned in Section 3.2, we only care about the zero sensor output, so this assumption is not necessary for our proposed model. However, we could provide the assumption in case that some possible applications need it. And also it makes 90˚ inclination a specific example and easy to extend. In our model, it is beneficial to know that the maximum amplitude of the sine curve in Figure 3(b) is also an approximate sine function to inclination when the inclination is . Therefore, we give another assumption below.
(a) Sensor output versus roll angle
(b) Difference of sensor output versus roll angle
Assumption 2. The difference of sensor output under different temperature could be approximate sine functions to roll angle, and the amplitude of the sine curve is also an approximate sine function to inclination.
3. Temperature Compensation Method
Inspired by abovementioned observations, we propose a novel temperature compensation procedure, which can simultaneously reduce the negative effect brought by temperature and roll angle. The proposed method is composed of the following four steps.
3.1. Modeling Temperature Curves of Maximum and Minimum
The bias and scale factor may suffer deteriorated performance when the accelerometer is used in high temperature and high accuracy required measurement system, probably because of different nonlinear temperature drift of the complex peripheral interface circuit. Therefore, for better expressing the temperature drift, we adopt a polynomial model here to obtain the temperature curve of maximum and minimum output of the sensor, which can be written aswhere is the maximum or minimum output, is the corresponding nonlinear coefficient, and is the temperature. is empirically set to 4 in our system due to the trade-off between the effective repeatability of the maximum and minimum output and overfitting. In fact, the sensor is not sensitive at the position of the maximum or minimum output, so it is easy to find the proper position for maximum or minimum output in the actual operation (for convenience, the compensated accelerometer is mounted into a cylinder structure, and the direction of the maximum sensor output is the same with the gravity).
3.2. Modeling Temperature Curve of Zero Sensor Output
The key step, also the significant difference of our proposed procedure, is to obtain the temperature curve of the zero sensor output. Due to the manufacturing technique or internal structure, the temperature curve of the zero sensor output changes with the roll angle at a fixed inclination, which makes the usual temperature compensation method less effective. The direct result is that the normalized positive and negative parts are asymmetrical at high temperature. Formally, that is, , where is the actual output of the sensor and , , and are maximum, minimum, and zero sensor output, respectively. To handle this problem, we need to obtain a more accurate temperature curve of zero sensor output.
Assume that the function between the zero sensor output and the temperature is linear (Assumption 1), which could be defined aswhere is the slope and is the offset. Suppose that two points and can determine this linear function, in which is a point at room temperature , so could be calculated by and could be . According to Assumption 2, if we set , then could satisfy the functionwhere is the amplitude, is the phase, and is the offset. In addition, since we only need to obtain the temperature curve of zero sensor output, the inclination is always about ; that is, . Here we assume that , and we experimentally found that this assumption has no obvious influence on the whole compensation procedure. is set to 1 in our model, and it is reasonable to consider it as a complete sine period when the sensor turns over from roll angle to . The experimental data in Figure 3 could also prove that. Therefore, (3) could be simplified asTherefore, as long as two points and are obtained, then and . Substitute , , and (4) into (2); then temperature compensation function of zero sensor output could be represented by
3.3. Solving Model Parameters
Since we have obtained three expressions of maximum, minimum, and zero sensor output changing with temperature, we need to get all involved parameters, that is, in (1) and , , and in (4). The temperature is required from a temperature sensor, and the roll angle could be controlled manually in laboratory tests or approximately calculated by another two accelerometers in the three-axis tilt sensing system. For simplicity, we adopt the nonlinear least squares regression method [17] to optimize these parameters by collecting some necessary data. To obtain the coefficient estimates, the least squares method minimizes the summed square of residuals. The residual for the th data point is defined as the difference between the observed response value and the fitted response value and is identified as the error associated with the data; that is, . The summed square of residuals is given bywhere is the number of data points included in the fit and is the sum of squares error estimate. For example, (5) can be rewritten as , where is an -by-1 vector of responses, is a function of and , is an -by-1 vector of coefficients, is the -by- design matrix for the model, is an -by-1 vector of errors, is the number of the parameter, and is the number of data points. We choose trust-region fitting algorithm [17] in the iterative procedure until the fit reaches the specified convergence criteria.
3.4. Normalizing Sensor Output
Since the output of accelerometer is often considered as symmetrical even at high temperature, that is, , classical normalization method [18] is often defined aswhere is the normalized sensor output. As mentioned above, due to the possible asymmetrical performance of low cost quartz accelerometer, classical normalization method may be not accurate, so that the usage is limited. Assuming that we have obtained a feasible temperature curve of zero sensor output, then the normalization of our model could beSince in static condition should be equal to 1 g and should be equal to −1 g, the unit of the normalized sensor output could be . Note that our proposed normalization is in fact an approximation to a nonlinear characteristic of the accelerometer, which may not be strictly realistic. However, it is reasonable to make such an approximation based on the observations mentioned above which is also proved to be effective by the following experimental results.
4. Experimental Results and Validation in Tilt Sensing
In order to test the proposed temperature compensation method, we first make temperature compensation to two low cost quartz accelerometers and then introduce the model into the three-axis tilt sensing system.
4.1. Experimental Results in Sensors Compensation
First we place each sensor (mounted into a cylinder structure) at the maximum and minimum positions in the thermal chamber, respectively, and collect the sensor output when temperature rises from room temperature to the upper limit temperature, that is, , and then these recorded data are used to solve the parameters in (1) by using the least squares regression method. We also collect some data points from two quartz accelerometers to obtain the parameters in (5), which is shown in Table 2. For simplicity, the data points should be as few as possible, because less data collection can simplify the production process. Besides, the data in our experiments are often recorded at high temperature, so less data points could also largely reduce the heating time, which is helpful for the experiments or production. In our experiments, we found that the readouts at 4 roll angles and 2 temperature points are enough to solve the parameters in (5). As shown in Table 2, we choose , , , and as 4 tested roll angles and a room temperature point (about ) and a high temperature point (close to ) as 2 tested temperature points. It is worth noting that any 4 roll angles increasing by about are acceptable for training the model. In addition, the accuracy of roll angles and temperatures are unnecessarily too high, because we found that these factors do not affect the result too much. In our actual operation, we first record 4 output values with 4 different roll angles at room temperature before heating, then heat up the accelerometers to the upper limit temperature under the heating rate per minute, and rotate the sensor to record another 4 values with corresponding 4 roll angles when the temperature is almost stable (after staying at least 20 minutes at ). The heating rate is set empirically by considering the operating environment of MWD or LWD systems. Besides, since the accelerometer cannot be very accurately mounted along the -axis of the cylinder structure, the data points at room temperature are not identical, as shown in Table 2, which means that the inclination of the accelerometer is not accurately while recording the data points. However, assisted by a standard measurement apparatus, which can provide the indication of the inclination with an accuracy of about 0.01, we found that at room temperature the zero sensor values at different roll angles are nearly equal. It means that the normalized positive and negative parts are symmetrical at room temperature, so the first point in (5) can be obtained by , where and are the maximum and minimum at room temperature , respectively. In this experiment, the first point is (19, 16386) for sensor #1 and (21, 16373) for sensor #2.
Utilizing the collected data, we solve the model parameters by the least squares regression method. As a consequence, the model parameters (, , ) of sensors #1 and #2 are (38.55, 0.8778, −33.75) and (36.88, 0.7726, −5.25), respectively. As shown in Figure 4, the proposed model could be fitted well. Since we have simplified the model, it is easy to determine the whole sine curve by using only 4 data points.
(a) Sensor #1
(b) Sensor #2
Once we have obtained the parameters in (1) and (5), we can get the temperature curves of the maximum, minimum, and zero sensor output, respectively, and then the normalized sensor output could be calculated. The comparison of the normalized sensor output by using (7) and (8) is illustrated in Figure 5, in which “classical” and “proposed” represent the results obtained by using (7) and (8), respectively. We randomly place the sensor at different roll angles where the inclination is close to . Figure 5 shows some results when the roll angles are about , , and for both sensors #1 and #2. We can see that our normalized sensor output changes much less than the classical method when the temperature rises from room temperature to even . Comparing with the classical method, the proposed temperature method could reduce the temperature drift error typically from 0.01 g (at least 0.005 g) to 0.001 g ( is the unit of the normalized sensor output obtained by (8)), which demonstrates that our method is more stable and less sensitive to high temperature environment. The compensated accelerometers could be easily adopted in tilt sensing systems in high temperature environment.
(a) Sensor #1. The roll angles are about 30°, 180°, and 330° and the inclination is close to 90°
(b) Sensor #2. The roll angles are about 30°, 180°, and 330° and the inclination is close to 90°
4.2. Experimental Validation in Tilt Sensing
To further verify the feasibility of the proposed temperature compensation method, we apply it to tilt sensing at high temperature. We comply with the basic assumption that an inertial sensor, in static condition, is subjected only to the gravity force. As a consequence, the module of the acceleration vector measured by the sensor has to be equal to , independent of the sensor orientation [20]. As shown in Figure 6, in a 3D Cartesian coordinate system, has the same direction with the gravity , and plane is perpendicular to . In a cylindrical inclinometer instrument, three quartz accelerometers are mounted along three axes of coordinate system. Due to the limitation of the installation process, sensitive axes of accelerometers are often not accurately identical with the corresponding direction. This results in the misalignment. Usually, the misalignment should be calibrated to obtain a satisfying accuracy. Since in this paper we do not focus on the misalignment calibration algorithm, for simplicity, we use the model in [19] to obtain the final inclination. After the temperature compensation procedure, we denote the normalized actual accelerometer output as (temperature compensation and normalization steps are introduced in Section 3) and the acceleration vector in the ideal coordinate system as .
The mathematical model describes the accelerometer output in matrix form aswhere is the sensitivity matrix and is the offset. The diagonal elements of represent the scale factors along the three axes, whereas the other elements of are called cross-axis factors. These terms allow describing the axes misalignment, the cross-axis sensitivity, and the mechanical imperfection effect between different channels caused by the sensor electronics. Finally, the inclination angle could be computed by [20, 21]
We test the proposed temperature compensation method on two inclinometer devices, which are noted as devices #1 and #2. A two-step procedure is adopted in our experiments, in which the temperature performance of each sensor is first compensated, and then a calibration method for three accelerometers is used for tilt sensing. We first make temperature compensation for each quartz accelerometer as mentioned in Section 3 and then place them in each inclinometer device as their own -, -, and -axes. After obtaining the model parameters of each accelerometer, each compensated sensor output of , , and could be calculated, respectively, where the roll angle of each sensor could be estimated by another two axes; for example, for -axis, the roll angle could be . To evaluate the effect of our proposed temperature compensation method, we fix each device into a thermal chamber with different inclination angles (here we randomly choose , , and for evaluation) and then heat it up from room temperature to (staying at least 20 minutes to keep the temperature at stable ) under the heating rate of per minute. The heating rate of is set empirically by considering the operating environment of MWD or LWD systems. We record all inclination angles with and without our proposed temperature compensation method during the heating process, and the experimental results are illustrated in Figure 7. Inclination angles (-axes) change as the temperature (-axes) rises from the room temperature to the upper limit temperature of . Similarly, “classical” and “proposed” represent the results obtained by using classical normalized method (see (7)) and ours (see (8)), respectively. To distinguish the influence brought by temperature and other factors (e.g., scale factor, bias, and misalignment) in tilt sensing, both “classical” and “proposed” methods calibrate all parameters by the model in [19] except the temperature compensation coefficients, so that the differences of evaluation results are only caused by different temperature compensation methods. Typically the temperature drift error of the inclination calculated by our method is only about , which is more robust than the drift error computed by using the classical temperature method. Furthermore, according to (10), when the inclination angle is near , plays the most important role and and have little or even no effect on the tilt angle. In fact, we found in our experiments that at least up to we can only compensate by our new method and solve and using the classical compensation method for simplicity, and the results are very close.
(a) Device #1. The three-axis inclinometer is placed with 90°, 50°, and 30° inclination with arbitrary roll angles
(b) Device #2. The three-axis inclinometer is placed with 90°, 50°, and 30° inclination with arbitrary roll angles
5. Conclusions
This paper proposes a mathematical temperature compensation procedure for low cost quartz accelerometer based on the observation that zero sensor output changes with both roll angle and temperature. By modeling maximum, minimum, and zero sensor output, respectively, we solve the problem that the normalized positive and negative parts of sensor output under high temperature environment are asymmetrical. To make it easy to operate, we only choose several representative data points (responding to 4 roll angles at low and high temperature and a common point at room temperature) to fit the zero sensor output curve. Although the model is not very accurate because some approximations and assumptions are adopted, the experimental results demonstrate that it can distinctly reduce the temperature drift error of the sensor output. Moreover, we apply the proposed temperature compensation method in the tilt sensing system, which shows the actual efficiency at high temperature environment.
In addition, some points should be noticed. First, we ignore the long-term drifts of the tested accelerometers in our experiments because the duration of the experiments is short enough. Second, 8 accelerometers are tested in our experiments (2 in Section 4.1 and 6 in Section 4.2), and all the tested sensors have similar temperature characteristics. By using our proposed temperature compensation model, the temperature performances of all the tested sensors are greatly improved. Third, low cost quartz accelerometers are usually several times cheaper than well-produced ones; however, by using our proposed temperature compensation model, they can be also adopted for accurate tilt sensing even in harsh temperature environment. Fourth, it is important to obtain the data points accurately, for example, and , for a small mistake in those readings can cause a significant change of the unknown parameters. Fortunately these values are quite easy to obtain in the process and should not be mistaken easily, because the temperature readout and the sensor output are the base of the whole measurement system. Meanwhile, if the data points are wrongly collected, a simple verification at high temperature would help. Finally, this paper gives some experimental data and proposes an effective compensation model to improve the temperature performance of the low cost quartz accelerometer, but, limited to the existing test condition, it is hard to compare with many types of low cost accelerometers to build a more general compensation model. The proposed mathematical model is a data-driven method essentially, which tries to compensate all possible errors (misalignment, temperature drift, etc.) of low cost accelerometers. Therefore, low cost quartz accelerometers with these similar characteristics could use our proposed model to measure tilt in static situation. Future work may lie in comparing with different types of accelerometers and proposing a more general model.
Competing Interests
The authors declare that there are no competing interests regarding the publication of this paper.
Acknowledgments
This research was supported in part by National Natural Science Foundation of China (61472053, 91420102), China Postdoctoral Science Foundation (2014M550456), Chongqing Postdoctoral Special Funding Project (Xm2014087), and Basic Science and Advanced Technology Research Project of Chongqing (cstc2015jcyjA40041).