CN115296752A - Mud pulse data coding and transmitting method, device and equipment - Google Patents
Mud pulse data coding and transmitting method, device and equipment Download PDFInfo
- Publication number
- CN115296752A CN115296752A CN202210922404.0A CN202210922404A CN115296752A CN 115296752 A CN115296752 A CN 115296752A CN 202210922404 A CN202210922404 A CN 202210922404A CN 115296752 A CN115296752 A CN 115296752A
- Authority
- CN
- China
- Prior art keywords
- dimension
- data
- similarity
- time slots
- coordinate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 83
- 238000005259 measurement Methods 0.000 claims abstract description 108
- 230000005540 biological transmission Effects 0.000 claims abstract description 93
- 238000004364 calculation method Methods 0.000 claims description 30
- 238000007405 data analysis Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 11
- 238000005553 drilling Methods 0.000 claims description 4
- 108091026890 Coding region Proteins 0.000 abstract description 8
- 238000007906 compression Methods 0.000 abstract description 2
- 230000006835 compression Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 10
- 229910002056 binary alloy Inorganic materials 0.000 description 5
- 238000013459 approach Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 210000002445 nipple Anatomy 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000002002 slurry Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000035485 pulse pressure Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B14/00—Transmission systems not characterised by the medium used for transmission
- H04B14/02—Transmission systems not characterised by the medium used for transmission characterised by the use of pulse modulation
- H04B14/04—Transmission systems not characterised by the medium used for transmission characterised by the use of pulse modulation using pulse code modulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/138—Devices entrained in the flow of well-bore fluid for transmitting data, control or actuation signals
Landscapes
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Signal Processing (AREA)
- Remote Sensing (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method, a device and equipment for encoding and transmitting mud pulse data, wherein the encoding method comprises the following steps: determining the measurement value and the total number to be coded; determining a decimal number range corresponding to each dimension; arranging the measurement values on each dimension in sequence from small to large; determining the number x of first dimension time slots theoretically required for each dimension 0 (ii) a Determining the similarity H and the number y of similarity time slots; the number x of second dimension time slots required by a certain dimension i with dimension coordinate value difference in two adjacent transmission data is determined i Determining the number x of time slots in the first dimension 0 (ii) a The coding sequence is determined as: number of similarity time slots y plus number of second dimension time slots x i In combination with each other. And the transmission method completes data transmission after coding according to the coding method. The coding mode of the invention is simple, and the similarity index is introduced to the transmission numberAccording to the compression, the transmission rate is improved.
Description
Technical Field
The invention relates to the technical field of underground data transmission in the field of oil exploration, in particular to a mud pulse data encoding method based on a multi-dimensional space coordinate, a mud pulse data encoding device based on the multi-dimensional space coordinate, a mud pulse data transmission method based on the multi-dimensional space coordinate, a mud pulse data transmission device based on the multi-dimensional space coordinate, computer equipment for realizing the encoding method and/or the transmission method and a computer readable storage medium for realizing the encoding method and/or the transmission method.
Background
The common MWD and LWD instrument coding modes at home and abroad mainly comprise three modes of pulse position modulation coding, manchester coding and optimized combination code. The pulse position modulation coding takes time intervals as data stream transmission information, the transmission time of the coding mode can be increased along with the increase of a measured value, the required time period can be greatly increased when more parameters need to be measured and the measured value is larger, the transmission rate is changed constantly, and when interference occurs in the transmission process, the position where data dislocation occurs is not easy to check, so that the pulse position modulation coding is not suitable for the transmission of a large amount of data. The invention discloses a method for encoding and transmitting slurry pulse data in 1 month and 7 days in 2020, and Chinese patent document with publication number CN110661580A describes a method for encoding slurry pulse data, which comprises the following steps: determining the number of pressure pulse combinations and the number of binary digits according to the range of data to be transmitted and the required precision; carrying out binary system conversion on the binary system, and segmenting the converted binary system into time segments according to digits, wherein each segment comprises time slots with different numbers; and respectively representing the numerical value of each digit by utilizing pressure pulse combinations with different amplitudes, carrying out permutation and combination coding on the amplitudes, the pulse pressure number and the positions corresponding to the pulses, establishing a coding protocol for the used numerical value and the time slot position according to different amplitudes, and finishing mud pulse data coding.
The Manchester code coding is a synchronous clock coding technology, uses level jump to represent the coding of 1 or 0, and has the advantages of simple change rule, high signal transmission rate and large transmitted data volume. However, the encoding method can only be used in the electromagnetic wave or continuous wave method, and the current domestic electromagnetic wave and continuous wave technology is not mature and is greatly limited in the field application process.
The optimized combined code is the most common mud MWD coding mode at present, and the coding mode has the advantages that the time length of data transmission does not change along with the change of binary values, so that whether signal pulses are lost or not is convenient to detect, but each numerical value of the coding mode needs to be combined by applying one pulse, the combination of a coding library is complex, the transmission rate is only 1-2bps, the problem of data repeated transmission is not considered, and the waste of a transmission channel is caused.
The stable transmission of MWD signal can be realized through analyzing above-mentioned three kinds of data encoding modes, but receive mud pulse generator and ground decoding system's restriction during the transmission, transmission rate can't satisfy data transmission demands such as engineering parameter measurement nipple joint.
Disclosure of Invention
The present invention aims to address at least one of the above-mentioned deficiencies of the prior art. For example, one of the purposes of the invention is to provide a mud pulse data encoding method with a transmission rate capable of meeting data transmission requirements of engineering parameter measurement nipples and the like.
In order to achieve the above object, an aspect of the present invention provides a mud pulse data encoding method based on multidimensional space coordinates, including: determining the measurement value and the total number M of the value of the underground measurement data to be coded according to the measurement range S and the precision K of the underground measurement data; determining a decimal number range G corresponding to each dimensionality according to an arithmetic n-th square root of the total number M of the numerical values, wherein n is the total number of the dimensionalities of the multidimensional space coordinate; sequentially arranging all measurement values of the underground measurement data on each dimension from small to large so that each multi-dimensional space coordinate corresponds to one measurement value; converting the decimal number range corresponding to each dimension into binary, and determining the number x of the first dimension time slots theoretically required by each dimension 0 (ii) a Comparing the coordinate values of all dimensions of the data transmitted twice, determining the similarity H, and determining the number y of similarity time slots according to the similarity; based on the similarity H, the number x of second dimension time slots required by a certain dimension i with dimension coordinate value difference in two adjacent transmission data is determined i Determining the number x of time slots in the first dimension 0 The number x of the third dimension time slots needed by a certain dimension j without dimension coordinate value difference in the two adjacent times of transmission data j Is determined to be 0; determining an encoding sequence of downhole measurement data as: number of similarity time slots y plus number of second dimension time slots x i The combination form of the mud pulse data coding is completed; wherein the value of H is 1-2 n I ranges from 0 to n, and j ranges from 0 to n.
In one exemplary embodiment of the multi-dimensional spatial coordinate-based mud pulse data encoding method of the present invention, the downhole measurement data may include directional parameters including well deviation, azimuth, and toolface, and drilling engineering parameters including weight on bit, torque, temperature, pressure, vibration, and torque.
In an exemplary embodiment of the method for encoding mud pulse data based on multidimensional space coordinates of the present invention, the total number of values may be calculated by: m = S/K, wherein M is the total number of numerical values of the downhole measurement data to be encoded, S is the measurement range of the downhole measurement data, and K is the precision of the downhole measurement data.
In an exemplary embodiment of the mud pulse data encoding method based on the multidimensional space coordinate, a value obtained by rounding the nth root of the total number M of values up may be determined as the decimal range G corresponding to each dimension.
In one exemplary embodiment of the multi-dimensional spatial coordinate-based mud pulse data encoding method of the present invention, n is 3 ≦ n ≦ 5.
In an exemplary embodiment of the method for encoding mud pulse data based on multi-dimensional space coordinates of the present invention, the calculation formula of the first-dimension time slot number may be:
wherein x is 0 And G is the decimal number range corresponding to each dimension.
In an exemplary embodiment of the mud pulse data encoding method based on the multidimensional space coordinate, the number y of the similarity time slots may be determined as a total number n of dimensions of the multidimensional space coordinate.
The invention also provides a mud pulse data transmission method based on the multidimensional space coordinate, which is used for transmitting the mud pulse data after coding the downhole measurement data according to the coding method.
The invention further provides a mud pulse data coding device based on multi-dimensional space coordinates, which comprises a measured data analysis unit, a dimension range determination unit, a coordinate matching unit, a first calculation unit, a similarity determination unit, a second calculation unit, a third calculation unit and a coding unit, wherein the measured data analysis unit is configured to determine a measured value and a total number M of the measured values of downhole measured data to be coded according to a measured range S and a precision K of the downhole measured data; dimension range determination sheetThe element is configured to determine a decimal number range corresponding to each dimensionality according to an arithmetic n-th-order root of the total number M of the numerical values, wherein n is the total number of the dimensionalities of the multidimensional space coordinate; the coordinate matching unit is respectively connected with the measured data analysis unit and the dimension range determination unit, and is configured to establish multi-dimensional space coordinates, and arrange the measured values of the underground measured data on each dimension in sequence from small to large so that each multi-dimensional space coordinate corresponds to one measured value; the first calculation unit is connected with the dimension range determination unit and is configured to convert the decimal number range corresponding to each dimension into binary, and determine the number x of the first dimension time slots theoretically required by each dimension 0 (ii) a The similarity determining unit is configured to compare coordinate values of each dimension of the two adjacent transmission data to determine similarity H; the second calculation unit is connected with the similarity determination unit and is configured to determine the number y of the similarity time slots according to the similarity; the third calculating unit is respectively connected with the similarity determining unit and the first calculating unit and is configured to calculate the number x of second dimension time slots required by a certain dimension i with dimension coordinate value difference in two adjacent transmission data based on the similarity H i Determining as a first dimension time slot number x 0 (ii) a The coding unit is respectively connected with the second calculation unit and the third calculation unit and is configured to add the second dimension time slot number x to the similarity time slot number y i The combined form of (a) forms a coded sequence of downhole measurement data.
The invention further provides a mud pulse data transmission device based on multi-dimensional space coordinates, which comprises the coding device and a transmission unit, wherein the transmission unit is connected with the coding unit and is configured to carry out mud pulse data transmission after coding the downhole measurement data according to the coding sequence.
Yet another aspect of the present invention provides a computer apparatus, the apparatus comprising: a processor; a memory storing a computer program which, when executed by the processor, implements at least one of the encoding method as described above, the transmission method as described above.
A further aspect of the invention provides a computer-readable storage medium storing a computer program which, when executed by a processor, implements at least one of the encoding method as described above, the transmission method as described above.
Compared with the prior art, the beneficial effects of the invention comprise at least one of the following:
(1) The coding mode of the invention is simple, all measurement parameter values only need to be numbered from small to large according to the coordinates in sequence, and the precision of different measurement data can be displayed by controlling the number of dimensions and the range of each dimension;
(2) According to the method, the measurement value is divided into a plurality of areas by utilizing dimensionality, each dimensionality can determine a measurement value range, and the measurement value range gradually approaches to a real pressure value along with the increase of the dimensionality;
(3) The data similarity index is introduced, and in the measuring process, because the data generally cannot be suddenly changed, the similarity of the data transmitted for the second time and the data transmitted for the first time is generally high, and the data transmitted by utilizing the characteristic is compressed, so that the transmission rate is obviously improved compared with the data which is not compressed.
Drawings
The above and other objects and/or features of the present invention will become more apparent from the following description taken in conjunction with the accompanying drawings, in which:
fig. 1 shows a flow chart of an encoding method of an exemplary embodiment of the mud pulse data encoding method based on multidimensional space coordinates of the present invention.
FIG. 2A illustrates a 3-dimensional spatial coordinate diagram of an exemplary embodiment of a multi-dimensional spatial coordinate based mud pulse data encoding method of the present invention; figure 2B shows an X =1 region coordinate diagram of an exemplary embodiment of a mud pulse data encoding method based on multidimensional spatial coordinates of the present invention;
fig. 2C shows a schematic diagram of X =1,y =8 region coordinates of an exemplary embodiment of the mud pulse data encoding method based on multidimensional space coordinates.
FIG. 3A is a schematic diagram of a corresponding code sequence when the similarity of an exemplary embodiment of the mud pulse data coding method based on multi-dimensional space coordinates is 0; FIG. 3B is a schematic diagram of a corresponding code sequence when the similarity of an exemplary embodiment of the method for encoding mud pulse data based on multidimensional space coordinates is 1; FIG. 3C is a schematic diagram of a corresponding code sequence with a similarity of 2 according to an exemplary embodiment of the method for encoding mud pulse data based on multidimensional space coordinates; fig. 3D shows a schematic diagram of a corresponding code sequence when the similarity of an exemplary embodiment of the mud pulse data coding method based on multidimensional space coordinates is 3.
Fig. 4 shows a schematic structural diagram of an encoding device of an exemplary embodiment of the mud pulse data encoding device based on multidimensional space coordinates.
Fig. 5 is a schematic structural diagram of a computer device of an exemplary embodiment of the multi-dimensional space coordinate-based mud pulse data encoding device of the present invention.
Description of the reference numerals:
100-mud pulse data coding device, 101-measured data analysis unit, 102-dimension range determination unit, 103-coordinate matching unit, 104-first calculation unit, 105-similarity determination unit, 106-second calculation unit, 107-third calculation unit, 108-coding unit, 200-computer device, 201-memory, 202-processor.
Detailed Description
Hereinafter, the mud pulse data encoding and transmitting method, apparatus and device of the present invention will be described in detail with reference to exemplary embodiments.
It should be noted that "first," "second," "third," and the like are merely for convenience of description and for ease of distinction, and are not to be construed as indicating or implying relative importance.
The invention provides a mud pulse data coding method based on multi-dimensional space coordinates.
In an exemplary embodiment of the mud pulse data encoding method based on multi-dimensional space coordinates of the present invention, the mud pulse data encoding method includes the steps of:
s1, determining a measurement numerical value and a total numerical value M of the downhole measurement data to be coded according to the measurement range S and the precision K of the downhole measurement data.
Specifically, the total number of numerical values may be calculated by: m = S/K, wherein M is the total number of numerical values of the downhole measurement data to be encoded, S is the measurement range of the downhole measurement data, and K is the precision of the downhole measurement data.
S2, determining a decimal number range G corresponding to each dimensionality according to an arithmetic n-th square root of the total number M of the numerical values.
Specifically, the decimal range G corresponding to each dimension may be determined as the value obtained by rounding the nth root of the total number M of the numerical values.
Where n is the total number of dimensions of the multi-dimensional space coordinate, 1 ≦ n ≦ 11, e.g., n may be 1, 2, 3, 4, 5, 6, etc.
And S3, sequentially arranging all measurement values of the underground measurement data on each dimension from small to large so that each multi-dimensional space coordinate corresponds to one measurement value.
S4, converting the decimal number range corresponding to each dimension into a binary system, and determining the number x of the first-dimension time slots theoretically required by each dimension 0 。
For example, the first dimension time slot number can be calculated as:
wherein x is 0 And G is the decimal number range corresponding to each dimension.
And S5, comparing the coordinate values of all dimensions of the two adjacent transmission data, determining the similarity H, and determining the number y of similarity time slots according to the similarity.
For example, the number of similarity time slots y may be determined as the total number of dimensions n of the multi-dimensional space coordinate.
S6, based on the similarity H, the number x of second dimension time slots needed by a certain dimension i with dimension coordinate value difference in two adjacent transmission data is determined i Determining the number x of time slots in the first dimension 0 The number x of the third dimension time slots needed by a certain dimension j without dimension coordinate value difference in the two adjacent transmission data j Is determined to be 0.
S7, determining the coding sequence of the underground measurement data as follows: number of similarity time slots y plus number of second dimension time slots x i And completing mud pulse data encoding.
The value of H is 1 to 2 n The value of i is 0-n, and the value of j is 0-n.
When H is 1, the coordinates of the underground measurement data transmitted twice in the adjacent directions on all dimensions are completely different, namely the value of i corresponding to a certain dimension i with dimension coordinate value difference in the data transmitted twice in the adjacent directions comprises 1, 2, 3, 4 \8230 \8230n, and the value of j corresponding to a certain dimension j without dimension coordinate value difference in the data transmitted twice in the adjacent directions comprises 0.
When H is 2, the 1 st dimensional coordinate of the underground measurement data which are transmitted in two adjacent times is the same, and the 2 nd to nth dimensional coordinates are different, namely the value of i corresponding to a certain dimension i in the data transmitted in two adjacent times, in which the dimensional coordinate value difference occurs, comprises 2, 3, 4, \ 8230 \ 8230n, and the value of j corresponding to a certain dimension j in the data transmitted in two adjacent times, in which the dimensional coordinate value difference does not occur, comprises 1.
When H is 3, the 2 nd dimensional coordinate of the underground measurement data which are transmitted in two adjacent times is the same, and the 1 st dimensional coordinate and the 3 rd to nth dimensional coordinates are different, namely the value of i corresponding to a certain dimension i in the data transmitted in two adjacent times, in which the dimensional coordinate value difference occurs, comprises 1, 3, 4 \8230 \82308230n, and the value of j corresponding to a certain dimension j in the data transmitted in two adjacent times, in which the dimensional coordinate value difference does not occur, comprises 2.
By analogy, when H is 1+ n, the nth dimensional coordinate of the underground measurement data transmitted for two adjacent times is the same, and the 1 st to the n-1 st dimensional coordinates are different, that is, the value of i corresponding to a certain dimension i with dimensional coordinate value difference in the data transmitted for two adjacent times includes 1, 2, 3, \ 8230, n-1, and the value of j corresponding to a certain dimension j without dimensional coordinate value difference in the data transmitted for two adjacent times includes n.
When H is 2+ n, the 1 st and 2 nd dimensional coordinates of the underground measurement data transmitted twice adjacently are the same, and the 3 rd to nth dimensional coordinates are different, that is, the value of i corresponding to a certain dimension i with dimension coordinate value difference in the data transmitted twice adjacently includes 3, 4, 5 \8230 \ 8230n, and the value of j corresponding to a certain dimension j without dimension coordinate value difference in the data transmitted twice adjacently includes 1 and 2.
When H is 3+ n, the 1 st and 3 rd dimensional coordinates of the underground measurement data transmitted twice adjacently are the same, and the 2 nd and 4 th to nth dimensional coordinates are different, that is, the value of i corresponding to a certain dimension i with dimension coordinate value difference in the data transmitted twice adjacently includes 2, 4, 5 \8230 \8230n, and the value of j corresponding to a certain dimension j without dimension coordinate value difference in the data transmitted twice adjacently includes 1 and 3.
By analogy, when H takes 2 n In the process, the coordinates of the underground measurement data transmitted twice in the two adjacent directions are completely the same in all dimensions, namely the value of i corresponding to a certain dimension i with dimension coordinate value difference in the data transmitted twice in the two adjacent directions comprises 0, and the value of j corresponding to a certain dimension j without dimension coordinate value difference in the data transmitted twice in the two adjacent directions comprises 1, 2, 3, 4 \8230, 8230n.
In the data coding process, y similarity time slots and the existence of pulses are combined to represent that j-dimension coordinate values are the same and corresponding measurement values of underground measurement data are also the same in adjacent two-time transmission data; using x i And combining the second dimension time slot with the existence of the pulse to represent the measurement value of the underground measurement data corresponding to each i-dimension coordinate value.
In this embodiment, the value of n is set to 3-5 for the underground engineering parameter measurement, which is most beneficial to ensuring the data measurement precision, and can also save the time slot number and ensure the transmission efficiency. Preferably, n may take 4.
For example, assuming that the measurement range and the accuracy requirement can be achieved only when the total number of values of the measurement data to be represented downhole is M = S/K =4096, the number of time slots required in each dimensional space when M =4096 can be analyzed to determine the most preferable dimensional number.
For convenience of calculation and analysis, it is assumed that when each dimension is subjected to mud pulse coding, a plurality of time slots are occupied by a plurality of numerical values, for example, an X coordinate has 64 numerical values, so that the number of the time slots required by the X dimension is 64, and the number of pulses at the position of the pulse represents that X is several; for another example, if the similarity H has 4 values, the number of time slots required for the similarity H is 4.
Then, when the dimension n =2, i.e. 4096 to the power of 2:
the formed coordinates are plane coordinates (X, Y), the values of the data ranges corresponding to two dimensions in the coordinates are all 1-64, and if the coordinate transmission is completed by using the assumed encoding method, 64+64=128 time slots are required. n =2 in this case, the 2 nd data (X) 2 ，Y 2 ) And data (X) at time 1 1 ，Y 1 ) The value of the similarity H between the two groups is 2 n In this case, 2 n =4 ((1)) 1 ＝X 2 And Y is 1 ＝Y 2 ；②X 1 ＝X 2 And Y is 1 ≠Y 2 ；③X 1 ≠X 2 And Y is 1 ＝Y 2 ；④X 1 ≠X 1 And Y is 1 ≠Y 2 ). Therefore, the similarity H needs 4 time slots, and therefore 128+4=132 time slots are needed in total after all data are transmitted.
When the dimension n =3, i.e. 4096 to the 3 th power:
the coordinates formed are three-dimensional (X, Y, Z) coordinates, in-coordinateThe values of the data ranges corresponding to the three dimensions are all 1-16, and 16+16=48 time slots are needed to complete the coordinate transmission by using the assumed coding mode. n =3 in this case, the value of the similarity H is 2 n In this case, 2 n =8 ((1)) 1 ＝X 2 ，Y 1 ＝Y 2 And Z is 1 ＝Z 2 ；②X 1 ＝X 2 ，Y 1 ＝Y 2 And Z is 1 ≠Z 2 ；③X 1 ＝X 2 ，Y 1 ≠Y 2 And Z is 1 ≠Z 2 ；④X 1 ＝X 2 ，Y 1 ≠Y 2 And Z is 1 ＝Z 2 ；⑤X 1 ≠X 2 ，Y 1 ＝Y 2 And Z is 1 ＝Z 2 ；⑥X 1 ≠X 1 ，Y 1 ＝Y 2 And Z is 1 ≠Z 2 ；⑦X 1 ≠X 1 ，Y 1 ≠Y 2 And Z is 1 ≠Z 2 ；⑧X 1 ≠X 1 ，Y 1 ≠Y 2 And Z is 1 ＝Z 2 ). Therefore, the similarity H needs 8 time slots, so that 48+8=56 time slots are needed in total after all data are transmitted.
When the dimension n =4, i.e. 4096 to the 4 th power:
the formed coordinates are four-dimensional coordinates (X, Y, Z, W), the values of the data ranges corresponding to the four dimensions in the coordinates are all 1-8, and if the coordinates are transmitted by using the assumed encoding method, 8+8=32 time slots are required. In the case where n =4, the similarity H takes a value of 2 n In this case, 2 n And (5) = 16. Therefore, the similarity H needs 16 time slots, so that a total of 32+16=48 time slots are needed after all data are transmitted.
When dimension n =5, i.e. 4096 to the power of 5:
the value is 6 after rounding up, and the formed coordinates are five-dimensional coordinates (X, Y, Z,w, U), the values of the data ranges corresponding to the five dimensions in the coordinate are all 1-6, and if the transmission of the coordinate is completed by using the assumed coding method, 6+6=30 time slots are needed. n =5 in this case, the similarity H takes a value of 2 n In this case, 2 n =32 kinds. Therefore, the similarity H needs 32 time slots, so that 30+32=62 time slots are needed in total after all data are transmitted.
When the dimension n =6, i.e. 4096 to the power of 6 is:
the formed coordinates are six-dimensional coordinates (X, Y, Z, W, U, V), the values of the data ranges corresponding to the six dimensions in the coordinates are all 1-4, and if the transmission of the coordinates is completed by using the assumed encoding method, 4+4 =24time slots are needed. n =6 in this case, the value of the similarity H is 2 n In this case, 2 n =64 kinds. Therefore, the similarity H needs 32 time slots, so that a total of 24+64=88 time slots is needed after all data is transmitted.
Thus, several conclusions can be drawn from the comparative analysis above:
(1) the larger the dimension n is, the less time slots occupied by coordinate coding will be, but the larger the dimension n is, the exponential increase of the number of permutation and combination of the similarity H will be caused, and the similarity class 2 n The occupied time slots are obviously increased, so that the value of n is not suitable to be overlarge, and 2 is the best mode n The value is as much as the value of each dimension (for example, when n = 4).
(2) Aiming at underground engineering parameter measurement, the occupied time slot is the minimum when n is taken as 4 dimensions, and the transmission efficiency is the highest. Therefore, n is preferably 3 to 5.
In this embodiment, the downhole measurement data may include directional parameters including well deviation, azimuth, and toolface, and drilling engineering parameters, which may include weight-on-bit, torque, temperature, pressure, vibration, and torque.
The invention further provides a mud pulse data transmission method based on the multi-dimensional space coordinates.
In an exemplary embodiment of the mud pulse data transmission method based on multi-dimensional space coordinates of the present invention, the mud pulse data transmission method includes: and coding the downhole measurement data according to the coding method and then transmitting mud pulse data. In the data transmission process, y similarity time slots and the existence of pulses are combined to represent that j-dimension coordinate values are the same and corresponding measurement values of underground measurement data are the same in the two adjacent times of data transmission; by x i And combining the second dimension time slot with the existence of the pulse to represent the measurement value of the underground measurement data corresponding to each i-dimension coordinate value.
In another aspect, the invention provides a mud pulse data coding device based on multi-dimensional space coordinates.
In an exemplary embodiment of the invention, the mud pulse data coding device based on multidimensional space coordinates comprises a measurement data analysis unit, a dimension range determination unit, a coordinate matching unit, a first calculation unit, a similarity determination unit, a second calculation unit, a third calculation unit and a coding unit.
The measuring data analysis unit is configured to determine a measuring numerical value and a total numerical value M of the downhole measuring data needing to be coded according to the measuring range S and the precision K of the downhole measuring data.
And the dimension range determining unit is configured to determine a decimal number range corresponding to each dimension according to an arithmetic n-th-order root of the total number M of the numerical values, wherein n is the total number of the dimensions of the multidimensional space coordinate.
The coordinate matching unit is respectively connected with the measurement data analysis unit and the dimension range determination unit and is configured to establish multi-dimensional space coordinates, and the measurement values of the underground measurement data are sequentially arranged on each dimension from small to large, so that each multi-dimensional space coordinate corresponds to one measurement value.
The first calculating unit is connected with the dimension range determining unit and is configured to convert the decimal number range corresponding to each dimension into binary number and determine each dimensionNumber x of first dimension time slots theoretically required by dimension 0 。
And the similarity determining unit is configured to compare the coordinate values of the dimensions of the two adjacent transmission data to determine the similarity H.
The second calculating unit is connected with the similarity determining unit and is configured to determine the number y of the similarity time slots according to the similarity.
The third calculating unit is respectively connected with the similarity determining unit and the first calculating unit and is configured to calculate the number x of second dimension time slots required by a certain dimension i with dimension coordinate value difference in two adjacent transmission data based on the similarity H i Determining as a first dimension time slot number x 0 。
The coding unit is respectively connected with the second calculation unit and the third calculation unit and is configured to add the second dimension time slot number x to the similarity time slot number y i The combined form of (a) forms a coded sequence of downhole measurement data.
In another aspect, the invention provides a mud pulse data transmission device based on multi-dimensional space coordinates.
In an exemplary embodiment of the mud pulse data transmission device based on multidimensional space coordinates, the mud pulse data transmission device comprises the encoding device and a transmission unit, wherein the transmission unit is connected with the encoding unit and is configured to encode the downhole measurement data according to the encoding sequence and then perform mud pulse data transmission.
The mud pulse data encoding method and the mud pulse data transmission method according to the present invention may be programmed as a computer program and corresponding program code or instructions may be stored in a computer readable storage medium, which when executed by a processor causes the processor to perform at least one of the above-described encoding method and transmission method, the processor and memory may be included in a computer device.
Exemplary embodiments according to still another aspect of the present invention also provide a computer-readable storage medium storing a computer program. The computer-readable storage medium stores a computer program that, when executed by a processor, causes the processor to perform at least one of an encoding method and a transmission method according to the present invention. The computer readable recording medium is any data storage device that can store data read by a computer system. Examples of the computer-readable recording medium include: read-only memory, random access memory, read-only optical disks, magnetic tapes, floppy disks, optical data storage devices, and carrier waves (such as data transmission through the internet via wired or wireless transmission paths).
Exemplary embodiments according to still another aspect of the present invention also provide a computer apparatus. The computer device includes a processor and a memory. The memory is for storing a computer program. The computer program is executed by a processor, which causes the processor to execute at least one of the encoding method and the transmission method according to the present invention.
For a better understanding of the above-described exemplary embodiments of the present invention, reference is made to the following description taken in conjunction with the accompanying drawings and specific examples.
Example 1
As shown in fig. 1, the mud pulse data encoding method based on multidimensional space coordinates is implemented by adopting the following technical scheme:
Taking the measurement of the pressure in the well as an example, the pressure measurement range is 0-100 MPa, the measurement precision K is 0.2%, the numerical values required to be coded are 0.2MPa, 0.4MPa, 0.6MPa \8230 \ 8230, 99.6MPa, 99.8MPa, and 100MPa, respectively, and the number of the numerical values required to be coded is M = S/K =100/0.2= 500.
And 2, opening the total number of each kind of measurement data by the power of n to determine a decimal number range corresponding to each dimension.
Taking three-dimensional space coordinates of three dimensions (i.e., n = 3) as an example, M values to be encoded are decomposed into 3-dimensional space coordinates, and the range of each dimension (i.e., X-axis, Y-axis, and Z-axis) is 500 to the power of 3 and is equal to 7.94 (rounded to 8), that is, the three coordinate ranges of X, Y, and Z are all from 1 to 8.
And 3, arranging the specific numerical values corresponding to each kind of measurement data from small to large in sequence according to n-dimensional coordinates, wherein each coordinate corresponds to 1 specific numerical value.
Namely, the coordinate ranges of the three axes are 1 to 8, 8 × 8=512 points can be formed, the number of the coordinate points is greater than the number of the pressure values, and each pressure value can correspond to a unique three-dimensional coordinate.
When coding is carried out, the numerical value of the X axis is fixed to be 1, the coding of the pressure values of the Y and the Z is finished sequentially from small to large, and the pressure value corresponding to the small coordinate is small. When the coordinate of the X axis is 1, the expressed pressure value is between 0.2MPa and 12.6 MPa; when the X-axis coordinate is 1 and the Y-axis coordinate is 1, the expressed pressure value is between 0.2MPa and 1.4 MPa; when all three axis coordinate positions are determined, the pressure value is also determined. When the X coordinate is 1, the coordinate table composed of the Y and Z coordinates is shown in table 1.
TABLE 1 coordinate Table of Y and Z coordinates when the X coordinate is 1
| Coordinate values | (1,1) | (1,2) | (1,3) | (1,4) | (1,5) | (1,6) | (1,7) | (1,8) |
| Pressure value | 0.2 | 0.4 | 0.4 | 0.6 | 0.8 | 1 | 1.2 | 1.4 |
| Coordinate values | (2,1) | (2,2) | (2,3) | (2,4) | (2,5) | (2,6) | (2,7) | (2,8) |
| Pressure value | 1.6 | 1.8 | 2 | 2.2 | 2.4 | 2.6 | 2.8 | 3 |
| Coordinate values | (3,1) | (3,2) | (3,3) | (3,4) | (3,5) | (3,6) | (3,7) | (3,8) |
| Pressure value | 3.2 | 3.4 | 3.6 | 3.8 | 4 | 4.2 | 4.4 | 4.6 |
| Coordinate values | (4,1) | (4,2) | (4,3) | (4,4) | (4,5) | (4,6) | (4,7) | (4,8) |
| Pressure value | 4.8 | 5 | 5.2 | 5.4 | 5.6 | 5.8 | 6 | 6.2 |
| Coordinate values | (5,1) | (5,2) | (5,3) | (5,4) | (5,5) | (5,6) | (5,7) | (5,8) |
| Pressure value | 6.4 | 6.6 | 6.8 | 7 | 7.2 | 7.4 | 7.6 | 7.8 |
| Coordinate values | (6,1) | (6,2) | (6,3) | (6,4) | (6,5) | (6,6) | (6,7) | (6,8) |
| Pressure value | 8 | 8.2 | 8.4 | 8.6 | 8.8 | 9 | 9.2 | 9.4 |
| Coordinate values | (7,1) | (7,2) | (7,3) | (7,4) | (7,5) | (7,6) | (7,7) | (7,8) |
| Pressure value | 9.6 | 9.8 | 10 | 10.2 | 10.4 | 10.6 | 10.8 | 11 |
| Coordinate values | (8,1) | (8,2) | (8,3) | (8,4) | (8,5) | (8,6) | (8,7) | (8,8) |
| Pressure value | 11.2 | 11.4 | 11.6 | 11.8 | 12 | 12.2 | 12.4 | 12.6 |
By utilizing the dimension-based coding method, the pressure value can be divided into a plurality of regions, each dimension can determine a pressure value range, and the pressure range gradually approaches to the real pressure value along with the increase of the dimension.
And 4, converting the decimal number range of each dimension into a binary system, and representing each dimension by using a certain number of time slots.
In the three-dimensional space coordinate described above, the maximum number of possible occurrences for each axis is 8 (i.e., 2) 3 ) Only binary codes of 0 and 1 can be identified in mud pulse transmission, and 8 combinations can be obtained by only using 3 time slots and combining the existence of pulses through PPM pulse position modulation, wherein 1-8 are shown. The maximum possible three-dimensional coordinates is 8, that is, each coordinate occupies 3 time slots, and any one coordinate in the space needs at least 9 time slots, that is, all 512 positions can be represented. The pulse sequence corresponding to coordinates (3, 8, 5) is shown in FIG. 2A; x =1 regional coordinates are shown in fig. 2B; x =1,y =8 region coordinates are shown in fig. 2C.
And 5, determining the similarity of the two adjacent transmission data, and representing the similarity by using a certain number of time slots.
Here, the data similarity H means: the similarity degree of the pressure data transmitted for the first time and the pressure data transmitted for the second time is larger, and the two pressure data are closer.
There are 8 cases where the value of the data similarity H is 1, 2, 3, 4, 5, 6, 7, and 8, respectively. The correspondence between different values of the similarity H and the degree of similarity of the dimensional pressure data is shown in table 2.
TABLE 2 correspondence of different values of similarity H to the degree of similarity of each dimension of pressure data
It should be noted that, when the dimension n =1, there are only 2 similarity cases (i.e., the data is the same and the data is different) between the pressure data transmitted for the first time and the pressure data transmitted for the second time, and the similarity H has 2 values. If the mud pulse is adopted to represent the two conditions, only 1 time slot needs to be occupied, namely the pressure data which represents two times of transmission are the same when the time slot has the pulse, and the pressure data which represents two times of transmission are different when the time slot does not have the pulse.
When the dimension n =2, there are 4 similarity conditions (4 types can be combined according to the same or different X and Y) between the pressure data transmitted for the first time and the pressure data transmitted for the second time, and the similarity H has 4 numerical values. If the mud pulse is used to represent the 4 cases, only 2 time slots need to be occupied, specifically representing the case: (1) the X-axis coordinate and the Y-axis coordinate of the pressure data transmitted for the second time and the pressure data transmitted for the first time are the same, and the pressure data are empty in 2 time slots; (2) when the X-axis coordinates of the pressure data transmitted for the second time are the same as those of the pressure data transmitted for the first time and the Y-axis coordinates are different, the 1 st time slot is empty, and the 2 nd time slot is provided with a pulse; (3) the X-axis coordinate of the pressure data transmitted for the second time is different from that of the pressure data transmitted for the first time, the Y-axis coordinate is the same, the 1 st time slot has pulses, and the 2 nd time slot is empty; (4) when the X-axis coordinate and the Y-axis coordinate of the pressure data transmitted for the second time are different from those of the pressure data transmitted for the first time, pulses exist on the 1 st time slot and the 2 nd time slot.
When the dimension n =3, there are 8 similarity conditions (8 types can be combined according to the identity and the non-identity of X, Y and Z) between the pressure data transmitted for the first time and the pressure data transmitted for the second time, and then the similarity H has 8 numerical values. If the mud pulse is used to represent these 8 cases, only 3 time slots need to be occupied, specifically represented as follows: (1) the X-axis coordinates, the Y-axis coordinates and the Z-axis coordinates of the pressure data transmitted for the second time and the pressure data transmitted for the first time are the same, and all the time slots are empty in 3 time slots; (2) when the X-axis coordinates of the pressure data transmitted for the second time are the same as those of the pressure data transmitted for the first time, and the Y-axis coordinates are different from those of the pressure data transmitted for the first time, the 1 st time slot is empty, and pulses are arranged on the 2 nd time slot and the 3 rd time slot; (3) when the X-axis coordinate of the pressure data transmitted for the second time is different from that of the pressure data transmitted for the first time, the Y-axis coordinate is the same, and the Z-axis coordinate is different, a pulse exists in the 1 st time slot, the 2 nd time slot is empty, and a pulse exists in the 3 rd time slot; (4) the pressure data transmitted for the second time is different from the pressure data transmitted for the first time in X-axis coordinates, the Y-axis coordinates are different, the Z-axis coordinates are the same, the 1 st time slot and the 2 nd time slot have pulses, and the 3 rd time slot is empty; (5) the X-axis coordinate of the pressure data transmitted for the second time is the same as that of the pressure data transmitted for the first time, the Y-axis coordinate is the same, the Z-axis coordinate is different, the 1 st time slot and the 2 nd time slot are empty, and the 3 rd time slot is provided with a pulse; (6) the pressure data transmitted for the second time and the pressure data transmitted for the first time have the same X-axis coordinate, different Y-axis coordinates and the same Z-axis coordinate, the 1 st time slot and the 3 rd time slot are empty, and the 2 nd time slot has pulses; (7) the pressure data transmitted for the second time is different from the pressure data transmitted for the first time in X-axis coordinates, the Y-axis coordinates are the same, the Z-axis coordinates are the same, a pulse exists on the 1 st time slot, and the 2 nd time slot and the 3 rd time slot are empty; (8) when the X-axis, Y-axis and Z-axis coordinates of the second transmitted pressure data and the first transmitted pressure data are different, there are pulses in all 3 time slots.
Therefore, when the data similarity H has 8 values, it can be represented by a combination of 3 time slots and the presence or absence of a pulse.
And 6, each measured data coding sequence is a data similarity H + each dimension coordinate mark, and when the previous dimension values are the same, only the dimension values with different values are transmitted.
(1) When the data similarity H =1, that is, the coordinates of the three dimensions X (i.e., the X axis, the Y axis, and the Z axis) are all different, the pressure value sequence transmitted for the second time is: similarity 3 time slots + X-axis coordinate 3 time slots + Y-axis coordinate 3 time slots + Z-axis coordinate 3 time slots, for a total of 12 time slots. The corresponding coding sequence with similarity of 1 is schematically shown in FIG. 3A.
(2) When the data similarity H =2, that is, the X-axis coordinate is the same, and the Y-axis coordinate and the Z-axis coordinate are different, the pressure value sequence transmitted for the second time is: similarity 3 time slots + Y-axis coordinate 3 time slots + Z-axis coordinate 3 time slots, for a total of 9 time slots. That is, the repeated X-axis coordinates are not transmitted, and only the Y-axis and Z-axis coordinates are transmitted. The corresponding coding sequence with similarity of 2 is schematically shown in FIG. 3B.
(3) When the data similarity H =3, that is, the Y-axis coordinate is the same, and the X-axis coordinate and the Z-axis coordinate are different, the pressure value sequence transmitted for the second time is: similarity 3 time slots + X-axis coordinate 3 time slots + Z-axis coordinate 3 time slots, for a total of 9 time slots. That is, the repeated Y-axis coordinates are not transmitted, only the X-axis and Z-axis coordinates are transmitted.
(4) When the data similarity H =4, namely the Z-axis coordinate is the same, and the X-axis coordinate and the Y-axis coordinate are different, the pressure value sequence transmitted for the second time is as follows: similarity 3 time slots + X-axis coordinate 3 time slots + Y-axis coordinate 3 time slots, for a total of 9 time slots. That is, the repeated Z-axis coordinates are not transmitted, only the X-axis and Y-axis coordinates are transmitted.
(5) When the data similarity H =5, that is, the X-axis and Y-axis coordinates are the same, and the Z-axis coordinate is different, the pressure value sequence transmitted for the second time is: similarity 3 time slots + Z-axis coordinates 3 time slots, for a total of 6 time slots. That is, the repeated X-axis and Y-axis coordinates are not transmitted, only the Z-axis coordinate is transmitted. A schematic diagram of the corresponding coding sequence with a similarity of 5 is shown in fig. 3C.
(6) When the data similarity H =6, that is, the X-axis and Z-axis coordinates are the same, and the Y-axis coordinate is different, the pressure value sequence transmitted for the second time is: similarity 3 time slots + Y-axis coordinates 3 time slots, for a total of 6 time slots. That is, the repeated X-axis and Z-axis coordinates are not transmitted, and only the Y-axis coordinate is transmitted.
(7) When the data similarity H =7, that is, the Y-axis and Z-axis coordinates are the same, and the X-axis coordinate is different, the pressure value sequence transmitted for the second time is: similarity 3 time slots + X axis coordinate 3 time slots for a total of 6 time slots. That is, the repeated Y-axis and Z-axis coordinates are not transmitted, only the X-axis coordinate is transmitted.
(8) When the data similarity H =8, the two transmission values are completely the same, namely the coordinates of the X axis, the Y axis and the Z axis are completely the same, and the pressure value sequence transmitted for the second time is as follows: similarity is 3 time slots, and 3 time slots are calculated. That is, the repeated X-axis, Y-axis and Z-axis coordinates are not transmitted, and only expressed by 3 time slots of the data similarity H. The corresponding coding sequence with similarity of 8 is schematically shown in FIG. 3D. It can be seen that, the conventional mud pulse data coding algorithm does not adopt a similarity compression time slot, 8 data are transmitted, the data of the X axis, the Y axis and the Z axis need to be transmitted once every time of transmission, 24 data need to be transmitted, each data occupies 3 time slots, and 72 time slots are occupied. And the mud pulse data coding algorithm of the above example only occupies 9 time slots at most when transmitting data, and the transmission rate can be improved by 12%.
Example 2
As shown in fig. 4, the mud pulse data encoding device 100 includes a measurement data analysis unit 101, a dimensional range determination unit 102, a coordinate matching unit 103, a first calculation unit 104, a similarity determination unit 105, a second calculation unit 106, a third calculation unit 107, and an encoding unit 108.
The measurement data analysis unit 101 is configured to determine a measurement numerical value and a total number M of numerical values of the downhole measurement data to be encoded according to a measurement range S and precision K of the downhole measurement data.
The dimension range determining unit 102 is configured to determine a decimal number range corresponding to each dimension according to an arithmetic n-th-order root of the total number M of the numerical values, where n is the total number of dimensions of the multidimensional space coordinate.
The coordinate matching unit 103 is connected to the measurement data analysis unit 101 and the dimension range determination unit 102, and is configured to establish multidimensional space coordinates, and sequentially arrange the measurement values of the downhole measurement data on each dimension in the order from small to large, so that each multidimensional space coordinate corresponds to one measurement value.
The first computing unit 104 is connected to the dimension range determining unit 102, and is configured to convert the decimal number range corresponding to each dimension into binary, and determine what each dimension theoretically needsNumber of first dimension time slots x 0 。
The similarity determining unit 105 is configured to compare coordinate values of each dimension of two adjacent transmission data to determine a similarity H.
The second calculating unit 106 is connected to the similarity determining unit 105 and configured to determine the number y of similarity time slots according to the similarity.
The third calculating unit 107 is respectively connected to the similarity determining unit 105 and the first calculating unit 104, and configured to, based on the similarity H, count the number x of second dimension time slots required for a certain dimension i with a difference in dimension coordinate values occurring in two adjacent transmission data i Determining as a first dimension time slot number x 0 。
The encoding unit 108 is connected to the second and third calculation units 106 and 107, respectively, and is configured to add the second dimension time slot number x to the similarity time slot number y i Forms a coded sequence of downhole measurement data.
Example 3
As shown in fig. 5, a computer device 200 includes a memory 201 and a processor 202. The memory is for storing a computer program.
The computer program is executed by a processor, causing the processor to execute at least one of the encoding method and the transmission method according to the present invention.
In summary, the beneficial effects of the invention include:
(1) The coding mode is simple, all measurement parameter values only need to be numbered from small to large in sequence according to the coordinates, and the precision of different measurement data can be displayed by controlling the number of dimensions and the range of each dimension.
(2) The encoding of the measurement data divides the pressure value into a plurality of regions by utilizing dimension division, each dimension can determine a pressure value range, and the pressure range gradually approaches to the real pressure value along with the increase of the dimension.
(3) The data similarity index is introduced, and in the measuring process, because data generally cannot be suddenly changed, the similarity of the data transmitted for the second time and the data transmitted for the first time is generally high, and the transmitted data is compressed by utilizing the characteristic, so that the transmission rate is obviously improved compared with the transmission rate without data compression.
Although the present invention has been described above in connection with the exemplary embodiments and the accompanying drawings, it will be apparent to those of ordinary skill in the art that various modifications may be made to the above-described embodiments without departing from the spirit and scope of the claims.
Claims (12)
1. A mud pulse data coding method based on multi-dimensional space coordinates is characterized by comprising the following steps:
determining the measurement value and the total number M of the value of the underground measurement data to be coded according to the measurement range S and the precision K of the underground measurement data;
determining a decimal number range G corresponding to each dimensionality according to an arithmetic n-th square root of the total number M of the numerical values, wherein n is the total number of the dimensionalities of the multidimensional space coordinate;
sequentially arranging all measurement values of the underground measurement data on each dimension from small to large so that each multi-dimensional space coordinate corresponds to one measurement value;
converting the decimal number range corresponding to each dimension into binary, and determining the number x of the first dimension time slots theoretically required by each dimension 0 ；
Comparing coordinate values of each dimension of the two adjacent times of transmission data, determining similarity H, and determining the number y of similarity time slots according to the similarity;
based on the similarity H, the number x of second dimension time slots required by a certain dimension i with dimension coordinate value difference in two adjacent times of transmission data is determined i Determining the number x of time slots in the first dimension 0 The number x of the third dimension time slots needed by a certain dimension j without dimension coordinate value difference in the two adjacent times of transmission data j Is determined to be 0; and
determining an encoding sequence of downhole measurement data as: number of similarity time slots y plus number of second dimension time slots x i In combination with each other, andcompleting mud pulse data coding;
wherein the value of H is 1-2 n I ranges from 0 to n, and j ranges from 0 to n.
2. The method of claim 1, wherein the downhole measurement data comprises directional parameters and drilling engineering parameters, the directional parameters comprise well deviation, azimuth and toolface, and the drilling engineering parameters comprise weight on bit, torque, temperature, pressure, vibration and torque.
3. The method for encoding mud pulse data according to claim 1, wherein the total number of values is calculated by: m = S/K, wherein M is the total number of numerical values of the downhole measurement data to be encoded, S is the measurement range of the downhole measurement data, and K is the precision of the downhole measurement data.
4. The method for encoding mud pulse data based on multidimensional space coordinates of claim 1, wherein the integer n of the total number M of numerical values is determined as the decimal range G corresponding to each dimension.
5. The method for encoding mud pulse data based on multidimensional space coordinates of claim 1, wherein n is greater than or equal to 3 and less than or equal to 5.
6. The method of claim 1, wherein the number of first dimension time slots is calculated as:
wherein x is 0 And G is the decimal number range corresponding to each dimension.
7. The method of claim 1, wherein the number y of similarity time slots is determined as a total number n of dimensions of the multidimensional space coordinate.
8. A mud pulse data transmission method based on multi-dimensional space coordinates is characterized in that the transmission method encodes downhole measurement data according to the encoding method of any one of claims 1 to 7 and then carries out mud pulse data transmission.
9. The mud pulse data coding device based on the multidimensional space coordinates is characterized by comprising a measured data analysis unit, a dimension range determination unit, a coordinate matching unit, a first calculation unit, a similarity determination unit, a second calculation unit, a third calculation unit and a coding unit,
the measurement data analysis unit is configured to determine a measurement numerical value and a total numerical value M of the downhole measurement data to be coded according to the measurement range S and the precision K of the downhole measurement data;
the dimension range determining unit is configured to determine a decimal number range corresponding to each dimension according to an arithmetic n-th-order root of the total number M of the numerical values, wherein n is the total number of the dimensions of the multidimensional space coordinate;
the coordinate matching unit is respectively connected with the measured data analysis unit and the dimension range determination unit, and is configured to establish multi-dimensional space coordinates, and arrange the measured values of the underground measured data on each dimension in sequence from small to large so that each multi-dimensional space coordinate corresponds to one measured value;
the first calculation unit is connected with the dimension range determination unit and is configured to convert the decimal number range corresponding to each dimension into binary, and determine the number x of the first dimension time slots theoretically required by each dimension 0 ；
The similarity determining unit is configured to compare coordinate values of each dimension of the two adjacent transmission data to determine similarity H;
the second calculation unit is connected with the similarity determination unit and is configured to determine the number y of similarity time slots according to the similarity;
the third calculation unit is respectively connected with the similarity determination unit and the first calculation unit and is configured to calculate the number x of second dimension time slots required by a certain dimension i with dimension coordinate value difference in two adjacent transmission data based on the similarity H i Determining the number x of time slots in the first dimension 0 (ii) a And
the coding unit is respectively connected with the second calculation unit and the third calculation unit and is configured to add the second dimension time slot number x to the similarity time slot number y i The combined form of (a) forms a coded sequence of downhole measurement data.
10. A mud pulse data transmission device based on multidimensional space coordinates, wherein the mud pulse data transmission device comprises the encoding device and a transmission unit according to claim 9, and the transmission unit is connected with the encoding unit and is configured to encode the downhole measurement data according to the encoding sequence and then perform mud pulse data transmission.
11. A computer device, the device comprising:
a processor; and
memory storing a computer program which, when executed by a processor, implements at least one of the encoding method of any one of claims 1 to 7, the transmission method of claim 8.
12. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out at least one of the encoding method of any one of claims 1 to 7, the transmission method of claim 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210922404.0A CN115296752B (en) | 2022-08-02 | 2022-08-02 | Mud pulse data coding and transmitting method, device and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210922404.0A CN115296752B (en) | 2022-08-02 | 2022-08-02 | Mud pulse data coding and transmitting method, device and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115296752A true CN115296752A (en) | 2022-11-04 |
| CN115296752B CN115296752B (en) | 2024-04-30 |
Family
ID=83826083
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210922404.0A Active CN115296752B (en) | 2022-08-02 | 2022-08-02 | Mud pulse data coding and transmitting method, device and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115296752B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115664601A (en) * | 2022-12-28 | 2023-01-31 | 合力（天津）能源科技股份有限公司 | Wireless measurement while drilling information transmission method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120039151A1 (en) * | 2010-08-12 | 2012-02-16 | Precision Energy Services, Inc. | Mud pulse telemetry synchronous time averaging system |
| CN108343429A (en) * | 2018-02-01 | 2018-07-31 | 中国石油大学(华东) | A kind of mud signal recognition methods based on Analysis on confidence |
| CN110661580A (en) * | 2019-11-04 | 2020-01-07 | 中国石油集团川庆钻探工程有限公司 | Slurry pulse data coding method and transmission method |
| CN110753224A (en) * | 2018-07-23 | 2020-02-04 | 中国石油化工股份有限公司 | Data rearrangement method and system for measurement while drilling data compression |
| CN110792429A (en) * | 2019-11-04 | 2020-02-14 | 中国石油集团川庆钻探工程有限公司 | Method for encoding downhole data by simultaneously utilizing positive and negative pressure pulses and transmission method |
| CN114650103A (en) * | 2020-12-21 | 2022-06-21 | 航天科工惯性技术有限公司 | Slurry pulse data transmission method, device, equipment and storage medium |
-
2022
- 2022-08-02 CN CN202210922404.0A patent/CN115296752B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120039151A1 (en) * | 2010-08-12 | 2012-02-16 | Precision Energy Services, Inc. | Mud pulse telemetry synchronous time averaging system |
| CN108343429A (en) * | 2018-02-01 | 2018-07-31 | 中国石油大学(华东) | A kind of mud signal recognition methods based on Analysis on confidence |
| CN110753224A (en) * | 2018-07-23 | 2020-02-04 | 中国石油化工股份有限公司 | Data rearrangement method and system for measurement while drilling data compression |
| CN110661580A (en) * | 2019-11-04 | 2020-01-07 | 中国石油集团川庆钻探工程有限公司 | Slurry pulse data coding method and transmission method |
| CN110792429A (en) * | 2019-11-04 | 2020-02-14 | 中国石油集团川庆钻探工程有限公司 | Method for encoding downhole data by simultaneously utilizing positive and negative pressure pulses and transmission method |
| CN114650103A (en) * | 2020-12-21 | 2022-06-21 | 航天科工惯性技术有限公司 | Slurry pulse data transmission method, device, equipment and storage medium |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115664601A (en) * | 2022-12-28 | 2023-01-31 | 合力（天津）能源科技股份有限公司 | Wireless measurement while drilling information transmission method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115296752B (en) | 2024-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112218079B (en) | Point cloud layering method based on spatial sequence, point cloud prediction method and point cloud prediction equipment | |
| CN102804230B (en) | Use repeat patterns compression 3D grid | |
| CN115296752A (en) | Mud pulse data coding and transmitting method, device and equipment | |
| JP2004532458A (en) | Image retrieval using distance measurement | |
| Sleator et al. | Short encodings of evolving structures | |
| CN104881449A (en) | Image retrieval method based on manifold learning data compression hash | |
| CN115913247A (en) | Deep lossless compression method and system for high-frequency power data | |
| CN113806601B (en) | Peripheral interest point retrieval method and storage medium | |
| CN112100308B (en) | Similarity query method for Beidou space-time trajectory | |
| CN113312809B (en) | Spacecraft telemetry data multi-parameter anomaly detection method based on related cluster division | |
| CN107327292B (en) | Coding method of logging-while-drilling signal | |
| CN110489606B (en) | Packet Hilbert coding and decoding method | |
| CN115882867B (en) | Data compression storage method based on big data | |
| Chung et al. | Space-filling approach for fast window query on compressed images | |
| CN102084595B (en) | Method for counting vectors in regular point networks | |
| DeVore et al. | Processing terrain point cloud data | |
| CN106373189A (en) | Space-time body data coding method based on linear hexadecimal tree | |
| CN106572093B (en) | A kind of wireless sensor array data compression method and system | |
| Hasan et al. | Piecemeal Clustering: a Self-Driven Data Clustering Algorithm | |
| CN115879020A (en) | Reservoir classification depicting method and device, electronic equipment and storage medium | |
| CN113961752A (en) | Method and related device for analyzing basic reachability graph of information physical system | |
| CN110792429A (en) | Method for encoding downhole data by simultaneously utilizing positive and negative pressure pulses and transmission method | |
| Kulyk et al. | Development of binary information compression methods based on the binomial numerical function | |
| Gamble et al. | Adaptive tracking of representative cycles in regular and zigzag persistent homology | |
| Hsu et al. | Bandwidth-efficient distributed k-nearest-neighbor search with dynamic time warping |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |