CN111611422A - SVC-based method and system for automatically generating qualitative graph in earthquake disaster risk assessment - Google Patents
SVC-based method and system for automatically generating qualitative graph in earthquake disaster risk assessment Download PDFInfo
- Publication number
- CN111611422A CN111611422A CN202010435703.2A CN202010435703A CN111611422A CN 111611422 A CN111611422 A CN 111611422A CN 202010435703 A CN202010435703 A CN 202010435703A CN 111611422 A CN111611422 A CN 111611422A
- Authority
- CN
- China
- Prior art keywords
- svc
- data
- imaging data
- risk assessment
- qualitative
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/56—Information retrieval; Database structures therefor; File system structures therefor of still image data having vectorial format
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/55—Clustering; Classification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/58—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually
- G06F16/583—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually using metadata automatically derived from the content
- G06F16/5846—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually using metadata automatically derived from the content using extracted text
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/58—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually
- G06F16/587—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually using geographical or spatial information, e.g. location
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
- G06N20/10—Machine learning using kernel methods, e.g. support vector machines [SVM]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Library & Information Science (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method for automatically generating a qualitative graph in earthquake disaster risk assessment based on SVC, which comprises the following steps of mapping: s1: first imaging data is acquired based on the system. S2: and normalizing the first imaging data to obtain second imaging data. S3: and constructing matrix points by using the second imaging data to obtain third imaging data. S4: and taking the second imaging data as a training sample, and training the third imaging data to obtain fourth imaging data. S5: the grid is constructed with the fourth imaging data. S6: and carrying out grid vectorization to generate a vector distribution map. The steps are integrated in a programming mode, and an automatic mapping function module in the system is formed. The invention also provides a system for automatically generating the qualitative map in the earthquake disaster risk assessment based on the SVC, and the method is adopted to automatically generate the qualitative map in the earthquake disaster risk assessment so as to assess the earthquake disaster risk. The method of the invention can carry out more reasonable classification on qualitative discrete points, rather than simply carrying out classification by interpolation by using specific numerical values. The mapping process can correct defects of interpolation classification and provide accurate evaluation maps for engineers and other professionals. In addition, the mapping method does not use manual drawing any more, thereby improving the mapping efficiency and precision and reducing the labor cost.
Description
Technical Field
The invention relates to the technical field of earthquake disaster risk assessment, in particular to a method and a system for automatically generating a qualitative map in earthquake disaster risk assessment based on SVC.
Background
China is a country with serious earthquake disaster loss in the world. Once an earthquake happens, very serious casualties and huge economic losses are caused. People's correct assessment and scientific management of earthquake disaster risks become the best way to achieve disaster reduction. For risk assessment of earthquake disasters, the chart can provide visual assessment results for people. The map in the earthquake disaster risk assessment is a field soil type distribution map, a field type distribution map, a sandy soil liquefaction distribution map, a soft soil seismic subsidence distribution map, and the like, and further includes various map of urban land generated from the maps and disaster maps similar to the map. The field soil type distribution map is generated according to the type of soil drilled at each discrete point of the research area, such as rock, hard soil or soft rock, medium hard soil, medium soft soil and soft soil types. The field category distribution map is generated based on the field category of the borehole at each discrete point in the area of interest, e.g. I0、Ⅰ1II, III, IV. The sand liquefaction profile is generated based on whether the borehole is liquefied (i.e., liquefied, not liquefied) or slightly liquefied, moderately liquefied, or heavily liquefied at each discrete point in the study area. The distribution diagram of soft soil seismic subsidence isAnd (4) generating a seismic trap (namely seismic trap and non-seismic trap) according to whether drilling holes at each discrete point of the research area are generated. The above figures are all of the planar type.
The commonly used mapping method at present is manually sketched according to the result of discrete points of a research area, and is time-consuming and labor-consuming. When the area of a research area is large, such as thousands of square kilometers, and the discrete points are unevenly distributed, the defects are more obvious. The other method is an evaluation method which forcibly gives 0 and 1 numerical values to the discrete point attribute and then carries out interpolation gridding, namely, a grid is generated by interpolation (adopting inverse distance ratio) of the discrete points, and the grid is classified and output as a vector diagram according to the grid.
Disclosure of Invention
In order to overcome the defects of the existing analysis and evaluation method and solve the defects in the prior art, the invention provides an automatic mapping method and system based on SVC (support vector classification), which has a perfect mapping process taking SVC as a classification tool and embodies the technical advantages of less sample points and no density of SVC.
The invention provides a method for automatically generating a qualitative graph in earthquake disaster risk assessment based on SVC, which comprises the following steps of mapping:
s1: first imaging data is acquired based on the system. The first mapping data comprises original point positions of the drill hole and corresponding attribute information. Preferably, the first imaging data is saved in a text file in the format of a record X, Y, Z for each line, where the X and Y axes represent position coordinates and Z is a classification and must be an integer value (INT).
S2: and normalizing the first imaging data to obtain second imaging data. Preferably, the X, Y axis data of the original point in the first mapping data is normalized.
S3: and constructing matrix points by using the second imaging data to obtain third imaging data. Preferably, interpolation is performed by taking an outsourcing rectangle of the original point location of the second mapping data as a boundary and a set interval after normalization to construct matrix points, so as to obtain third mapping data, where the third mapping data are nxm matrix points, and are the SVC training model. The ranges of N and M need to be determined empirically, typically within 2000 to balance efficiency with resultant fidelity. M, N, the larger the value, the denser the corresponding grid, the closer the SVC prediction is to the trusted data, but the more computationally intensive it is. Preferably, N is 2000 and M is 2000.
S4: and taking the second imaging data as a training sample, and training the third imaging data to obtain fourth imaging data. Preferably, the fourth graph data is an SVC prediction model, and the SVC prediction model is saved as a text file. Preferably, the normalized drilling point information is used as a training sample, the SVC training model is trained to obtain an attribute value of the N × M matrix prediction point, the SVC prediction model is generated, and the SVC prediction model is stored in a text file. Preferably, the kernel function K in the training SVC model is a gaussian kernel function.
S5: the grid is constructed with the fourth imaging data. Preferably, the grid is constructed with X, Y, Z data of the fourth graphing data.
S6: and carrying out grid vectorization to generate a vector distribution map. Preferably, the estimation range is used for clipping, and finally, a distribution diagram in the designated range is generated.
The steps are integrated in a programming mode, and an automatic mapping function module in the system is formed.
The invention also provides a system for automatically generating the qualitative map in the earthquake disaster risk assessment based on the SVC, and the method is adopted to automatically generate the qualitative map in the earthquake disaster risk assessment so as to assess the earthquake disaster risk.
The SVC-based method for automatically generating the qualitative map in the earthquake disaster risk assessment can classify qualitative discrete points more reasonably, and does not simply use specific numerical values for interpolation for classification. The interpolation classification method depends on the density of discrete points, is very hard, easily weakens some non-prominent points, and classifies the points classified as A into B. The SVC-based method of the invention obtains an accurate evaluation graph by constructing a matrix, training an SVC model to correct attribute values of matrix points, and vectorizing a corrected matrix grid. The mapping process can correct defects of interpolation classification and provide accurate evaluation maps for engineers and other professionals. In addition, the mapping method does not use manual drawing any more, thereby improving the mapping efficiency and precision and reducing the labor cost.
Drawings
FIG. 1 is a flow chart illustrating a method for automatically generating a qualitative map in earthquake disaster risk assessment based on SVC according to the present invention;
FIG. 2 is a flow chart of a method for automatically generating a site soil type qualitative map in earthquake disaster risk assessment based on SVC in embodiment 1;
FIG. 3-1 shows a prior art generated plot of field soil type (with point location data) for the study area of example 1;
3-2 show a prior art generated plot field soil type profile (without point location data) for the study area of example 1;
3-3 show the field soil type profile (local magnification, with point location data) for the study area in example 1 generated by the prior art;
FIG. 4-1 shows the plot soil type distribution map (with point location data) of the study area in example 1 automatically generated based on SVC;
4-2 show the plot soil type profile (without point location data) for the study area in example 1 automatically generated based on SVC;
4-3 show the SVC-based automatically generated plot soil type distribution map (local magnification, with point location data) for the study area in example 1;
FIG. 5 is a view showing a program interface for automatically generating a field soil type distribution map based on SVC in example 1;
fig. 6 shows a flowchart of a method for automatically generating a soft soil seismic trapping qualitative map in earthquake disaster risk assessment based on SVC in embodiment 2;
FIG. 7-1 shows a prior art generated soft soil sag profile (with point location data) for the study area in example 2;
FIG. 7-2 shows a prior art generated soft soil sag profile (without point location data) for the study area of example 2;
7-3 show the soft soil sag profiles (local amplification with point location data) for the study area in example 2 generated by the prior art;
FIG. 8-1 shows the soft soil collapse distribution map (with point location data) of the research area in example 2 based on SVC automatic generation;
fig. 8-2 shows the soft soil collapse distribution map (without point location data) of the research area in example 2 based on SVC automatic generation;
8-3 show the soft soil sag distribution map (local magnification, with point location data) of the study area in example 2 based on SVC automatic generation;
fig. 9 shows a program interface for automatically generating a soft soil collapse distribution map based on SVC in example 2.
Detailed Description
SVM (Support Vector Machine) is a supervised learning model related to the associated learning algorithm, and can analyze data, recognize models, and use them for classification and regression analysis. Different models can be made according to different input data. When SVM is used for classification analysis, its input label is a classification value, which is also referred to as SVC at this time. SVC can be used not only for general linear classification, but also for some linearly indivisible problems, since some may be non-linearly separable. When using SVC for non-linear classification, it is often necessary to use a kernel function to assist in the classification. The kernel function can be constructed according to actual conditions, and the existing kernel function can also be selected. However, for a constructed kernel function, it is difficult to verify that it is a semi-positive definite matrix for any gram matrix in the input space, so a common choice is to use an off-the-shelf kernel function. When SVC classification analysis is carried out, the construction and selection of the kernel function can greatly simplify the complexity of calculation.
FIG. 1 shows a flow chart of the SVC-based method for automatically generating a qualitative graph in earthquake disaster risk assessment, which specifically comprises the following steps:
s1: acquiring first imaging data based on a system; the first mapping data comprises original point positions of the drill hole and corresponding attributesAnd (4) information. Preferably, the first imaging data is saved in a text file in the format of a record X, Y, Z for each line, where the X and Y axes represent position coordinates and Z is a classification and must be an integer value (INT). The Z value is obtained by the system through calculation according to a series of national standards and industrial specifications such as building earthquake-resistant design specification GB50011-2010, geotechnical engineering survey specification GB50021-2001(2009 edition) and geotechnical engineering survey specification JGJ83-2011 by utilizing actually measured information of the drilling hole, such as shear wave speed, standard penetration or stratum information. This original information is saved in a text file in the format of a record per line X, Y, Z. The classification information Z correspondingly generates attribute information of the map pieces such as a field soil type distribution map, a field type distribution map, a sandy soil liquefaction distribution map, a soft soil seismic subsidence distribution map and the like, for example, (1) the field soil type distribution map corresponds to the field soil type of each discrete point drilled hole in the research area: rock, hard or soft soil, medium hard soil, medium soft soil, soft soil; (2) the field category distribution map corresponds to the field category of each discrete point drilling in the research area: i0、Ⅰ1II, III, IV; (3) whether the sand liquefaction distribution diagram corresponds to each discrete point drilling hole of the research area is liquefied: liquefaction, no liquefaction, or liquefaction grade: light liquefaction, medium liquefaction, severe liquefaction; and (4) whether the drilling hole of each discrete point of the corresponding research area of the soft soil seismic subsidence distribution map is seismic subsided or not: collapse due to vibration and collapse due to vibration.
S2: and normalizing the first imaging data to obtain second imaging data. Preferably, the X, Y axis data of the original point in the first mapping data is normalized.
Because the coordinate values of the original data collected in step S1 are actual coordinates and have large values, the later SVC training model is abnormal and the classification effect is poor, so normalization is performed here, the default range is [ -1,1], and the grid spacing is more conveniently adjusted by a program.
S3: the second imaging data is used to construct matrix points. Preferably, interpolation is performed by taking an outsourcing rectangle of the original point location in the second mapping data as a boundary and a set distance after normalization to construct matrix points, so as to obtain third mapping data, where the third mapping data is N × M matrix prediction points, and is an SVC training model.
Interpolation is carried out according to the outsourcing range of the original coordinate point and the distance set after normalization in the step S2, an N multiplied by M matrix prediction point is constructed, the attribute value of the N multiplied by M matrix prediction point needs to be predicted, namely the range of the predicted value Z, N and M needs to be determined empirically, and balance can be obtained on efficiency and result fidelity within 2000 generally. N, M, the larger the mesh, the closer the SVC prediction is to the trusted data, but the more computationally intensive and time consuming.
S4: and taking the second imaging data as a training sample, and training the third imaging data to obtain fourth imaging data. Preferably, the fourth graph data is an SVC prediction model. Preferably, the normalized drilling point information is used as a training sample, the SVC training model is trained to obtain an attribute value of the N × M matrix prediction point, the SVC prediction model is generated, and the SVC prediction model is saved in a text file (txt format).
In SVC, a certain rule is used to map a sample that cannot be linearly segmented into a space with a higher latitude, and then find out a hyperplane, so that the core is to find out appropriate parameters, so that the distance between segmented hyperplanes is maximized, and data can be correctly classified. The kernel function K is K (x, y) < f (x), f (y) >, where x and y are input values for n dimensions, and f (x) and f (y) are n-dimensional to m-dimensional mappings (typically m > > n), < f (x), and f (y) > is the inner product of x and y. The kernel function includes a linear kernel function, a polynomial kernel function, a gaussian kernel function, and the like. In this step, python script code is used, which is packaged into exe by pyinstteller, invoked in a procedural manner in the. NET environment, which greatly simplifies the complexity of system integration.
In the case of dense and disordered points, a linear kernel function cannot be used, but a gaussian kernel function (kernel) is selected because training takes time and time is too expensive when a polynomial kernel function (kernel) is used. The gaussian kernel Function is also called Radial Basis Function (RBF) in SVM, and is the most dominant kernel Function of non-linear classification SVM, and the formula is
Due to K (x)(i),x(j))＝Ф(x(i))TФ(x(i)) Through proper mathematical transformation, the characteristic transformation function corresponding to the Gaussian kernel function can be obtained as
The accumulator physical meaning of the foregoing infinite polynomial is to convert the feature vector into an infinite multidimensional vector space, i.e., the gaussian function can expand the input feature into an infinite multidimensional space. The derivation of the formula would be to use taylor's formula.
Wherein K (x)(i),x(j)) Is a Gaussian kernel function, aiOnly at the samples to which the support vector corresponds is not 0. It can be seen that, in predicting the linear combination of the gaussian functions with the central point at the support vector machine, the coefficient of the linear combination is aiy(i). Thus, the gaussian kernel function is also referred to as RBF kernel function, i.e. a linear combination of inverse bell-shaped functions.
When training an SVC with a Radial Basis (RBF) kernel, two parameters have to be considered: c penalty factor and gamma. The parameter C is used for balancing model accuracy and complexity, is in conflict with the simplicity of a decision surface, and can perform valuable conversion on the misclassification of the training sample. A smaller C makes the decision surface smoother, while a higher C aims at correctly classifying all training samples. Gamma defines how much a single training sample can play. Larger gamma will allow other samples to be affected even more. Intuitively, the gamma parameter defines how far the effect of a single training example is reached, with a low value meaning "far" and a high value meaning "close". The gamma parameter can be seen as the inverse of the radius of influence of the sample selected by the model support vector. This C-parameter corresponds the misclassification of the training examples to the simplicity of the decision surface. A low value of C smoothes the decision surface, while a high value of C aims at correctly classifying all training samples by giving the model freedom to choose more samples as support vectors.
S5: the grid is constructed with the fourth imaging data. Preferably, the grid is constructed with the X, Y, Z axis data of the fourth graphical data.
S6: and (4) carrying out grid vectorization on the grid constructed in the step (S5) to generate a vector distribution map. Preferably, the estimation range is used for clipping, and finally, a distribution diagram in the designated range is generated.
Because the matrix coordinates and the corresponding attribute value data of step S4 are trained by the SVC model, and the data is accurate, vectorizing the grid in step S5 is equivalent to including points of the entire evaluation range, thereby obtaining an evaluation map of the corresponding attribute of the entire evaluation range.
The steps are integrated in a programming mode, and an automatic mapping function module in the system is formed. When the graph is formed, the qualitative evaluation graph can be generated only by clicking the graph forming menu.
The method for automatically generating different qualitative graphs in earthquake disaster risk assessment based on SVC according to the present invention will be described in detail below.
Example 1: field soil type distribution map
518km of a part of an urban area in Guangzhou City2The method is specifically described by taking an embodiment as an example, wherein the range is a research area, 517 drill holes are discrete points, and a field soil type distribution diagram in earthquake disaster assessment is automatically generated, and the specific steps are shown in FIG. 2.
S1: and in the mapping system, adding a drilling layer and obtaining drilling point data. According to a related algorithm in building earthquake resistance design specification GB50011-2010 which is already programmed into the system, the system calculates and obtains site soil type attribute information (rock, hard soil or soft rock, medium hard soil, medium soft soil and soft soil) of each drilling point by using drilling actual wave velocity test data, and the site soil type attribute information is shown in table 1. Because the data is too huge, only the fields and parts of the data obtained by calculation are displayed in the table for reference.
TABLE 1 field soil type Attribute Table
The obtained data information is saved in a text file (txt format), and each line records X, Y and Z. Where X, Y represent the specific location coordinates of the drilling point, respectively, and Z is the classification, which must be a shaping (INT), such as 1/2/3/4/5 for rock/hard earth or soft rock/medium hard earth/medium soft earth/weak earth, respectively.
S2: the original X, Y axis data in the step S1 is normalized to avoid that the interpolation value in the large data range is too intensive to exceed the calculation time.
S3: and (4) carrying out interpolation by taking the outsourcing rectangle of the original point location as a boundary and the interval set after normalization to construct N multiplied by M matrix points, namely the SVC training model. The ranges of N and M need to be determined empirically, typically within 2000 to balance efficiency with resultant fidelity. M, N, the larger the value, the denser the corresponding grid, the closer the SVC prediction is to the trusted data, but the more computationally intensive it is. In this embodiment, N, M takes on the value of 2000.
S4: and taking the drilling point data of the field soil type result obtained after normalization in the second step as training samples, training the SVC training model to obtain the attribute values of the N multiplied by M matrix prediction points, generating the SVC prediction model, and storing the SVC prediction model in a text file. For data of a field soil type, a linear kernel function cannot be used due to dense points and uneven distribution, and a polynomial kernel function (when kernel) is too time-consuming and time-consuming for sample data training. Therefore, when a gaussian kernel function (kernel) is selected, and parameters are repeatedly adjusted to determine that gamma is 50, the effect and the time cost are relatively optimal. Since the samples are not uniformly distributed, in order to avoid that part of the data is weakened as a noise point, an equalization weight (class _ weight) is used. Other parameters use default values, and the sample data SVC training code is as follows:
s5: the lattice is constructed from the X, Y, Z data in the text file at step S4, resulting in a 2000 × 2000 lattice.
S6: vectorizing the 2000 x 2000 grid and clipping the results with the evaluation range to generate a field soil type map, see fig. 4-2.
And programming and integrating the steps to form a functional module for automatically generating the field soil type distribution map in the system. And selecting the drill holes participating in the mapping during mapping, clicking the 'field soil type' in a program interface (shown in an attached drawing 5), entering a 'field soil distribution diagram' interface, and clicking a 'generation partition diagram' menu at the lower right corner in the diagram to generate the field soil type distribution diagram. Figures 3-1 (with point location data), 3-2 (without point location data) and 3-3 (with partial magnification, with point location data) show field soil type profiles for the study area of example 1 generated using the prior art. Fig. 4-1 (with point location data), fig. 4-2 (without point location data) and fig. 4-3 (partially enlarged, with point location data) show the field soil type profiles of the study area in example 1 automatically generated based on SVC using the present invention. The two methods of local magnification range are the same. By comparing the graphs automatically generated by the two methods, it can be obviously found that in the research area, the graph 3-1 generated by the prior art has more errors when different soils are classified into graphs, and if drilling points of soft soil and drilling points of medium soft soil are distributed in the range of the soft soil at the two positions at the bottom, the division is obviously unreasonable; the middle hard soil in the upper left part is distributed with middle soft soil drilling points with larger occupied range, which is not reasonable … …, the phenomenon can be more clearly shown in figures 3-3, and meanwhile, the middle hard soil and the middle soft soil with small ranges are not distributed with drilling holes in subareas. By analyzing the generated figure 4-1 and combining the figure 4-3, the invention can find that the regions of various soils correspond well to the distribution of the types of the soil at the drilling points, because the method of the invention carries out more reasonable classification on qualitative discrete points instead of simply carrying out interpolation classification by using specific numerical values, the analysis result is more accurate. It can be seen that the two mapping methods to obtain the map provide the technician with different analysis results.
Data such as terrain, landform and water system are superposed with results generated by the two methods, and the results are found by combining with expert experience, so that a result graph generated based on the SVC method is more reliable.
Example 2: soft soil seismic subsidence distribution map
638km of urban area in Guangzhou city2The method is specifically described for an embodiment in which 5766 boreholes are discrete points in a research area and a soft soil seismic subsidence distribution map in seismic disaster evaluation is automatically generated, and the specific steps are shown in fig. 6.
S1: and in the mapping system, adding a drilling layer to obtain drilling point position data. According to relevant judgment methods and algorithms in specifications of geotechnical engineering survey specification GB50021-2001(2009 edition), soft soil area geotechnical engineering survey specification JGJ83-2011, building earthquake-resistant design specification GB50011-2010 and the like which are already programmed into the system, the system calculates and obtains soft soil earthquake-subsidence attribute information (earthquake subsidence/non-earthquake subsidence) of each drill hole according to actual wave speed test data of the drill hole and the thickness of soft soil in stratum information, and the like, and the soft soil earthquake-subsidence attribute information (earthquake subsidence/non-earthquake subsidence) is shown in table 2. Since the data is too voluminous, only a portion of the data is displayed for reference.
TABLE 2 Soft soil seismic subsidence attribute table
And storing the obtained drilling point location information in a txt format, and recording X, Y and Z in each row. Where X, Y represent the specific location coordinates of the drilling point, respectively, and Z is the classification, which must be the shaping (INT), e.g. 2, 5 represent seismic/non-seismic.
S2: the original X, Y axis data in the step S1 is normalized to avoid that the interpolation value in the large data range is too intensive to exceed the calculation time.
S3: and (4) carrying out interpolation by taking the outsourcing rectangle of the original point location as a boundary and the interval set after normalization to construct N multiplied by M matrix points, namely the SVC training model. The ranges of N and M need to be determined empirically, typically within 2000 to balance efficiency with resultant fidelity. M, N, the larger the value, the denser the corresponding grid, the closer the SVC prediction is to the trusted data, but the more computationally intensive it is. In this embodiment, N, M takes on the value of 2000.
S4: and taking the drill hole point data of the normalized soft soil seismic subsidence result obtained in the second step as a training sample, training the SVC training model to obtain the attribute values of the N multiplied by M matrix prediction points, generating the SVC prediction model, and storing the SVC prediction model in a text file. For the data of soft soil seismic subsidence, because points are dense and are unevenly distributed, a linear kernel function cannot be used, and a polynomial kernel function (when kernel is 'poly'), sample data training is time-consuming and time-consuming. Therefore, when a gaussian kernel function (kernel) is selected, through repeated parameter adjustment, the effect of gamma taking 50 and the time cost are determined to be relatively optimal. Since the samples are not uniformly distributed, in order to avoid that part of the data is weakened as a noise point, an equalization weight (class _ weight) is used. Other parameters use default values, and the sample data SVC training code is as follows:
s5: the lattice is constructed from the X, Y, Z data in the text file at step S4, resulting in a 2000 × 2000 lattice.
S6: vectorizing the 2000 x 2000 grid and clipping the results with the evaluation range to generate a soft soil subsidence distribution map, see fig. 8-2.
And programming and integrating the steps to form a functional module for automatically generating the soft soil seismic subsidence distribution diagram in the system. And selecting the drill holes participating in the image forming during the image forming, clicking a 'soft soil seismic subsidence' menu in a program interface (see figure 9), entering a 'soft soil seismic subsidence distribution map interface', clicking a 'seismic subsidence judgment' menu at the lower right corner in the image to obtain a judgment result of each drill hole point, and clicking a 'partition map generation' menu at the lowest right corner in the image to generate a soft soil seismic subsidence distribution map.
Fig. 7-1 (with point location data), fig. 7-2 (without point location data), and fig. 7-3 (with local magnification, with point location data) show soft soil subsidence profiles for the study area of example 2 generated using the prior art. Fig. 8-1 (with point location data), fig. 8-2 (without point location data), and fig. 8-3 (partially enlarged, with point location data) show soft soil sag profiles for the study area in example 2 that were automatically generated based on SVC using the present invention. The two methods of local magnification range are the same. Comparing fig. 7-1 with fig. 8-1, it can be obviously found that in the research area, many non-seismic-subsidence drilling points are still distributed in the seismic subsidence range divided in the graph generated by the prior art, that is, the non-seismic-subsidence drilling points which can be divided are not divided into the non-seismic subsidence area, but the seismic subsidence drilling points which are originally divided into the seismic subsidence range are not divided. The generated graph reasonably classifies the results of the existing drilling points, and divides a proper region into the seismic subsidence range, wherein the seismic subsidence drilling points are reasonably distributed; different from the graph generated by the prior art, a partial range at the right lower corner is divided into non-seismic subsidence areas, so that more ranges which cannot generate seismic subsidence risks are eliminated. This error is more clearly shown by comparing fig. 7-3 (partially enlarged, with point location data) with fig. 8-3 (partially enlarged, with point location data). The method of the invention carries out more reasonable classification on qualitative discrete points, rather than simply carrying out classification by interpolation by using specific numerical values, and the analysis result is more accurate. It can be seen that the two mapping methods to obtain the map provide the technician with different analysis results.
The results generated by the two methods are superimposed on the image layers of the terrain, the landform, the water system and the like, so that the result image generated based on the SVC method is more reliable.
The present invention is not limited to the above-described embodiments, and various changes and modifications of the present invention are intended to be included within the scope of the claims and the equivalent technology of the present invention if they do not depart from the spirit and scope of the present invention.
Claims (10)
1. A method for automatically generating a qualitative graph in earthquake disaster risk assessment based on SVC comprises the following steps of mapping:
s1: acquiring first imaging data based on a system;
s2: normalizing the first imaging data to obtain second imaging data;
s3: constructing matrix points by using the second imaging data to obtain third imaging data;
s4: training the third imaging data by taking the second imaging data as a training sample to obtain fourth imaging data;
s5: constructing a grid with the fourth mapping data;
s6: carrying out grid vectorization to generate a vector distribution map;
the steps are integrated in a programming mode, and an automatic mapping function module in the system is formed.
2. The method for SVC automatic generation of qualitative charts in seismic disaster risk assessment according to claim 1, characterized in that the first chart data comprises borehole origin point locations and corresponding property information.
3. The method for SVC automatic generation of qualitative graphs in earthquake disaster risk assessment according to claim 1 is characterized in that the X, Y axis data of the original point in the first graph data is normalized.
4. The method of claim 1, wherein a third mapping data is obtained by interpolating a construction matrix with an outer rectangle of an original point of the second mapping data as a boundary and a distance set after normalization, and the third mapping data is nxm matrix points, which is an SVC training model.
5. The method for automatically generating a qualitative graph in an SVC disaster risk assessment according to claim 1, wherein said fourth graph data is an SVC predictive model.
6. The method of claim 1, wherein the normalized drill point information is used as a training sample, the third imaging data is trained to generate an SVC prediction model, and the SVC prediction model is saved in a text file.
7. The method for automatically generating a qualitative graph in an SVC risk assessment according to claim 1, wherein the kernel function K in the training of the third graph data is selected from Gaussian kernel function and gamma is 50.
8. The method for SVC based automatic generation of qualitative graphs in earthquake disaster risk assessment according to claim 4 characterized in that N2000, M2000.
9. A qualitative graph automatically generated using the method of automatically generating a qualitative graph in an SVC risk assessment based on any of claims 1-8.
10. A system for automatically generating a qualitative graph using the method in SVC automatic generation of seismic disaster risk assessment based on any of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010435703.2A CN111611422B (en) | 2020-05-21 | 2020-05-21 | SVC-based method and system for automatically generating qualitative graph in earthquake disaster risk assessment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010435703.2A CN111611422B (en) | 2020-05-21 | 2020-05-21 | SVC-based method and system for automatically generating qualitative graph in earthquake disaster risk assessment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111611422A true CN111611422A (en) | 2020-09-01 |
| CN111611422B CN111611422B (en) | 2020-12-01 |
Family
ID=72205223
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010435703.2A Active CN111611422B (en) | 2020-05-21 | 2020-05-21 | SVC-based method and system for automatically generating qualitative graph in earthquake disaster risk assessment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111611422B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113269852A (en) * | 2021-05-27 | 2021-08-17 | 广东省地震局 | Stratum partition map and stratum isobath mapping method and device |
| CN113419271A (en) * | 2021-05-07 | 2021-09-21 | 中铁二院工程集团有限责任公司 | Earthquake magnitude prediction method, device, equipment and readable storage medium |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103308942A (en) * | 2012-03-12 | 2013-09-18 | 中国石油天然气股份有限公司 | Method and system for visualizing seismic data |
| CN105093285A (en) * | 2014-05-16 | 2015-11-25 | 许丁友 | Artificial intelligence earthquake judgment method and earthquake detection system |
| KR20160024232A (en) * | 2014-08-25 | 2016-03-04 | 서울대학교산학협력단 | 3-Dimensional Space Modeling Method based on the Geotechnical Information |
| US20170193176A1 (en) * | 2015-12-30 | 2017-07-06 | Insilico Medicine, Inc. | System, Method, and Software for Improved Drug Efficacy and Safety in a Patient |
| CN107085662A (en) * | 2017-05-12 | 2017-08-22 | 首都师范大学 | A kind of method that earthquake intensity is extracted based on nighttime light data |
| CN108390380A (en) * | 2018-05-16 | 2018-08-10 | 国网山东省电力公司青岛供电公司 | A kind of transformer state parameter trend forecasting method and system |
| CN109507724A (en) * | 2018-11-01 | 2019-03-22 | 中国石油大学(华东) | A kind of earthquake TEC geochemical anomalies studying method based on non-vibration state ambient field |
| CN109709612A (en) * | 2018-11-07 | 2019-05-03 | 中国石油大学(北京) | Sand-body distribution prediction technique, device, equipment and readable storage medium storing program for executing |
| CN110501743A (en) * | 2019-08-05 | 2019-11-26 | 中海石油(中国)有限公司 | A kind of sandstone thickness isogram preparation method of reservoir architecture constraint |
| CN110837111A (en) * | 2018-08-15 | 2020-02-25 | 中国石油化工股份有限公司 | Seismic data interpolation method and system |
-
2020
- 2020-05-21 CN CN202010435703.2A patent/CN111611422B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103308942A (en) * | 2012-03-12 | 2013-09-18 | 中国石油天然气股份有限公司 | Method and system for visualizing seismic data |
| CN105093285A (en) * | 2014-05-16 | 2015-11-25 | 许丁友 | Artificial intelligence earthquake judgment method and earthquake detection system |
| KR20160024232A (en) * | 2014-08-25 | 2016-03-04 | 서울대학교산학협력단 | 3-Dimensional Space Modeling Method based on the Geotechnical Information |
| US20170193176A1 (en) * | 2015-12-30 | 2017-07-06 | Insilico Medicine, Inc. | System, Method, and Software for Improved Drug Efficacy and Safety in a Patient |
| CN107085662A (en) * | 2017-05-12 | 2017-08-22 | 首都师范大学 | A kind of method that earthquake intensity is extracted based on nighttime light data |
| CN108390380A (en) * | 2018-05-16 | 2018-08-10 | 国网山东省电力公司青岛供电公司 | A kind of transformer state parameter trend forecasting method and system |
| CN110837111A (en) * | 2018-08-15 | 2020-02-25 | 中国石油化工股份有限公司 | Seismic data interpolation method and system |
| CN109507724A (en) * | 2018-11-01 | 2019-03-22 | 中国石油大学(华东) | A kind of earthquake TEC geochemical anomalies studying method based on non-vibration state ambient field |
| CN109709612A (en) * | 2018-11-07 | 2019-05-03 | 中国石油大学(北京) | Sand-body distribution prediction technique, device, equipment and readable storage medium storing program for executing |
| CN110501743A (en) * | 2019-08-05 | 2019-11-26 | 中海石油(中国)有限公司 | A kind of sandstone thickness isogram preparation method of reservoir architecture constraint |
Non-Patent Citations (2)
| Title |
|---|
| 毕丽思 等: "基于高密度钻孔揭露的广州城区软土分布特征及震陷预测", 《地震工程学报》 * |
| 黄元敏 等: "广东省县(市)地震应急能力评价指标体系的构建与分析", 《华南地震》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113419271A (en) * | 2021-05-07 | 2021-09-21 | 中铁二院工程集团有限责任公司 | Earthquake magnitude prediction method, device, equipment and readable storage medium |
| CN113419271B (en) * | 2021-05-07 | 2023-08-29 | 中铁二院工程集团有限责任公司 | Earthquake magnitude prediction method, device, equipment and readable storage medium |
| CN113269852A (en) * | 2021-05-27 | 2021-08-17 | 广东省地震局 | Stratum partition map and stratum isobath mapping method and device |
| CN113269852B (en) * | 2021-05-27 | 2021-12-07 | 广东省地震局 | Stratum partition map and stratum isobath mapping method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111611422B (en) | 2020-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3535607B1 (en) | Seismic data processing artificial intelligence | |
| McGarigal | FRAGSTATS help | |
| CN111611422B (en) | SVC-based method and system for automatically generating qualitative graph in earthquake disaster risk assessment | |
| Wan et al. | A landslide expert system: image classification through integration of data mining approaches for multi-category analysis | |
| Wan | Entropy-based particle swarm optimization with clustering analysis on landslide susceptibility mapping | |
| Yarus et al. | Practical geostatistics-an armchair overview for petroleum reservoir engineers | |
| CN110441820B (en) | Intelligent interpretation method of geological structure | |
| Olalotiti-Lawal et al. | Post-combustion carbon dioxide enhanced-oil-recovery development in a mature oil field: model calibration using a hierarchical approach | |
| CN111709072B (en) | Underground vibration amplitude parameter prediction method | |
| CN115859094A (en) | Method for automatically identifying homodyne landslide of neural network of coupling convolution and entropy importance graph | |
| Hatampour et al. | Hydraulic flow units' estimation from seismic data using artificial intelligence systems, an example from a gas reservoir in the Persian Gulf | |
| Han et al. | Comprehensive analysis for production prediction of hydraulic fractured shale reservoirs using proxy model based on deep neural network | |
| Lock et al. | Visualizing spatial data: The importance of geographic information systems | |
| Waggoner et al. | Improved reservoir modeling with time-lapse seismic in a Gulf of Mexico gas condensate reservoir | |
| CN114402233A (en) | Automatic calibration of forward deposition model | |
| CN117332640A (en) | Earthquake motion amplification factor prediction method and prompting system based on BP neural network | |
| CN114036829B (en) | Geological profile generation method, system, equipment and storage medium | |
| CN113946691A (en) | Foundation soil layering system and method | |
| CN112016048B (en) | Geological attribute prediction method based on multi-point geological statistics | |
| Karimi Ghalehjough et al. | Peak ground acceleration prediction by fuzzy logic modeling for Iranian plateau | |
| Zare et al. | Assessment of land cover changes using Taguchi-based optimized SVM classification approach | |
| WO2021206755A1 (en) | Systems and methods for evaluating a simulation model of a hydrocarbon field | |
| CN117094234B (en) | Landslide vulnerability evaluation method integrating convolutional neural network and self-attention model | |
| Berry et al. | Integrated Assisted History Matching and Forecast Optimisation Under Uncertainty for More Robust Mature Field Redevelopment Project | |
| CN118797764A (en) | Foundation pit deformation prediction method, medium and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |