Design of a Marine Sediments Resistivity Measurement System Based on a Circular Permutation Electrode

: Marine sediments are rich in mineral resources, organic resources, and microbial life. The study of marine sediments is of great signiﬁcance for the development and utilization of marine resources and understanding the life process. Resistivity is the overall characteristic of the electrical conductivity of marine sediments. Measuring the resistivity of marine sediments is helpful to ascertain the marine geological structure, study the distribution of marine mineral resources, and evaluate the marine soil environment. Measuring the resistivity of marine sediments is of great signiﬁcance to promote marine exploration. At present, the resistivity measurement device on the market can be directly used to measure soil and water on land, but if used to measure marine sediments, it will be disturbed by seawater temperature and pressure, resulting in large errors. In this paper, a high-precision pressure-maintaining transfer system of marine sediment resistivity measurement instrument based on circular permutation electrode is designed, which can measure the resistivity of marine sediment samples after pressure-maintaining transfer. At the same time, a new type of circular permutation electrode measurement method is proposed, which makes the resistivity value more accurate, reduces the length of the probe appropriately, and saves the cost. By measuring the resistivity of marine sediments, the type of sediments can be inverted, which provides a way of thinking about the promotion of the research and development and utilization of marine resources. Author Contributions: Conceptualization, P.Z., J.C. and D.R.; Methodology, P.Z., D.R., X.P.; software, X.W.; validation, D.R., X.P. and Q.G.; formal analysis, Z.R.; investigation, X.W.; resources, Z.R.; data curation, D.R.; writing—original draft preparation, P.Z., D.R. and X.P.; writing—review and editing, J.C.; visualization, Q.G.; supervision, J.C.; project administration,


Introduction
Marine sediments [1] are composed of insoluble materials (rocks, soil particles, etc.) transported from land to sea by wind, ice, and rivers. They also contain chemical sediments in seawater, products of submarine volcanoes, and debris of marine organisms. Marine sediments contain abundant mineral resources and organic resources, including offshore oilfields [2], sulfur, phosphate rocks, polymetallic crust, and nodules. Marine sediments are also rich in microbial life, which promotes global biomass to a large extent and is also a key component of the earth system. The study of marine sediments is of great significance in finding marine minerals, developing and utilizing marine microbial resources, studying paleoclimate evolution, understanding the formation and evolution of the ocean, and understanding the evolution and origin of life.
Resistivity is the overall characteristic of the electrical conductivity of marine sediments. Measuring the resistivity of marine sediments is of great significance for exploring marine geological structures [3], studying the distribution of marine mineral resources, and evaluating the corrosiveness of a marine soil environment. The resistivity characteristics of marine sediments are mainly affected by porosity [4], water content, hydrate saturation and sediment types, and the properties of sediments can be inverted by resistivity values. At present, the equipment for measuring the resistivity of marine sediments has not been 2 of 21 used for commercialization and is still in its infancy. Some resistivity measurement devices are mainly used for land, such as soil and water measurement. If used to measure marine sediments, due to the interference of seawater temperature and pressure, there will be large errors. Therefore, it is not only of great scientific significance, but also of important application value to develop a marine sediment resistivity measuring instrument suitable for deep sea.
The commonly used methods for measuring the resistivity of marine sediments are mainly drag method and penetration method. Jackson [5] first developed a device to measure the resistivity of marine surface sediments. The electrode was installed on the plate, and the resistivity of marine sediments was directly measured by the dragging method, but the measurement depth was only 0.5 m. Dale F. Rucker et al. [6] used the drag method to investigate the water resistivity of sediments in the Panama Canal. In order to adapt to the required depth of investigation, the electrode spacing along the cable should reach 15 m and the total length of the cable should reach 170 m. The penetration method mainly includes the ring electrode measurement method and the point electrode measurement method. Won [7] designed a four-ring Wenner electrode resistivity measurement probe. A Rosenberger [8] designed a four-point electrode resistivity measurement probe. According to the principle of ring electrode resistivity measurement, Fu et al. [9] in 2015 designed a ring copper electrode probe with as many as 94 components of equidistant (1.0 cm apart). Three people adopted the ring electrode measurement method. The continuous ring electrode for resistivity measurement can improve the accuracy of the resistivity value by increasing the number of measurements, but it also increases the length of the probe. Therefore, this paper combines the measurement number and the length of the probe, and proposes a new circular permutation electrode method.
The resistivity meter for monitoring water quality and measuring soil has become mature, but most of the devices used in marine sediments do not consider the seawater temperature or pressure. Therefore, it is very important to design a resistivity meter for marine sediments considering seawater temperature and pressure. The accuracy of the resistivity meter for measuring soil is ±5%, and the error of the resistivity meter for measuring marine sediments will be larger, so it is also very important to improve the accuracy of the resistivity meter. In view of the lack of consideration of seawater temperature and pressure in the current resistivity meter for marine sediments, this paper designed a new type of resistivity measurement device for marine sediment samples after pressure-maintaining transfer. At the same time, the original ring electrode measurement method was improved, and a new type of circular permutation electrode measurement method was proposed. The corresponding device and probe rod were designed to avoid the influence of temperature and pressure on the resistivity value of sediments.

Seabed Sediment Resistivity Measuring Instrument and Marine Sediment Resistivity Measuring Instrument Pressure-Maintaining Transfer System Design
The resistivity measurement system designed in this paper is mainly used for the measurement of marine sediments after pressure-maintaining transfer. Therefore, it is necessary to consider the influence of pressure and temperature changes of sediments on the measurement results in the process of measurement. In view of this, the marine sediment resistivity measurement instrument and the pressure-maintaining transfer system of marine sediment resistivity measurement instrument are designed in this paper.
The marine sediment resistivity measuring device is shown in Figure 1. The cabin is made of 316 stainless steel. The wall thickness is 12.5 mm and the inner diameter is 65 mm, which can withstand the pressure of 36 MPa. The two ends can be connected with the pressure-maintaining transfer system of the marine sediment resistivity measuring device through the hoop. The cabin has a thread hole, which can maintain environmental pressure when measuring the resistivity of the sediment. tool, the plastic pipe outside the sediment sample is opened. The probe enters the sediment sample through the opening for data collection.
The probe consists of a probe rod and probe head. A POM (Polyoxymethylene) tube was selected as the shell material of probe rod. The probe rod is 360 mm in length and 10 mm in diameter. Six groups of ring electrode grooves were set on the POM tube with equal spacing. Brass was selected as the electrode material to facilitate electrode placement. A lead hole was set in each groove to facilitate electrode connection to the probe. A principle block diagram and device diagram of pressure-maintaining transfer system of the marine sediment resistivity measuring instrument are shown in Figures 2 and 3 as follows. After the pressure-maintaining transfer device obtains the marine sediment core, it is necessary to transfer the sediment core into the marine sediment resistivity measuring instrument without releasing the sediment pressure through the transfer system. The sediment is pushed backward by the long-stroke push device of the pressure maintaining core tube in the pressure-maintaining transfer system, the core samples are cut by the high-fidelity cutting device, the samples are scanned by the CT scanning device to obtain the preliminary detection results, and then the results can be compared with the measured resistivity. After the detection and processing, the sediment core is pushed into The detecting rod is composed of the counter-wire joint, ball valve, O-ring, thread guide rod, thread cap, thrust ball bearing, flange bearing, and elastic retaining ring. The detecting rod is connected to the cabin by the wire joint, and the ball valve control can open and close to control the penetration of the detecting rod and the cabin. O-ring is used for sealing between screw guide rod and ball valve. The thread guide rod is nested with the thread cap, and the tool or probe is rotated through the thread to enter the cabin. The upper part of the detecting rod is equipped with thrust ball bearings, flange bearings, and elastic retaining rings.
The cutting tool is made of stainless steel and adopts an integrated structure. The lower end face is machined with a cutter head. During the drilling process of the cutting tool, the plastic pipe outside the sediment sample is opened. The probe enters the sediment sample through the opening for data collection.
The probe consists of a probe rod and probe head. A POM (Polyoxymethylene) tube was selected as the shell material of probe rod. The probe rod is 360 mm in length and 10 mm in diameter. Six groups of ring electrode grooves were set on the POM tube with equal spacing. Brass was selected as the electrode material to facilitate electrode placement. A lead hole was set in each groove to facilitate electrode connection to the probe.
A principle block diagram and device diagram of pressure-maintaining transfer system of the marine sediment resistivity measuring instrument are shown in Figures 2 and 3 as follows. After the pressure-maintaining transfer device obtains the marine sediment core, it is necessary to transfer the sediment core into the marine sediment resistivity measuring instrument without releasing the sediment pressure through the transfer system. The sediment is pushed backward by the long-stroke push device of the pressure maintaining core tube in the pressure-maintaining transfer system, the core samples are cut by the high-fidelity cutting device, the samples are scanned by the CT scanning device to obtain the preliminary detection results, and then the results can be compared with the measured resistivity. After the detection and processing, the sediment core is pushed into the marine sediment resistivity measuring instrument controlled by the ball valve. Firstly, the drilling device composed of the detecting rod and the cutting tool is drilled. Three test holes are drilled above the plastic tube outside the sediment core sample. The opening position of the test hole is consistent with the probe position of the three test rod devices. The resistivity test rod is moved to the pressure chamber by manual spiral propulsion, so that the probe enters the plastic tube outside the sediment core sample and extends into the drilled test hole opening to directly contact the core of the sediment for resistivity testing and data acquisition. J. Mar. Sci. Eng. 2021, 9, x FOR PEER REVIEW 4 of 24 the marine sediment resistivity measuring instrument controlled by the ball valve. Firstly, the drilling device composed of the detecting rod and the cutting tool is drilled. Three test holes are drilled above the plastic tube outside the sediment core sample. The opening position of the test hole is consistent with the probe position of the three test rod devices. The resistivity test rod is moved to the pressure chamber by manual spiral propulsion, so that the probe enters the plastic tube outside the sediment core sample and extends into the drilled test hole opening to directly contact the core of the sediment for resistivity testing and data acquisition.

Structure Design of the Probe
POM (Polyoxymethylene) tube is selected as the probe shell material. POM is an engineering plastic with excellent performance, which has the hardness, strength, and steel properties similar to those of metals. It has good self-lubrication, fatigue resistance, and chemical resistance in a wide range of temperature and humidity. The probe consists of a probe rod and a probe head. The structure of the probe is shown in Figure 4, and the profile of the probe is shown in Figure 5. The probe rod is set to 360 mm in length and 10 mm in diameter. Six groups of ring electrode grooves are set on the probe rod with equal the marine sediment resistivity measuring instrument controlled by the ball valve. Firstly, the drilling device composed of the detecting rod and the cutting tool is drilled. Three test holes are drilled above the plastic tube outside the sediment core sample. The opening position of the test hole is consistent with the probe position of the three test rod devices. The resistivity test rod is moved to the pressure chamber by manual spiral propulsion, so that the probe enters the plastic tube outside the sediment core sample and extends into the drilled test hole opening to directly contact the core of the sediment for resistivity testing and data acquisition.

Structure Design of the Probe
POM (Polyoxymethylene) tube is selected as the probe shell material. POM is an engineering plastic with excellent performance, which has the hardness, strength, and steel properties similar to those of metals. It has good self-lubrication, fatigue resistance, and chemical resistance in a wide range of temperature and humidity. The probe consists of a probe rod and a probe head. The structure of the probe is shown in Figure 4, and the profile of the probe is shown in Figure 5. The probe rod is set to 360 mm in length and 10 mm in diameter. Six groups of ring electrode grooves are set on the probe rod with equal

Structure Design of the Probe
POM (Polyoxymethylene) tube is selected as the probe shell material. POM is an engineering plastic with excellent performance, which has the hardness, strength, and steel properties similar to those of metals. It has good self-lubrication, fatigue resistance, and chemical resistance in a wide range of temperature and humidity. The probe consists of a probe rod and a probe head. The structure of the probe is shown in Figure 4, and the profile of the probe is shown in Figure 5. The probe rod is set to 360 mm in length and 10 mm in diameter. Six groups of ring electrode grooves are set on the probe rod with equal spacing. Brass is selected as electrode material to facilitate electrode placement. A lead hole is set in each groove to facilitate electrode connection to the probe. At the same time, in order to obtain more accurate resistivity values, this paper set up a number of different electrode spacing and electrode diameter; the first probe rod diameter was fixed at 10 mm, and then changed to different electrode diameters-1 mm, 1.2 mm, 1.4 mm, 1.6 mm, 1.8 mm, and 2.0 mm, respectively, with 0.2 mm interval, a total of 6 groups, after selecting the best electrode diameter. The probe rod diameter were then changed to 10 mm, 11 mm, 12 mm, 13 mm, 14 mm, and 15 mm, respectively, with a 1 mm interval, a total of 6 groups, to determine the best electrode diameter.
spacing. Brass is selected as electrode material to facilitate electrode placement. A lead hole is set in each groove to facilitate electrode connection to the probe. At the same time, in order to obtain more accurate resistivity values, this paper set up a number of different electrode spacing and electrode diameter; the first probe rod diameter was fixed at 10 mm, and then changed to different electrode diameters-1 mm, 1.2 mm, 1.4 mm, 1.6 mm, 1.8 mm, and 2.0 mm, respectively, with 0.2 mm interval, a total of 6 groups, after selecting the best electrode diameter. The probe rod diameter were then changed to 10 mm, 11 mm, 12 mm, 13 mm, 14 mm, and 15 mm, respectively, with a 1 mm interval, a total of 6 groups, to determine the best electrode diameter.

Hardware Design of the Measurement System
The overall hardware structure of the resistivity measurement system is shown in Figure 6. The hardware circuit system mainly includes a master module, squeeze module, acquisition module, communication module, storage module, power module, and temperature module. The master module is mainly responsible for controlling and processing the relevant information of each sub-module and implementing thread management. The squeeze module generates sine wave signals with different frequencies according to the command information issued by the master module and provides them to the electrode through subsequent processing and conversion. The acquisition module is responsible for collecting the weak voltage signals generated at both ends of the electrode and reporting them to the master module after processing and conversion. The communication module is responsible for the information interaction between the master module and the upper machine. The power supply module provides the required waveform stable power for the whole hardware system. The temperature module is responsible for collecting the temperature of the system to provide temperature compensation and correct the resistivity value. The storage module stores the collected data in accordance with the preset format to the SD card. hole is set in each groove to facilitate electrode connection to the probe. At the same time, in order to obtain more accurate resistivity values, this paper set up a number of different electrode spacing and electrode diameter; the first probe rod diameter was fixed at 10 mm, and then changed to different electrode diameters-1 mm, 1.2 mm, 1.4 mm, 1.6 mm, 1.8 mm, and 2.0 mm, respectively, with 0.2 mm interval, a total of 6 groups, after selecting the best electrode diameter. The probe rod diameter were then changed to 10 mm, 11 mm, 12 mm, 13 mm, 14 mm, and 15 mm, respectively, with a 1 mm interval, a total of 6 groups, to determine the best electrode diameter.

Hardware Design of the Measurement System
The overall hardware structure of the resistivity measurement system is shown in Figure 6. The hardware circuit system mainly includes a master module, squeeze module, acquisition module, communication module, storage module, power module, and temperature module. The master module is mainly responsible for controlling and processing the relevant information of each sub-module and implementing thread management. The squeeze module generates sine wave signals with different frequencies according to the command information issued by the master module and provides them to the electrode through subsequent processing and conversion. The acquisition module is responsible for collecting the weak voltage signals generated at both ends of the electrode and reporting them to the master module after processing and conversion. The communication module is responsible for the information interaction between the master module and the upper machine. The power supply module provides the required waveform stable power for the whole hardware system. The temperature module is responsible for collecting the temperature of the system to provide temperature compensation and correct the resistivity value. The storage module stores the collected data in accordance with the preset format to the SD card.

Hardware Design of the Measurement System
The overall hardware structure of the resistivity measurement system is shown in Figure 6. The hardware circuit system mainly includes a master module, squeeze module, acquisition module, communication module, storage module, power module, and temperature module. The master module is mainly responsible for controlling and processing the relevant information of each sub-module and implementing thread management. The squeeze module generates sine wave signals with different frequencies according to the command information issued by the master module and provides them to the electrode through subsequent processing and conversion. The acquisition module is responsible for collecting the weak voltage signals generated at both ends of the electrode and reporting them to the master module after processing and conversion. The communication module is responsible for the information interaction between the master module and the upper machine. The power supply module provides the required waveform stable power for the whole hardware system. The temperature module is responsible for collecting the temperature of the system to provide temperature compensation and correct the resistivity value. The storage module stores the collected data in accordance with the preset format to the SD card.

Mechanical Mechanism Simulation
There are many widely accepted numerical simulation methods: Lagrangian method, Eulerian method, Computational Fluid Dynamics (CFD), and Arbitrary Lagrangian-Eulerian (ALE). The Lagrangian method can accurately describe the particle motion trajectory [10], but in the calculation of such large deformation problems as the penetration soil, the excessive deformation and distortion of the grid can easily lead to the deformation of the calculation grid, and eventually lead to the difficulty of convergence [11]. The Eulerian method overcomes the Lagrangian mesh distortion problem, but the material flow at the interface is difficult to calculate, so the pure Eulerian method is not generally used to analyze the geotechnical penetration process. The Arbitrary Lagrangian_Eulerian method (ALE) combines the advantages of the Lagrangian analysis method and the Eulerian analysis method. ALE allows the grid unit to move independently of the material without changing the grid unit and connection relationship [12]. A complete ALE analysis consists of two steps: first, create a new grid, and then transfer the answers to the old grid and state variables to the new grid. Through this method, the grid and the material point can be separated from each other, so even if the grid has a large distortion, the ALE method can also ensure high quality grid throughout the analysis process. Based on the above mechanism, the finite element simulation analysis of the whole penetration process can be carried out by using the ALE method to penetrate the probe into the core sample.
In order to reduce the calculation time and improve the simulation efficiency, the dynamic implicit analysis method is adopted and the axisymmetric model is used to simulate. Probe part: POM material strength is much greater than the strength of soil, so discrete rigid body simulation is used. Soil part: the soil cylinder with the probe penetration path as the axis, the radius of 0.5 m, and the depth of 1.5 m is selected as the analysis object, and the radius of the cylinder is 50 times the cone radius. It can be considered that the boundary does not affect the stress and deformation of the soil near the axis. The position of the soil surface at the axis is reserved, and the size is the same as the probe cone size. Before penetration, the probe is placed in the opening in advance, and the probe cone tip is flat with the soil surface.

Mechanical Mechanism Simulation
There are many widely accepted numerical simulation methods: Lagrangian method, Eulerian method, Computational Fluid Dynamics (CFD), and Arbitrary Lagrangian-Eulerian (ALE). The Lagrangian method can accurately describe the particle motion trajectory [10], but in the calculation of such large deformation problems as the penetration soil, the excessive deformation and distortion of the grid can easily lead to the deformation of the calculation grid, and eventually lead to the difficulty of convergence [11]. The Eulerian method overcomes the Lagrangian mesh distortion problem, but the material flow at the interface is difficult to calculate, so the pure Eulerian method is not generally used to analyze the geotechnical penetration process. The Arbitrary Lagrangian_Eulerian method (ALE) combines the advantages of the Lagrangian analysis method and the Eulerian analysis method. ALE allows the grid unit to move independently of the material without changing the grid unit and connection relationship [12]. A complete ALE analysis consists of two steps: first, create a new grid, and then transfer the answers to the old grid and state variables to the new grid. Through this method, the grid and the material point can be separated from each other, so even if the grid has a large distortion, the ALE method can also ensure high quality grid throughout the analysis process. Based on the above mechanism, the finite element simulation analysis of the whole penetration process can be carried out by using the ALE method to penetrate the probe into the core sample.
In order to reduce the calculation time and improve the simulation efficiency, the dynamic implicit analysis method is adopted and the axisymmetric model is used to simulate. Probe part: POM material strength is much greater than the strength of soil, so discrete rigid body simulation is used. Soil part: the soil cylinder with the probe penetration path as the axis, the radius of 0.5 m, and the depth of 1.5 m is selected as the analysis object, and the radius of the cylinder is 50 times the cone radius. It can be considered that the boundary does not affect the stress and deformation of the soil near the axis. The position of the soil surface at the axis is reserved, and the size is the same as the probe cone size. Before penetration, the probe is placed in the opening in advance, and the probe cone tip is flat with the soil surface.
Before meshing, in order to refine the vicinity of the cone tip, the soil is first separated and meshed, respectively. The soil element type is an eight-node axisymmetric four-order displacement bilinear quadrilateral pore pressure element (CAX8P). During the penetration process, the soil has a large deformation. In order to reduce the distortion of the soil after deformation, the grid is arranged in the form of upper left to lower right, as shown in Figure 7. The analysis of the whole penetration process is shown in Figure 8. Before meshing, in order to refine the vicinity of the cone tip, the soil is first separated and meshed, respectively. The soil element type is an eight-node axisymmetric four-order displacement bilinear quadrilateral pore pressure element (CAX8P). During the penetration process, the soil has a large deformation. In order to reduce the distortion of the soil after deformation, the grid is arranged in the form of upper left to lower right, as shown in Figure 7. The analysis of the whole penetration process is shown in Figure 8.

Design of Key Parameters for New Kind Circular Permutation Electrode
The ring electrode measurement method is one of the penetration measurement methods, which can directly deploy the electrode array to marine sediments. As shown in

Design of Key Parameters for New Kind Circular Permutation Electrode
The ring electrode measurement method is one of the penetration measurement methods, which can directly deploy the electrode array to marine sediments. As shown in Figure 9, the ring electrode installed at a certain distance is embedded on the surface of the vertical rod. After the probe is installed, it can be measured by selecting four continuous electrodes.

Design of Key Parameters for New Kind Circular Permutation Electrode
The ring electrode measurement method is one of the penetration measurement methods, which can directly deploy the electrode array to marine sediments. As shown in Figure 9, the ring electrode installed at a certain distance is embedded on the surface of the vertical rod. After the probe is installed, it can be measured by selecting four continuous electrodes. Figure 9. Ring electrode measuring probe. Where a is the distance between the two electrodes, and b is the radius of the probe rod.
Similar to the horizontal array four-electrode resistivity measurement technology, the ring-electrode resistivity measurement technology adopts the principle of the electric field of the uniform isotropic medium half-space point power supply [7]. The principle of electric field distribution is shown in Figure 10. Similar to the horizontal array four-electrode resistivity measurement technology, the ring-electrode resistivity measurement technology adopts the principle of the electric field of the uniform isotropic medium half-space point power supply [7]. The principle of electric field distribution is shown in Figure 10. It can be understood that a cylinder with radius of and infinite length is embedded in a homogeneous isotropic medium. The resistivity of the medium is . The current source connected to the annular electrode A2 injects uniform current in its radial direction. Assuming that A1 is at infinity, the current source will form an equipotential surface with the insulating probe rod as the boundary. It can be understood that a cylinder with radius of b and infinite length is embedded in a homogeneous isotropic medium. The resistivity of the medium is ρ. The current source connected to the annular electrode A 2 injects uniform current in its radial direction. Assuming that A 1 is at infinity, the current source will form an equipotential surface with the insulating probe rod as the boundary.
The current density at any point p(r, z) on the equipotential surface is: In the formula 4πR(R + πb 2 ) is the area of the circular equipotential surface, R is the distance from point P to the ring electrode, then: Formula (1) satisfies the boundary condition: J n ≈ I 4πR 2 at infinity. The current intensity on the cylinder surface is 0, according to Ohm law: where V is the P-point potential, the Formula (1) is combined with the Formula (3) to obtain: Then there are: lim b→0 V = Iρ 4πR is the electric potential of point source. The Wenner ring electrode arrangement is the four electrodes arranged according to Wenner, namely equidistant vertical installation on the insulation probe surface.
According to Wenner configuration of four electrodes, external A 1 and A 2 are current source electrodes, and internal V 1 and V 2 are potential difference measurement electrodes. All electrodes are arranged in equal spacing, and the spacing is a.
The potential of V 1 is: The potential of V 2 is: According to the Formulas (5)-(7), the potential difference between V 1 and V 2 is: Let resistivity ρ be: The electrode geometry configuration factor G of the vertical four-ring Wenner arrangement is: After calculating the electrode geometry configuration factor, the resistivity ρ can be calculated.
The traditional Wiener ring electrode method only uses four electrodes, which is not conducive to obtaining accurate values. Therefore, we propose a new kind circu-lar permutation electrode measurement method, and the developed probe is shown in Figure 11. At present, there are three methods of resistivity measurement using the C1C2C3C4, C2C3C4C5, C3C4C5C6 electrode combination. The advantage of this method is to increase the number of measurements and make the calculated value more accurate, but this will also make the probe length too long, which is not conducive to the later design. Therefore, considering the measurement times and the length of the probe, we designed a circular permutation electrode probe, that is, we can use the similar measurement of C1C2C5C6, that is, the current is injected into C1 and C6 to measure the potential difference between C2 and C5. If we design x (x > = 4) electrodes, according to the traditional method, we only have the x−3 group measurement method. If we use the new kind ring circuit measurement method, we can have 4 i C kinds of measurement method, which reduces the length of the probe, saves cost, and facilitates subsequent calculation. At present, there are three methods of resistivity measurement using the C 1 C 2 C 3 C 4 , C 2 C 3 C 4 C 5 , C 3 C 4 C 5 C 6 electrode combination. The advantage of this method is to increase the number of measurements and make the calculated value more accurate, but this will also make the probe length too long, which is not conducive to the later design. Therefore, considering the measurement times and the length of the probe, we designed a circular permutation electrode probe, that is, we can use the similar measurement of C 1 C 2 C 5 C 6 , that is, the current is injected into C 1 and C 6 to measure the potential difference between C 2 and C 5 . If we design x (x > = 4) electrodes, according to the traditional method, we only have the x−3 group measurement method. If we use the new kind ring circuit measurement method, we can have C 4 i kinds of measurement method, which reduces the length of the probe, saves cost, and facilitates subsequent calculation.
For example, for C 1 , C 2 , C 3 , C 5 four electrodes have: According to Formulas (5), (11), and (12), the potential difference between C 1 and C 2 is: For C 1 , C 2 , C 3 , C 6 four electrodes have: According to Formulas (5), (14), and (15), the potential difference between C 1 and C 2 is: The electrode geometry configuration factor G of the vertical four-ring Wenner arrangement is: After calculating the electrode geometry configuration factor, the resistivity ρ can be calculated.
After calculating the resistivity of each electrode combination, the average value is: In this way, more accurate resistivity value can be obtained, and the length of the probe can be appropriately reduced to facilitate the subsequent design.

Key Parameter Design of the Filter and Amplification Module
Because DDS [13] adopts an all-digital structure, spurious is inevitably introduced. There are three main sources of error: spurious caused by phase round-off error of phase accumulator, spurious caused by amplitude quantization error as a result of limited word length of memory, and spurious caused by non-ideal characteristics of a D/A converter. Therefore, in order to obtain a better quality output signal, a low-pass filter is added to the signal output end. The main function of the filter is to filter out the high frequency noise in the output spectrum. Common types of low-pass filters include the Butterworth filter, Chebyshev filter, and elliptical filter. Among them, the transition bands of the Butterworth filter and Chebyshev filter are relatively flat, while the elliptical filter has a steep transition band, which is more suitable for filtering the noise near the passband frequency. Therefore, the elliptical filter is selected in this paper. Considering the chip output signal cut-off frequency of 70 MHZ, a low-pass elliptic filter with cut-off frequency of 70 MHZ is designed.
The order of the elliptic filter is calculated according to Matlab [14], and the calling format is as follows: [n, W p] = ellipord(W p, Ws, Rp, Rs) where n is the minimum order of elliptic filter, Wp is the cut-off angular frequency of the elliptic filter, Ws is the stopband starting angular frequency of the elliptic filter, Rp is the passband ripple, and Rs is the minimum stopband attenuation. The Filter Solution software is used to calculate the parameters of the internal components of the filter, and a 7-order elliptic low-pass filter with a cut-off frequency of 70 MHz is obtained. At the same time, due to the small output signal of AD9851 [15], the signal is amplified to the required range through the amplitude, and the amplification factor is (R 8 + R 9 )/R 7 . Here, OPA454 is used as the amplification chip. OPA454 is a high-voltage operational amplifier produced by TI company. Its dynamic range is 100 V, and the maximum output current is 50 mA. The filter and amplification circuit are shown in Figure 12. The order of the elliptic filter is calculated according to Matlab [14], and the calling format is as follows: [ , ] ( , , , ) n Wp ellipord Wp Ws Rp Rs = where n is the minimum order of elliptic filter, Wp is the cut-off angular frequency of the elliptic filter, Ws is the stopband starting angular frequency of the elliptic filter, Rp is the passband ripple, and Rs is the minimum stopband attenuation. The Filter Solution software is used to calculate the parameters of the internal components of the filter, and a 7order elliptic low-pass filter with a cut-off frequency of 70 MHz is obtained. At the same time, due to the small output signal of AD9851 [15], the signal is amplified to the required range through the amplitude, and the amplification factor is (R8 + R9)/R7. Here, OPA454 is used as the amplification chip. OPA454 is a high-voltage operational amplifier produced by TI company. Its dynamic range is 100 V, and the maximum output current is 50 mA. The filter and amplification circuit are shown in Figure 12.

Key Parameters Design of the Constant Current Source Module
Since AD9851 outputs voltage signal, if voltage power supply is adopted, the device may be damaged due to excessive current, so the power supply circuit of this system adopts a constant current source power supply. Therefore, it is necessary to convert the

Key Parameters Design of the Constant Current Source Module
Since AD9851 outputs voltage signal, if voltage power supply is adopted, the device may be damaged due to excessive current, so the power supply circuit of this system adopts a constant current source power supply. Therefore, it is necessary to convert the AC voltage signal generated by the signal generation circuit to the AC current signal, and the conversion of voltage and current is realized by the voltage control current source. Therefore, as shown in Figure 13, a voltage-current conversion circuit is designed, which is composed of a voltage follower and an in-phase summation operation circuit.

Key Parameters Design of the Constant Current Source Module
Since AD9851 outputs voltage signal, if voltage power supply is adopted, the device may be damaged due to excessive current, so the power supply circuit of this system adopts a constant current source power supply. Therefore, it is necessary to convert the AC voltage signal generated by the signal generation circuit to the AC current signal, and the conversion of voltage and current is realized by the voltage control current source. Therefore, as shown in Figure 13, a voltage-current conversion circuit is designed, which is composed of a voltage follower and an in-phase summation operation circuit.  For the voltage follower: For the in-phase summation operation circuit: The output voltage U 7 is: Substituting Formula (22) into Formula (21), we get: The voltage on R 15 is: So the output current I OUTC is: Therefore, the output current can be adjusted by changing the value of R 15 , and finally a sinusoidal current signal with stable waveform and constant amplitude is output.

System Software Design
The system software design is mainly divided into the lower machine program design and the upper machine program design; the lower machine software design flow chart is shown in Figure 14, mainly through the Stm32Cube to build the corresponding frame-work, and the initialization configuration, the use of multi-threaded processing method in the system of DDS driver module, analog-to-digital conversion module, temperature acquisition module and communication module processing; at the same time, the collected voltage and temperature values are stored in the SD card after the initial Kalman filter processing and the communication module uses the serial port to receive an idle interrupt way to monitor the data sent by the upper machine and respond according to the analysis. The upper machine selects the integrated environment Visual Studio 2019 [16], based on NET Framework 4.8, and uses multi-threaded processing to process the communication module, data analysis module, and data display module in the system. The software design flow chart of the upper machine is shown in Figure 15. After the system is initialized, the man-machine interface instruction waits to be received, and the inquiry instruction is sent to the lower machine after receiving the inquiry instruction. The data is received by the serial port in real time and judges whether it enters the interrupt function. The correct data are preprocessed, displayed, and stored, and the wrong data are discarded. The software design flow chart of the upper machine is shown in Figure 15. After the system is initialized, the man-machine interface instruction waits to be received, and the inquiry instruction is sent to the lower machine after receiving the inquiry instruction. The data is received by the serial port in real time and judges whether it enters the interrupt function. The correct data are preprocessed, displayed, and stored, and the wrong data are discarded. The upper machine software interface of the marine sediment resistivity measuring instrument is shown in Figure 16. After the software is opened, the circuit is connected and the serial port is set to the corresponding port number, and then clicked to start; the resistivity value of the resistivity measuring instrument can be directly measured.

Resistivity Measuring Instrument Experiment
Firstly, the distilled water with standard value was measured, and the appropriate electrode diameter and electrode length of the detecting rod were determined by comparing with the standard value. Secondly, the conductivity of the standard liquid was measured, and the error range was determined by comparing with the standard value. Finally, the resistivity of the marine sediments was measured. The upper machine software interface of the marine sediment resistivity measuring instrument is shown in Figure 16. After the software is opened, the circuit is connected and the serial port is set to the corresponding port number, and then clicked to start; the resistivity value of the resistivity measuring instrument can be directly measured. The upper machine software interface of the marine sediment resistivity measuring instrument is shown in Figure 16. After the software is opened, the circuit is connected and the serial port is set to the corresponding port number, and then clicked to start; the resistivity value of the resistivity measuring instrument can be directly measured.

Resistivity Measuring Instrument Experiment
Firstly, the distilled water with standard value was measured, and the appropriate electrode diameter and electrode length of the detecting rod were determined by comparing with the standard value. Secondly, the conductivity of the standard liquid was measured, and the error range was determined by comparing with the standard value. Finally, the resistivity of the marine sediments was measured.

Resistivity Measuring Instrument Experiment
Firstly, the distilled water with standard value was measured, and the appropriate electrode diameter and electrode length of the detecting rod were determined by comparing with the standard value. Secondly, the conductivity of the standard liquid was measured, and the error range was determined by comparing with the standard value. Finally, the resistivity of the marine sediments was measured.

Electrical Resistivity Measurement of Distilled Water in a Laboratory
The resistivity of the purchased industrial distilled water is 18 MΩ.cm, the electrode diameter is fixed to 1 mm, and different detecting rod diameters are changed to 10 mm, 11 mm, 12 mm, 13 mm, 14 mm, and 15 mm, respectively. The measurement results are shown in Figure 17.

Electrical Resistivity Measurement of Distilled Water in a Laboratory
The resistivity of the purchased industrial distilled water is 18 MΩ.cm, the electrode diameter is fixed to 1 mm, and different detecting rod diameters are changed to 10 mm, 11 mm, 12 mm, 13 mm, 14 mm, and 15 mm, respectively. The measurement results are shown in Figure 17. It can be seen that under the same electrode condition, when the electrode length is 10 mm, the resistivity value is closer to the accurate value, so the electrode length is selected as 10 mm. Calculation error is: A is the measured value, E is the standard value. Substituting the measured value and standard value into (26), the average error is: Keep the detecting rod diameter 10 mm, the electrode diameter is changed to 1.0 mm, 1.2 mm, 1.4 mm, 1.6 mm, 1.8 mm, and 2.0 mm, respectively.
The measured data are drawn by Origin software, as shown in Figure 18. It can be seen that under the same electrode condition, when the electrode length is 10 mm, the resistivity value is closer to the accurate value, so the electrode length is selected as 10 mm. Calculation error is: A is the measured value, E is the standard value. Substituting the measured value and standard value into (26), the average error is: Keep the detecting rod diameter 10 mm, the electrode diameter is changed to 1.0 mm, 1.2 mm, 1.4 mm, 1.6 mm, 1.8 mm, and 2.0 mm, respectively.
The measured data are drawn by Origin software, as shown in Figure 18. It can be seen that under the condition of equal electrode length, when the electrode diameter is 1 mm, the resistivity value is close to the real value, and it can be observed that the smaller the electrode diameter is, the smaller the error is. Due to the limitation of device and technology, it is impossible to make a smaller electrode diameter temporarily. Therefore, we choose the electrode diameter of 1 mm, and the average error is 2.7%. It can be seen that under the condition of equal electrode length, when the electrode diameter is 1 mm, the resistivity value is close to the real value, and it can be observed that the smaller the electrode diameter is, the smaller the error is. Due to the limitation of device and technology, it is impossible to make a smaller electrode diameter temporarily. Therefore, we choose the electrode diameter of 1 mm, and the average error is 2.7%.

Conductivity Measurement of Laboratory Standard Liquid
The EC conductivity buffer standard solution used in this experiment was purchased, and the standard values were 84 us/cm, 1413 us/cm, and 12.88 ms/cm, respectively. In order to prevent the artificial disturbance to the standard solution, the electrode of the resistivity measuring instrument was directly immersed in the standard solution, and the conductivity was measured. Since the resistivity value was measured, we only need to use 1/resistivity value, and then carry out unit conversion.
As shown in Table 1, the conductivity of three groups of different laboratory standard solutions was measured. The average conductivity values obtained from multiple measurements of 84 us/cm, 1413 us/cm and 12,880 us/cm conductivity solutions were compared with the standard Figure 18. Experimental plot of distilled water with equal electrode diameter and different electrode length.

Conductivity Measurement of Laboratory Standard Liquid
The EC conductivity buffer standard solution used in this experiment was purchased, and the standard values were 84 us/cm, 1413 us/cm, and 12.88 ms/cm, respectively. In order to prevent the artificial disturbance to the standard solution, the electrode of the resistivity measuring instrument was directly immersed in the standard solution, and the conductivity was measured. Since the resistivity value was measured, we only need to use 1/resistivity value, and then carry out unit conversion.
As shown in Table 1, the conductivity of three groups of different laboratory standard solutions was measured. The average conductivity values obtained from multiple measurements of 84 us/cm, 1413 us/cm and 12,880 us/cm conductivity solutions were compared with the standard values. The average errors calculated by Formula (26) were 2%, 2.8%, and 2.9%, respectively.

Soil Resistivity Measurement
According to the different resistivity of different soils, we selected pottery clay, peat soil, black calcium soil, kaolinite clay, and loess in the middle red to measure and compare the corresponding resistivity. The soil resistivity is measured using the standard soil resistivity instrument purchased, as shown in Table 2. The resistivity measured by the actual developed resistivity measuring instrument is shown in Table 3. The average resistivity values measured by pottery clay, peat soil (pond silt), black calcium soil (beach), kaolinite clay (lakeside) and loess in the middle red were compared with the standard resistivity values. The average errors calculated by Formula (26) were 2.9%, 2.2%, 2.7%, 0.9% and 1.8%, respectively. The average error obtained from the experiment is within the error range of device design.

Sediment Resistivity Measurement
Since no samples of deep-sea pressure-maintaining sediments were taken, two different sites on the shore of Qizhen Lake were selected in this experiment to measure the resistivity of the sediments on the shore. The field test diagrams of the development instrument and standard instrument at site 1 of Qizhen Lake are shown in Figure 19a The resistivity and average values measured by developed and standard measuring instruments at Site 1 of Qizhen Lake are shown in Table 4.  The resistivity and average values measured by developed and standard measuring instruments at Site 1 of Qizhen Lake are shown in Table 4. The average resistivity measured by the developed instrument is compared with that measured by the standard instrument. The average error calculated by Formula (26) is 3.9%.
The field test diagram of the developed instrument and standard instrument in site 2 of Qizhen Lake is shown in Figure 20a  The resistivity and average values measured by developed and standard measuring instruments at Site 1 of Qizhen Lake are shown in Table 4. The average resistivity measured by the developed instrument is compared with that measured by the standard instrument. The average error calculated by Formula (26) is 3.9%.
The field test diagram of the developed instrument and standard instrument in site 2 of Qizhen Lake is shown in Figure 20a The resistivity and average values measured by developed and standard measuring instruments at Site 2 of Qizhen Lake are shown in Table 5 below. The resistivity and average values measured by developed and standard measuring instruments at Site 2 of Qizhen Lake are shown in Table 5 below. The average resistivity measured by the developed instrument is compared with that measured by the standard instrument. The average error calculated by Formula (26) is 3.5%.

Overall Experiment
In order to verify the overall working performance of the marine sediment resistivity measuring instrument, the core management pressure-maintaining transfer device equipped with the system in the <Research and Application of Gas Hydrate Subsea Drilling and Shipborne Detection Technology> project was carried out in April 2021 with the No. 2 mother ship of the Guanghai Bureau of the China Geological Survey to the Shenhu area of the South China Sea for sea trial function verification.
After reaching the intended sea area, the seabed sediment was first drilled by the 'Hainiu' pressure-maintaining rig, and then the pressure-maintaining core tube was drilled on the shore by the rig. After the pressure monitoring met the pressure-maintaining index, it was docked with the pressure-maintaining transfer device of the shipborne core tube. As shown in Figure 21 below, after successful docking, the samples of the pressure-maintaining core tube were scanned by acoustic wave, CT scan, cutting and segmented pressure-maintaining encapsulation according to the preset control program of the pressuremaintaining transfer system. ing and Shipborne Detection Technology> project was carried out in April 2021 with the No. 2 mother ship of the Guanghai Bureau of the China Geological Survey to the Shenhu area of the South China Sea for sea trial function verification.
After reaching the intended sea area, the seabed sediment was first drilled by the 'Hainiu' pressure-maintaining rig, and then the pressure-maintaining core tube was drilled on the shore by the rig. After the pressure monitoring met the pressure-maintaining index, it was docked with the pressure-maintaining transfer device of the shipborne core tube. As shown in Figure 21 below, after successful docking, the samples of the pressure-maintaining core tube were scanned by acoustic wave, CT scan, cutting and segmented pressure-maintaining encapsulation according to the preset control program of the pressure-maintaining transfer system. The sample sediments containing natural gas hydrate were selected from the data results of high-fidelity cutting and CT scanning, and transferred to the marine sediment resistivity measuring instrument system connected with the pressure-maintaining transfer device. As shown in Figure 22, the resistivity measurement of transferred sediment samples was carried out by the resistivity measuring device inserted into the marine sediment resistivity measuring instrument system. The sample sediments containing natural gas hydrate were selected from the data results of high-fidelity cutting and CT scanning, and transferred to the marine sediment resistivity measuring instrument system connected with the pressure-maintaining transfer device. As shown in Figure 22, the resistivity measurement of transferred sediment samples was carried out by the resistivity measuring device inserted into the marine sediment resistivity measuring instrument system.
According to the preliminary scanning results, marine sediment samples containing natural gas hydrate were selected and measured by a self-made resistivity measuring instrument. A total of three sections of sediment sample tubes were selected for the measurement, and each section of sediment sample was measured five times. The data are shown in Table 6. After average calculation, and compared with the resistivity value of marine natural gas hydrate recorded in the previous literature, the CT scanner shows that the resistivity value of sediments in the area with high hydrate content is significantly lower than that in the area with low hydrate content, and it is obviously observed that the resistivity value of the area containing natural gas hydrate is lower than that of the area without natural gas hydrate.
The sample sediments containing natural gas hydrate were selected from the data results of high-fidelity cutting and CT scanning, and transferred to the marine sediment resistivity measuring instrument system connected with the pressure-maintaining transfer device. As shown in Figure 22, the resistivity measurement of transferred sediment samples was carried out by the resistivity measuring device inserted into the marine sediment resistivity measuring instrument system. According to the preliminary scanning results, marine sediment samples containing natural gas hydrate were selected and measured by a self-made resistivity measuring instrument. A total of three sections of sediment sample tubes were selected for the measurement, and each section of sediment sample was measured five times. The data are shown in Table 6. After average calculation, and compared with the resistivity value of marine natural gas hydrate recorded in the previous literature, the CT scanner shows that the resistivity value of sediments in the area with high hydrate content is significantly lower than that in the area with low hydrate content, and it is obviously observed that the

Conclusions and Foresight
In this paper, a new type of resistivity measurement method based on circular permutation electrode is proposed in the design of marine sediment resistivity measurement system. In the same detecting rod length, the measurement method can carry out more measurements and reduce the error to less than 4% through the selection of components and circuit design, which is better than 5% accuracy of some resistivity measuring instruments on the market when measuring soil. At the same time, the system also includes a set of instrument pressure test cabin suitable for deep-sea sediment resistivity measurement. The cabin takes into account the influence of temperature and pressure on the measurement results. When using this cabin to implement the resistivity measurement of marine sediment under high pressure, the pressure and temperature values of the measured sediment in situ are maintained, and the accuracy of the measurement is improved. By comparing the measured results with the CT scanning results, it is found that the content of hydrate in sediments is inversely proportional to the resistivity of sediments, and the measured resistivity of sediments without hydrate has little change, which is basically consistent with the existing literature data. The subsequent research will focus on the inversion of the composition and content of the sediment through the resistivity value. At the same time, the porosity sensor, penetration force sensor and shear force sensor are embedded in the measurement system to explore the relationship between the composition and content of the sediment and each parameter.