[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN103853882B - Drawing method and system for engineering geological profile - Google Patents

Drawing method and system for engineering geological profile Download PDF

Info

Publication number
CN103853882B
CN103853882B CN201410050271.8A CN201410050271A CN103853882B CN 103853882 B CN103853882 B CN 103853882B CN 201410050271 A CN201410050271 A CN 201410050271A CN 103853882 B CN103853882 B CN 103853882B
Authority
CN
China
Prior art keywords
cross
section
stratum
ground
exploration
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410050271.8A
Other languages
Chinese (zh)
Other versions
CN103853882A (en
Inventor
黄新文
钱国玉
易菊香
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Railway Engineering Consulting Group Co Ltd
Original Assignee
China Railway Engineering Consulting Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Railway Engineering Consulting Group Co Ltd filed Critical China Railway Engineering Consulting Group Co Ltd
Priority to CN201410050271.8A priority Critical patent/CN103853882B/en
Publication of CN103853882A publication Critical patent/CN103853882A/en
Application granted granted Critical
Publication of CN103853882B publication Critical patent/CN103853882B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)

Abstract

The invention discloses a drawing method and a drawing system for an engineering geological profile. The drawing method comprises the following steps: drawing the engineering geological profile which extends along the direction of a road be to constructed; drawing ground lines in an engineering geological cross profile for a cross section of the road to be constructed according to imported elevations and horizontal coordinates of a plurality of ground test points on the ground of the cross section; determining an X-coordinate of the cross section in the engineering geological profile according to the imported geographic position data of the cross section; selecting an exploration body with the difference value between the X-coordinate of a ground exploration point and the X-coordinate of the cross section less than a preset threshold value; determining the depth of each stratum at the cross section according to the depth of each stratum in the selected exploration body; drawing each stratum line in the engineering geological cross profile according to the ground lines and each stratum depth at the cross section. According to the technical scheme of the invention, the stratum lines are not required to be drawn by technical staff, the workload is greatly reduced, and the work efficiency is improved.

Description

Method and system for drawing engineering geological profile
Technical Field
The invention relates to the field of engineering geological information processing, in particular to a method and a system for drawing an engineering geological profile.
Background
With the continuous development of infrastructure construction in China, a large number of roads are to be constructed. Prior to the construction of roads, the design units of roads are usually surveyed at a predetermined construction site of the road (e.g. railway, highway): the engineering personnel drill a plurality of exploration entities (such as static exploration holes, dynamic exploration holes, drilling holes or test pits and the like) in a predetermined construction site; each exploration entity extends vertically downward from the surface to a set depth. For each exploration entity, after measuring the ground elevation and horizontal plane coordinates of a ground exploration point of the exploration entity, engineering personnel acquire a rock-soil body within the volume range of the exploration entity; determining an interface between every two strata in the obtained rock-soil body so as to detect geological data of each stratum in the exploration entity; geological data of the formation may include: geological age, geological cause, rock name, humidity, plasticity, weathering degree, crushing degree and the like; regarding the interface between each stratum, taking the depth of the interface as the depth of the stratum with the interface as the bottom surface; thereby obtaining the depths of the various strata of the exploration entity.
The engineering personnel inputs the ground elevation and horizontal plane coordinates of the ground exploration points of each exploration entity, and the depth and geological data of each stratum in the exploration entity into engineering geological exploration software.
The engineering geological survey software summarizes geological data of each stratum in each exploration entity, divides the stratum with the same or similar geological data into the same stratum (for example, the same stratum number), and takes each divided stratum as each summarized stratum; for each exploration entity, if each stratum in the exploration entity is compared with each stratum summarized, and one or more strata are missing, one or more virtual strata are correspondingly supplemented for the exploration entity aiming at the missing stratum, the thickness of the supplemented virtual stratum is set to be 0, namely, the depth of the added virtual stratum is the same as the depth of the stratum which is the last layer of the virtual stratum.
Then, the engineering geological survey software draws a ground line in an engineering geological profile (as shown in figure 1) along the extending direction of the road to be built according to the horizontal plane coordinates and the ground elevation of ground exploration points of all exploration entities; according to the depth of each stratum (including virtual stratum) in each exploration entity, drawing a stratum line of each stratum in the engineering geological profile along the extending direction of the road to be built, wherein the stratum line represents the bottom surface of each stratum: in the engineering geological profile, an X axis in a rectangular coordinate system is parallel to the extending direction of a road to be built, and a Y axis represents the depth of a stratum.
However, in practical applications, it is often necessary to have engineering geological profiles of several cross sections of the road to be built, herein referred to as engineering geological cross sections for the sake of description.
At present, most exploration entities are generally distributed along the extending direction of a road to be built, namely few exploration entities are distributed in the direction vertical to the road to be built, so that engineering geological survey software cannot draw an engineering geological profile along the extending direction of the road to be built by a method for drawing the engineering geological profile; the engineering geological survey software cannot draw an engineering geological cross section according to an engineering geological section along the extending direction of the road to be built.
Currently, the process of drawing an engineered geological cross-section map of a determined cross-section by a skilled person typically comprises: for a cross section determined in a road to be built, a technician inputs horizontal plane coordinates and elevations of all ground test points of the cross section into engineering geological survey software; the engineering geological exploration software draws a ground line according to the input data; for each stratum of the cross section, a technician manually translates the copied ground line vertically downwards in engineering geological survey software according to the depth of the stratum to obtain a translated ground line; a technician manually smoothens the translated ground line in engineering geological survey software to obtain a stratum line of the stratum; thereby resulting in an engineered geological cross-section of this cross-section as shown in figure 2. In the engineering geological cross section, an X axis in a rectangular coordinate system is vertical to the extending direction of a road to be built, and a Y axis represents the depth of a stratum. And after the engineering geological cross section map of each preset cross section in the road to be built is obtained, printing and outputting the engineering geological cross section map and the engineering geological section map along the extending direction of the road to be built.
In practical applications, the number of the engineering geological cross sections of a road to be built is often large (for example, tens), and the existing method for drawing the engineering geological cross sections of the road to be built needs to manually draw a large number of formation lines (for example, hundreds), so that the workload is large, the time consumption is long, and the efficiency is low. Therefore, it is necessary to provide a method for drawing an engineering geological profile, so that technicians do not need to manually draw each horizon line in the drawing process of the engineering geological profile, the workload is reduced, and the working efficiency is improved.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a method and a system for drawing an engineering geological profile, so that technicians do not need to manually draw each horizon line in the drawing process of the engineering geological profile, the workload is reduced, and the working efficiency is improved.
According to one aspect, the invention provides a method for drawing an engineering geological profile, which comprises the following steps: drawing an engineering geological profile extending along a road to be built according to the recorded ground elevations and horizontal plane coordinates of ground exploration points of a plurality of exploration entities, and the depth and geological data of each stratum in each exploration entity; wherein the exploration entity is arranged on a construction ground of a road to be constructed; the drawing method further comprises the following steps:
and for one cross section of the road to be built, further drawing an engineering geological cross section of the cross section according to the drawn engineering geological profile:
drawing a ground line in the engineering geological cross section according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of the cross section; and are
Determining the X coordinate of the cross section in the engineering geological profile according to the recorded geographic position data of the cross section; selecting exploration entities of which the difference value between the X coordinate of the ground exploration point and the X coordinate of the cross section is smaller than a set threshold value; determining the depths of all stratums at the cross section according to the selected depths of all stratums in the exploration entity; and drawing each stratum line in the engineering geological cross section according to the ground line and each stratum depth at the cross section.
Preferably, the drawing a ground line in the engineering geological cross section according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of the cross section includes:
according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of the cross section, determining the coordinates of the rectangular coordinate system of each ground test point of the cross section in the engineering geological cross section; the X axis of a rectangular coordinate system in the engineering geological cross section is vertical to the extending direction of the road to be built, and the Y axis represents the stratum depth;
and drawing a ground line in the engineering geological profile according to the coordinates of the ground test points of the cross section in the rectangular coordinate system in the engineering geological cross section.
Preferably, the drawing each horizon in the engineering geological cross section according to the horizon and each horizon depth at the cross section specifically includes:
for each stratum at the cross section, subtracting the depth of the stratum at the cross section from the Y coordinate of the rectangular coordinate system of the ground test point in the engineering geological cross section for each ground test point of the cross section to obtain the coordinate of the virtual test point of the stratum at the cross section, which corresponds to the ground test point;
and for each stratum at the cross section, drawing a stratum line of the stratum in the engineering geological cross section after smoothing the coordinates of each virtual test point of the stratum at the cross section.
Preferably, the smoothing of the coordinates of each virtual test point of the formation at the cross section specifically includes:
fitting a curve equation with a set order according to the coordinates of each virtual test point of the stratum at the cross section;
and substituting the abscissa of each virtual test point of the stratum at the cross section into the curve equation to obtain the ordinate of the smoothed virtual test point.
Preferably, the set order is 4; and fitting a curve equation with a set order according to the coordinates of each virtual test point of the stratum at the cross section, wherein the fitting specifically comprises the following steps:
determining a constant coefficient a of 0 rd to 3 rd order unknowns in the curve equation according to the following formula 40、a1、a2、a3
Wherein p represents the pth formation from top to bottom at the cross-section, (X)01,Yp1)、…、(X0i,Ypi)、…、(X0n,Ypn) And respectively representing the coordinates of n virtual test points of the p-th stratum, wherein p is a natural number, n is the total number of the ground test points of the cross section, and i is a natural number not greater than n.
Further, after the step of drawing the stratigraphic line of the stratum in the engineering geological cross section, the method further comprises the following steps:
numbering each stratum at the cross section to obtain the stratum number of each stratum at the cross section;
and for each stratum at the cross section, correspondingly displaying the stratum number of the stratum near the stratum line of the stratum.
According to another aspect of the present invention, there is provided a system for mapping an engineering geological profile, comprising: the device comprises a survey data processing module and a CAD module; wherein,
the survey data processing module is used for determining the coordinates of the ground exploration points of the exploration entities in a rectangular coordinate system of an engineering geological profile along the extending direction of a road to be built according to the recorded ground elevations and horizontal plane coordinates of the ground exploration points of the exploration entities; summarizing the input geological data of each stratum in each exploration entity; for each stratum obtained by summarizing, determining the depth of the stratum in each exploration entity; then, sending the coordinates of the ground exploration points of the exploration entities in a rectangular coordinate system of the engineering geological profile map and the summarized depth of each stratum in the exploration entities to the CAD module;
the CAD module is used for drawing the engineering geological profile according to the received data; and
the reconnaissance data processing module is further used for determining coordinates of each ground test point in a rectangular coordinate system of the engineering geological cross section according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of one cross section of the road to be built; in a rectangular coordinate system of the engineering geological cross section, an X axis is vertical to the extending direction of the road to be built, and a Y axis represents the depth of the stratum; the survey data processing module is also used for determining the X coordinate of the cross section in the engineering geological profile according to the recorded geographic position data of the cross section; selecting exploration entities of which the difference value between the X coordinate of the ground exploration point and the X coordinate of the cross section is smaller than a set threshold value; determining the depths of all stratums at the cross section according to the selected depths of all stratums in the exploration entity; for each stratum at the cross section, subtracting the depth of the stratum at the cross section from the Y coordinate of the rectangular coordinate system of the ground test point in the engineering geological cross section for each ground test point of the cross section to obtain the coordinate of the virtual test point of the stratum at the cross section, which corresponds to the ground test point; then, sending the coordinates of the ground test points in a rectangular coordinate system of the engineering geological cross section and the coordinates of the virtual test points of the ground at the cross section to the CAD module;
the CAD module is also used for drawing a ground line in the engineering geological cross section according to the received coordinates of each ground test point in the rectangular coordinate system of the engineering geological cross section; and drawing each stratum line in the engineering geological cross section according to the received coordinates of the virtual test points of each stratum at the cross section.
Preferably, the reconnaissance data processing module is further configured to, before the sending the coordinates of the ground test points in the rectangular coordinate system of the engineering geological cross-section map and the coordinates of the virtual test points of the strata at the cross-section to the CAD module, smooth the coordinates of the virtual test points of each stratum for each stratum at the cross-section: fitting a curve equation with a set order according to the coordinates of each virtual test point of the stratum at the cross section; and substituting the abscissa of each virtual test point of the stratum at the cross section into the curve equation to obtain the ordinate of the smoothed virtual test point.
Preferably, the survey data processing module specifically includes:
the coordinate conversion unit is used for determining the coordinates of the ground exploration points of the exploration entities in a rectangular coordinate system of an engineering geological profile along the extending direction of a road to be built according to the recorded ground elevations and horizontal plane coordinates of the ground exploration points of the exploration entities; determining coordinates of each ground test point in a rectangular coordinate system of the engineering geological cross section according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of one cross section of the road to be built;
the stratum depth determining unit is used for summarizing the input geological data of each stratum in each exploration entity; for each stratum obtained by summarizing, determining the depth of the stratum in each exploration entity; determining the X coordinate of the cross section in the engineering geological profile according to the recorded geographic position data of the cross section; selecting exploration entities of which the difference value between the X coordinate of the ground exploration point and the X coordinate of the cross section is smaller than a set threshold value; determining the depths of all stratums at the cross section according to the selected depths of all stratums in the exploration entity;
a virtual test point coordinate determining unit, configured to, for each ground test point of the cross section, subtract, for each ground test point of the cross section, the Y coordinate of the rectangular coordinate system of the ground test point in the engineering geological cross-sectional view, which is determined by the coordinate converting unit, from the depth of the ground at the cross section, which is determined by the ground depth determining unit, and obtain a coordinate of the virtual test point of the ground at the cross section, which corresponds to the ground test point;
the curve equation fitting unit is used for fitting a curve equation with a set order according to the coordinates of each virtual test point of the stratum at the cross section, which are obtained by the virtual test point coordinate determination unit;
the coordinate smoothing processing unit is used for substituting the abscissa of each virtual test point determined by the virtual test point coordinate determining unit into the curve equation determined by the curve equation fitting unit to obtain the ordinate of each virtual test point after smoothing processing; obtaining the coordinates of each virtual test point of the stratum at the cross section after smoothing treatment;
the data output unit is used for sending the coordinates of the ground exploration points of the exploration entities determined by the coordinate conversion unit in a rectangular coordinate system of the engineering geological profile and the summarized depth of each stratum determined by the stratum depth determination unit in the exploration entities to the CAD module; and sending the coordinates of the ground test points determined by the coordinate conversion unit in a rectangular coordinate system of the engineering geological cross section and the coordinates of the virtual test points of the ground layers on the cross section after the coordinates are smoothed by the coordinate smoothing unit to the CAD module.
According to the technical scheme, after the engineering geological profile is drawn along the extending direction of the road to be built according to the exploration data of an exploration entity, for a cross section of the road to be built, the depth of each stratum at the cross section can be determined according to the position coordinates of the cross section in the engineering geological profile, and therefore the stratum lines of each stratum at the cross section can be drawn. That is to say, in the process of drawing the engineering geology cross section map of the cross section, technicians do not need to manually draw each horizon line, the workload is greatly reduced, and the working efficiency is improved.
Drawings
FIG. 1 is a schematic representation of a prior art engineered geological profile along the direction of elongation of a roadway to be constructed;
FIG. 2 is a schematic representation of an engineered geological cross-section of a prior art road to be constructed;
FIG. 3 is a block diagram of the internal structure of a system for mapping an engineered geological profile according to an embodiment of the present invention;
FIG. 4 is a flow chart of a method for drawing an engineering geological cross section of a road to be constructed by the system for drawing an engineering geological profile according to the embodiment of the invention;
FIG. 5 is a schematic flow chart of a method for smoothing the coordinates of the virtual test points of each stratum at the cross section by the system for drawing the engineering geological profile according to the embodiment of the present invention;
FIG. 6 is a block diagram of the internal structure of the survey data processing module according to the embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail below with reference to the accompanying drawings by way of examples of preferred embodiments. It should be noted, however, that the numerous details set forth in the description are merely for the purpose of providing the reader with a thorough understanding of one or more aspects of the present invention, which may be practiced without these specific details.
As used in this application, the terms "module," "system," and the like are intended to include a computer-related entity, such as but not limited to hardware, firmware, a combination of hardware and software, or software in execution. For example, a module may be, but is not limited to: a process running on a processor, an object, an executable, a thread of execution, a program, and/or a computer. For example, an application running on a computing device and the computing device may both be a module. One or more modules may reside within a process and/or thread of execution and a module may be localized on one computer and/or distributed between two or more computers.
The inventor of the invention considers that on the basis of the existing drawing method of the engineering geological profile, if the method exists, the depth of each stratum at the cross section can be automatically determined for the cross section of the profile to be drawn in the road to be built; outputting the determined depth of each stratum at the cross section to CAD software; CAD software can automatically draw each layer line of the cross section according to the received data; thereby reducing the workload of technical personnel and improving the working efficiency.
The technical solution of the embodiments of the present invention is described in detail below with reference to the accompanying drawings.
An internal structure block diagram of the system for drawing an engineering geological profile provided by the embodiment of the invention is shown in fig. 3, and includes: a survey data processing module 301 and a CAD module 302.
The method for drawing the engineering geological profile comprises the following steps: the method comprises the steps of drawing an engineering geological profile along the extending direction of a road to be built and drawing engineering geological cross sections of a plurality of cross sections in the road to be built.
Before drawing an engineering geological profile along the extending direction of a road to be built, an engineer inputs the ground elevations and horizontal plane coordinates of ground exploration points of a plurality of exploration entities along the road to be built, and the depth and geological data of each stratum in each exploration entity into an exploration data processing module 301 in a drawing system of the engineering geological profile.
The survey data processing module 301 determines the coordinates of the ground exploration points of each exploration entity in a rectangular coordinate system of an engineering geological profile extending along a road to be built according to the recorded ground elevations and horizontal plane coordinates of the ground exploration points of each exploration entity; the X axis (namely the transverse axis) of the rectangular coordinate system of the engineering geological profile is parallel to the extending direction of the road to be built, and the Y axis (namely the longitudinal axis) represents the depth of the stratum.
The survey data processing module 301 summarizes the geological data of each stratum in each exploration entity, which is entered: the stratum with the same or similar geological data in each exploration entity is divided into the same stratum (for example, the stratum is numbered as the same stratum number), and each divided stratum is used as each stratum for summarizing, so that each stratum is summarized. For each stratum obtained in the summary, determining the depth of the stratum in each exploration entity: for example, for a exploration entity, if one or more strata are missing in each stratum of the exploration entity compared with the summarized strata, one or more virtual strata are correspondingly supplemented to the exploration entity for the missing strata, and the thickness of the supplemented virtual strata is set to be 0, that is, the depth of the added virtual strata is the same as the depth of the stratum which is the last of the virtual strata. Thereby obtaining the aggregated depth of each formation in each exploration entity.
Then, the survey data processing module 301 sends the coordinates of the ground exploration points of the exploration entities in the rectangular coordinate system of the engineering geological profile and the summarized depth of each stratum in the exploration entities to the CAD module 302 in the drawing system of the engineering geological profile.
The CAD module 302 performs engineering geological profile mapping based on the received data: the CAD module 302 draws a ground line along the extending direction of the road to be built according to the received coordinates of the ground exploration points of the exploration entities in the rectangular coordinate system of the engineering geological profile; for each stratum in the summary, a stratum line between the stratum and the next stratum, namely the stratum bottom line of the stratum, is drawn according to the ground line and the depth of the stratum in each exploration entity. So as to obtain an engineering geological profile along the extension direction of the road to be built: in the engineering geological profile, an X axis in a rectangular coordinate system is parallel to the extending direction of a road to be built, and a Y axis represents the depth of a stratum.
Based on an engineering geological profile along the extending direction of a road to be built, a method for drawing an engineering geological cross section of a cross section for one cross section in a plurality of cross sections of the road to be built by a drawing system of the engineering geological profile is shown in a flow chart of fig. 4, and comprises the following steps:
s401: technical personnel input the geographic position data of a cross section of a road to be built, and the horizontal plane coordinates and elevations of each ground test point of the cross section into a drawing system of an engineering geological profile.
Specifically, for a cross section of an engineering geological cross section to be drawn of a road to be built, a technician inputs geographic position data (for example, mileage along the extending direction of the road to be built) of the cross section into a survey data processing module 301 in a drawing system of the engineering geological cross section; horizontal plane coordinates and elevations of a plurality of ground test points selected on the ground of the cross section in a Beijing 54 coordinate system are recorded into the survey data processing module 301.
Preferably, the technician may enter the horizontal plane coordinates of the intersection point of the cross section of the road to be constructed and the road center line in the beijing 54 coordinate system into the reconnaissance data processing module 301.
S402: and the drawing system of the engineering geological profile map determines the coordinates of each ground test point of the cross section in the rectangular coordinate system of the cross section for the horizontal plane coordinates and elevations of each ground test point of the cross section of the recorded road to be built.
Specifically, for the cross section of the road to be built, which is related in step S401, the survey data processing module 301 determines, for each ground test point of the cross section, according to the horizontal plane coordinate and the elevation of the recorded ground test point in the beijing 54 coordinate system, the abscissa and the ordinate of the ground test point in the rectangular coordinate system of the cross section respectively; in a rectangular coordinate system of the cross section, an X axis (namely a horizontal axis) is vertical to the extending direction of the road to be built, and a Y axis (namely a vertical axis) represents the depth of the stratum.
For example, the number of the ground test points of the cross section of the road to be built related in the step S401 is n, where n is a natural number; the coordinates of the n ground test points determined by the survey data processing module 301 can be respectively expressed as (X)01,Y01)、…、(X0i,Y0i)、…、(X0n,Y0n) (ii) a 0 represents a ground layer, and i is a natural number not greater than n.
Preferably, for the cross section of the road to be built involved in step S401, the survey data processing module 301 determines the abscissa of the intersection point in the rectangular coordinate system of the cross section according to the horizontal plane coordinates of the intersection point of the cross section and the center line of the road in the coordinate system of beijing 54.
S403: and the system for drawing the engineering geological profile map determines the X coordinate (namely the abscissa) of the cross section in the engineering geological profile map along the extending direction of the road to be built according to the recorded geographic position data of the cross section, and selects the exploration entity of which the X coordinate is adjacent to the X coordinate of the cross section.
Specifically, for the cross section of the road to be built involved in step S401, the survey data processing module 301 determines, according to the entered geographic position data of the cross section (for example, the mileage along the extending direction of the road to be built), the X coordinate of the cross section in the engineering geological profile along the extending direction of the road to be built; for each exploration entity, calculating an absolute value obtained by subtracting the X coordinate of the cross section in the engineering geological profile from the X coordinate of the ground exploration point of the exploration entity in the engineering geological profile, and taking the calculated absolute value as a difference value of the X coordinate of the exploration entity and the X coordinate of the cross section; and selecting exploration entities of which the difference value between the X coordinate and the X coordinate of the cross section is smaller than a set threshold value. The set threshold may be a difference in the corresponding X-coordinates of a geographic threshold distance (e.g., a distance of 150 meters) in the engineered geological profile along the direction of travel of the roadway to be constructed.
Preferably, if the survey data processing module 301 determines that the number of the selected exploration entities is not greater than two, the step S404 is continuously executed; if the number of the selected exploration entities is judged to be more than two, after the process can be finished, the method for drawing the engineering geological profile along the extending direction of the road to be built in the embodiment of the invention is adopted to draw the engineering geological cross section of the cross section, and the specific method is not repeated here.
S404: and the drawing system of the engineering geological profile map determines the depth of each stratum at the cross section involved in the step S401 according to the selected depth of each stratum in the exploration entity.
Specifically, for the cross section of the road to be built related to step S401, the survey data processing module 301 may determine the depth of each stratum in the exploration entity selected in step S403 from the depths of each stratum in each exploration entity determined by the survey data processing module 301; the depth of each formation at the cross-section may then be determined in a number of ways.
For example, for the cross section of the road to be built involved in step S401, for one exploration entity in the exploration entities selected in step S403, the depth of each stratum in the exploration entity may be respectively used as the depth of each stratum at the cross section.
For another example, the exploration data processing module 301 obtains each stratigraphic line between the exploration entities selected in step S403 in the engineering geological profile along the extending direction of the road to be built; for the cross section of the road to be built related to the step S401, after coordinates of intersection points between the cross section and each obtained floor line are respectively determined, depths of each floor at the cross section are determined according to the determined coordinates of each intersection point.
S405: and the system for drawing the engineering geological profile map determines the coordinates of the virtual test points of each stratum at the cross section according to the coordinates of each ground test point of the cross section involved in the step S401 and the depth of each stratum at the cross section.
In this step, for each stratum at the cross section of the road to be built, which is referred to in step S401, the reconnaissance data processing module 301 subtracts, for each ground test point of the cross section, the Y coordinate (i.e., the vertical coordinate) of the rectangular coordinate system of the ground test point in the engineering geological cross section of the cross section from the depth of the stratum at the cross section (i.e., the depth of the ground test point translating downward into the stratum), so as to obtain the coordinates of the virtual test point of the stratum at the cross section, which correspond to the ground test point.
For example, the survey data processing module 301 determines, for the cross section of the road to be constructed involved in step S401, the p-th stratum from top to bottom at the cross section, where p is a natural number; if the determined depth of the p-th stratum is hpThen the coordinates (X) of n ground test points can be determined01,Y01)、…、(X0i,Y0i)、…、(X0n,Y0n) The depth h of the p-th stratum is subtracted from the Y coordinate of (1)pObtaining the coordinates of each virtual test point of the p-th stratum at the cross section as (X)01,Y01-hp)、…、(X0i,Y0i-hp)、…、(X0n,Y0n-hp) (ii) a The coordinates of these virtual test points can be respectively noted as (X)01,Yp1)、…、(X0i,Ypi)、…、(X0n,Ypn)。
Preferably, the survey data processing module 301 performs a smoothing process on the coordinates of each virtual test point of each stratum at the cross section of the road to be built, which is referred to in step S401, and a specific method will be described in detail later.
S406: and the system for drawing the engineering geological profile maps draws the engineering geological cross-section maps for the cross sections of the roads to be built, which are related in the step S401, according to the coordinates of the ground test points in the rectangular coordinate system of the cross section and the coordinates of the virtual test points of the ground layers at the cross section.
Specifically, for the cross section of the road to be built, which is referred to in step S401, the reconnaissance data processing module 301 outputs the coordinates of the ground test points in the rectangular coordinate system of the cross section and the coordinates of the virtual test points of the ground layers at the cross section to the CAD module 302 in the drawing system of the engineering geological profile; drawing an engineering geological cross section of the cross section by the CAD module 302 according to the received data; in the engineering geological cross section, the X axis is vertical to the extending direction of the road to be built, and the Y axis represents the depth of the stratum.
The process of drawing the engineering geological cross section by the CAD module 302 may specifically include: for the cross section of the road to be built, which is related in the step S401, after the ground test points are generated according to the coordinates of the ground test points in the rectangular coordinate system of the cross section, every two adjacent ground test points are connected by adopting a straight line segment to form a ground line of the cross section; for each stratum at the cross section, generating each virtual test point of the stratum according to the coordinate of each virtual test point of the stratum, and connecting every two adjacent virtual test points by adopting a straight line segment to form a stratum line of the stratum at the cross section; thereby obtaining an engineered geological cross section of the cross section.
Preferably, the CAD module 302 may further draw a stratigraphic line of each stratigraphic layer at the cross section of the road to be built, which is related in step S401, according to the coordinates of each virtual test point of the stratigraphic layer after smoothing processing; compared with the stratum line drawn according to the coordinates of the virtual test points without smoothing, the method is smoother and smoother, and better conforms to the variation trend of the actual stratum.
In fact, the reconnaissance data processing module 301 may also output the coordinates of the intersection point of the cross section of the road to be created, which is referred to in step S401, and the center line of the road to be created, to the CAD module 302; the CAD module 302 can draw a plumb line passing through a point in the engineering geological cross section map of the cross section after drawing the point according to the received coordinates of the intersection point. The reconnaissance data processing module 301 may also number each stratum at the cross section of the road to be built, which is related to step S401, to obtain stratum numbers of each stratum at the cross section, and output the stratum numbers to the CAD module 302; the CAD module 302 can correspondingly display the received stratigraphic number of the stratigraphic layer near the stratigraphic line of the stratigraphic layer for each stratigraphic layer at the cross-section in the engineered geological cross-section map of the cross-section.
Describing the method for the reconnaissance data processing module 301 in the mapping system of the engineering geological profile in step S405 in detail, for each stratum at the cross section of the road to be built involved in step S401, smoothing the coordinates of each virtual test point of the stratum, as shown in fig. 5, the flow may include the following steps:
s501: the survey data processing module 301 fits a curve equation of a set order to the coordinates of each virtual test point of the formation at the cross section involved in step S401 for each formation at the cross section.
Specifically, the survey data processing module 301 may assume a curve equation with a set order; for each stratum at the cross section of the road to be built related to the step S401, calculating the sum of distances from each virtual test point of the stratum to a curve represented by an assumed curve equation, namely the sum of squares of deviations, according to the coordinates of each virtual test point of the stratum; and determining the value of the constant coefficient of each unknown number in each order in the assumed curve equation by solving the extreme value of the deviation sum of squares, and then determining the curve equation with the set order.
For example, the survey data processing module 301 assumes a 4 th order curve equation, which can be expressed by the following equation 1:
Y=a0+a1X+a2X2+a3X3(formula 1)
In formula 1, a0、a1、a2、a3Respectively assuming constant coefficients of 0 th order unknown number to 3 rd order unknown number in a curve equation; x and Y represent the abscissa and ordinate, respectively, of any one point on the curve.
The survey data processing module 301 calculates, for the p-th formation at the cross section of the road to be built involved in step S401, a deviation sum of squares from each virtual test point of the p-th formation to the curve represented by the curve equation of the above formula 1 according to the coordinates of each virtual test point of the formation, where the deviation sum of squares can be represented by the following expression 2:
when the reconnaissance data processing module 301 makes expression 2 reach an extreme value (e.g., a minimum value), the following formula 3 may be determined:
in formula 3, j is a natural number not greater than n, and k is a non-negative integer from 0 to 3.
The survey data processing module 301 may rewrite equation 3 in a matrix form as equation 4 below:
the reconnaissance data processing module 301 may calculate a constant coefficient a of each unknown order in the assumed curve equation according to formula 40、a1、a2、a3A value of (d); the 4 th order curve equation can thus be determined.
S502: for each stratum at the cross section involved in step S401, the reconnaissance data processing module 301 substitutes the abscissa of each virtual test point of the stratum into the fitted curve equation, and then obtains the ordinate of each virtual test point after smoothing.
Specifically, for each stratum at the cross section involved in step S401, the reconnaissance data processing module 301 substitutes the abscissa of each virtual test point of the stratum at the cross section into the fitted curve equation to obtain the smoothed ordinate of each virtual test point; and obtaining the coordinates of each virtual test point of the stratum at the cross section after smoothing treatment.
For example, the survey data processing module 301 maps the abscissa X of the n virtual test points of the p-th formation with respect to the p-th formation at the cross section involved in step S40101To X0nRespectively substituting the n-th order curve equation determined in the step S501 into the 4-th order curve equation to respectively obtain the vertical coordinates of the n virtual test points of the p-th stratum after smoothing; and obtaining the coordinates of the n virtual test points of the p-th stratum after smoothing treatment.
Further, as shown in fig. 6, the internal structure block diagram of the survey data processing module 301 specifically includes: the system comprises a coordinate conversion unit 601, a formation depth determination unit 602, a virtual test point coordinate determination unit 603, a curve equation fitting unit 604, a coordinate smoothing processing unit 605 and a data output unit 606.
The coordinate conversion unit 601 is used for determining the coordinates of the ground exploration points of each exploration entity in a rectangular coordinate system of an engineering geological profile along the extending direction of a road to be built according to the recorded ground elevation and horizontal plane coordinates of the ground exploration points of each exploration entity; and determining the coordinates of each ground test point in the rectangular coordinate system of the engineering geological cross section according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of one cross section of the road to be built. And then, outputting a notice of finishing the coordinate conversion.
The stratum depth determining unit 602 is configured to summarize the entered geological data of each stratum in each exploration entity; for each stratum obtained by summarizing, determining the depth of the stratum in each exploration entity; determining the X coordinate of the cross section in the engineering geological profile according to the recorded geographic position data of the cross section; selecting exploration entities of which the difference value between the X coordinate of the ground exploration point and the X coordinate of the cross section is smaller than a set threshold value; and determining the depths of all the stratums at the cross section according to the selected depths of all the stratums in the exploration entity. And then, outputting a notice of completion of formation depth determination.
The virtual test point coordinate determining unit 603 is configured to, after receiving the coordinate conversion completion notification and the formation depth determination completion notification, subtract, for each ground test point of the cross section, the Y coordinate of the ground test point determined by the coordinate converting unit 601 in the orthogonal coordinate system in the engineering geological cross-sectional view from the depth of the ground test point at the cross section determined by the formation depth determining unit 602, and obtain the coordinate of the virtual test point of the ground test point of the cross section corresponding to the ground test point. After that, the virtual test point coordinate determination unit 603 sends a smoothing process notification to the curve equation fitting unit 604.
The curve equation fitting unit 604 is configured to fit a curve equation of a set order according to the coordinates of the virtual test points of the formation at the cross section determined by the virtual test point coordinate determining unit 603 after receiving the smoothing notification.
The coordinate smoothing unit 605 is configured to substitute, for each virtual test point of the formation at the cross section, the abscissa of the virtual test point determined by the virtual test point coordinate determining unit 603 into the curve equation determined by the curve equation fitting unit 604, and obtain the ordinate of the virtual test point after smoothing; and obtaining the coordinates of each virtual test point of the stratum at the cross section after the smoothing treatment. After that, a data output notification is sent to the data output unit 606.
The data output unit 606 is configured to, after receiving the data output notification, send the coordinates of the ground exploration points of the exploration entities determined by the coordinate conversion unit 601 in the rectangular coordinate system of the engineering geological profile and the summarized depth of each stratum in the exploration entities determined by the stratum depth determination unit 602 to the CAD module 302; the coordinates of the ground test points determined by the coordinate conversion unit 601 in the rectangular coordinate system of the engineering geological cross section and the coordinates of the virtual test points of the ground layers at the cross section determined by the virtual test point coordinate determination unit 603 are sent to the CAD module 302.
The specific implementation method of the functions of the coordinate conversion unit 601, the formation depth determination unit 602, the virtual test point coordinate determination unit 603, the curve equation fitting unit 604, the coordinate smoothing processing unit 605, and the data output unit 606 may refer to the specific content of the method flow shown in fig. 5, and is not described herein again.
According to the technical scheme, after the engineering geological profile is drawn along the extending direction of the road to be built according to the exploration data of an exploration entity, for a cross section of the road to be built, the depth of each stratum at the cross section can be determined according to the position coordinates of the cross section in the engineering geological profile, and therefore the stratum lines of each stratum at the cross section can be drawn. That is to say, in the process of drawing the engineering geology cross section map of the cross section, technicians do not need to manually draw each horizon line, the workload is greatly reduced, and the working efficiency is improved.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that those skilled in the art can make various improvements and modifications without departing from the principle of the present invention, and these improvements and modifications should also be construed as the protection scope of the present invention.

Claims (9)

1. A method of mapping an engineered geological profile, comprising: drawing an engineering geological profile extending along a road to be built according to the recorded ground elevations and horizontal plane coordinates of ground exploration points of a plurality of exploration entities, and the depth and geological data of each stratum in each exploration entity; wherein the exploration entity is arranged on a construction ground of a road to be constructed; it is characterized by also comprising:
and for one cross section of the road to be built, further drawing an engineering geological cross section of the cross section according to the drawn engineering geological profile:
drawing a ground line in the engineering geological cross section according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of the cross section; and are
Determining the X coordinate of the cross section in the engineering geological profile according to the recorded geographic position data of the cross section; selecting exploration entities of which the difference value between the X coordinate of the ground exploration point and the X coordinate of the cross section is smaller than a set threshold value; determining the depths of all stratums at the cross section according to the selected depths of all stratums in the exploration entity; and drawing each stratum line in the engineering geological cross section according to the ground line and each stratum depth at the cross section.
2. The method according to claim 1, wherein the step of drawing the ground line in the cross-sectional profile of the engineering geology according to the recorded elevation and horizontal plane coordinates of a plurality of selected ground test points on the ground of the cross-section comprises the following steps:
according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of the cross section, determining the coordinates of the rectangular coordinate system of each ground test point of the cross section in the engineering geological cross section; the X axis of a rectangular coordinate system in the engineering geological cross section is vertical to the extending direction of the road to be built, and the Y axis represents the stratum depth;
and drawing a ground line in the engineering geological profile according to the coordinates of the ground test points of the cross section in the rectangular coordinate system in the engineering geological cross section.
3. The method of claim 2, wherein said mapping each horizon in the engineered geological cross-sectional profile based on the horizon and each depth of horizon at the cross-section comprises:
for each stratum at the cross section, subtracting the depth of the stratum at the cross section from the Y coordinate of the rectangular coordinate system of the ground test point in the engineering geological cross section for each ground test point of the cross section to obtain the coordinate of the virtual test point of the stratum at the cross section, which corresponds to the ground test point;
and for each stratum at the cross section, drawing a stratum line of the stratum in the engineering geological cross section after smoothing the coordinates of each virtual test point of the stratum at the cross section.
4. The method of claim 3, wherein smoothing the coordinates of each virtual test point of the formation at the cross-section comprises:
fitting a curve equation with a set order according to the coordinates of each virtual test point of the stratum at the cross section;
and substituting the abscissa of each virtual test point of the stratum at the cross section into the curve equation to obtain the ordinate of the smoothed virtual test point.
5. The method of claim 4, wherein the set order is 4; and fitting a curve equation with a set order according to the coordinates of each virtual test point of the stratum at the cross section, wherein the fitting specifically comprises the following steps:
determining a constant coefficient a of 0 rd to 3 rd order unknowns in the curve equation according to the following formula 40、a1、a2、a3
Wherein p represents the pth formation from top to bottom at the cross-section, (X)01,Yp1)、…、(X0i,Ypi)、…、(X0n,Ypn) Respectively representing the coordinates of n virtual test points of the p-th stratum, wherein p is a natural number, and n isAnd the total number i of the ground test points of the cross section is a natural number not greater than n.
6. The method of any one of claims 1-5, further comprising, after said mapping each horizon in said engineered geological cross-section:
numbering each stratum at the cross section to obtain the stratum number of each stratum at the cross section;
and for each stratum at the cross section, correspondingly displaying the stratum number of the stratum near the stratum line of the stratum.
7. A system for mapping an engineered geological profile, comprising: the device comprises a survey data processing module and a CAD module; wherein,
the survey data processing module is used for determining the coordinates of the ground exploration points of the exploration entities in a rectangular coordinate system of an engineering geological profile along the extending direction of a road to be built according to the recorded ground elevations and horizontal plane coordinates of the ground exploration points of the exploration entities; summarizing the input geological data of each stratum in each exploration entity; for each stratum obtained by summarizing, determining the depth of the stratum in each exploration entity; then, sending the coordinates of the ground exploration points of the exploration entities in a rectangular coordinate system of the engineering geological profile and the summarized depth of each stratum in the exploration entities to the CAD module;
the CAD module is used for drawing the engineering geological profile according to the received data; and
the reconnaissance data processing module is further used for determining coordinates of each ground test point in a rectangular coordinate system of the engineering geological cross section according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of one cross section of the road to be built; in a rectangular coordinate system of the engineering geological cross section, an X axis is vertical to the extending direction of the road to be built, and a Y axis represents the depth of the stratum; the survey data processing module is also used for determining the X coordinate of the cross section in the engineering geological profile according to the recorded geographic position data of the cross section; selecting exploration entities of which the difference value between the X coordinate of the ground exploration point and the X coordinate of the cross section is smaller than a set threshold value; determining the depths of all stratums at the cross section according to the selected depths of all stratums in the exploration entity; for each stratum at the cross section, subtracting the depth of the stratum at the cross section from the Y coordinate of the rectangular coordinate system of the ground test point in the engineering geological cross section for each ground test point of the cross section to obtain the coordinate of the virtual test point of the stratum at the cross section, which corresponds to the ground test point; then, sending the coordinates of the ground test points in the rectangular coordinate system of the engineering geological cross section and the coordinates of the virtual test points of the ground at the cross section to the CAD module;
the CAD module is also used for drawing a ground line in the engineering geological cross section according to the received coordinates of each ground test point in the rectangular coordinate system of the engineering geological cross section; and drawing each stratum line in the engineering geological cross section according to the received coordinates of the virtual test points of each stratum at the cross section.
8. The system of claim 7,
the reconnaissance data processing module is further configured to, before the coordinates of the ground test points in the rectangular coordinate system of the engineering geological cross section and the coordinates of the virtual test points of the strata at the cross section are sent to the CAD module, smooth the coordinates of the virtual test points of each stratum for each stratum at the cross section: fitting a curve equation with a set order according to the coordinates of each virtual test point of the stratum at the cross section; and substituting the abscissa of each virtual test point of the stratum at the cross section into the curve equation to obtain the ordinate of the smoothed virtual test point.
9. The system of claim 8, wherein the survey data processing module specifically comprises:
the coordinate conversion unit is used for determining the coordinates of the ground exploration points of the exploration entities in a rectangular coordinate system of an engineering geological profile along the extending direction of a road to be built according to the recorded ground elevations and horizontal plane coordinates of the ground exploration points of the exploration entities; determining coordinates of each ground test point in a rectangular coordinate system of the engineering geological cross section according to the recorded elevation and horizontal plane coordinates of a plurality of ground test points selected on the ground of one cross section of the road to be built;
the stratum depth determining unit is used for summarizing the input geological data of each stratum in each exploration entity; for each stratum obtained by summarizing, determining the depth of the stratum in each exploration entity; determining the X coordinate of the cross section in the engineering geological profile according to the recorded geographic position data of the cross section; selecting exploration entities of which the difference value between the X coordinate of the ground exploration point and the X coordinate of the cross section is smaller than a set threshold value; determining the depths of all stratums at the cross section according to the selected depths of all stratums in the exploration entity;
a virtual test point coordinate determining unit, configured to, for each ground test point of the cross section, subtract, for each ground test point of the cross section, the Y coordinate of the rectangular coordinate system of the ground test point in the engineering geological cross-sectional view, which is determined by the coordinate converting unit, from the depth of the ground at the cross section, which is determined by the ground depth determining unit, and obtain a coordinate of the virtual test point of the ground at the cross section, which corresponds to the ground test point;
the curve equation fitting unit is used for fitting a curve equation with a set order according to the coordinates of each virtual test point of the stratum at the cross section, which are obtained by the virtual test point coordinate determination unit;
the coordinate smoothing processing unit is used for substituting the abscissa of each virtual test point determined by the virtual test point coordinate determining unit into the curve equation determined by the curve equation fitting unit to obtain the ordinate of each virtual test point after smoothing processing; obtaining the coordinates of each virtual test point of the stratum at the cross section after smoothing treatment;
the data output unit is used for sending the coordinates of the ground exploration points of the exploration entities determined by the coordinate conversion unit in a rectangular coordinate system of the engineering geological profile and the summarized depth of each stratum in the exploration entities determined by the stratum depth determination unit to the CAD module; and sending the coordinates of the ground test points determined by the coordinate conversion unit in a rectangular coordinate system of the engineering geological cross section and the coordinates of the virtual test points of the ground layers on the cross section after the coordinates are smoothed by the coordinate smoothing unit to the CAD module.
CN201410050271.8A 2014-02-13 2014-02-13 Drawing method and system for engineering geological profile Active CN103853882B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410050271.8A CN103853882B (en) 2014-02-13 2014-02-13 Drawing method and system for engineering geological profile

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410050271.8A CN103853882B (en) 2014-02-13 2014-02-13 Drawing method and system for engineering geological profile

Publications (2)

Publication Number Publication Date
CN103853882A CN103853882A (en) 2014-06-11
CN103853882B true CN103853882B (en) 2017-02-22

Family

ID=50861533

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410050271.8A Active CN103853882B (en) 2014-02-13 2014-02-13 Drawing method and system for engineering geological profile

Country Status (1)

Country Link
CN (1) CN103853882B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105740574B (en) * 2016-03-03 2019-02-01 中国石油大学(北京) A kind of ancient landform section restoration methods
CN112819919B (en) * 2021-01-29 2022-03-18 北京城建勘测设计研究院有限责任公司 Profile map full-flow drawing method applied to geotechnical engineering investigation industry
CN115984497B (en) * 2022-12-30 2024-02-23 中国铁路设计集团有限公司 Automatic filling method for geological cross section

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05212999A (en) * 1992-02-05 1993-08-24 Japan Servo Co Ltd Recorder
CN102103758A (en) * 2011-03-21 2011-06-22 四川大学 Computer graph generating method for engineering geological section map
CN102360404A (en) * 2011-11-01 2012-02-22 中国水电顾问集团华东勘测设计研究院 Method for compiling geological display map of chamber

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05212999A (en) * 1992-02-05 1993-08-24 Japan Servo Co Ltd Recorder
CN102103758A (en) * 2011-03-21 2011-06-22 四川大学 Computer graph generating method for engineering geological section map
CN102360404A (en) * 2011-11-01 2012-02-22 中国水电顾问集团华东勘测设计研究院 Method for compiling geological display map of chamber

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于AutoCAD开发纵横断面辅助生成软件及其应用;袁兆奎;《北京测绘》;20120530(第5期);第34--39,103页 *
工程地质岩性剖面图绘图软件的研制;欧湘萍;《武汉水运工程学院学报》;19900630;第14卷(第2期);第199-204页 *

