Abstract

This paper presents a multiparameter water quality sensor chip with temperature compensation. The sensor chip was manufactured using microelectromechanical system (MEMS) technology. The surface of the chip integrated pH, dissolved oxygen (DO), ammonia nitrogen, and temperature sensors. To compensate for the solution temperature, the chip was also designed with a sandwich, plate-type serpentine Pt resistance heater. The experimental results showed that the pH sensor had a high sensitivity of 0.288 mA/pH with good linearity (), the sensitivity of the temperature sensor was 0.949 Ω/°C, the sensitivity of the ammonia nitrogen sensor was 0.1139 mA/ppm, the sensitivity of the dissolved oxygen (DO) sensor was 2.22 μA/ppm, and the sensitivity of the temperature changes with respect to the heater power was 0.3126°C/mW. Compared with a single water quality parameter sensor, the as-prepared sensor chip could simultaneously detect multiple parameters of a water sample and had a good temperature compensation effect. Moreover, the sensor chip was small in size, rugged, and highly accurate.

1. Introduction

Water is closely related to human life and is an important factor affecting human health. However, with the rapid development of industry and agriculture, water pollution has become increasingly serious. As water sources continue to deteriorate, the number of water-related diseases increases every year. Therefore, it is of great significance to develop a miniature sensor that can simultaneously detect multiple water quality parameters [1].

The pH is a critical evaluation parameter of aqueous solutions, which affects the solution properties and chemical reaction rate [2]. The World Health Organization recommends that, for residents, the pH of drinking water be maintained between 6.5 and 8.5 [3]. If the pH of drinking water exceeds the specified range, a few problems may arise. For example, if the pH is extremely low, water pipelines will corrode. If the pH value is extremely high, calcium and magnesium ions in the water will precipitate and form scale [4]. Traditional glass pH electrodes have the advantages of good stability, an ideal Nernst slope, good reproducibility, and long life. However, glass electrodes also have some disadvantages, such as being bulky and fragile. To overcome these shortcomings, several scientific research teams have conducted research and developed new pH sensors, for example, fiber-optic pH sensors [5] and ion-selective field-effect transistor (ISFET) pH sensors [6], which have received attention in recent years. Based on the characteristics of mechanical stability, potential metal oxide pH sensors have been extensively studied, such as Ta2O5 [7], WO3 [8], TiO2 [9], and CuO [10]. In recent years, pH sensors consisting of a conductive polyaniline (PANI) polymer have also been reported. In particular, PANI most often exists in the form of an intermediate oxidation state (emeraldine base (EB)). Through a reduction reaction, PANI can be converted from an intermediate state to a full benzene structure (leucoemeraldine base (LEB)). In contrast, through an oxidation reaction, PANI can also be converted from an intermediate state to a full quinone structure (pernigraniline base). As a conductive polymer, the conductivity and optical properties of PANI vary with pH. Therefore, based on these characteristics, PANI can be used to manufacture electrical and optical pH sensors [11].

Ammonia nitrogen (NH3-N) is also an important water quality parameter. Ammonia nitrogen (NH3-N) includes ammonium ions () and free ammonia (NH3) [12]. Such pollutants often appear in liquid organic waste such as sewage and manure [13]. Excessive ammonia nitrogen can cause eutrophication of the water environment. Therefore, to prevent water pollution, it is very important to monitor the ammonia nitrogen content. The method for determining ammonia nitrogen has been studied worldwide. The traditional method for determining the NH3-N concentration is to use a colorimetric method. However, the colorimetric test process is tedious, and it easily causes secondary pollution; thus, real-time monitoring cannot be achieved. Among the various methods for ammonia nitrogen measurement, electrochemical detection exhibits a low cost, high sensitivity, good selectivity, and the possibility of real-time monitoring. Therefore, various ammonium ion-selective electrodes have been developed [14]. However, due to the inherent sensing characteristics of ion-selective electrodes, the potential sensor has a low response at low target concentrations and is susceptible to interference from other ions [15]. To overcome the shortcomings of the potentiometric method, many laboratory models of amperometric ammonium ion sensors have been developed [16]. Compared with traditional colorimetric and potentiometry methods, amperometric sensors generally have higher sensitivity and lower detection limits while also being easy to implement in real-time monitoring.

Dissolved oxygen is also an important water pollution indicator. Dissolved oxygen is to molecular oxygen dissolved in water, which is denoted by DO. The concentration of DO is the number of milligrams of oxygen per liter of water, expressed in mg/L or ppm. Normally, the DO in water is close to saturation. When the oxygen in the water is consumed by organic or reducing substances, the DO concentration in the water will decrease. When the DO concentration is lower than a certain level (<4 mg/L), anaerobic bacteria in the water can multiply and give the water an odor. Therefore, the detection of DO is a necessary condition for judging the quality of water. The traditional detection method is an iodometric method. The iodometric method is an internationally recognized reference method for measuring DO in water [17]. Although the iodometric method is accurate, it requires chemical reagents. Furthermore, the steps are time consuming, and the system is bulky. Therefore, real-time monitoring is not possible [18, 19]. Electrochemical measurements are a commonly used method and include potentiometric methods [20, 21] and amperometric methods [22]. The electrochemical measurement method has the advantages of a simple operation and a fast and efficient response, which are convenient for real-time monitoring. Therefore, it has great research value.

On the other hand, water temperature is also an important factor that affects the pH, DO concentration, and ammonia nitrogen concentration [2326]. To measure each parameter accurately, the solution temperature needs to be accounted for. At present, there are many studies on microscale temperature sensors [27, 28]. The most commonly used thermistors are metal [27], polymer [28], etc. However, to integrate a temperature sensor with other water quality parameter sensors on the same silicon chip, thus, manufacturing an integrated multiparameter sensor chip has yet to be explored.

This paper proposes a multiparameter water quality sensing system that can simultaneously measure the following four parameters: pH, DO concentration, ammonia nitrogen concentration, and temperature of the water sample. The electrode array is made of microelectromechanical system (MEMS) technology and includes three three-electrode arrays and a resistive temperature sensor. PANI/CuO is used as the pH-sensitive material and the ammonia nitrogen-sensitive material. High-purity Pt resistance sensors are used for the temperature sensor. Considering the influence of temperature on various parameters, an integrated serpentine Pt resistance heater is used to compensate for temperature. Experimental results demonstrate that the chip exhibits good sensing performance. The chip is small and structurally stable; therefore, it is suitable for mass production as well as for use in remote and distributed measurements.

2. Materials and Methods

2.1. Reagents and Instruments

The reagents used in the experiment were of analytical grade. Reagents were purchased from Sinopharm Group Beijing Chemical Reagent Company. Deionized water was used for dilution and cleaning. Ethanol, aniline, oxalic acid, ammonium persulfate (APS), copper chloride (CuCl2), nanocopper powder, and ammonia (NH3•H2O) were used to prepare polymers, such as PANI, PANI/CuO, and PANI/Cu. Acetic acid (HOAc), phosphoric acid (H3PO4), and NaOH were used to prepare a pH buffer solution. Sodium sulfite (Na2SO3) particles and potassium chloride (KCl) were used to prepare a DO measurement solution. Ammonium chloride (NH4Cl) and potassium chloride were used to prepare an ammonia nitrogen measurement solution. Silver (Ag) powder, silver chloride (AgCl) powder, and binder (mixture of polyvinyl butyral, ethyl cellulose, ethyl acetate, and terpineol) were used to prepare silver/silver chloride (Ag/AgCl) paste.

The instruments used in the experiment are as follows: an electronic balance (Sartorius, Germany), an ultrapure water machine (Beijing Yingan Meicheng Scientific Instrument Co., Ltd.), a laboratory-developed potentiostat, a data acquisition system (PCI-1716/1716L, Advantech, China), a signal generator (DG1032, RIGOL, China), a peristaltic pump (BT100-2J, Dichuang, China), a magnetic pH meter (pHS-3E, Thunder magnetic, China), a thermal imaging analyzer (Fluke, PTi120, USA), beakers, and a silicone tube.

2.2. Microelectrode Array (MEA) Sensor Design

The sensor chip contains three three-electrode arrays (ammonia nitrogen sensor, pH sensor, and DO sensor) and a Pt resistance temperature sensor (Figure 1(a)). The three electrodes of the electrode array include a working electrode, a counter electrode (Au), and a reference electrode (Ag/AgCl) [29, 30]. The working electrodes of the three sensors are modified with different sensitive membranes. The working electrode of the DO sensor was unmodified and based on Au. The working electrode of the ammonia nitrogen sensor was modified using PANI/Cu [31], and the working electrode of the pH sensor was modified using PANI/CuO [32, 33]. Considering the influence of temperature on the measurement of various parameters, a flat sandwich sensor structure with a serpentine heater was prepared (Figure 1(b)) [34, 35].

When the electrochemical reaction occurs, if the solution temperature is excessively low, the power of the heating track can be adjusted according to the real-time monitored temperature, thus allowing for the surface of the electrochemical sensor to reach the ideal temperature, which ensures the stability and reliability of the various parameter measurements. To obtain an optimal heating effect, Pt was used for the heater design, and four types of heater structures were designed, as shown in Figure 2.

A finite element steady-state thermal analysis was conducted to investigate the heat transfer process of the four heater structures. Under the assumption that the ambient temperature was 0°C, the heating generation rate was considered to be . To simplify the simulation process, the changes in the thermal conductivity due to the variation in temperature were neglected, in addition to the thermal contact resistances among the sensor layers. Figure 3 presents the thermal simulation results of the four heater types, which were subject to the same boundary conditions and heating power.

Based on the simulation results, the first type of heater structure (Figure 3(a)) demonstrated the most significant heat transfer effect, and the temperature distribution was relatively uniform. In addition, it can cover the entire electrochemical reaction region. Therefore, in this study, the first heater structure was adopted.

2.3. Fabrication of the Sensor
2.3.1. Fabrication of the MEA

The detailed process of the sensor fabrication is shown in Figure 4. A single 520 μm crystal silicon wafer was employed as the substrate (Figure 4(a)). The substrate was ultrasonically cleaned with acetone and absolute ethanol. After cleaning, a 300-nm-thick SiO2 isolation layer was grown via a thermal oxidation method (Figure 4(b)). A serpentine heating electrode with a thickness of 200 nm was prepared by photolithography, and the target material was Pt (Figure 4(c)). An SiO2 insulating layer with a thickness of 500 nm was grown on the heating layer via a thermal oxidation method, and the insulating layer was reduced to a level between the metal layer and the insulating layer (200 nm). A heating electrode pad with a thickness of 500 nm was prepared via photolithography. Thereafter, a SiO2 layer with a thickness of 300 nm was grown via the oxidation method (Figure 4(d)). Using a photolithographic technique, a 200-nm-thick electrochemical sensor working electrode (Figure 4(e)) and a counter electrode (Figure 4(f)) were prepared, and the target material was Au. A reference electrode with a thickness of 200 nm was then prepared (Figure 4(g)), and the target material was Ag. Finally, a Pt temperature sensor with a thickness of 200 nm was prepared via photolithography (Figure 4(h)).

2.3.2. Fabrication of the Ag/AgCl Electrode

Silver powder, silver chloride powder, and binder were mixed according to a given ratio and then ground using a ball mill for 10 h to obtain an Ag/AgCl slurry. The obtained slurry was coated on the Ag electrodes of the three three-electrode arrays. Then, the slurry was cured and dried at 100°C for 1 h, which completed the preparation of the Ag/AgCl electrodes (Figure 4(i)).

2.3.3. Fabrication of the PANI/Cu Electrode

PANI/Cu was prepared using a slight modification of a method reported in the literature [36]. In particular, 11.4 g of APS was weighed and dissolved in 25 mL of water for later use. Thereafter, 100 mL of a 0.5 M oxalic acid solution was prepared. Next, 0.265 g of nanocopper was ultrasonically dispersed in 100 mL of 0.5 M oxalic acid and sonicated for 30 min. Then, 2.5 mL of aniline solution was prepared. The prepared APS, aniline, and oxalic acid solutions were mixed and stirred at room temperature for 6 h. A dark green suspension was obtained and then suction filtered. The PANI/Cu precipitate was collected on filter paper after suction filtration and then washed using deionized water and absolute ethanol. After washing, the PANI/Cu was placed in a vacuum oven and dried at 65°C for 24 h. The obtained PANI material was ground using a ball mill for 1 h to obtain a PANI/Cu powder. PANI/Cu powder and binder were mixed at a specific ratio, and the mixed material was evenly coated on the ammonia nitrogen working electrode. Finally, it was placed in a vacuum oven and dried at 80°C for 6 h (Figure 4(j)).

2.3.4. Fabrication of the PANI/CuO Electrode

Nano-CuO was mixed with APS solution, oxalic acid solution, and aniline solution, followed by the preparation of a PANI/CuO powder in accordance with the PANI/Cu synthesis method. PANI/CuO powder and binder were mixed in a specific ratio, and the mixed material was evenly coated on the pH sensor working electrode. Finally, it was placed in a vacuum oven and dried at 80°C for 6 h (Figure 4(j)).

2.4. Experimental Setup and Measurement
2.4.1. Test Container

Combined with polymethylmethacrylate (PMMA) microfluidic technology, a microfluidic test container was prepared. The test container, which was fabricated using plexiglass processing, is shown in Figure 5(c). Figure 5(a) is a schematic diagram of the internal structure of the test container. Figure 5(b) is a three-dimensional schematic diagram of the test container. The test container was adhered to the printed circuit board (PCB) of the sensor using insulating glue. The test container and sensor package are shown in Figure 5(d).

2.4.2. Experimental Device

The schematic of the entire experimental device is shown in Figure 6(a). A peristaltic pump transferred the test liquid into the test container at a speed of 50-60 rpm. After the test chamber was filled, the temperature, pH, DO concentration, and ammonia nitrogen concentration of the solution were successively measured. After the test was completed, the peristaltic pump was switched on, and the liquid in the test container was discharged at the same speed. The test circuit is shown in Figure 6(b).

3. Results and Discussion

3.1. Temperature Measurement

The microtemperature sensor manufactured in this paper is a Pt resistance temperature detector (RTD). Due to the need for accurate temperature measurements, the purity of the Pt resistance is 99.9999%. A water bath was used to heat the temperature sensor from 0°C to 50°C (the water temperature was monitored by a commercial temperature sensor (TES1310)), and a multimeter (17B+, Fluke) was used to measure the resistance of the Pt RTD. The regression analysis results of the Pt RTD resistance Rt and temperature T are shown in Figure 7(a). The sensitivity is 0.949 Ω/°C, and the linearity is 0.9994, which shows that the Pt RTD can be used to detect temperature. The stability measurement of the temperature sensor was performed for 10 consecutive days, and the measurement results are shown in Figure 7(b). The stability test shows that the temperature is relatively stable, and even though the output fluctuates slightly, the change rate is less than 1%.

3.2. DO Measurement

Regarding the measurements of the DO sensor, the three electrodes of the DO sensor (Figure 1) were connected to the electrodes at the corresponding positions of the potentiostat (Figure 6(b)). Sodium sulfite particles (Na2SO3) were slowly dissolved in a 0.1 M KCl solution to prepare six solutions with different DO concentrations (1.5 mg/L, 3 mg/L, 4.5 mg/L, 6 mg/L, 7.5 mg/L, and 9 mg/L). Scanning was performed using cyclic voltammetry (CV) at a scan rate of 0.1 V/s and a scan range of -1 V to 0 V. The measurement results were calibrated by colorimetry.

The measurement principle of the DO sensor is to reduce the oxygen dissolved in the solution through an electrochemical reaction and measure the current generated by it. The chemical reaction on the electrode surface can be expressed as [20]:

The reaction on the working electrode surface of the DO sensor causes a change in current/potential. The CV model of the DO sensor is described below [22]: where is the alternating current, is the number of electrons, is the overpotential, is the transfer coefficient, is the ratio of the Faraday constant to the ideal gas constant and temperature, and and are the diffusion-limited cathode current and anode current, respectively. Figure 8(a) shows the CV model of the sensor at different DO concentrations. The model describes two regions: a low-deviation kinetic control region (region 1) and a high-deviation diffusion control region (region 2). We used the current value at a potential of 0.8 V in region 2 as the concentration analysis current. Figure 8(b) shows that the current is proportional to the DO concentration in region 2. The sensitivity of the sensor is 2.22 μA/ppm, and the linearity is also relatively high (). The stability test of the DO sensor was performed for 10 consecutive days. The measurement temperature was . The measurement results are shown in Figure 8(c). The stability test shows that the DO concentration value is generally stable, and even though the output fluctuates slightly, the change rate is less than 4%.

3.3. pH Measurement

For the pH sensor, the three electrodes of the pH sensor (Figure 1) were connected to the interface of the potentiostat (Figure 6(b)). Moreover, 0.1 M H3PO4, 0.1 M HOAc, and 0.1 M NaOH solutions were mixed to prepare pH (, 5, 6, 7, 8, and 9) buffer solutions. Additionally, a commercial pH meter (pHS-3E, Thunder Magnetic, China) was used to monitor the pH of the buffer solution. Each sample was measured at a positive potential of 1 V using the ampere method, the duration of the measurement was 5 s, and the data acquisition system was used to monitor the variations in current.

The pH is detected by measuring the concentration of hydrogen ions. Among the PANI composites used in this study, the most common form of polyaniline is referred to as EB. Moreover, PANI is a remarkable pH sensing material. Under low pH conditions, there are more H+ ions present in the solution, and the EB of PANI is converted into a pernigraniline base (PN) via an oxidation process. When the pH is greater than 7, there is an excess of OH- ions in the sample. When these ions approach the nanocomposite, EB is converted to LEB through a reduction reaction [11, 37], as shown in Figure 9. Considering that polyaniline and copper oxide composites have a better sensing effect on pH [32, 33], this paper used PANI/CuO composites as the pH-sensitive material.

The pH sensing mechanism of CuO can be described by the following redox reactions [32]:

Scanning electron microscopy (SEM; Zeiss Supra 55) was used to analyze the morphology of the PANI and PANI/CuO composites. Figures 10(a) and 10(b) present the morphology of PANI, which reveals the porous structure of the macromolecular chain. Figures 10(c) and 10(d) present the morphology of the PANI/CuO composite. As seen from the figure, compared with those of PANI, the agglomeration effect of PANI/CuO was more significant, and the interparticle gap was larger.

Figure 11(a) presents the current response of the pH sensor at different pH values. The measurement conditions were . The fitted curve is shown in Figure 11(b), and the sensitivity of the sensor is 0.288 mA/pH. Based on the test results, the pH sensor demonstrates a high sensitivity and linearity (). The stability test of the pH sensor was performed for 10 consecutive days, and the measurement conditions were . The measurement results are shown in Figure 11(c). The stability test shows that the pH value is generally stable, and even though the output fluctuates slightly, the change rate is less than 6%.

3.4. Ammonia Nitrogen Measurement

A total of 382.14 mg of NH4Cl powder was dissolved in 1 L of water and configured as a 100 mg/L (calculated based on N content) ammonia nitrogen solution. A 0.1 M KCl solution was prepared as the supporting electrolyte. Different volumes of NH4Cl solution were dissolved in the KCl solution and diluted to different concentrations (4 mg/L, 3 mg/L, 2 mg/L, 1 mg/L, 0.5 mg/L, 0.25 mg/L) of ammonia nitrogen solution. The concentration of the solution was calibrated by the ammonia nitrogen standard color chart. During the test, the three electrodes of the ammonia nitrogen sensor (Figure 1) were connected to the corresponding positions of the potentiostat (Figure 6(b)). Each sample was measured at a positive potential of 1 V using the ampere method, the duration of the measurement was 3 s, and the data acquisition system was used to monitor the variations in current.

Considering the excellent conductivity of PANI, it can be used for the detection of ammonia nitrogen (NH3-N) concentration, but undoped PANI has a low sensitivity to the concentration of NH3-N [38]. Nanocopper-doped PANI leads to significant changes in polymer properties. The SEM images of PANI/Cu are shown in Figures 12(a) and 12(b). Compared with that of the undoped PANI, the agglomeration effect of nanocopper-doped PANI is more obvious.

Furthermore, the ammonia nitrogen in the solution mainly exists in the form of ammonium ions () and free ammonia (NH3). A modification of PANI/Cu on the electrode surface will lead to the formation of a complex between Cu2+ and free ammonia. This complex can be easily reduced to Cu+ and subsequently oxidized to Cu2+ [39]. When the cathode current increases, it indicates that the formation of the ammonia complex promotes the occurrence of redox reactions. The redox reaction of the ammonia complex is as follows:

Figure 13(a) shows the current response of the ammonia nitrogen sensor at different concentrations. The measurement conditions were . The fitted curve is shown in Figure 13(b), and the sensitivity of the sensor is 0.1139 mA/ppm. The test results show that the ammonia nitrogen sensor has high sensitivity and good linearity (). The stability test of the ammonia nitrogen sensor was performed for 10 consecutive days. The measurement temperature was , and the measurement results are shown in Figure 13(c). The stability test shows that the ammonia nitrogen concentration is relatively stable, and even though the output fluctuates slightly, the change rate is less than 5%.

3.5. Temperature Compensation

A thermal analysis was conducted to investigate the heater using a thermal image analyzer (C3, FlIR). The initial sensor temperature was 22°C before heating. Different power levels were applied to the heater to obtain the temperature after heating. The temperature variation can be expressed as follows:

Table 1 summarizes the relationship between heating power and temperature. Figure 14 shows a linear relationship between the heating power and temperature variation.

From Figure 14, the linear relationship between the heating power and the temperature variation can be expressed as follows:

According to the related literature [2326], the pH, ammonia nitrogen concentration, and DO concentration in an aqueous solution are affected by the solution temperature. Therefore, keeping the solution temperature constant (or within a small temperature range) is of great significance to improve the measurement accuracy. When the measured ambient temperature is low, the heater can be applied to reach 25-28°C according to the temperature value measured by the temperature sensor. When the solution temperature is high, external circulation of cold water can be applied to cool the entire sensor device. According to the temperature value monitored in real time, the measurement will be performed when the temperature is within 25-28°C. Therefore, the influence of the solution temperature on the measurement accuracy is reduced. Based on this feature, the sensor detection system is suitable for real-time measurement of remote waters.

3.6. Monitoring Real Water Samples

We used our developed multiparameter sensor to measure the parameters of tap water and Songhua River water. For pH detection, commercial glass electrodes (pHS-3E, Thunder magnetic, China) were used as a reference. For temperature detection, a commercial temperature sensor was used as a reference. Regarding the detection of the DO and ammonia nitrogen concentrations, a colorimetric test box was used as a reference. The values of the temperature, pH, DO concentration, and ammonia nitrogen concentration from the multifunctional sensor were compared with the results from the reference methods, and the results are shown in Table 2. From Table 2, the results measured by our system are close to those measured by the reference method.

4. Conclusions

We proposed an integrated multiparameter sensor chip for monitoring the following four parameters of water samples: temperature, pH, DO concentration, and ammonia nitrogen concentration. To achieve temperature compensation, a serpentine heater was designed on the chip to ensure the ideal temperature for sensor measurements. The sensor chip size was approximately . The sensitivity of the temperature sensor was 0.949 Ω/°C, and the sensitivity of the DO sensor was 2.22 μA/ppm. Additionally, PANI/CuO was modified on the electrode surface as a pH-sensitive material, and the sensor showed a good sensitivity of 0.288 mA/pH. PANI/Cu was modified on the electrode as an ammonia nitrogen-sensitive material and showed a good sensitivity of 0.1139 mA/ppm. The power and temperature sensitivity of the heater was 0.3126°C/mW. The advantages of this sensor chip and detection system are as follows: (i)The following four parameters of a solution can be measured simultaneously: temperature, pH, DO concentration, and ammonia nitrogen concentration(ii)PANI/Cu as a sensitive material for ammonia nitrogen and PANI/CuO as a pH-sensitive material both show good sensitivity(iii)The sensor integrated multiple identical three-electrode arrays, which can be used as a chip laboratory to compare and measure a single parameter and find the best measurement solution(iv)The sensor array can be modified by different sensitive membranes, which allows it to measure other water quality parameters such as nitrate, heavy metal ions, and chloride ions(v)The sensor has a heating element, which can heat the solution to be measured to an ideal measurement temperature, thus improving measurement accuracy(vi)The sensor detection system is small in size, economical, structurally stable, and durable. Moreover, a microperistaltic pump with a control valve can realize automatic control of the detection system

In this paper, the corresponding electrodes were modified by using an adhesive coating method. This method may cause the thickness of the modified sensitive film to be slightly different, which will slightly affect the measurement results. Therefore, the next research direction will use other methods to modify the electrode to ensure that the thickness of the modified sensitive film is consistent. In summary, the multiparameter sensor chip studied in this paper was suitable for drinking water detection. The sensor chip was small and rugged and suitable for portable or online water quality detection systems.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This study has been supported by a Project of Defense Industrial Technology Development Program No. JCKY2017412C003.