A method and system for estimating human activity intensity in primitive settlements

Abstract

The present invention relates to a method and system for estimating the intensity of human activities in a primitive settlement. The seismic profile of simulated seismic wave field data and the geological model are compared and analyzed to obtain the seismic waveform characteristics and amplitude threshold value of the pottery accumulation area; according to the seismic waveform characteristics, the amplitude values ​​that meet the amplitude threshold value are extracted from the real seismic wave field data to form a sorted amplitude data body that only meets the characteristics of pottery fragments; the pottery fragment content and the amplitude average value in the seismic body corresponding to the soil sample are calculated, and the correlation model between the two is fitted; the pottery fragment content at different points in the site area is obtained according to the correlation model and the sorted amplitude data body; the pottery fragment content at different points in the site area is calculated according to the pottery fragment content at different points in the site area The cumulative value of the pottery fragment content within the depth range is calculated; the cumulative value of the pottery fragment content is projected onto a plane to obtain a pottery fragment content plane map; the intensity of human activities in the site area is estimated according to the pottery fragment content plane map. The present invention can improve the accuracy of estimating the intensity of human activities.

Description

Technical Field

[0001] The present invention relates to the field of human activity intensity estimation, and in particular to a method and system for estimating human activity intensity in primitive settlements. Background Art

[0002] Generally speaking, researchers use the method of excavation and discovery, and cannot make inferences about the intensity of human activities in the site area before excavation. In addition, most site excavations only select some samples, ignoring some settlement groups in the site area. This poses a great challenge to the estimation of the intensity of human activities, and new methods are urgently needed. Previous researchers have introduced seismic exploration methods into the study of archaeological sites, but they were all used to detect the specific location of the site, and no quantitative calculation research was introduced. Summary of the invention

[0003] The purpose of the present invention is to provide a method and system for estimating the intensity of human activities in primitive settlements, which uses seismic wave exploration methods to quantitatively calculate the content of pottery fragments buried underground, thereby improving the accuracy of estimating the intensity of human activities.

[0004] To achieve the above object, the present invention provides the following solutions:

[0005] A method for estimating the intensity of human activities in primitive settlements includes:

[0006] Based on the geological data, a geological model consistent with the site study area is constructed;

[0007] The software is used for forward modeling to obtain simulated seismic wave field data after stacking and migration;

[0008] Comparative analysis is performed on the seismic profile of the simulated seismic wave field data and the geological model to obtain seismic waveform characteristics and amplitude threshold values ​​of the pottery fragment accumulation area;

[0009] Conduct three-dimensional seismic wave exploration and data processing in the actual archaeological site study area to obtain real seismic wave field data after superposition and migration;

[0010] According to the seismic waveform characteristics, the amplitude values ​​that meet the amplitude threshold value are extracted from the real seismic wave field data to form a sorted amplitude data body that only meets the characteristics of the pottery fragments;

[0011] Selecting a set number of soil samples from the pottery sherd accumulation area, screening the pottery sherds in the soil samples, and calculating the pottery sherd content;

[0012] Calculating the average value of the seismic amplitude in the body corresponding to the soil sample, and fitting and generating a correlation model between the pottery fragment content and the average value of the seismic amplitude in the body;

[0013] According to the correlation model and the sorted amplitude data body, the pottery content at different points in the site area is obtained;

[0014] Select a depth range, and calculate the cumulative value of the pottery content within the depth range according to the pottery content at different points in the site area;

[0015] Projecting the accumulated value of the pottery sherd content onto a plane to obtain a pottery sherd content plane diagram;

[0016] The intensity of human activities in the site area was estimated based on the plan view of pottery content.

[0017] Optionally, the forward modeling using software to obtain simulated seismic wave field data after stacking and migration specifically includes:

[0018] The simulation software SOFI2D is used for forward modeling to obtain the original single-channel seismic data;

[0019] Seismic data processing is performed on the original single-channel seismic data to obtain simulated seismic wave field data after stacking and migration.

[0020] Optionally, selecting a set number or more of soil samples from the pottery sherd accumulation area and screening the pottery sherds in the soil samples to calculate the pottery sherd content specifically includes:

[0021] Select more than 20 soil samples from the pottery sherd accumulation area and screen the pottery sherds in the samples;

[0022] The pottery fragments are calculated by using the formula Ri=Mi / Vi;

[0023] Among them, R i is the content of pottery fragments unearthed in the ith sample plot, Mi is the total mass of pottery fragments unearthed in the ith sample plot, and Vi is the volume of the ith soil sample plot; V i It is the product of the length, width and height of the soil sample.

[0024] Optionally, the calculating the average value of the seismic in-vivo amplitude corresponding to the soil sample plot and fitting a correlation model between the pottery fragment content and the average value of the seismic in-vivo amplitude specifically includes:

[0025] The average amplitude value corresponding to each soil sample in the seismic data volume is calculated using Petrel software;

[0026] According to the pottery fragment content and average amplitude value corresponding to each soil sample, the correlation model between the two was fitted using statistical software:

[0027] R i =a+b*H i ;

[0028] Among them, a, b are constants, R i is the content of pottery fragments unearthed in the i-th sample plot, H i is the average amplitude value corresponding to the i-th sample in the seismic data volume.

[0029] A system for estimating the intensity of human activities in primitive settlements includes:

[0030] The geological model building module is used to build a geological model that is consistent with the site study area based on geological data;

[0031] A simulated seismic wave field data determination module is used to perform forward simulation using software to obtain simulated seismic wave field data after stacking and migration;

[0032] A seismic waveform feature and amplitude threshold value determination module, used to compare and analyze the seismic profile of the simulated seismic wave field data with the geological model to obtain the seismic waveform feature and amplitude threshold value of the pottery fragment accumulation area;

[0033] The real seismic wave field data determination module is used to conduct three-dimensional seismic wave exploration and data processing on the actual archaeological site research area to obtain the real seismic wave field data after superposition and migration;

[0034] A sorted amplitude data volume determination module is used to extract amplitude values ​​that meet the amplitude threshold value from the real seismic wave field data according to the seismic waveform characteristics, so as to form a sorted amplitude data volume that only meets the characteristics of the pottery fragments;

[0035] A pottery sherd content determination module is used to select a set number of soil samples from the pottery sherd accumulation area, screen the pottery sherds in the soil samples, and calculate the pottery sherd content;

[0036] A correlation model fitting module is used to calculate the average value of the seismic amplitude in the body corresponding to the soil sample, and to fit and generate a correlation model between the pottery fragment content and the average value of the seismic amplitude in the body;

[0037] A module for determining the content of pottery fragments at different points, used to obtain the content of pottery fragments at different points in the site area according to the correlation model and the amplitude data body after sorting;

[0038] A pottery fragment content cumulative value determination module is used to select a depth range and calculate the pottery fragment content cumulative value within the depth range according to the pottery fragment content at different points in the site area;

[0039] A sherd content plane map determination module is used to project the accumulated value of the sherd content onto a plane to obtain a sherd content plane map;

[0040] The module for estimating the intensity of human activities is used to estimate the intensity of human activities in the site area according to the plan view of the pottery fragment content.

[0041] Optionally, the simulated seismic wave field data determination module specifically includes:

[0042] The original single-channel seismic data determination unit is used to perform forward modeling using simulation software SOFI2D to obtain the original single-channel seismic data;

[0043] The simulated seismic wave field data determining unit is used to perform seismic data processing on the original single-channel seismic data to obtain simulated seismic wave field data after stacking and migration.

[0044] Optionally, the pottery fragment content determination module specifically includes:

[0045] The pottery sherd screening unit is used to select more than 20 soil samples from the pottery sherd accumulation area and screen the pottery sherds in the samples;

[0046] A pottery fragment content determination unit, used to calculate the pottery fragment content by using the formula Ri=Mi / Vi for the pottery fragment;

[0047] Among them, R i is the content of pottery fragments unearthed in the ith sample plot, Mi is the total mass of pottery fragments unearthed in the ith sample plot, and Vi is the volume of the ith soil sample plot; V i It is the product of the length, width and height of the soil sample.

[0048] Optionally, the correlation model fitting module specifically includes:

[0049] The average amplitude value determination unit is used to calculate the average amplitude value corresponding to each soil sample in the seismic data volume according to the Petrel software;

[0050] The correlation model fitting unit is used to fit the correlation model between the pottery fragment content and the average amplitude value corresponding to each soil sample using statistical software:

[0051] R i =a+b*H i ;

[0052] Among them, a, b are constants, R i is the content of pottery fragments unearthed in the i-th sample plot, H i is the average amplitude value corresponding to the i-th sample in the seismic data volume.

[0053] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0054] The present invention provides a method for estimating the intensity of human activities in a primitive settlement. The method comprises the following steps: comparing and analyzing a seismic profile of simulated seismic wave field data with a geological model to obtain seismic waveform characteristics and an amplitude threshold value of a pottery fragment accumulation area; extracting amplitude values ​​that meet the amplitude threshold value from real seismic wave field data according to the seismic waveform characteristics to form a sorted amplitude data body that only meets the characteristics of pottery fragments; calculating the amplitude average values ​​in the seismic body corresponding to the pottery fragment content and the soil sample quadrat, and fitting a correlation model between the two; obtaining the pottery fragment content at different points in the site area according to the correlation model and the sorted amplitude data body; calculating the accumulated value of the pottery fragment content within a depth range according to the pottery fragment content at different points in the site area; projecting the accumulated value of the pottery fragment content onto a plane to obtain a pottery fragment content plane map; estimating the intensity of human activities in the site area according to the pottery fragment content plane map, thereby improving the accuracy of estimating the intensity of human activities. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0056] Figure 1 This is a flow chart of the method for estimating the intensity of human activities in primitive settlements according to the present invention;

[0057] Figure 2 This is a system structure diagram for estimating the intensity of human activities in primitive settlements for the present invention. DETAILED DESCRIPTION

[0058] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0059] The purpose of the present invention is to provide a method and system for estimating the intensity of human activities in primitive settlements, which can improve the accuracy of estimating the intensity of human activities by using seismic wave exploration methods to quantitatively calculate the content of pottery fragments.

[0060] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0061] Figure 1 This is a flow chart of the method for estimating the intensity of human activities in primitive settlements in the present invention. Figure 1As shown in Figure 1, a method for estimating the intensity of human activities in primitive settlements includes:

[0062] Step 101: Based on the geological data, construct a geological model that is consistent with the site study area.

[0063] The geological model parameters include model length, depth, longitudinal wave velocity of soil, longitudinal wave velocity of pottery fragments, thickness of pottery fragments accumulation, etc.

[0064] Step 102: Use software to perform forward simulation to obtain simulated seismic wave field data after stacking and migration, specifically including:

[0065] The simulation software SOFI2D is used to perform forward modeling to obtain original single-channel seismic data. The simulation software used includes but is not limited to SOFI2D.

[0066] Seismic data processing is performed on the original single-channel seismic data to obtain simulated seismic wave field data after stacking and migration.

[0067] Step 103: Compare and analyze the seismic profile of the simulated seismic wave field data and the geological model to obtain the seismic waveform characteristics and amplitude threshold value of the pottery fragment accumulation area.

[0068] Step 104: Conduct three-dimensional seismic wave exploration and data processing on the actual archaeological site study area to obtain real seismic wave field data after superposition and migration.

[0069] Step 105: According to the seismic waveform characteristics, the amplitude values ​​that meet the amplitude threshold value are extracted from the real seismic wave field data to form a sorted amplitude data volume that only meets the characteristics of the pottery fragments.

[0070] Step 106: selecting a set number or more of soil samples from the pottery sherd accumulation area, screening the pottery sherds in the soil samples, and calculating the pottery sherd content, specifically including:

[0071] Select more than 20 soil samples from the pottery fragment accumulation area and screen the pottery fragments in the samples.

[0072] The pottery content is calculated using the formula Ri=Mi / Vi.

[0073] Among them, R i is the content of pottery fragments unearthed in the ith sample plot, Mi is the total mass of pottery fragments unearthed in the ith sample plot, and Vi is the volume of the ith soil sample plot; V i It is the product of the length, width and height of the soil sample.

[0074] Step 107: Calculate the average value of the seismic amplitude in the soil sample, and fit the correlation model between the pottery content and the average value of the seismic amplitude in the soil sample, which specifically includes:

[0075] The average amplitude value corresponding to each soil sample in the seismic data volume is calculated according to Petrel software, and the software here includes but is not limited to Petrel software.

[0076] According to the pottery fragment content and average amplitude value corresponding to each soil sample, the correlation model between the two was fitted using statistical software:

[0077] R i =a+b*H i ;

[0078] Among them, a, b are constants, R i is the content of pottery fragments unearthed in the i-th sample plot, H i is the average amplitude value corresponding to the i-th sample in the seismic data volume.

[0079] Step 108: According to the correlation model and the sorted amplitude data body, the pottery content at different points in the site area is obtained.

[0080] H is a variable about the coordinates (x, y, z), so R is also a variable about (x, y, z).

[0081] According to the correlation model, the pottery content R(x,y,z) of each coordinate point is obtained:

[0082] R(x,y,z)=a+b*H(x,y,z)

[0083] Among them, R(x,y,z) is the pottery content at the coordinate point (x,y,z); H(x,y,z) is the amplitude value at the coordinate point (x,y,z).

[0084] Step 109: Select a depth range, and calculate the cumulative value of the pottery content within the depth range according to the pottery content at different points in the site area. (A is the cumulative value of the pottery content at different depths at a certain point, and R is the pottery content at a certain point)

[0085] The accumulated value A of the pottery content is a variable related to (x, y, z). When calculating the accumulated content of a point, x and y are fixed values, and the depth value z changes. Therefore, when calculating the accumulated value A' of a certain point A, the formula can be used:

[0086]

[0087] Step 110: Project the accumulated value of the pottery content onto a plane to obtain a pottery content plane diagram.

[0088] The accumulation of pottery sherd contents is to accumulate the pottery sherd contents corresponding to different depths under the same coordinate; finally, it is projected onto a plane to form a pottery sherd content plane within a specified depth range.

[0089] Step 111: Estimate the intensity of human activities in the site area based on the pottery fragment content plan.

[0090] According to the plan view of pottery content, it can be inferred that the human activity intensity is high in areas with high pottery content; conversely, the human activity intensity is low.

[0091] The changes in the scale of human land use in different eras are evaluated by the planar distribution map of pottery content between different depth ranges (representing different eras). Among them, different depth ranges represent different periods of human activity. Different depth layers are selected, and the changes in the intensity of human activities at the same place in different eras are explained according to steps 110-111. At the same time, the changes in the scale of human land use can be estimated based on the change results.

[0092] The present invention uses seismic wave exploration methods to quantitatively calculate the content of pottery fragments, supplementing scientific and technological archaeological research methods. This technology is aimed at shallow ancient human sites. Due to the needs of scientific research, it is necessary to estimate the intensity of human activities to characterize the changes in the scale of human land use.

[0093] Figure 2 This is the system structure diagram for estimating the human activity intensity in the primitive settlements in this invention. Figure 2 As shown in Figure 1, a system for estimating the intensity of human activities in primitive settlements includes:

[0094] The geological model building module 201 is used to build a geological model that is consistent with the site study area based on geological data;

[0095] A simulated seismic wave field data determination module 202 is used to perform forward modeling using software to obtain simulated seismic wave field data after stacking and migration;

[0096] A seismic waveform feature and amplitude threshold value determination module 203 is used to compare and analyze the seismic profile of the simulated seismic wave field data with the geological model to obtain the seismic waveform feature and amplitude threshold value of the pottery fragment accumulation area;

[0097] The real seismic wave field data determination module 204 is used to perform three-dimensional seismic wave exploration and data processing on the actual archaeological site research area to obtain real seismic wave field data after superposition and migration;

[0098] A sorted amplitude data volume determination module 205 is used to extract amplitude values ​​that meet the amplitude threshold value from the real seismic wave field data according to the seismic waveform characteristics, and form a sorted amplitude data volume that only meets the characteristics of the pottery fragments;

[0099] The pottery sherd content determination module 206 is used to select a set number of soil samples from the pottery sherd accumulation area, screen the pottery sherds in the soil samples, and calculate the pottery sherd content;

[0100] A correlation model fitting module 207 is used to calculate the average value of the seismic amplitude in the body corresponding to the soil sample, and to fit and generate a correlation model between the pottery fragment content and the average value of the seismic amplitude in the body;

[0101] The pottery content determination module 208 at different points is used to obtain the pottery content at different points in the site area according to the correlation model and the sorted amplitude data body;

[0102] The pottery fragment content cumulative value determination module 209 is used to select a depth range and calculate the pottery fragment content cumulative value within the depth range according to the pottery fragment content at different points in the ruins area;

[0103] A sherd content plane map determination module 210 is used to project the accumulated sherd content values ​​onto a plane to obtain a sherd content plane map;

[0104] The human activity intensity estimation module 211 is used to estimate the intensity of human activities in the site area according to the pottery fragment content plan.

[0105] The simulated seismic wave field data determination module 202 specifically includes:

[0106] The original single-channel seismic data determination unit is used to perform forward modeling using simulation software SOFI2D to obtain the original single-channel seismic data;

[0107] The simulated seismic wave field data determination unit is used to perform seismic data processing on the original single-channel seismic data to obtain simulated seismic wave field data after stacking and migration.

[0108] The pottery fragment content determination module 206 specifically includes:

[0109] The pottery sherd screening unit is used to select more than 20 soil samples from the pottery sherd accumulation area and screen the pottery sherds in the samples;

[0110] A pottery fragment content determination unit, used to calculate the pottery fragment content by using the formula Ri=Mi / Vi for the pottery fragment;

[0111] Among them, R i is the content of pottery fragments unearthed in the ith sample plot, Mi is the total mass of pottery fragments unearthed in the ith sample plot, and Vi is the volume of the ith soil sample plot; V i It is the product of the length, width and height of the soil sample.

[0112] The correlation model fitting module 207 specifically includes:

[0113] The average amplitude value determination unit is used to calculate the average amplitude value corresponding to each soil sample in the seismic data volume according to the Petrel software;

[0114] The correlation model fitting unit is used to fit the correlation model between the pottery fragment content and the average amplitude value corresponding to each soil sample using statistical software:

[0115] R i =a+b*H i ;

[0116] Among them, a, b are constants, R i is the content of pottery fragments unearthed in the i-th sample plot, H i is the average amplitude value corresponding to the i-th sample in the seismic data volume.

[0117] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0118] The principles and implementation methods of the present invention are described in this article using specific examples. The description of the above embodiments is only used to help understand the method and core idea of ​​the present invention. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for estimating the intensity of human activities in primitive settlements, characterized in that: include: Based on the geological data, a geological model consistent with the site study area is constructed; The software is used for forward modeling to obtain simulated seismic wave field data after stacking and migration; Comparative analysis is performed on the seismic profile of the simulated seismic wave field data and the geological model to obtain seismic waveform characteristics and amplitude threshold values ​​of the pottery fragment accumulation area; Conduct three-dimensional seismic wave exploration and data processing in the actual archaeological site study area to obtain real seismic wave field data after superposition and migration; According to the seismic waveform characteristics, the amplitude values ​​that meet the amplitude threshold value are extracted from the real seismic wave field data to form a sorted amplitude data body that only meets the characteristics of the pottery fragments; Selecting a set number of soil samples from the pottery sherd accumulation area, screening the pottery sherds in the soil samples, and calculating the pottery sherd content; Calculating the average value of the seismic amplitude in the body corresponding to the soil sample, and fitting and generating a correlation model between the pottery fragment content and the average value of the seismic amplitude in the body; According to the correlation model and the sorted amplitude data body, the pottery content at different points in the site area is obtained; Select a depth range, and calculate the cumulative value of the pottery content within the depth range according to the pottery content at different points in the site area; Projecting the accumulated value of the pottery sherd content onto a plane to obtain a pottery sherd content plane diagram; The intensity of human activities in the site area was estimated based on the plan view of pottery content.

2. The method for estimating the intensity of human activities in primitive settlements according to claim 1, characterized in that: The forward modeling using software to obtain simulated seismic wave field data after stacking and migration specifically includes: The simulation software SOFI2D is used for forward modeling to obtain the original single-channel seismic data; Seismic data processing is performed on the original single-channel seismic data to obtain simulated seismic wave field data after stacking and migration.

3. The method for estimating human activity intensity in primitive settlements according to claim 1, characterized in that: The method of selecting a set number or more of soil samples from the pottery sherd accumulation area and screening the pottery sherds in the soil samples to calculate the pottery sherd content specifically includes: Select more than 20 soil samples from the pottery sherd accumulation area and screen the pottery sherds in the samples; The pottery fragments are calculated by using the formula Ri=Mi / Vi; Among them, R i is the content of pottery fragments unearthed in the ith sample plot, Mi is the total mass of pottery fragments unearthed in the ith sample plot, and Vi is the volume of the ith soil sample plot; V i It is the product of the length, width and height of the soil sample.

4. The method for estimating human activity intensity in primitive settlements according to claim 1, characterized in that: The calculating of the mean value of the seismic amplitude in the body corresponding to the soil sample and fitting the correlation model between the pottery fragment content and the mean value of the seismic amplitude in the body specifically includes: The average amplitude value corresponding to each soil sample in the seismic data volume is calculated using Petrel software; According to the pottery fragment content and average amplitude value corresponding to each soil sample, the correlation model between the two was fitted using statistical software: R i =a+b*H i ; Among them, a, b are constants, R i is the content of pottery fragments unearthed in the i-th sample plot, H i is the average amplitude value corresponding to the i-th sample in the seismic data volume.

5. A system for estimating the intensity of human activities in primitive settlements, characterized in that: include: The geological model building module is used to build a geological model that is consistent with the site study area based on geological data; A simulated seismic wave field data determination module is used to perform forward modeling using software to obtain simulated seismic wave field data after stacking and migration; A seismic waveform feature and amplitude threshold value determination module, used to compare and analyze the seismic profile of the simulated seismic wave field data with the geological model to obtain the seismic waveform feature and amplitude threshold value of the pottery fragment accumulation area; The real seismic wave field data determination module is used to conduct three-dimensional seismic wave exploration and data processing on the actual archaeological site research area to obtain the real seismic wave field data after superposition and migration; A sorted amplitude data volume determination module is used to extract amplitude values ​​that meet the amplitude threshold value from the real seismic wave field data according to the seismic waveform characteristics, so as to form a sorted amplitude data volume that only meets the characteristics of the pottery fragments; A pottery sherd content determination module is used to select a set number of soil samples from the pottery sherd accumulation area, screen the pottery sherds in the soil samples, and calculate the pottery sherd content; A correlation model fitting module is used to calculate the average value of the seismic amplitude in the body corresponding to the soil sample, and to fit and generate a correlation model between the pottery fragment content and the average value of the seismic amplitude in the body; A module for determining the content of pottery fragments at different points, used to obtain the content of pottery fragments at different points in the site area according to the correlation model and the amplitude data body after sorting; A pottery fragment content cumulative value determination module is used to select a depth range and calculate the pottery fragment content cumulative value within the depth range according to the pottery fragment content at different points in the site area; A sherd content plane map determination module is used to project the accumulated value of the sherd content onto a plane to obtain a sherd content plane map; The module for estimating the intensity of human activities is used to estimate the intensity of human activities in the site area according to the plan view of the pottery fragment content.

6. The system for estimating human activity intensity in primitive settlements according to claim 5, characterized in that: The simulated seismic wave field data determination module specifically includes: The original single-channel seismic data determination unit is used to perform forward modeling using simulation software SOFI2D to obtain the original single-channel seismic data; The simulated seismic wave field data determining unit is used to perform seismic data processing on the original single-channel seismic data to obtain simulated seismic wave field data after stacking and migration.

7. The system for estimating human activity intensity in primitive settlements according to claim 5, characterized in that: The pottery fragment content determination module specifically includes: The pottery sherd screening unit is used to select more than 20 soil samples from the pottery sherd accumulation area and screen the pottery sherds in the samples; A pottery fragment content determination unit, used to calculate the pottery fragment content by using the formula Ri=Mi / Vi for the pottery fragment; Among them, R i is the content of pottery fragments unearthed in the ith sample plot, Mi is the total mass of pottery fragments unearthed in the ith sample plot, and Vi is the volume of the ith soil sample plot; V i It is the product of the length, width and height of the soil sample.

8. The system for estimating human activity intensity in primitive settlements according to claim 5, characterized in that: The correlation model fitting module specifically includes: The average amplitude value determination unit is used to calculate the average amplitude value corresponding to each soil sample in the seismic data volume according to the Petrel software; The correlation model fitting unit is used to fit the correlation model between the pottery fragment content and the average amplitude value corresponding to each soil sample using statistical software: R i =a+b*H i ; Among them, a, b are constants, R i is the content of pottery fragments unearthed in the i-th sample plot, H i is the average amplitude value corresponding to the i-th sample in the seismic data volume.

Images