Method for quickly calculating floor position seismic intensity and seismic sensation based on source depth

Abstract

The application discloses a method for quickly calculating floor position seismic intensity and seismic feeling based on a seismic source depth, which comprises the following steps: firstly, reading the information of the time, place, magnitude, source depth and rupture direction of an earthquake through a webpage; secondly, obtaining the position of a target building and the epicentral distance, introducing the source depth to correct the traditional seismic intensity attenuation relation model, and obtaining a source depth simplified reduction model; thirdly, quickly calculating the surface seismic intensity of the target building according to the seismic intensity attenuation relation coefficient; fourthly, deducing the calculation formula of the seismic intensity of different floors by combining the floor amplification coefficient formula and the instrument intensity calculation formula, and calculating the seismic intensity of different floors of the target building by the formula; and finally, determining the seismic feeling and possible damage of different floors by referring to the contrast table of the seismic intensity, seismic feeling and building influence. The application considers the influence of the source depth and floor position on the intensity, improves the accuracy of the seismic intensity, and overcomes the shortcomings of the traditional calculation.

Description

Technical Field

[0001] This invention relates to the field of earthquake monitoring technology, and in particular to a method for rapidly calculating the earthquake intensity and sensation at different floor locations based on the focal depth. Background Technology

[0002] Seismic intensity reflects the degree of damage an earthquake causes to the Earth's surface. It is typically manifested in the actual impact of the earthquake on the ground, such as the degree of human sensation, the reaction of objects, the extent of damage or destruction to buildings, and changes in the landscape. While a single earthquake has only one magnitude, its intensity varies from place to place within the same earthquake, depending on factors such as the depth of the earthquake's focal point and its distance from the epicenter (epicentral distance). Commonly used seismic intensity attenuation relationships generally do not take into account the influence of focal depth and building height.

[0003] Currently, most methods for determining earthquake intensity and sensation at different floors, both domestically and internationally, involve deploying strong-motion seismometers or performing complex numerical simulations. Instrumental observation methods require extensive deployment, are costly, and are difficult to implement across large areas. They also struggle to provide guidance and a basis for post-earthquake behavior for the general public in different regions of the country. Existing calculation methods are numerically complex and involve a large workload of simulation calculations; extending them to systems with thousands of units would require enormous computer resources. Summary of the Invention

[0004] The purpose of this invention is to provide a method for rapidly calculating seismic intensity and seismic sensation at different floor locations based on focal depth. This method considers the influence of focal depth and floor location on intensity, corrects the traditional seismic intensity attenuation relationship model, improves the accuracy of seismic intensity at different focal depths, and derives a calculation formula for rapidly calculating seismic intensity at different floors. This method overcomes the shortcomings of traditional numerical observation and numerical calculation methods, such as "high cost, difficult deployment" and "high computational complexity, making it difficult to promote".

[0005] To achieve the above objectives, this invention provides a method for rapidly calculating seismic intensity and perceived tremors at different floor levels based on focal depth, comprising the following steps:

[0006] S1 acquires earthquake information: It compares the earthquake rapid report catalog information released by China Earthquake Networks Center with the time and determines whether to add it to the database, and reads the earthquake's occurrence time, occurrence location (including latitude and longitude), earthquake magnitude, focal depth, and focal rupture direction;

[0007] S2 calculates surface intensity: It obtains the latitude and longitude coordinates of the target building and determines the epicentral distance by combining them with the coordinates of the earthquake epicenter; it introduces the focal depth to correct the traditional seismic intensity attenuation model, deriving a simplified focal depth reduction model; then, based on the seismic intensity attenuation coefficients of different tectonic zones, it quickly calculates the surface seismic intensity of the target building. The simplified focal depth reduction model is as follows:

[0008] PGI = A + B M+C lg(R+R0 (L+D) / E)

[0009] Where: PGI is the surface seismic intensity, M, R, and L are the magnitude, epicentral distance, and focal depth, and A, B, C, D, E, and R0 are the seismic intensity attenuation coefficients;

[0010] S3 Determines seismic intensity of different floors: Utilizing the relationship between peak ground acceleration (PGA) and ground surface acceleration (GFA) in seismic-resistant building design, and the quantitative relationship between instrumental intensity and GFA in the Chinese Seismic Intensity Scale, combined with the floor amplification factor formula and the instrumental intensity calculation formula, a calculation formula for the seismic intensity of different floors is derived. This formula is then used to calculate the seismic intensity of different floors in the target building. The derived formula for calculating the seismic intensity of different floors is as follows:

[0011] PFI = PGI + 3.17lg(1 + (f-1) / F)

[0012] Where: PFI is the seismic intensity of the floor, f is the location of the floor, and F is the total number of floors in the target building;

[0013] S4 Determine the floor level vibration sensation: Refer to the comparison table of earthquake intensity, earthquake vibration sensation, and building impact to determine the floor level vibration sensation and possible damage of different floors.

[0014] Furthermore, step S1 also includes the following steps:

[0015] S11. Access the earthquake rapid report catalog information released by China Earthquake Networks Center through web page mode, and read information including the earthquake's occurrence time, location (including latitude and longitude), and magnitude.

[0016] S12. Obtain earthquake focal depth information from the official earthquake rapid report. For automatic earthquake rapid reports that do not contain focal depth information, select the focal depths of earthquakes within 10 kilometers of the epicenter and with a magnitude difference of about 1 that have been saved in the local earthquake catalog database. Calculate the average value of the selected focal depths as the reference focal depth for this earthquake.

[0017] S13. Select the dominant source rupture direction of the potential source zone where the earthquake source is located, and determine it as the tectonic rupture direction of this earthquake, so as to achieve rapid determination of the tectonic rupture direction of the earthquake.

[0018] Furthermore, in step S1, if the catalog information released by the China Earthquake Networks Center is newly occurring earthquake information, all information of the newly occurring earthquake is entered into the database after obtaining the required earthquake information; if the catalog information released by the China Earthquake Networks Center is historical earthquake information, the required earthquake information is obtained directly.

[0019] Furthermore, in step S2, the latitude and longitude coordinates of the target building are obtained through the smartphone LBS system.

[0020] Furthermore, in step S2, the tectonic zoning method of the fifth-generation seismic zoning map is used to determine the seismic intensity attenuation parameters of different zones.

[0021] Furthermore, in step S3, the floor response amplification factor is the ratio of the floor peak acceleration to the ground peak acceleration, used to characterize the amplification effect of each floor on the ground motion acceleration. The derivation process of the floor seismic intensity calculation formula is as follows:

[0022] Referencing my country's seismic design code (GB50011) The formula for the floor magnification factor provided in 2010 is as follows:

[0023] PFA / PGA = 1 + z / h

[0024] Where: PFA is the peak floor acceleration, PGA is the peak ground acceleration, z is the floor height, and h is the total building height;

[0025] Referring to the instrumental intensity calculation method specified in the "China Seismic Intensity Scale" (GB / T 17742—2020), seismic intensity I is calculated by the following quantitative relationship:

[0026] I=3.17 lgPGA+6.59

[0027] The seismic intensity PFI and surface seismic intensity PGI at different floor locations were obtained:

[0028] PFI=3.17 lgPFA+6.59

[0029] PGI=3.17 lgPGA+6.59

[0030] This leads to the conclusion that:

[0031] PFI-PGI=3.17 lg(PFA / PGA) = 3.17 lg(1+z / h)

[0032] Considering that the height of each floor in the building is relatively uniform, the proportion of each floor's height to the total building height is calculated using the number of floors, leading to the formula for calculating the seismic intensity (PFI) of different floors:

[0033] PFI = PGI + 3.17lg(1 + (f-1) / F)

[0034] Where: f is the floor location, and F is the total number of floors in the target building.

[0035] After adopting the above solution, the beneficial effects of the present invention are as follows:

[0036] 1. This invention comprehensively considers the influence of focal depth and floor location. By introducing focal depth correction to the traditional seismic intensity attenuation relationship model, a simplified focal depth reduction model is derived, improving the accuracy of seismic intensity at different focal depths. This invention also considers the amplification effect of floor locations, deriving a method that can quickly calculate seismic intensity at different floor locations. The method is simple, highly efficient, and easy to perform real-time calculations, providing a faster calculation method for determining seismic intensity and sensation in multi-story or high-rise buildings.

[0037] 2. This invention can calculate the surface seismic intensity of the target building and the seismic intensity of different floors by obtaining earthquake information and target building information through a webpage, without the need to deploy strong motion meters for observation, thus overcoming the disadvantages of traditional numerical observation that are "high cost and difficult to deploy".

[0038] 3. This invention can quickly determine the seismic intensity at the floor location by calculating the seismic intensity at the ground surface, meeting the need for rapid calculation, reducing unnecessary complex numerical simulation processes, and overcoming the shortcomings of traditional numerical calculation methods, such as "high computational complexity and difficulty in promotion".

[0039] 4. The description of the intensity of the earthquake and the impact on buildings adopted in this invention can better help the public understand the situation and provide a basis for subsequent post-earthquake behavior guidance. Attached Figure Description

[0040] Figure 1 This is a flowchart of the calculation method of the present invention. Detailed Implementation

[0041] The present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.

[0042] This invention provides a method for rapidly calculating seismic intensity and sensation at floor locations based on focal depth, such as... Figure 1 As shown, it includes the following steps:

[0043] S1 acquires earthquake information: It compares the earthquake rapid report catalog information published by the China Earthquake Networks Center with the time and determines whether to add it to the database, and reads the earthquake's occurrence time, location (including latitude and longitude), magnitude, focal depth, and focal rupture direction; the specific information acquisition process is as follows:

[0044] S11. Access the earthquake rapid report catalog information released by China Earthquake Networks Center through web page mode, and read information including the earthquake's occurrence time, location (including latitude and longitude), and magnitude.

[0045] S12. Earthquake focal depth information can be obtained from the official earthquake rapid reports. Focal depth information is available in all official earthquake rapid reports. However, there are some automatic earthquake rapid reports that do not report focal depth information. If only such an automatic earthquake rapid report is available, the focal depth can be estimated. The estimation process is as follows: select the focal depths of earthquakes within 10 kilometers of the earthquake origin that have been saved in the local earthquake catalog database and whose magnitudes differ from the earthquakes by about 1. Calculate the average of the multiple selected focal depths as the reference focal depth of this earthquake.

[0046] S13. Select the dominant source rupture direction of the potential source area where the earthquake source is located. The dominant source rupture direction is the tectonic rupture direction of the earthquake. If the metronome zone of the earthquake is elliptical, the dominant source rupture direction is the major axis direction of the ellipse. Selecting the major axis direction of the ellipse can quickly determine the tectonic rupture direction of the earthquake.

[0047] S2 calculates surface intensity by obtaining the latitude and longitude coordinates of the target building and the distance between the target building and the epicenter (i.e., epicentral distance) using a smartphone LBS (Location Based Service) system or other location-based instruments. It then incorporates focal depth to correct the traditional seismic intensity attenuation model, using the seismic intensity attenuation model used in the compilation of the fifth-generation seismic ground motion zoning map as an example.

[0048] The fifth-generation seismic motion attenuation relationship is divided into the following regions: 1) Qinghai-Tibet region: the entire Qinghai-Tibet seismic zone, including the West Kunlun-Pamir seismic belt, Longmenshan seismic belt, Liupanshan-Qilianshan seismic belt, Qaidam-Altun seismic belt, Bayan Har Mountain seismic belt, Xianshuihe-Eastern Yunnan seismic belt, Himalayas seismic belt, Southwestern Yunnan seismic belt, and Central Tibet seismic belt; 2) Xinjiang region: the Xinjiang seismic zone excluding the Tarim-Alashan seismic belt, including the Altai Mountain seismic belt, North Tianshan seismic belt, Central Tianshan seismic belt, and South Tianshan seismic belt. Seismic zones; 3) Strong seismic activity zone in the east, including other areas of the North China seismic zone except the Ordos seismic zone (Yinchuan-Hetao seismic zone, Fenwei seismic zone, North China Plain seismic zone, Tanlu seismic zone, lower Yangtze River-South Yellow Sea seismic zone) and the coastal areas of South China; 4) Moderate to strong seismic zones, including the Northeast seismic zone, the middle and lower Yangtze River seismic zone, the Youjiang seismic zone, the Ordos seismic zone, and the Tarim-Alashan seismic zone; The coefficients of the seismic intensity attenuation relationship (ellipse model) of the different zones obtained through the above partitioning are shown in Table 1:

[0049] Table 1. Coefficients of Seismic Intensity Attenuation Relationship (Elliptic Model) in Different Zones

[0050]

[0051] The seismic intensity attenuation model used in the compilation of the fifth-generation seismic ground motion zoning map is as follows:

[0052] I = A + B M+C lg(R+R0)

[0053] Where: I is the surface seismic intensity, M and R are the magnitude and epicentral distance, and A, B, C, and R0 are regression parameters, also known as seismic intensity attenuation coefficients, which can be obtained from Table 1;

[0054] However, the aforementioned seismic intensity attenuation model does not consider the influence of focal depth on seismic intensity, resulting in inaccurate surface seismic intensity readings. Therefore, this invention introduces focal depth to correct the traditional seismic intensity attenuation model, deriving a simplified focal depth reduction model for rapid calculation of the surface seismic intensity of the target building. The simplified focal depth reduction model is as follows:

[0055] PGI = A + B M+C lg(R+R0 (L+D) / E)

[0056] Where: PGI is the surface seismic intensity, M, R, and L are the magnitude, epicentral distance, and focal depth, A, B, C, and R0 are the seismic intensity attenuation coefficients, and D and E are the seismic intensity attenuation coefficients related to the focal depth, which can be obtained through regression calculation. Based on the research of actual earthquake calculated seismic intensity and surveyed intensity, D is taken as the optimal value of 25, and E is taken as the optimal value of 40.

[0057] S3 Determines the seismic intensity of each floor: Floors amplify seismic intensity. A floor response amplification factor can be introduced to calculate the seismic intensity of a floor. This factor is the ratio of peak ground acceleration (PGA) to ground acceleration (PGA), representing the amplification effect of each floor on ground motion acceleration. This invention utilizes the relationship between the design PGA of seismic-resistant floors and the ground PGA, as well as the quantitative relationship between instrumental intensity and PGA in the Chinese Seismic Intensity Scale. Combining the floor amplification factor formula and the instrumental intensity calculation formula, a calculation formula for the seismic intensity of different floors can be derived. This formula is then used to calculate the seismic intensity of different floors in the target building. The derivation process is as follows:

[0058] Referencing my country's seismic design code (GB50011) The formula for the floor magnification factor provided in 2010 is as follows:

[0059] PFA / PGA = 1 + z / h

[0060] Where: PFA is the peak floor acceleration, PGA is the peak ground acceleration, z is the floor height, and h is the total building height;

[0061] Referring to the instrumental intensity calculation method specified in the "China Seismic Intensity Scale" (GB / T 17742—2020), seismic intensity can be calculated using the following quantitative relationships:

[0062] I=3.17 lgPGA+6.59

[0063] The seismic intensity PFI and surface seismic intensity PGI at different floor locations can be obtained:

[0064] PFI=3.17 lgPFA+6.59

[0065] PGI=3.17 lgPGA+6.59

[0066] This leads to the conclusion that:

[0067] PFI-PGI=3.17 lg(PFA / PGA) = 3.17 lg(1+z / h)

[0068] Considering that the height of each floor in the building is roughly the same, the formula for calculating the seismic intensity PFI of different floors is derived by calculating the percentage of each floor's height to the total building height using the number of floors:

[0069] PFI = PGI + 3.17lg(1 + (f-1) / F)

[0070] Where: PFI is the seismic intensity of the floor, f is the location of the floor, and F is the total number of floors in the target building;

[0071] S4 Determine the Earthquake Sensation on Each Floor: Step S3 calculates the seismic intensity for different floors of the target building. Then, referring to the table comparing seismic intensity, earthquake sensation, and building impact, the earthquake sensation and potential damage on different floors can be determined, completing all steps. The comparison table of seismic intensity, earthquake sensation, and building impact is shown in Table 2:

[0072] Table 2 Comparison of Earthquake Intensity, Sensation, and Building Impact

[0073]

[0074] Furthermore, it should be noted that in step S1, if the catalog information released by the China Earthquake Networks Center is newly occurring earthquake information, all information on the newly occurring earthquake is entered into the database after obtaining the required earthquake information to provide reference data for subsequent earthquake intensity calculations; if the catalog information released by the China Earthquake Networks Center is historical earthquake information, the required earthquake information is obtained directly. In step S12, this invention considers the case where earthquake rapid reports do not have focal depth information. The average of multiple focal depths that are not significantly different from the magnitude and location of this earthquake is calculated as the focal depth of this earthquake, ensuring that the estimated focal depth is basically the same as the actual focal depth. In step S2, it is preferable to use a smartphone LBS system to obtain the latitude and longitude coordinates of the target building. This method provides fast location information acquisition, and smartphones are portable and widely applicable.

[0075] The following specific example further illustrates the calculation of the seismic intensity at the location of the 10th floor of a 20-story building on Xinghai Ninth Road (Xinghai Bay) in Xiamen, taking the magnitude 4.7 earthquake (focal depth 100 km) that occurred in Yilan County, Taiwan Province, on November 1, 2022 at 16:30 as an example:

[0076] S1 obtains earthquake information: Earthquake information is obtained from the China Earthquake Networks Center rapid reporting system.

[0077] S11. Access the earthquake bulletin information released by China Earthquake Networks Center on November 1, 2022 at 16:30, regarding the 4.7-magnitude earthquake in Yilan County, Taiwan Province, my country. The information includes the earthquake's occurrence time (November 1, 2022, 16:30), the latitude and longitude of the epicenter (122.46 degrees east longitude, 24.74 degrees north latitude), and the earthquake magnitude (4.7).

[0078] S12. The earthquake focal depth information is included in the official rapid earthquake report, which is 100 kilometers.

[0079] S13. Using the dominant rupture direction of the potential source area at the location of the earthquake source, the magnitude zone of this earthquake is an ellipse in the northwest direction, and the major axis of the ellipse is 20 degrees west of north. That is, the dominant rupture direction of this earthquake source is 20 degrees west of north.

[0080] S2 calculates surface intensity: The latitude and longitude coordinates (118.131311°E, 24.615695°N) of the 20-story building on Xinghai Ninth Road in Xiamen were determined using a smartphone LBS system. This coordinates were then compared with the latitude and longitude coordinates of the epicenter of this earthquake to determine the epicentral distance, which was 438 kilometers. The surface intensity was calculated using the seismic intensity attenuation model employed in the compilation of the fifth-generation seismic ground motion zoning map.

[0081] This earthquake was located in the southeastern coastal seismic belt, belonging to the eastern strong earthquake zone. Based on the latitude and longitude coordinates of the target location and the epicenter location, trigonometric functions were used to calculate the angle between this earthquake and the major axis of the ellipse to be 21.64 degrees. The ellipse model requires calculation of seismic intensity along both the major and minor axes; the maximum value was selected as the final reference value. The formula for calculating the intensity along the major axis is as follows:

[0082] I 1长 =A 1长 +B 1长 M1+C 1长 lg(R 1长 +R 0长 )

[0083] Among them: I 1长 The earthquake intensity along the long axis of the 10th floor of a 10-story building on Xinghai Ninth Road, Xiamen is M1=4.7, and the projection of the epicentral distance along the long axis is R. 1长 =438 cos21.64° = 407 km. Both the location of this earthquake and the target building are located along the coast of South China, a region prone to strong earthquakes. Option A can be selected from Table 1. 1长 =5.7123、B 1长 =1.3626、C 1长 =-4.2903, and R0长 =25;

[0084] Substituting the above parameters into the calculation formula for the major axis direction, we get:

[0085] I 1长 =5.7123+1.3626 4.7+ (-4.2903) lg(407+25))

[0086] Finally, I was obtained. 1长 =0.8

[0087] Similarly, the formula for calculating the minor axis direction is as follows:

[0088] I 1短 =A 1短 +B 1短 M1+C 1短 lg(R 1短 +R 0短 )

[0089] Among them: I 1短 The earthquake intensity along the minor axis of the 10th floor of a 10-story building on Xinghai Ninth Road, Xiamen, is given by the projection R of the epicentral distance along the minor axis. 1短 =438 sin21.64° = 162 km, select A from Table 1. 1短 =3.6588、B 1短 =1.3626、C 1短 =-3.5406、R 0短 =13;

[0090] Substituting the above parameters into the calculation formula for the minor axis direction, we get:

[0091] I 1短 =3.6588+1.3626 4.7+ (-3.506) lg(162+13))

[0092] Finally, I was obtained. 1短 =2.2

[0093] By selecting the maximum value, i.e. the surface intensity along the minor axis, the surface seismic intensity at Xinghai Ninth Road in this earthquake is 2.2 degrees. Using the latitude and longitude coordinates of the location of Xinghai Ninth Road in Xiamen, the seismic intensity value of the location is read as 1 degree from the map overlaid with the typical seismic intensity attenuation relationship curve. The error is not large. However, the seismic intensity calculated by this seismic intensity attenuation relationship model does not take into account the influence of the focal depth, so the obtained intensity value is not accurate.

[0094] This invention introduces the concept of focal depth and proposes the following simplified reduction model for focal depth:

[0095] PGI = A + B M+C lg(R+R0 (L+D) / E)

[0096] Where M, R, and L represent magnitude, epicentral distance, and focal depth, respectively, with L = 100 km. A, B, C, D, E, and R0 are regression parameters, obtained from Table 1. Based on the study of earthquake intensity calculated from actual earthquakes and the intensity surveyed, D was set to the optimal value of 25, and E was set to the optimal value of 40.

[0097] Referring to the calculation parameters and process described above, the seismic intensity along the major and minor axes was calculated similarly, and the values ​​were substituted to obtain:

[0098] PGI 1长 =5.7123+1.3626 4.7+ (-4.2903) lg(407+25 (100+25) / 40)

[0099] PGI 1短 =3.6588+1.3626 4.7+ (-3.506) lg(162+13 (100+25) / 40)

[0100] PGI was obtained 1长 =0.6, PGI 1短 =1.7;

[0101] Selecting the maximum value, i.e., the surface intensity along the minor axis, we obtain a surface seismic intensity of 1.7 degrees at Xinghai Ninth Road. However, when referring to Table 1 to determine the perceived intensity, the seismic intensity needs to be rounded down to the nearest integer. Following the method of rounding up to the nearest 7 for seismic intensity (i.e., no rounding up for decimals less than or equal to 7, and rounding up for decimals greater than or equal to 8), the surface seismic intensity at Xinghai Ninth Road is 1 degree. Using the latitude and longitude coordinates of Xinghai Ninth Road in Xiamen, the seismic intensity value at this location, obtained from the seismic intensity attenuation curve after focal depth correction overlaid on the map, is also 1 degree. The calculation is correct.

[0102] S3 Determine the seismic intensity of a floor: The seismic intensity of a floor is quickly calculated using the formula derived in this invention. The formula is as follows:

[0103] PFI = PGI + 3.17lg(1 + (f-1) / F)

[0104] Where: PGI is the surface seismic intensity calculated in step S2, PGI=1.7; f is the floor location, which is selected as the 10th floor in this embodiment; F is the total number of floors in the building, which is selected as 20 floors in this embodiment.

[0105] Substituting PGI=1.7, f=10, and F=20 into the calculation formula, we get:

[0106] PFI = 1.7 + 3.17lg(1 + (10-1) / 20) = 2.2 degrees

[0107] According to the method of taking the earthquake intensity as 7-8, the intensity of the earthquake on the 10th floor of the 20-story building on Xinghai Ninth Road in Xiamen is 2 degrees.

[0108] S4 determines the floor level of the earthquake: Referring to Table 2, without considering the floor location, the ground intensity on Xinghai Ninth Road in Xiamen is 1 degree, and the earthquake is imperceptible; considering the floor location, the earthquake intensity at the 10th floor of a 20-story building on Xinghai Ninth Road in Xiamen is 2 degrees, and the earthquake is slightly felt – a particularly sensitive person would feel it even when completely still. The earthquake intensity and the earthquake sensation are highly consistent with the actual situation, conform to the objective physical laws of earthquake intensity attenuation, and meet the needs of rapid calculation.

[0109] The above description is only a preferred embodiment of the present invention and is not intended to limit the design of this case. All equivalent changes made based on the key design features of this case shall fall within the protection scope of this case.

Claims

1. A method for rapidly calculating seismic intensity and sensation at floor locations based on focal depth, characterized in that: The method comprises the following steps: S1: obtaining earthquake information: comparing the time of earthquake information published by China Earthquake Network and judging the storage, reading the earthquake time, earthquake location, earthquake magnitude, focal depth and focal rupture direction, wherein the earthquake location comprises longitude and latitude; S2: calculating the ground intensity: obtaining the longitude and latitude coordinates of the target building and determining the epicenter distance with the epicenter position coordinates; introducing the focal depth to correct the traditional earthquake intensity attenuation relationship model to obtain a simplified focal depth reduction model, and then quickly calculating the ground seismic intensity of the target building according to the seismic intensity attenuation relationship coefficients of different structural partitions, wherein the simplified focal depth reduction model is: PGI = A + B M+C *lg(R+R0* (L+D) / E) Wherein: PGI is the ground seismic intensity, M, R and L are the magnitude, epicentral distance and focal depth, and A, B, C, D, E and R0 are the seismic intensity attenuation relationship coefficients; S3: determining the floor seismic intensity: using the relationship between the building engineering anti-seismic floor design peak acceleration and the ground peak acceleration, and the quantitative relationship between the Chinese seismic intensity table instrument intensity and the peak acceleration, combining the floor amplification coefficient formula and the instrument intensity calculation formula, deriving the calculation formula of the seismic intensity of different floors and calculating the seismic intensity of different floors of the target building through the calculation formula, and the derived calculation formula of the floor seismic intensity is: PFI=PGI+3.17lg(1+(f-1) / F) Wherein: PFI is the floor seismic intensity, f is the floor position, and F is the total number of floors of the target building; S4: determining the floor vibration: referring to the comparison table of seismic intensity, seismic sensation and building influence to determine the floor vibration and possible damage of different floors.

2. The method for quickly calculating floor position seismic intensity and seismic feeling based on source depth according to claim 1, characterized in that: Step S1 further comprises the following steps: S11: accessing the earthquake information published by China Earthquake Network through a web page mode, reading the earthquake time, earthquake location and earthquake magnitude information, wherein the earthquake location comprises longitude and latitude; S12: obtaining the earthquake focal depth information in the formal rapid report earthquake information, for the automatic earthquake rapid report without focal depth information, selecting the focal depth of the earthquake within 10 kilometers of the earthquake location and the difference in magnitude of about 1 level saved in the local earthquake catalog database, calculating the average value of the selected focal depth as the reference focal depth of this earthquake; S13: selecting the dominant focal rupture direction of the potential focal region where the focal source is located to determine the tectonic rupture direction of this earthquake, realizing the rapid determination of the earthquake tectonic rupture direction.

3. The method for quickly calculating floor position seismic intensity and seismic feeling based on source depth according to claim 1, characterized in that: In step S1, if the catalog information published by China Earthquake Network is new earthquake information, after obtaining the required earthquake information, all information of the new earthquake is entered into the database; if the catalog information published by China Earthquake Network is historical earthquake information, the required earthquake information is directly obtained.

4. The method for quickly calculating floor position seismic intensity and seismic feeling based on source depth according to claim 1, characterized in that: In step S2, the longitude and latitude coordinates of the target building are obtained through the LBS system of the smart phone.

5. The method for quickly calculating floor position seismic intensity and seismic feeling based on source depth according to claim 1, characterized in that: In step S2, the tectonic zoning method of the fifth generation earthquake zoning map is used to determine the seismic intensity attenuation parameters of different partitions.

6. The method for quickly calculating floor position seismic intensity and seismic feeling based on source depth according to claim 1, characterized in that: In step S3, the floor amplification coefficient is the ratio of the floor peak acceleration to the ground peak acceleration, which is used to represent the amplification effect of each floor on the ground motion acceleration. The derivation process of the floor seismic intensity calculation formula is as follows: Referring to the Chinese building seismic code GB50011 The floor amplification factor formula provided in 2010 is as follows: PFA / PGA =1+z / h Where: PFA is the floor peak acceleration, PGA is the ground peak acceleration, z is the height of the floor, and h is the total height of the building. Referring to the instrument intensity calculation method specified in the Chinese Seismic Intensity Table GB / T 17742-2020, the seismic intensity I is calculated by the following quantitative relationship: I=3.17*lgPGA+6.59 The seismic intensity PFI at different floor positions and the ground seismic intensity PGI are obtained: PFI=3.17*lgPFA+6.59 PGI=3.17*lgPGA+6.59 Further derivation is as follows: PFI-PGI=3.17*lg(PFA / PGA)=3.17*lg(1+z / h) Considering that the height of each floor of the building is equivalent, the floor height ratio of the total height of the building is calculated by the number of floors, and the final formula for calculating the seismic intensity PFI at different floor positions is obtained: PFI=PGI+3.17lg(1+(f-1) / F) Where: f is the floor position, and F is the total number of floors of the target building.

Images