A Technical Note on the Application of Interpolation Method on Modeling of 3D Surface Rock Layer

+e planning stages of dredging require comprehensive and detailed analyses. Identifying the dredging environment is one of the important points. A three-dimensional (3D) geological modeling technology has been shown to be a robust tool for representing and analyzing the conditions of geology. From a general perspective, a 3Dmodel is established by some spatial surfaces. Based on a dredging project, this study investigated the estimation capability of an interpolation method of triangulation combined with BP neural network, for modeling a rock layer surface. +e performance of the proposed model is compared with some conventional methods in the literature.+e results showed that this interpolation method is effective to be employed for surface modeling of the rock layer.


Introduction
A reasonable judgment about the location of all soil types in the dredging area is one of the most important elements for the planning of maritime dredging operations, as the material to be dredged determines the selection of dredging equipment and drives the productivity computations [1][2][3]. With the advance of science, a 3D geological modeling technology has been shown to be a robust tool for representing and analyzing the conditions of geology [4,5]. From a general perspective, a 3D model is established by some spatial surfaces. In other words, if surfaces are determined, the 3D model can be obtained easily among these surfaces.
Surface modeling needs some necessary information, such as real data at sampled sites and predicted values at unsampled sites. Interpolation is the process of predicting the values of attributes at unsampled sites. Many researchers have carried out studies on interpolation methods for the purpose of providing an optimum model structure [6][7][8][9][10]. For an unsampled point in the spatial position, the closer it is to the sampled point, the more likely the attribute value is similar, and this is the most basic assumption of spatial interpolation methods [11]. On this basis, various interpolation methods have been formed for surface estimation, such as triangulation interpolation, inverse distance weighting, nearest neighbor, and splines methods [12,13].
In addition, some scholars try to apply artificial intelligence technology to the study of spatial interpolation methods. Lin and Chen proposed an attribute interpolation method based on radial basis function networks and variograms [14]. Samanta and Bandopadhyay also used radial basis function network for attribute interpolation, and a genetic algorithm was introduced in their study [15]. Tutmez carried out the research of interpolation methods based on fuzzy mathematical algorithms [16]. In the present study, a triangulation interpolation method combined with BP neural network is used for surface modeling. According to the real borehole data information in a dredging project, the availability for rock formation identification based on this interpolation method is examined. Moreover, a 3D model of rock stratum distribution in the construction area is established, and the amount of rock is estimated.

Description of Problem.
ere are three types of 3D models which are wireframe, surface, and solid models, and these models completely and unambiguously define the stratigraphy for the site being modeled [17]. As it is known, a solid model can be created by extending the top surface down to the bottom surface of geologic units. Surface modeling is crucial as it affects the accuracy of solid modeling. us, the utility of different methods over conventional methods should be investigated.

Triangulation Interpolation.
Triangulation interpolation is the most common and simplest interpolation method, especially in geographic information systems. For this method, the equation of a simple planar surface is where Z(x, y) is the attribute value at the specified location (x, y) and a 0 , a 1 , and a 2 are unknown coefficients that can be found by solving three normal equations of the known coordinates (Z(x, y), x, y) of three points. e triangular region surrounded by the three points is called the natural neighbors [18]. Except the three points, any sampling points are not allowed to exist in this region. By assuming that the calculated coefficients are constant within the natural neighbors, the attribute value Z(x, y) of any location (x, y) within the region can be estimated.

Inverse Distance Weighting Interpolation.
Inverse distance weighting interpolation is a method in which the distance between the interpolation point and the sample point is used as the weight [19]. at means the closer the sample point is to the interpolation point, the greater the weight assigned in the calculation. In this method, a simple planar surface can be expressed as where Z is the attribute value of the interpolation point, is the attribute value of the sampling point, n is the number of the sampling points, and D i is the distance between the interpolation point and the ith sampling point.

Interpolation Method of Triangulation Combined with BP Neural Network.
A neural network, constructed under the inspiration of biological neural networks, is a data processing model with a large amount of interconnecting artificial neuronst [20]. Researchers perform input data modeling by adjusting the weights and thresholds among neurons, thereby ultimately achieving the capability to solve real problems [21]. e BP neural network is a type of forward neural network, which consists of an input layer, an output layer, and several hidden layers. e weight coefficient between layers is adjusted by network learning, and only when the coefficient is determined can the correct output value be given according to the input. e BP network learning belongs to the error minimization type of supervised learning, i.e., it requires a learning sample set with known output and minimizes the error between the expected output and actual output through training. Figure 1 shows a three-layer BP network structure with a vector of length M as the input layer and a vector of length J as the output layer. If u i and v i are used to represent the input and output of the ith layer, respectively, n as the number of iterations, the actual network output is , the error signal of the nth iteration can be defined as e j (n) � d j (n) − Y j (n).
us, the error energy can be defined as e interpolation method of triangulation combined with BP neural network, which is used for surface modeling in this study, actually contains two interpolation processes. Figure 2 shows the steps of this method. Firstly, the data points obtained by the actual measurement are generated (Figure 2(a)). Secondly, using these points as learning sample set, virtual data points are randomly generated through BP neural network (Figure 2(b)). Finally, based on all data points (measured and virtual), the natural neighbors are determined by triangulation (Figure 2(c)). It can be found that the vertices of each triangle region contain at least one virtual data point.

Project Description and Borehole Data Information.
e interpolation method of triangulation combined with BP neural network is applied in the 3D geological modeling of a dredging area. is area, which has a design excavation depth of 13.5 m, is located at one coast of China, and the

Input layer
Output layer Hidden layer main dredged materials are soil and rock. Compared with soil, rock excavation is more difficult and requires higher construction costs. erefore, the 3D model research mainly focuses on the rock layer, and the amount of excavation is paid special attention. A total of 12 boreholes can be used for model establishment. To make the method more general, a reference point (B01) is determined, and the coordinate of which is defined as (X, Y) � (0, 0). Once the reference point is determined, other borehole data can be easily converted. Some important information are listed in Figure 3.

Network Establishment and Evaluation.
e establishment and training of BP neural networks are implemented using computer software. Quasi-Newton algorithm is used for network training, and some relative parameters are shown in Table 1. Figure 4 shows the distribution of boreholes. It can be seen that 10 of the 12 boreholes are used for network establishment, and the other two are used to test the network. e borehole information should be preprocessed before training. Each borehole is divided into several data points from the bottom up, with an interval of 0.2 m. e output of the data points above the elevation of boundary between the soil layer and rock layer is counted as 0, and others as 1. erefore, the input layer of the BP neural network consists of 3 factors which are the location (x, y) and depth, and the output layer only consists of a signal sequence of 0 and 1. After preprocessing, there are totally 670 data points for network establishment and 134 data points for the test. Random assignment is made among the data points used for network establishment, which is 70% for training and 30% for verification. For a certain where n is the number of elements with an output of 1 in the depth direction. e operation results of BP neural network are shown in Figures 5 and 6. Figure 5 reflects the trend of network error performance under different iterations. It can be concluded that the error of training samples, verification samples, and test samples decreases as the number of iterations increases. e minimum errors of the training sample and the test sample are 0.0298 and 0.0238, respectively, and the correct rate of the network is 98.5%. Figure 6 is the regression analysis diagram of BP neural network sample data. As can be seen from the figure, the linear correlation coefficients of the training sample, the verification sample, the test sample, and all the samples are 0.897, 0.929, 0.918, and 0.906, respectively. According to the operation results, the trained values of the network roughly match the actual values.

Comparison of Results.
As mentioned above, 2 boreholes are used for testing. Table 2 is a comparison of the 3 interpolation methods described in the present study. For both boreholes, the predicted value obtained with the BP neural network method is much closer to the actual value. Compared with the inverse distance weighting, the means for adjusting the weight by multiple iterations in the BP neural network method seems to be more reasonable than simply using the square of the distance as the weight, especially in dealing with boundary line interpolation. Moreover, all the borehole information provided is utilized in the BP neural network method, which is different from the triangulation interpolation.
In addition, the BP neural network also provides a way to estimate the maximum error of the elevation of boundary between the soil layer and rock layer. If η,m, and k represent the correct rate of the network, the sampling interval, and the count of sample points, respectively, the maximum error can be expressed as follows:    Advances in Civil Engineering 5

Advances in Civil Engineering
In this study, η is 98.5%, m is 0.2 m, and k is 67, so that the maximum error E is 0.2 m.

3D Modeling.
In the modeling based on triangulation interpolation, the number and distribution of boreholes have a great impact on the partitioning of natural neighbors. erefore, it can be inferred that the virtual drilling arrangement is an important part in the modeling using the interpolation method of triangulation combined with BP neural network. 4 steps of virtual drilling arrangements are made in this paper, aiming at obtaining a more accurate 3D model. Figure 7 shows these steps. As can be seen, the first step mainly focuses at the arrangement of key nodes on the boundary lines, while the other steps are the uniform arrangement of the entire dredging area. Figure 8 shows the amount of rocks at each stage, calculated by the model. It can be seen that as the drilling arrangement is refined, the change of rock amount gradually reduced. When step 4 of virtual drilling arrangement is completed, the change of rock amount is less than 5,000 m 3 , compared to step 3.
is result has reached the owner's request, and the final statistic of rock is 362,131 m 3 .
After completing the arrangement of virtual drilling, the 3D model can be established, as shown in Figure 9. It can be seen from the figure that the distribution of the rock layer is between −4 and −13.5 in depth, the trend is gradually increasing from west to east, and the highest part is located in the upper right corner of the dredging area.
is model provides strong support for the dredging operation planning of the project.

Conclusions
Sensitively generated surface models have a crucial importance for the optimum planning and 3D modeling. In this study, a new interpolation method is used for modeling a surface of the rock layer. Based on a dredging project, the result is compared with 2 conventional interpolation methods.
As a result, some merits of the interpolation method of triangulation combined with BP neural network have been recorded. First, this method has provided more accurate results. Second, it gives a way to calculate the maximum error of interpolation points, and the training accuracy rate of BP neural network is considered to be closely related to the error.
Finally, the 3D model of the dredging area is established using this new interpolation method, and the statistic of rock is determined. is case of application aims to provide a new example for material reserve estimation.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request. Disclosure e BP neural network calculation has been done by using the MATLAB2016a software, e relative code is wrote by the corresponding author of this paper, and the name of code is sci-code. e repository link is https://share.weiyun. com/DP6mGX5q.

Conflicts of Interest
e authors declare that they have no conflicts of interest.