CN104933117B - Processing method for water pumping test data - Google Patents

Abstract

The present invention relates to a processing method of pumping test data, mainly including pumping test data conversion and processing technology of .mdf format file, non-target layer segment cutting technology of hydrological comprehensive histogram, used for drawing Q and S-t curves of pumping test And the multi-period variable scale and data positioning linkage technology required to restore the water level curve, the graphic element attribute identification and mouse dynamic positioning technology, dynamic data management technology, etc. . The present invention develops a set of data collection, data processing, data management and computer graphics from the perspective of data collection—>data conversion—>data processing and processing—>data processing<——>data application (drawing) integrated software module. The invention generally solves the problem of integrated application of data collection, data processing and automatic charting in a pumping test.

Description

A Processing Method of Pumping Test Data

technical field

The invention relates to a processing method for pumping test data, which belongs to the field of computer-aided drawing of geological and mineral maps in geological information engineering.

Background technique

Pumping test is an important means of obtaining hydrological information in resource exploration. The comprehensive result map of hydrological pumping test is the most important hydrological map file in mine exploration and mining, and it is the most direct embodiment of pumping test data visualization. It carries important hydrological information. , It plays a decisive role in obtaining the hydrogeological parameters of the aquifer, the hydraulic connection with other aquifers and the surface, timely evaluating the quality of pumping work, and judging the hydrological situation in the research area. Driven by the introduction of GIS visualization technology, database technology and computer automatic mapping technology into multi-disciplinary and multi-channel geological scientific research and the development trend of geosciences in actual production work, how to use large-scale pumping test data for batch data mapping to obtain a complete The result map is the most concerned issue of computer cartography experts and geoscientists, and it has also become an important indicator to measure the technical level of the geological field.

At present, the preprocessing and mapping of collected data are mostly decentralized. The drawing technology for hydrological pumping tests is basically based on semi-manual mode. For example, Liu Yeqing (2010) proposed the method of using Excel to draw pumping test curves, but this method is only limited to generating curves, which is far from the completion of thematic maps. Far. The basic drawing mode of some software is: draw a single chart by Excel, import it into MapGis, and vectorize the raster graphics; use manual map and integrated drawing for the wiring of the curve, and then vectorize it by MapGis, and calculate the hydrological parameters Manually calculated, stored in an Excel table, etc. The software of the domestic geology and mining industry basically does not have the automatic compilation function of the comprehensive result map of the pumping test based on database management, and most of the manual complex operations are still used from data collection to effective data extraction, and the drawing mode of the comprehensive result map of the pumping test is still semi-manual The process of drawing is complex, the work efficiency is low, the amount of calculation is large, the intermediate and result data cannot be automatically stored and managed, and the integrated management, processing and application of data cannot be formed, and finally a set of data collection and data processing cannot be formed. , data management and computer graphics in one integrated software module. It will cause certain difficulties for geological personnel in data storage, analysis, calculation, management and automatic drawing of thematic maps. The drawing automation degree of the comprehensive results map of the pumping test is not high, and the comprehensive automatic drawing technology is not mature enough, and there is a contradiction between efficiency and quality. Therefore, it is very necessary to establish a mechanism from data collection to data conversion and processing, and then to automatic drawing of pumping test data based on database management.

Contents of the invention

In order to solve the deficiencies of the prior art, the present invention provides a method for processing pumping test data, which quickly and efficiently realizes the computer automatic drawing and storage of the comprehensive result graph of the pumping test.

The technical scheme that the present invention adopts for solving its technical problem is: provide a kind of processing method of pumping test data, comprise the following steps:

(1) collect the pumping test data, obtain the pumping test data file produced by the .mdf format of the flow water level logging instrument in the well; the pumping test data file includes static water level observation data file, pumping test data file and recovery water level data file III Class, all kinds of pumping test data files include the basic information of the borehole pumping test and the result data, each line of result data includes time, temperature, water level, flow velocity and depth;

(2) the pumping test data file is preprocessed, and the pumping test data file after the pretreatment is written into the database to form the pumping test data;

(3) The pumping test data in the database and other relevant data of the stored hydrological drilling are retrieved, and the following subgraphs and data are respectively drawn according to the rules of longitudinal distribution under the same project according to the pumping test data:

(a) Set the basic parameters of the histogram and the range of non-target layers, and use the clipping technique to draw a comprehensive hydrological histogram omitting the non-target layers;

(b) Construct the initial date time axis according to the time data of each record in the pumping test data; construct the initial cumulative time axis according to the time interval between each record in the pumping test data and the first record; the time of the first record The world coordinates of the data location are (0, 0);

Using multi-period variable scale and data linkage positioning technology, according to the date at the variable scale, calculate the world coordinates of each variable scale boundary between the date time axis and the cumulative time axis to re-plan the scale distribution of the date time axis and the cumulative time axis, Then dynamically calculate the world coordinates of each data point in the pumping test data according to the size of the variable scale in each period;

Take the date and time axis and the cumulative time axis as the abscissa axis, and use the Q and S ordinate axes as the ordinate axis to draw the Q, S-t and restoration water level duration curves of the pumping test;

(c) Calculate and draw Q-S curve, S-lgt curve, S'-lgt' curve and lgS-lgt curve by computer according to pumping test data; Under the dual coordinate system of the lgS-lgt curve, combined with the performance characteristics of the pumping test data in the curve under the coordinate system, use the graphic element attribute recognition technology to obtain the graphic element attribute value of the first and last scale stored in the curve coordinate system, and then use the mouse to dynamically locate The technology dynamically calculates the actual coordinate values of the current mouse movement position in each coordinate system, and displays them on the screen in real time. After drawing, obtain the characteristic hydrological points, and choose whether to write the coordinate values of the characteristic hydrological points into the database;

(d) According to the pumping test data in the database and the characteristic hydrological points obtained during the mapping process, the general scientific calculation method and iterative recursive algorithm are used to automatically analyze and calculate the relevant hydrogeological parameters and other comprehensive results data, and the hydrogeological parameters and other comprehensive results Data is written to the database;

(4) Move each sub-graph drawn in step (3) to obtain a complete single-hole pumping test comprehensive result map and multi-hole pumping test comprehensive result map.

The preprocessing described in step (2) includes obtaining the pumping test data file, checking the validity of the file and content, analyzing, processing and extracting the sampled data, and automatically writing the pumping test data obtained by sampling into the database.

As described in step (3), use the multi-period variable scale and data linkage positioning technology to calculate the world coordinates of each variable scale boundary between the date time axis and the cumulative time axis according to the date at the variable scale to replan the date time axis and Accumulate the scale distribution of the time axis, and then dynamically calculate the world coordinates of each data point in the pumping test data according to the size of the variable scale in each period. The world coordinates at the boundary of the ith variable scale are calculated using the following formula:

When i=0,

DX _{points 0} = 0

When i>=1,

DX _{points i} = DX _{points i-1} + (T _i -T _i-1 ) × d _i-1 ;

Among them, DX _{point i} is the world coordinate of the current variable scale boundary, DX _{point i-1} is the world coordinate of the previous variable scale boundary, T _i is the date and time of the current variable scale boundary, T _i-1 is the previous The date and time at the boundary of the variable scale, d _i-1 is the distance on the map of each unit hour before the current variable scale boundary; T _i , T _i-1 and d _i-1 are all set values;

The use of data linkage positioning technology to calculate the world coordinates of each data point in the pumping test data, wherein the world coordinates of the i-th data point is calculated using the following formula:

DX _i ＝DX _{min i} +(T _cur -T _i )×d _i

Where i>=1, DX _i is the world coordinate of the current data point; T _cur is the date and time of the current data point; d _i is the distance on the map of each unit hour after the boundary of the current variable scale; T _cur and d _i are Settings.

Combined with the performance characteristics of the pumping test data in the curve under the coordinate system described in step (3), use the graphic element attribute recognition technology to obtain the graphic element attribute value of the first and last scale stored in the curve coordinate system, and then use the mouse dynamic positioning technology to dynamically calculate The actual coordinate value of the current mouse movement position in each coordinate system, including the following process:

Use the primitive attribute recognition technology to obtain the attribute values of the first and last scales of the horizontal and vertical coordinate axes of the coordinate system to determine the positioning range of the feature points. Write attribute values at the leftmost and rightmost ends of the abscissa axis and the uppermost and lower ends of the ordinate axis. The attribute values include S _max , S _min , t _max , t _min , dLX, dRX, dBY and dTY, and then use The mouse dynamic positioning technology dynamically calculates the actual coordinate values of the current mouse movement position in each coordinate system, and draws the coordinate values of the characteristic hydrological points in the S-lgt and S'-lgt single coordinate systems, among which the single coordinate system S-lgt The coordinates (S, t) of the characteristic hydrological points are calculated by the following formula:

The coordinates (S’, t’) of the characteristic hydrological points in the single coordinate system S’-lgt’ are calculated by the following formula:

Among them, dY is the world ordinate of the view area where the mouse is currently located, S _max is the maximum scale of the ordinate axis S of the S-lgt or S'-lgt' coordinate system, and S _min is the S-lgt or S'-lgt' coordinate system The minimum scale of the ordinate axis S, t _max is the maximum scale of the abscissa axis t in the S-lgt or S'-lgt' coordinate system, and t _min is the minimum scale of the abscissa axis t in the S-lgt or S'-lgt' coordinate system Scale, dLX is the world abscissa coordinate at the leftmost end of the t abscissa axis in the S-lgt or S'-lgt' coordinate system, dRX is the world abscissa coordinate at the rightmost end of the t abscissa axis in the S-lgt or S'-lgt' coordinate system, dBY is the world ordinate at the bottom of the S ordinate axis in the S-lgt or S'-lgt' coordinate system, and dTY is the world ordinate at the top of the S ordinate axis in the S-lgt or S'-lgt' coordinate system.

Combined with the performance characteristics of the pumping test data in the curve under the coordinate system described in step (3), use the graphic element attribute recognition technology to obtain the graphic element attribute value of the first and last scale stored in the curve coordinate system, and then use the mouse dynamic positioning technology to dynamically calculate The actual coordinate value of the current mouse movement position in each coordinate system, including the following process:

Use the primitive attribute identification technology to obtain the attribute values of the first and last scales of the horizontal and vertical coordinate axes of the coordinate system to determine the positioning range of the feature points, and draw the lgS-lgt curve and lgW(u)-lg(1/u) curve and their coordinates In the dual-coordinate system, write attribute values to the leftmost and rightmost ends of the abscissa axis and the uppermost and lowermost ends of the ordinate axis respectively in the dual coordinate system. The attribute values include S _max , S _min , t _max , t _min , dLX, dRX, dBY, dTY, W(u) _max , W(u) _min , u _max , u _min , dLX', dRX', dBY' and dTY', and then use the mouse dynamic positioning technology to dynamically calculate the current mouse movement position in each The actual coordinate values in the coordinate system, draw the coordinate values of the characteristic hydrological points in the dual coordinate system of lgS-lgt and standard curve, in which the lgS-lgt coordinate system is combined with the lgW(u)-lg(1/u) coordinate system The coordinates (S, t, W(u), 1/u) of the characteristic hydrological points in the resulting dual coordinate system are calculated by the following formula:

Among them, S _max is the maximum scale of the ordinate axis S of the lgS-lgt coordinate system; S _min is the minimum scale of the ordinate axis S of the lgS-lgt coordinate system; dY is the world ordinate of the view area of the lgS-lgt coordinate system where the current mouse is located ;dTY is the world ordinate at the top of the ordinate axis of the lgS-lgt coordinate system; dBY is the world ordinate at the bottom of the ordinate axis of the lgS-lgt coordinate system;

Among them, t _max is the maximum scale of the abscissa axis t in the lgS-lgt coordinate system; t _min is the minimum scale of the abscissa axis t in the lgS-lgt coordinate system; dX is the world abscissa coordinate of the view area of the lgS-lgt coordinate system where the current mouse is located ;dLX is the world abscissa coordinate at the leftmost end of the abscissa axis in the lgS-lgt coordinate system; dRX is the world abscissa coordinate at the rightmost end of the abscissa axis in the lgS-lgt coordinate system;

W(u) _max is the maximum scale of the ordinate axis W(u) of the lgW(u)-lg(1/u) coordinate system; W(u) _min is the ordinate of the lgW(u)-lg(1/u) coordinate system The minimum scale of the coordinate axis W(u); dY' is the world ordinate of the view area of the lgW(u)-lg(1/u) coordinate system where the current mouse is located; dTY' is lgW(u)-lg(1/u) The world ordinate at the top of the ordinate axis of the coordinate system; dBY' is the world ordinate at the bottom of the ordinate axis of the lgW(u)-lg(1/u) coordinate system;

1/u is the 1/u value displayed by the real-time positioning of the mouse; (1/u) _max is the maximum scale of the abscissa axis 1/u of the lg W(u)-lg(1/u) coordinate system; (1/u) _min is the minimum scale of the abscissa axis 1/u of the lg W(u)-lg(1/u) coordinate system; dX' is the world of the view area of the lg W(u)-lg(1/u) coordinate system where the current mouse is located Abscissa; dLX' is the world abscissa at the leftmost end of the abscissa axis of the lg W(u)-lg(1/u) coordinate system; dRX' is the abscissa axis of the lg W(u)-lg(1/u) coordinate system The rightmost world abscissa.

Relevant hydrogeological parameters and other comprehensive result data described in step (3), including S-lgt line diagram method data table, S'-lgt' line diagram method data table, lgS-lgt line diagram method data table, single and double hole Comprehensive result table of pumping test, water chemical analysis table and construction technical parameter table.

The beneficial effect that the present invention has based on its technical scheme is:

(1) Utilize the processing method of pumping test data provided by the present invention, can utilize computer to realize data automatic processing and redrawing, real-time data storage;

(2) The present invention overcomes a number of key technologies, uses cutting technology to draw a comprehensive hydrological histogram that omits non-target layers, uses multi-period variable scale technology and data linkage positioning technology to unify the time axis, and obtains accurate and easy-to-store data. Timestamp, using primitive attribute identification and mouse dynamic positioning technology to obtain the coordinates of characteristic hydrological points for storage, which can realize the computer intelligently according to the characteristics and constraints of the pumping experiment data, simulate the traditional manual drawing work, and carry out the comprehensive result map of the pumping experiment Automatic drawing by computer, so as to quickly and accurately obtain the final comprehensive result map, which is used to guide the work of hydrogeological personnel;

(3) By the processing method of the pumping test data provided by the present invention, it is possible to implement the pattern data storage of drawing-real-time storage, replace the interaction between manual labor and database, and improve the production accuracy and efficiency;

(4) The operation of redrawing the subgraphs in the comprehensive results graph of the pumping test is flexible, and the new data generated by the redrawing can be dynamically updated to the database in real time.

Description of drawings

Fig. 1 is a schematic diagram of multi-period variable scale and data linkage positioning technology processing in the present invention.

Fig. 2 is a simulation diagram of graphic entity attribute identification and mouse dynamic positioning in the present invention.

Fig. 3 is a schematic flow chart of the present invention.

detailed description

The present invention will be further described below in conjunction with drawings and embodiments.

The present invention provides a kind of processing method of pumping test data, with reference to Fig. 3, comprises the following steps:

(1) collect the pumping test data, obtain the pumping test data file produced by the .mdf format of the flow water level logging instrument in the well; the pumping test data file includes static water level observation data file, pumping test data file and recovery water level data file III Class, all kinds of pumping test data files include the basic information of the borehole pumping test and the result data, each line of result data includes time, temperature, water level, flow velocity and depth;

(2) Through the visualization of the interface, the pumping test data file and file type are selected by the user, and the computer automatically preprocesses the pumping test data file, and writes the preprocessed pumping test data file into the database to form the pumping test data ;Including the acquisition of pumping test data files, checking the validity of files and contents, sampling data analysis, processing and extraction, and automatically writing the pumping test data obtained by sampling into standardized and coded databases with different data models;

The traversal method of the file content adopts line-by-line traversal, and the line content is read in the form of a string. According to the data characteristics of the .mdf file, check the module identification "[Basic Information]", "[Result Data Start]", and the content " Measurement date = XX-XX-XX", "Start time XX: XX: XX", "time, temperature, water level, flow rate", for the format of "measurement date = XX-XX-XX", take an item-by-item inspection method .

The raw time data format read in the basic information is:

Measurement Date = 2013-08-04

start_time = 18:53:23

The data format of the data area is:

Time, temperature, water level, velocity, (depth)

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

99.50, 20.06, 169.29, 32.69, 355.76

104.44, 20.06, 169.31, 32.64, 355.76

109.34, 20.06, 169.33, 32.58, 355.76

114.50, 20.06, 169.36, 32.51, 355.76

The observation time of the target database is a time type, and its format is:

Year-Month-Day Hour:Minute:Second

The time in the data area is a number type, and the unit is second (s), which needs to be calculated and converted. According to user needs, the general measurement start time starts at 30 or 00 minutes. Taking the above data as an example, the first time record after processing is:

2013-08-04 19:00:23

The sampling data interval rule is: 5 pieces of 1 minute, 5 pieces of 2 minutes, 9 pieces of 5 minutes, and some 30 minutes. If the processing time is longer than 2 hours, the records of the same row will be kept. The time of each record finally processed is Data type, which is year-month-day hour:minute:second.

Write the preprocessed pumping test data files into the database;

(3) Retrieve the pumping test data in the database and other relevant data of the stored hydrological boreholes. Other relevant data of hydrological boreholes are the original data of hydrological boreholes collected in the field, and may also include some preprocessed data, including borehole General situation information, geological stratification information, aquifer information, and well logging information, according to the pumping test data, draw the following sub-graphs and data according to the rules of longitudinal distribution under the same project:

(a) Set the basic parameters of the histogram and the range of non-target layers, and use the clipping technique to draw a comprehensive hydrological histogram omitting the non-target layers;

For example, if the data within 100-300m is omitted, the height of the omitted part on the map shall be the value represented by a scale scale on the map (the same as the treatment method for the water level above the surface). The omitted part is marked with a tilde "≈≈" in the middle, and its depth is dSignDepth. The starting depth of the columnar area to be drawn is dbeginDepth, and the ending depth is dendDepth, so the processing method for the omitted part is as follows:

1. When denddepth<100m, draw normally;

2. Otherwise, when denddepth>100m, there are five situations,

2-1. If dbeginDepth<100m and denddepth<300m, when the current columnar area is drawn on the drawing, take the depth dbeginDepth, stop depth dSignDepth, and the pattern code is the pattern code of the current area;

2-2. dbeginDepth<100m and denddepth>300m, the omitted part is distributed in the same columnar area, and the omitted thickness in this area is directly subtracted;

2-3. dbeginDepth>100m and denddepth<300m, this columnar area is distributed in the omitted part and is not drawn;

2-4, 100m<dbeginDepth<300m and denddepth>300m, then the current columnar area is drawn on the drawing with the depth dSignDepth, the stop depth dendDepth, and the pattern code is the pattern code of the current area;

2-5. dbeginDepth>300m and denddepth>300m, the starting and ending depths remain unchanged, and the thickness of the omitted part is subtracted at the same time.

(b) Construct the initial date time axis according to the time data of each record in the pumping test data; construct the initial cumulative time axis according to the time interval between each record in the pumping test data and the first record; the time of the first record The world coordinates of the data location are (0, 0);

Using multi-period variable scale and data linkage positioning technology, according to the date at the variable scale, calculate the world coordinates of each variable scale boundary between the date time axis and the cumulative time axis to re-plan the scale distribution of the date time axis and the cumulative time axis, Then dynamically calculate the world coordinates of each data point in the pumping test data according to the size of the variable scale in each period;

Referring to Figure 1, assume that the first variable scale of the cumulative time axis is a unit of 4h, and the second variable scale is a unit of 12h; for the drawing of the variable scale of the date and time axis, the first variable scale uses a unit of 1h, and the second paragraph The variable scale is 24h one unit, and the third section variable scale is 48h one unit. If the scale is changed on the 2nd, the cumulative time of the scale change corresponding to the cumulative time axis is 18, and the first scale change is 4h before the 18th point, and the second scale change is 12h after 18;

Using multi-period variable scale and data linkage positioning technology, according to the date at the variable scale, calculate the world coordinates of each variable scale boundary between the date time axis and the cumulative time axis to re-plan the scale distribution of the date time axis and the cumulative time axis, Then dynamically calculate the world coordinates of each data point in the pumping test data according to the size of the variable scale in each period, and the world coordinates at the boundary of the ith variable scale are calculated using the following formula:

When i=0,

DX _{points 0} = 0

When i>=1,

DX _{points i} = DX _{points i-1} + (T _i -T _i-1 ) × d _i-1 ;

Among them, DX _{point i} is the world coordinate of the current variable scale boundary, DX _{point i-1} is the world coordinate of the previous variable scale boundary, T _i is the date and time of the current variable scale boundary, T _i-1 is the previous The date and time at the boundary of the variable scale, d _i-1 is the distance on the map of each unit hour before the current variable scale boundary; T _i , T _i-1 and d _i-1 are all set values;

The use of data linkage positioning technology to calculate the world coordinates of each data point in the pumping test data, wherein the world coordinates of the i-th data point is calculated using the following formula:

DX _i ＝DX _{min i} +(T _cur -T _i )×d _i

Where i>=1, DX _i is the world coordinate of the current data point; T _cur is the date and time of the current data point; d _i is the distance on the map of each unit hour after the boundary of the current variable scale; T _cur and d _i are Settings;

Take the date and time axis and the cumulative time axis as the abscissa axis, and use the Q and S ordinate axes as the ordinate axis to draw the Q, S-t and restoration water level duration curves of the pumping test;

(c) Calculate and draw Q-S curve, S-lgt curve, S'-lgt' curve and lgS-lgt curve by computer according to pumping test data; Under the dual coordinate system of the lgS-lgt curve, combined with the performance characteristics of the pumping test data in the curve under the coordinate system, use the graphic element attribute recognition technology to obtain the graphic element attribute value of the first and last scale stored in the curve coordinate system, and then use the mouse to dynamically locate The technology dynamically calculates the actual coordinate values of the current mouse movement position in each coordinate system, and displays them on the screen in real time. After drawing, obtain the characteristic hydrological points, and choose whether to write the coordinate values of the characteristic hydrological points into the database;

Among them, the S-lgt and S'-lgt' coordinate systems are both single coordinate systems. In the process of drawing the curve coordinate system, record the maximum value of the scale of the horizontal and vertical axes of the curve and its world coordinates in the view window as the coordinates Store the attribute values of the system, which are respectively represented by S _max , S _min , t _max , t _min , dLX, dRX, dBY and dTY (generally S _min = 0, t _min = 1), and are dynamically displayed and positioned by moving the mouse (S,t) coordinate value. Mouse positioning requires real-time conversion between mouse view logical coordinates and world coordinates.

Use the primitive attribute recognition technology to obtain the attribute values of the first and last scales of the horizontal and vertical coordinate axes of the coordinate system to determine the positioning range of the feature points. Write attribute values at the leftmost and rightmost ends of the abscissa axis and the uppermost and lower ends of the ordinate axis. The attribute values include S _max , S _min , t _max , t _min , dLX, dRX, dBY and dTY, and then use The mouse dynamic positioning technology dynamically calculates the actual coordinate values of the current mouse movement position in each coordinate system, and draws the coordinate values of the characteristic hydrological points in the S-lgt and S'-lgt single coordinate systems, among which the single coordinate system S-lgt The coordinates (S, t) of the characteristic hydrological points are calculated by the following formula:

The coordinates (S’, t’) of the characteristic hydrological points in the single coordinate system S’-lgt’ are calculated by the following formula:

Among them, dY is the world ordinate of the view area where the mouse is currently located, S _max is the maximum scale of the ordinate axis S of the S-lgt or S'-lgt' coordinate system, and S _min is the S-lgt or S'-lgt' coordinate system The minimum scale of the ordinate axis S, t _max is the maximum scale of the abscissa axis t in the S-lgt or S'-lgt' coordinate system, and t _min is the minimum scale of the abscissa axis t in the S-lgt or S'-lgt' coordinate system Scale, dLX is the world abscissa coordinate at the leftmost end of the t abscissa axis in the S-lgt or S'-lgt' coordinate system, dRX is the world abscissa coordinate at the rightmost end of the t abscissa axis in the S-lgt or S'-lgt' coordinate system, dBY is the world ordinate at the bottom of the S ordinate axis in the S-lgt or S'-lgt' coordinate system, and dTY is the world ordinate at the top of the S ordinate axis in the S-lgt or S'-lgt' coordinate system.

Use the primitive attribute identification technology to obtain the attribute values of the first and last scales of the horizontal and vertical coordinate axes of the coordinate system to determine the positioning range of the feature points, and draw the lgS-lgt curve and lgW(u)-lg(1/u) curve and their coordinates In the dual-coordinate system, write attribute values to the leftmost and rightmost ends of the abscissa axis and the uppermost and lowermost ends of the ordinate axis respectively in the dual coordinate system. The attribute values include S _max , S _min , t _max , t _min , dLX, dRX, dBY, dTY, W(u) _max , W(u) _min , u _max , u _min , dLX', dRX', dBY' and dTY', and then use the mouse dynamic positioning technology to dynamically calculate the current mouse movement position in each The actual coordinate values in the coordinate system, draw the coordinate values of the characteristic hydrological points in the dual coordinate system of lgS-lgt and standard curve, in which the lgS-lgt coordinate system is combined with the lgW(u)-lg(1/u) coordinate system The coordinates (S, t, W(u), 1/u) of the characteristic hydrological points in the resulting dual coordinate system are calculated by the following formula:

Among them, S _max is the maximum scale of the ordinate axis S of the lgS-lgt coordinate system; S _min is the minimum scale of the ordinate axis S of the lgS-lgt coordinate system; dY is the world ordinate of the view area of the lgS-lgt coordinate system where the current mouse is located ;dTY is the world ordinate at the top of the ordinate axis of the lgS-lgt coordinate system; dBY is the world ordinate at the bottom of the ordinate axis of the lgS-lgt coordinate system;

Among them, t _max is the maximum scale of the abscissa axis t in the lgS-lgt coordinate system; t _min is the minimum scale of the abscissa axis t in the lgS-lgt coordinate system; dX is the world abscissa coordinate of the view area of the lgS-lgt coordinate system where the current mouse is located ;dLX is the world abscissa coordinate at the leftmost end of the abscissa axis in the lgS-lgt coordinate system; dRX is the world abscissa coordinate at the rightmost end of the abscissa axis in the lgS-lgt coordinate system;

W(u) _max is the maximum scale of the ordinate axis W(u) of the lgW(u)-lg(1/u) coordinate system; W(u) _min is the ordinate of the lgW(u)-lg(1/u) coordinate system The minimum scale of the coordinate axis W(u); dY' is the world ordinate of the view area of the lgW(u)-lg(1/u) coordinate system where the current mouse is located; dTY' is lgW(u)-lg(1/u) The world ordinate at the top of the ordinate axis of the coordinate system; dBY' is the world ordinate at the bottom of the ordinate axis of the lgW(u)-lg(1/u) coordinate system;

1/u is the 1/u value displayed by the real-time positioning of the mouse; (1/u) _max is the maximum scale of the abscissa axis 1/u of the lg W(u)-lg(1/u) coordinate system; (1/u) _min is the minimum scale of the abscissa axis 1/u of the lg W(u)-lg(1/u) coordinate system; dX' is the world of the view area of the lg W(u)-lg(1/u) coordinate system where the current mouse is located Abscissa; dLX' is the world abscissa at the leftmost end of the abscissa axis of the lg W(u)-lg(1/u) coordinate system; dRX' is the abscissa axis of the lg W(u)-lg(1/u) coordinate system The rightmost world abscissa.

Figure 2 shows the simulation diagram of primitive attribute recognition and mouse dynamic positioning.

(d) According to the pumping test data in the database and the characteristic hydrological points obtained during the mapping process, use general scientific calculation methods and iterative recursive algorithms to automatically analyze and calculate relevant hydrogeological parameters and other comprehensive result data, including S-lgt line diagram method data table, S'-lgt' straight line graphical method data table, lgS-lgt straight line graphical method data table, single and double hole pumping test comprehensive result table, water chemical analysis table and construction technical parameter table, hydrogeological parameters and other comprehensive result data write to the database;

(4) Since each sub-graph is automatically arranged and generated according to the vertical and vertical rules in the engineering view, each sub-graph drawn in step (3) is moved to the frame, that is, the fine-tuning of the position of the graph in the engineering view is obtained to meet the industry standard The complete single-hole pumping test comprehensive result graph and multi-hole pumping test comprehensive result graph.

Claims (3)

1. a processing method of pumping test data, is characterized in that comprising the following steps: (1) collect the pumping test data, obtain the pumping test data file produced by the .mdf format of the flow water level logging instrument in the well; the pumping test data file includes static water level observation data file, pumping test data file and recovery water level data file III Class, all kinds of pumping test data files include the basic information of the borehole pumping test and the result data, each line of result data includes time, temperature, water level, flow velocity and depth; (2) the pumping test data file is preprocessed, and the pumping test data file after the pretreatment is written into the database to form the pumping test data; (3) The pumping test data in the database and other relevant data of the stored hydrological drilling are retrieved, and the following subgraphs and data are respectively drawn according to the rules of longitudinal distribution under the same project according to the pumping test data: (a) Set the basic parameters of the histogram and the range of non-target layers, and use the clipping technique to draw a comprehensive hydrological histogram omitting the non-target layers; (b) Construct the initial date time axis according to the time data of each record in the pumping test data; construct the initial cumulative time axis according to the time interval between each record in the pumping test data and the first record; the time of the first record The world coordinates of the data location are (0, 0); Using multi-period variable scale and data linkage positioning technology, according to the date at the variable scale, calculate the world coordinates of each variable scale boundary between the date time axis and the cumulative time axis to re-plan the scale distribution of the date time axis and the cumulative time axis, Then dynamically calculate the world coordinates of each data point in the pumping test data according to the size of the variable scale in each period, and the world coordinates at the boundary of the ith variable scale are calculated using the following formula: When i=0, DX _{points 0} = 0 When i>=1, DX _{points i} = DX _{points i-1} + (T _i -T _i-1 ) × d _i-1 ; Among them, DX _{point i} is the world coordinate of the current variable scale boundary, DX _{point i-1} is the world coordinate of the previous variable scale boundary, T _i is the date and time of the current variable scale boundary, T _i-1 is the previous The date and time at the boundary of the variable scale, d _i-1 is the distance on the map of each unit hour before the current variable scale boundary; T _i , T _i-1 and d _i-1 are all set values; The use of data linkage positioning technology to calculate the world coordinates of each data point in the pumping test data, wherein the world coordinates of the i-th data point is calculated using the following formula: DX _i ＝DX _{min i} +(T _cur -T _i )×d _i Where i>=1, DX _i is the world coordinate of the current data point; T _cur is the date and time of the current data point; d _i is the distance on the map of each unit hour after the boundary of the current variable scale; T _cur and d _i are Settings; Take the date and time axis and the cumulative time axis as the abscissa axis, and use the Q and S ordinate axes as the ordinate axis to draw the Q, S-t and restoration water level duration curves of the pumping test; (c) Calculate and draw Q-S curve, S-lgt curve, S'-lgt' curve and lgS-lgt curve by computer according to pumping test data; Under the dual coordinate system of the lgS-lgt curve, combined with the performance characteristics of the pumping test data in the curve under the coordinate system, use the graphic element attribute recognition technology to obtain the graphic element attribute value of the first and last scale stored in the curve coordinate system, and then use the mouse to dynamically locate The technology dynamically calculates the actual coordinate values of the current mouse movement position in each coordinate system, and displays them on the screen in real time. After drawing, obtain the characteristic hydrological points, and choose whether to write the coordinate values of the characteristic hydrological points into the database; Use the primitive attribute recognition technology to obtain the attribute values of the first and last scales of the horizontal and vertical coordinate axes of the coordinate system to determine the positioning range of the characteristic hydrological points. When drawing the S-lgt, S'-lgt' curves and their coordinate systems, respectively Write attribute values at the leftmost and rightmost ends of the abscissa axis and the uppermost and lower ends of the ordinate axis. The attribute values include S _max , S _min , t _max , t _min , dLX, dRX, dBY and dTY, and then Use the mouse dynamic positioning technology to dynamically calculate the actual coordinate values of the current mouse movement position in each coordinate system, and draw the coordinate values of the characteristic hydrological points in the S-lgt, S'-lgt' single coordinate system, among which in the single coordinate system S- The coordinates (S, t) of the characteristic hydrological points in lgt are calculated by the following formula: $\mathrm{<span\; class="notranslate">\; S</span>}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; T</span>}\mathrm{<span\; class="notranslate">\; Y</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; Y</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; T</span>}\mathrm{<span\; class="notranslate">\; Y</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; B</span>}\mathrm{<span\; class="notranslate">\; Y</span>}}<span\; class="notranslate">\; \&times;</span><span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; a</span>}\mathrm{<span\; class="notranslate">\; x</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}}$ $\mathrm{<span\; class="notranslate">\; t</span>}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; 10</span>}}^{<span\; class="notranslate">\; (</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; x</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; L</span>}\mathrm{<span\; class="notranslate">\; x</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; R</span>}\mathrm{<span\; class="notranslate">\; x</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; L</span>}\mathrm{<span\; class="notranslate">\; x</span>}}<span\; class="notranslate">\; \&times;</span><span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; lgt</span>}}_{\mathrm{<span\; class="notranslate">\; max</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; lgt</span>}}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; lgt</span>}}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span>}$ The coordinates (S’, t’) of the characteristic hydrological points in the single coordinate system S’-lgt’ are calculated by the following formula: ${\mathrm{<span\; class="notranslate">\; S</span>}}^{<span\; class="notranslate">\; ,</span>}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; max</span>}}<span\; class="notranslate">\; -</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; T</span>}\mathrm{<span\; class="notranslate">\; Y</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; Y</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; T</span>}\mathrm{<span\; class="notranslate">\; Y</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; B</span>}\mathrm{<span\; class="notranslate">\; Y</span>}}<span\; class="notranslate">\; \&times;</span><span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; max</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span>$ ${\mathrm{<span\; class="notranslate">\; t</span>}}^{<span\; class="notranslate">\; ,</span>}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; 10</span>}}^{<span\; class="notranslate">\; (</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; x</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; L</span>}\mathrm{<span\; class="notranslate">\; x</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; R</span>}\mathrm{<span\; class="notranslate">\; x</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; L</span>}\mathrm{<span\; class="notranslate">\; x</span>}}<span\; class="notranslate">\; \&times;</span><span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; lgt</span>}}_{\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; a</span>}\mathrm{<span\; class="notranslate">\; x</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; lgt</span>}}_{\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; lgt</span>}}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span>}$ Among them, dY is the world ordinate of the view area where the mouse is currently located, S _max is the maximum scale of the ordinate axis S of the S-lgt or S'-lgt' coordinate system, and S _min is the S-lgt or S'-lgt' coordinate system The minimum scale of the ordinate axis S, t _max is the maximum scale of the abscissa axis t in the S-lgt or S'-lgt' coordinate system, and t _min is the minimum scale of the abscissa axis t in the S-lgt or S'-lgt' coordinate system Scale, dLX is the world abscissa coordinate at the leftmost end of the t abscissa axis in the S-lgt or S'-lgt' coordinate system, dRX is the world abscissa coordinate at the rightmost end of the t abscissa axis in the S-lgt or S'-lgt' coordinate system, dBY is the world ordinate at the bottom of the S ordinate axis in the S-lgt or S'-lgt' coordinate system, dTY is the world ordinate at the top of the S ordinate axis in the S-lgt or S'-lgt' coordinate system; Utilize the primitive attribute identification technology to obtain the attribute values of the first and last scales of the horizontal and vertical coordinate axes of the coordinate system to determine the positioning range of the characteristic hydrological points, and draw the lgS-lgt curve and lgW(u)-lg(1/u) curve and its In the coordinate system, write attribute values to the leftmost and rightmost ends of the abscissa axis and the uppermost and lowermost ends of the ordinate axis in the dual coordinate system respectively. The attribute values include S _max , S _min , t _max , t _min , dLX , dRX, dBY, dTY, W(u) _max , W(u) _min , u _max , u _min , dLX', dRX', dBY' and dTY', and then use the mouse dynamic positioning technology to dynamically calculate the current mouse movement position in The actual coordinate values in each coordinate system, draw the coordinate values of the characteristic hydrological points in the double coordinate system of lgS-lgt and standard curve, in which the lgS-lgt coordinate system and the lgW(u)-lg(1/u) coordinate system The coordinates (S,t,W(u),1/u) of the characteristic hydrological points in the combined dual coordinate system are calculated by the following formula: $\mathrm{<span\; class="notranslate">\; S</span>}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; 10</span>}}^{<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; lm\; w</span>}}_{\mathrm{<span\; class="notranslate">\; max</span>}}<span\; class="notranslate">\; -</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; T</span>}\mathrm{<span\; class="notranslate">\; Y</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; Y</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; T</span>}\mathrm{<span\; class="notranslate">\; Y</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; B</span>}\mathrm{<span\; class="notranslate">\; Y</span>}}<span\; class="notranslate">\; \&times;</span><span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; lm\; w</span>}}_{\mathrm{<span\; class="notranslate">\; max</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; lm\; w</span>}}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; )</span>}$ Among them, S _max is the maximum scale of the ordinate axis S of the lgS-lgt coordinate system; S _min is the minimum scale of the ordinate axis S of the lgS-lgt coordinate system; dY is the world ordinate of the view area of the lgS-lgt coordinate system where the current mouse is located ;dTY is the world ordinate at the top of the ordinate axis of the lgS-lgt coordinate system; dBY is the world ordinate at the bottom of the ordinate axis of the lgS-lgt coordinate system; $\mathrm{<span\; class="notranslate">\; t</span>}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; 10</span>}}^{<span\; class="notranslate">\; (</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; x</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; L</span>}\mathrm{<span\; class="notranslate">\; x</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; R</span>}\mathrm{<span\; class="notranslate">\; x</span>}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; L</span>}\mathrm{<span\; class="notranslate">\; x</span>}}<span\; class="notranslate">\; \&times;</span><span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; lgt</span>}}_{\mathrm{<span\; class="notranslate">\; max</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; lgt</span>}}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; lgt</span>}}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span>}$ Among them, t _max is the maximum scale of the abscissa axis t in the lgS-lgt coordinate system; t _min is the minimum scale of the abscissa axis t in the lgS-lgt coordinate system; dX is the world abscissa coordinate of the view area of the lgS-lgt coordinate system where the current mouse is located ;dLX is the world abscissa coordinate at the leftmost end of the abscissa axis in the lgS-lgt coordinate system; dRX is the world abscissa coordinate at the rightmost end of the abscissa axis in the lgS-lgt coordinate system; $\mathrm{<span\; class="notranslate">\; W</span>}<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; u</span>}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; 10</span>}}^{<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; lg</span>}\mathrm{<span\; class="notranslate">\; W</span>}{<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; u</span>}<span\; class="notranslate">\; )</span>}_{\mathrm{<span\; class="notranslate">\; max</span>}}<span\; class="notranslate">\; -</span>\frac{{\mathrm{<span\; class="notranslate">\; wxya</span>}}^{<span\; class="notranslate">\; ,</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; wxya</span>}}^{<span\; class="notranslate">\; ,</span>}}{{\mathrm{<span\; class="notranslate">\; wxya</span>}}^{<span\; class="notranslate">\; ,</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; dBY</span>}}^{<span\; class="notranslate">\; ,</span>}}<span\; class="notranslate">\; \&times;</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; lg</span>}\mathrm{<span\; class="notranslate">\; W</span>}{<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; u</span>}<span\; class="notranslate">\; )</span>}_{\mathrm{<span\; class="notranslate">\; max</span>}}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; lg</span>}\mathrm{<span\; class="notranslate">\; W</span>}{<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; u</span>}<span\; class="notranslate">\; )</span>}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; )</span>}$ W(u) _max is the maximum scale of the ordinate axis W(u) of the lgW(u)-lg(1/u) coordinate system; W(u) _min is the ordinate of the lgW(u)-lg(1/u) coordinate system The minimum scale of the coordinate axis W(u); dY' is the world ordinate of the view area of the lgW(u)-lg(1/u) coordinate system where the current mouse is located; dTY' is lgW(u)-lg(1/u) The world ordinate at the top of the ordinate axis of the coordinate system; dBY' is the world ordinate at the bottom of the ordinate axis of the lgW(u)-lg(1/u) coordinate system; $\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; /</span>\mathrm{<span\; class="notranslate">\; u</span>}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; 10</span>}}^{<span\; class="notranslate">\; (</span>\frac{{\mathrm{<span\; class="notranslate">\; wxya</span>}}^{<span\; class="notranslate">\; ,</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; wxya</span>}}^{<span\; class="notranslate">\; ,</span>}}{{\mathrm{<span\; class="notranslate">\; wxya</span>}}^{<span\; class="notranslate">\; ,</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; wxya</span>}}^{<span\; class="notranslate">\; ,</span>}}<span\; class="notranslate">\; \&times;</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; lg</span>}{<span\; class="notranslate">\; (</span>\frac{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; u</span>}}<span\; class="notranslate">\; )</span>}_{\mathrm{<span\; class="notranslate">\; max</span>}}<span\; class="notranslate">\; -</span>\mathrm{<span\; class="notranslate">\; lg</span>}{<span\; class="notranslate">\; (</span>\frac{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; u</span>}}<span\; class="notranslate">\; )</span>}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; lg</span>}{<span\; class="notranslate">\; (</span>\frac{\mathrm{<span\; class="notranslate">\; 1</span>}}{\mathrm{<span\; class="notranslate">\; u</span>}}<span\; class="notranslate">\; )</span>}_{\mathrm{<span\; class="notranslate">\; min</span>}}<span\; class="notranslate">\; )</span>}$ 1/u is the 1/u value displayed by the real-time positioning of the mouse; (1/u) _max is the maximum scale of the abscissa axis 1/u of the lg W(u)-lg(1/u) coordinate system; (1/u) _min is the minimum scale of the abscissa axis 1/u of the lg W(u)-lg(1/u) coordinate system; dX' is the world of the view area of the lg W(u)-lg(1/u) coordinate system where the current mouse is located Abscissa; dLX' is the world abscissa at the leftmost end of the abscissa axis of the lg W(u)-lg(1/u) coordinate system; dRX' is the abscissa axis of the lg W(u)-lg(1/u) coordinate system The world abscissa at the far right; (d) According to the pumping test data in the database and the characteristic hydrological points obtained during the mapping process, the general scientific calculation method and iterative recursive algorithm are used to automatically analyze and calculate the relevant hydrogeological parameters and other comprehensive results data, and the hydrogeological parameters and other comprehensive results Data is written to the database; (4) Move each sub-graph drawn in step (3) to obtain a complete single-hole pumping test comprehensive result map and multi-hole pumping test comprehensive result map. 2. the processing method of pumping test data according to claim 1, is characterized in that: the described pretreatment of step (2), comprises to pumping test data file acquisition, file and content validity inspection, sampling data analysis, processing And extract, and automatically write the pumping test data obtained by sampling into the database. 3. the processing method of pumping test data according to claim 1, is characterized in that: relevant hydrogeological parameter and other comprehensive achievement data described in step (3), comprise S-lgt linear graph method data table, S'- lgt' straight line graphical method data table, lgS-lgt straight line graphical method data table, single and double hole pumping test comprehensive results table, water chemical analysis table and construction technical parameter table.