Also Published As

Publication number Publication date
CN103853882A (en) 2014-06-11

Similar Documents

Publication Publication Date Title
CN107742015B (en) DC induced polarization method three-dimensional numerical simulation method based on arbitrary dipole-dipole device
US10921960B2 (en) Mapping geological features using facies-type data manipulation
AU2013402201B2 (en) In-situ wellbore, core and cuttings information system
CN105631935A (en) BIM-based three-dimensional geological automatic modeling method
CN108333627B (en) Igneous rock area is broken recognition methods and the device of the true and false
CN103853882B (en) Drawing method and system for engineering geological profile
CN105911592B (en) A kind of real three dimensional seismic data long wavelength's static correcting method for boring constraint
CN103527172B (en) Variable rock electric coupling index water saturation calculation method
CN104459800A (en) Sand body pinch-out prediction method and device
CN104570070A (en) Method and equipment for building two-dimensional near-surface geological model
WO2019056536A1 (en) Wellbore cross section shape measurement method, device and terminal
CN105507891B (en) Method and device for obtaining anisotropic coefficient of formation resistivity of highly deviated well
CN111208558B (en) Method and device for establishing ultra-deep low-amplitude three-dimensional geological structure
CN103615239A (en) Logging formation imaging system and method
CA2965867C (en) Determining non-linear petrofacies using cross-plot partitioning
US10180516B2 (en) Method of determining channelway trajectories
Esaki et al. A GIS-based prediction method to evaluate subsidence-induced damage from coal mining beneath a reservoir, Kyushu, Japan
CN108242069B (en) Oil reservoir profile compiling and drawing method based on seismic profile and time-depth relation
CN104408755A (en) Geologic plane graph drawing method based on structural plane attitude and geologic observation spot coordinate
CN115984501A (en) Geological three-dimensional model establishing method, device, equipment and storage medium
CN107942400B (en) Method and device for predicting drilling risk
CN113658335A (en) Method for creating complex geological model of karst development area
CN107015289A (en) Trap evaluation stock number determines method and device
CN100337125C (en) Method for processing data of well logging to calculate stratum attitude in drilling well
CN114114390B (en) Oil and gas exploration method, device, equipment and storage medium

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant