A method for detecting abnormality of suction sample

Abstract

The present invention is applicable to the field of in vitro detection technology, and provides a method for detecting abnormal sampling, including obtaining a hydraulic pressure curve to be tested when a sampling needle is sampling; obtaining a first sampling state through the negative pulse area of ​​the hydraulic pressure curve to be tested; obtaining a second sampling state through the jitter wave area and the number of jitter waves based on the hydraulic pressure curve to be tested; when the first sampling state or the second sampling state is abnormal, abnormal processing is performed and an alarm prompt is given; wherein, if the first sampling state is abnormal, it indicates that the sampling needle is in a blocked needle state; if the second sampling state is abnormal, it indicates that the sampling needle is in an empty suction state. The present invention can improve the accuracy of the test results of the analyzer.

Description

Technical Field

[0001] The present invention relates to the field of in vitro detection technology, and in particular to a method for detecting abnormality in sample aspiration. Background Art

[0002] The sampling needle of the in vitro detection analyzer is often used to extract samples and reaction reagents, but the inner diameter of the sampling needle is small, and improper operation will affect its sampling effect and thus affect the final test results.

[0003] For example, if the sample sucked by the sampling needle contains impurities or the sample forms a gel, it is easy to cause the needle to be blocked, resulting in inaccurate sample aspiration and inaccurate test results. For another example, generally speaking, in order to increase the stability of the reagent, the reagent bottle mouth will retain the sealing film, and only open a small hole for the sampling needle to enter. If condensed water gathers on the sealing film, the needle wall of the sampling needle is likely to touch the condensed water during the descent process, causing false detection of the liquid level detection module, which is intuitively manifested as the sampling needle starting to aspirate before it descends to the actual liquid level position, resulting in empty aspiration, which ultimately leads to inaccurate test results.

[0004] Therefore, there is an urgent need for an abnormality detection method based on a sampling needle to improve the accuracy of the test results of an in vitro detection analyzer. Summary of the invention

[0005] The main purpose of the present invention is to propose a method for detecting abnormal sample aspiration, so as to solve the problem in the prior art that when using a sampling needle to aspirate samples, needle blockage or empty aspiration is prone to occur, thereby affecting the accuracy of the test results.

[0006] To achieve the above-mentioned purpose, a first aspect of an embodiment of the present invention provides a method for detecting abnormality in aspiration of a sample, comprising:

[0007] Obtain the hydraulic pressure curve to be measured when the sampling needle aspirates the sample;

[0008] Calculating the negative pulse area of ​​the hydraulic pressure curve to be measured according to the hydraulic pressure curve to be measured includes:

[0009] Calculating a first mean value of the hydraulic pressure curve to be measured;

[0010] removing a curve portion whose hydraulic value is less than the first mean value in the hydraulic pressure curve to be measured, calculating the mean value of the remaining curve portion, and obtaining a second mean value;

[0011] Calculating and obtaining the negative pulse area of ​​the hydraulic pressure curve to be measured according to the pattern area formed by the hydraulic pressure curve to be measured and the second mean value;

[0012] Acquiring a first sampling state through the negative pulse area of ​​the hydraulic pressure curve to be measured;

[0013] The second sampling state is obtained based on the shaking wave area and the number of shaking waves of the hydraulic pressure curve to be measured, including:

[0014] Finding the minimum position of the hydraulic pressure value in the hydraulic pressure curve to be tested;

[0015] Starting from the minimum position of the hydraulic pressure value and heading leftward, 200 points in the hydraulic pressure curve to be measured are selected as target data, and curve smoothing is performed on the target data to obtain a smooth curve;

[0016] In the range of the minimum position of the hydraulic value and the displacement of 200 points of the minimum position of the hydraulic value, the pattern area formed by the smooth curve and the hydraulic pressure curve to be measured is calculated to obtain the jitter wave area of ​​the hydraulic pressure curve to be measured;

[0017] Calculating the number of jitter waves whose peak values ​​of the pattern formed by the smooth curve and the hydraulic pressure curve to be measured exceed a set threshold, and obtaining the number of jitter waves of the hydraulic pressure curve to be measured;

[0018] When the first sample suction state or the second sample suction state is abnormal, abnormal processing is performed and an alarm prompt is given;

[0019] If the first sample suction state is abnormal, it means that the sampling needle is in a blocked state; if the second sample suction state is abnormal, it means that the sampling needle is in an empty suction state.

[0020] Optionally, the first mean value of the hydraulic pressure curve to be measured is calculated using the following formula:

[0021] ;

[0022] Wherein, Pressure[i] is the hydraulic pressure curve to be tested, n is the data length of the hydraulic pressure curve to be tested, and average1 is the first average value of the hydraulic pressure curve to be tested;

[0023] The second mean is calculated as follows:

[0024] ;

[0025] in, is the number of data with hydraulic pressure values ​​greater than average1, and average2 is the second average value of the hydraulic pressure curve to be tested;

[0026] According to the pattern area formed by the hydraulic pressure curve to be measured and the second mean value, the negative pulse area of ​​the hydraulic pressure curve to be measured is calculated and obtained by the following formula:

[0027] ;

[0028] ;

[0029] in, Represents an array, is the negative pulse area of ​​the hydraulic curve; n is the data length of the hydraulic curve.

[0030] Optionally, the target data is subjected to curve smoothing to obtain a smooth curve, and the formula is as follows:

[0031] ;

[0032] Among them, p_min is the minimum position of the hydraulic value, and p_min-200 is the position after the minimum position of the hydraulic value is displaced by 200 points;

[0033] The pattern area formed by the smooth curve and the hydraulic pressure curve to be measured is calculated to obtain the jitter wave area of ​​the hydraulic pressure curve to be measured, and the formula is as follows:

[0034] ;

[0035] Among them, PeakSum is the jitter wave area, Pressure[i] is the hydraulic pressure curve to be tested, and Smooth[i] is the smooth curve.

[0036] Optionally, the first sample suction state is abnormal, including:

[0037] Obtaining the standard functional relationship between the sample aspiration volume and the sample negative pulse area;

[0038] Calculating the needle blocking calibration parameters through the standard functional relationship;

[0039] Acquire a standard hydraulic pressure curve, and calculate a needle blocking threshold value according to the standard hydraulic pressure curve and the needle blocking calibration parameter;

[0040] When the negative pulse area of ​​the hydraulic pressure curve to be measured is greater than the needle blocking threshold under the preset sensitivity, the first sampling state is abnormal;

[0041] The second sampling state is abnormal, including:

[0042] Obtain the standard hydraulic pressure curve and calculate the jitter wave area of ​​the air suction hydraulic pressure curve as the air suction calibration parameter;

[0043] Calculate and obtain the number of calibration jitter waves according to the air suction calibration parameters;

[0044] When the jitter wave area of ​​the hydraulic pressure curve to be measured is smaller than the calibrated jitter wave area under the preset sensitivity, or the number of jitter waves is smaller than a preset jitter wave number threshold, the second sampling state is abnormal.

[0045] Optionally, the standard functional relationship is a piecewise linear relationship;

[0046] Calculating the needle blocking threshold according to the needle blocking calibration parameter includes:

[0047] Through the standard functional relationship, the needle blocking calibration parameters are calculated, and the formula is:

[0048] ;

[0049] ;

[0050] ;

[0051] ;

[0052] Among them, Vol1 to Vol4, NegSum1 to NegSum4 represent the four known sample aspiration volumes and sample negative pulse areas respectively;

[0053] Using the needle blocking calibration parameters to describe the standard functional relationship, obtaining a first straight line equation and a second straight line equation;

[0054] The sample aspirated amount of the sampling needle during sample aspirating is substituted into the first linear equation or the second linear equation to obtain the needle blocking threshold.

[0055] Optionally, the number of calibration jitter waves is calculated according to the air suction calibration parameters, and the formula is:

[0056] ;

[0057] Among them, PeakSum_Th is the air suction calibration parameter obtained based on the air suction hydraulic pressure curve, and K is the adjustable coefficient.

[0058] Optionally, before obtaining the standard hydraulic pressure curve, the following steps may be performed:

[0059] Detecting whether the currently prepared hydraulic pressure curve is a hydraulic pressure curve used under standard conditions when the same sample volume as that of the sampling needle is aspirated;

[0060] If yes, use the currently prepared hydraulic pressure curve as the standard hydraulic pressure curve;

[0061] If not, the currently prepared hydraulic pressure curve is replaced.

[0062] The third aspect of an embodiment of the present invention provides a terminal device, including a memory, a processor, and a computer program stored in the above memory and executable on the above processor, and when the above processor executes the above computer program, the steps of the method provided in the first aspect are implemented.

[0063] A fourth aspect of an embodiment of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method provided in the first aspect are implemented.

[0064] An embodiment of the present invention proposes a method for detecting sample suction abnormality, by analyzing the hydraulic pressure curve to be tested when the sampling needle sucks samples, two values ​​of negative pulse area and jitter wave area are obtained, and the negative pulse area is used as a value to detect whether the first sample suction state of the sampling needle is normal, and the jitter wave area is used as a value to detect whether the second sample suction state of the sampling needle is normal, wherein the first sample suction state is abnormal, indicating that the sampling needle is in a blocked needle state; the second sample suction state is abnormal, indicating that the sampling needle is in an empty suction state. Therefore, the embodiment of the present invention uses a relatively small amount of calculation to achieve abnormality detection, is not limited to a specific analyzer, and has a wide range of applications, so that when the sampling needle is blocked or empty suction occurs, the system can promptly perform abnormal processing and give an alarm prompt, thereby reducing the probability of incorrect sample suction volume reported by the system due to blockage or empty suction, and improving the accuracy of the test results of the analyzer. BRIEF DESCRIPTION OF THE DRAWINGS

[0065] Figure 1 A schematic diagram of the implementation flow of the method for detecting abnormal sampling provided by an embodiment of the present invention;

[0066] Figure 2 A schematic diagram of the structure of a sampling system used in the method for detecting abnormal sampling provided by an embodiment of the present invention;

[0067] Figure 3 is a curve graph of the first mean;

[0068] Figure 4 is a curve graph of the second mean;

[0069] Figure 5 A curve diagram of the negative pulse area of ​​the hydraulic pressure curve to be measured;

[0070] Figure 6 It is a comparison diagram of the hydraulic pressure curve under normal suction and the hydraulic pressure curve under empty suction;

[0071] Figure 7 It is a schematic diagram of hydraulic pressure curve, smooth curve and jitter wave area based on normal sampling;

[0072] Figure 8 It is a schematic diagram of hydraulic pressure curve, smooth curve and jitter wave area based on the empty suction condition;

[0073] Fig. 9 It is a schematic diagram of the standard function relationship.

[0074] The realization of the purpose, functional features and advantages of the present invention will be further explained in conjunction with embodiments and with reference to the accompanying drawings. DETAILED DESCRIPTION

[0075] It should be understood that the specific embodiments described herein are only used to explain the present invention, and are not used to limit the present invention.

[0076] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element.

[0077] In this document, suffixes such as "module", "component" or "unit" used to represent elements are only used to facilitate the description of the present invention, and have no specific meanings. Therefore, "module" and "component" can be used interchangeably.

[0078] like Figure 1 As shown, the embodiment of the present invention provides a method for detecting abnormal sampling, including but not limited to the following steps:

[0079] S101, obtaining a hydraulic pressure curve to be measured when the sampling needle aspirates a sample.

[0080] In the embodiment of the present invention, by Figure 2 The sampling system 200 shown implements the sample suction abnormality detection method of the analyzer. The sampling system 200 includes a sampling needle 201, a hydraulic sensor 202, a detection plate 203 and a sample suction power source 204, wherein the sampling needle 201 is the sampling needle in the above step S101, and the hydraulic sensor 202 is used to provide the hydraulic curve to be tested in the above step S101.

[0081] S102, obtaining a first sampling state through the negative pulse area of ​​the hydraulic pressure curve to be measured; and obtaining a second sampling state based on the jitter wave area and the number of jitter waves of the hydraulic pressure curve to be measured.

[0082] In a specific application, when the sampling needle absorbs a normal sample, different sample aspiration amounts form different liquid pressure curves. The larger the sample aspiration amount, the larger the negative pulse area. In the embodiment of the present invention, this feature is used to determine whether the sampling needle is blocked, that is, whether the first sample aspiration state in the following step S103 is abnormal.

[0083] The embodiment of the present invention provides a method for calculating the negative pulse area of ​​the hydraulic pressure curve to be measured, that is, a method for calculating the negative pulse area of ​​the hydraulic pressure curve to be measured according to the hydraulic pressure curve to be measured, which, before the above step S102, includes:

[0084] Calculating a first mean value of the hydraulic pressure curve to be measured;

[0085] removing a curve portion whose hydraulic value is less than the first mean value in the hydraulic pressure curve to be measured, calculating the mean value of the remaining curve portion, and obtaining a second mean value;

[0086] Calculating and obtaining the negative pulse area of ​​the hydraulic pressure curve to be measured according to the pattern area formed by the hydraulic pressure curve to be measured and the second mean value;

[0087] The first mean value of the hydraulic pressure curve to be measured is calculated using the following formula:

[0088] ;

[0089] Wherein, Pressure[i] is the hydraulic pressure curve to be tested, n is the data length of the hydraulic pressure curve to be tested, and average1 is the first average value of the hydraulic pressure curve to be tested;

[0090] The second mean is calculated as follows:

[0091] ;

[0092] Among them, is the number of data with hydraulic pressure values ​​greater than average1, and average2 is the second mean value of the hydraulic pressure curve to be tested;

[0093] According to the pattern area formed by the hydraulic pressure curve to be measured and the second mean value, the negative pulse area of ​​the hydraulic pressure curve to be measured is calculated and obtained by the following formula:

[0094] ;

[0095] ;

[0096] in, Represents an array, is the negative pulse area of ​​the hydraulic curve; n is the data length of the hydraulic curve.

[0097] like Figures 3 to 5 As shown, the embodiment of the present invention also exemplarily shows a curve diagram of the first mean value, the second mean value and the negative pulse area of ​​the hydraulic pressure curve to be measured, Figure 3In the figure, the thin solid line represents the hydraulic pressure curve to be measured, the thick solid line represents the first mean value, the X-axis represents the sampling point, that is, the sequence number of the hydraulic point, which is equivalent to time, and the Y-axis represents the hydraulic pressure AD value, that is, the voltage AD value output by the analog-to-digital converter, the maximum value is 65535, and the minimum value is 0, which is equivalent to the hydraulic pressure value; Figure 4 In the figure, the thin solid line represents the hydraulic pressure curve to be measured, the thick solid line represents the second mean, the X-axis represents the sampling points, and the Y-axis represents the hydraulic pressure AD value; Figure 5 In the figure, the solid line represents the hydraulic pressure curve to be measured, the dotted line represents the second mean, the shaded part represents the negative pulse area, the X-axis represents the sampling point, and the Y-axis represents the hydraulic pressure AD value.

[0098] like Figure 6 As shown, an embodiment of the present invention exemplarily shows a comparison diagram of the hydraulic curve during normal sample suction and the hydraulic curve during empty suction. In a specific application, during normal sample suction, the liquid suction starting position of the hydraulic curve will have violent shaking, while during empty suction, the liquid suction starting position of the hydraulic curve will not have this feature. In an embodiment of the present invention, this feature is used to distinguish between empty suction and normal sample suction, that is, whether the second sample suction state in the following step S103 is abnormal.

[0099] An embodiment of the present invention provides a method for calculating the jitter wave area and the number of jitter waves of the hydraulic curve to be measured, that is, a method for calculating the jitter wave area and the number of jitter waves of the hydraulic curve to be measured according to the hydraulic curve to be measured, which, before the above-mentioned step S102, includes:

[0100] Finding the minimum position of the hydraulic pressure value in the hydraulic pressure curve to be tested;

[0101] Starting from the minimum position of the hydraulic pressure value and heading leftward, 200 points in the hydraulic pressure curve to be measured are selected as target data, and curve smoothing is performed on the target data to obtain a smooth curve;

[0102] In the range of the minimum position of the hydraulic value and the displacement of 200 points of the minimum position of the hydraulic value, the pattern area formed by the smooth curve and the hydraulic pressure curve to be measured is calculated to obtain the jitter wave area of ​​the hydraulic pressure curve to be measured;

[0103] The number of jitter waves whose peak values ​​of a pattern formed by the smooth curve and the hydraulic pressure curve to be measured exceed a set threshold is calculated to obtain the number of jitter waves of the hydraulic pressure curve to be measured.

[0104] The target data is subjected to curve smoothing to obtain a smooth curve, and the formula is as follows:

[0105] ;

[0106] Among them, p_min is the minimum position of the hydraulic value, and p_min-200 is the position after the minimum position of the hydraulic value is displaced by 200 points;

[0107] The pattern area formed by the smooth curve and the hydraulic pressure curve to be measured is calculated to obtain the jitter wave area of ​​the hydraulic pressure curve to be measured, and the formula is as follows:

[0108] ;

[0109] Among them, PeakSum is the jitter wave area, Pressure[i] is the hydraulic pressure curve to be tested, and Smooth[i] is the smooth curve.

[0110] like Figure 7 As shown, the embodiment of the present invention uses the hydraulic pressure curve under normal sampling conditions to exemplarily show the smooth curve, hydraulic pressure curve and jitter wave area schematic diagram obtained based on this; Figure 8 As shown, the embodiment of the present invention also uses the hydraulic pressure curve under the empty suction condition to exemplarily show the smooth curve, hydraulic pressure curve and jitter wave area schematic diagram obtained based on this. Figure 7 In the figure, the thin solid line represents the hydraulic pressure curve under normal conditions, and the thick solid line represents the smooth curve; the shaded part represents the jitter wave area, the X-axis represents the sampling point, that is, the sequence number of the hydraulic pressure point, which is equivalent to time, and the Y-axis represents the hydraulic pressure AD value, that is, the voltage AD value output by the analog-to-digital converter, the maximum value is 65535, and the minimum value is 0, which is equivalent to the hydraulic pressure value; Figure 8 In the figure, the thin solid line represents the hydraulic pressure curve under the empty suction condition, and the thick solid line represents the smooth curve; the shaded part represents the jitter wave area, the X-axis represents the sampling point, and the Y-axis represents the hydraulic pressure AD value.

[0111] S103: When the first sample aspirating state or the second sample aspirating state is abnormal, abnormal processing is performed and an alarm prompt is given.

[0112] In a specific application, the detailed implementation of step S103 may be:

[0113] Obtaining the first sampling state;

[0114] When the first sampling state is abnormal, abnormal processing is performed and an alarm prompt is given;

[0115] After the exception is handled, the first sampling state is reacquired until the first sampling state is normal;

[0116] When the first sample suction state is normal, the second sample suction state is obtained;

[0117] When the second sampling state is abnormal, abnormal processing is performed and an alarm prompt is given;

[0118] After the exception is handled, the second sampling state is acquired again until the second sampling state is normal.

[0119] It should be noted that if the alarm prompt is continuously given, it means that the sampling system cannot automatically complete the abnormal processing and needs to be processed by the user. Therefore, preferably, the alarm prompt can also include a first alarm message and a second alarm message. When the alarm prompt is the first alarm message, the abnormal processing is automatically performed. When the alarm prompt is the second alarm message, the user is notified to process. In the embodiment of the present invention, an alarm is given for abnormal sampling, which not only ensures the accuracy of sampling, but also protects the sampling pipeline and extends the service life of the equipment.

[0120] In the above step S103, if the first sample suction state is abnormal, it means that the sampling needle is in a needle blocking state, and the abnormal processing based on it includes but is not limited to:

[0121] (1) Spit the sample that has been sucked back into the sample tube;

[0122] (2) Adjust the needle insertion position, for example, to the left; adjust the needle insertion depth, for example, to deeper.

[0123] If the second sample suction state is abnormal, it means that the sampling needle is in an empty suction state, and the abnormal processing based on it includes but is not limited to: re-inserting the needle and adjusting the depth of the needle.

[0124] The embodiment of the present invention obtains the needle blocking state by analyzing the abnormal value of the negative pulse area of ​​the hydraulic pressure curve to be tested, and obtains the empty suction state by analyzing the abnormal value of the jitter wave area of ​​the hydraulic pressure curve to be tested.

[0125] Among them, in order to obtain the abnormality of the negative pulse area of ​​the hydraulic curve to be tested, first, based on the hydraulic curves with different sample suction amounts under normal sample suction conditions, the functional relationship between the sample suction amount and the negative pulse area is obtained, and then the hydraulic curve with the same sample suction amount as the hydraulic curve to be tested under normal sample suction conditions is taken as the standard hydraulic curve, and the negative pulse area is compared with the hydraulic curve to be tested. The comparison result is used to determine whether the first sample suction state is abnormal.

[0126] Therefore, in one embodiment, the first sample suction state is abnormal, including:

[0127] S01, obtaining a standard functional relationship between the sample aspiration amount and the sample negative pulse area;

[0128] S02, calculating the needle blocking calibration parameter through the standard function relationship;

[0129] S03, obtaining a standard hydraulic pressure curve, and calculating a needle blocking threshold value according to the standard hydraulic pressure curve and the needle blocking calibration parameter;

[0130] S04. When the negative pulse area of ​​the hydraulic pressure curve to be measured is greater than the needle blocking threshold under a preset sensitivity, the first sampling state is abnormal.

[0131] In the embodiment of the present invention, the standard hydraulic curve is a hydraulic curve with the same sample suction volume as the hydraulic curve to be measured under normal sample suction conditions. Therefore, before obtaining the standard hydraulic curve, the following steps are included:

[0132] Detecting whether the currently prepared hydraulic pressure curve is a hydraulic pressure curve used under standard conditions when the same sample volume as that of the sampling needle is aspirated;

[0133] If yes, use the currently prepared hydraulic pressure curve as the standard hydraulic pressure curve;

[0134] If not, the currently prepared hydraulic pressure curve is replaced.

[0135] In the embodiment of the present invention, the standard functional relationship is a piecewise linear relationship; therefore, the formula for calculating the needle blocking calibration parameter through the standard functional relationship is:

[0136] ;

[0137] ;

[0138] ;

[0139] ;

[0140] Among them, Vol1 to Vol4, NegSum1 to NegSum4 represent the four known sample aspiration volumes and sample negative pulse areas respectively;

[0141] The detailed implementation method of calculating the needle blocking threshold according to the standard hydraulic pressure curve and the needle blocking calibration parameter is as follows:

[0142] Using the needle blocking calibration parameters to describe the standard functional relationship, obtaining a first straight line equation and a second straight line equation;

[0143] The sample aspirated amount of the sampling needle during sample aspirating is substituted into the first linear equation or the second linear equation to obtain the needle blocking threshold.

[0144] In the embodiment of the present invention, the hydraulic curves with the sample suction volumes of Vol1=20uL, Vol2=50uL, Vol3=100uL, and Vol4=180uL are taken, and the negative pulse areas NegSum1, NegSum2, NegSum3, and NegSum4 of Vol1, Vol2, Vol3, and Vol4 are calculated, wherein Vol1, Vol2, Vol3, and Vol4 are used as the sample suction volumes, and NegSum1, NegSum2, NegSum3, and NegSum4 are used as the sample negative pulse areas, forming points P1 (Vol1, NegSum1), P2 (Vol2, NegSum2), P3 (Vol3, NegSum3), and P4 (Vol4, NegSum4). The piecewise linear relationship determined by points P1 to P4 is the standard functional relationship between the sample suction volume and the sample negative pulse area, and the needle blocking calibration parameters can be calculated.

[0145] For example, Fig. 9 It is a schematic diagram of the standard functional relationship of the sample aspiration volume, wherein the X-axis represents the aspiration volume, the unit is uL, the Y-axis represents the negative pulse area, which is the virtual area calculated according to the hydraulic curve, K1=8337.4, B1=141681, K2=10137, B2=212396, the rectangular line represents the first section of the functional relationship, that is, the first straight line equation, which is y=83337.4x+14681, which is used when the aspiration volume is 10~75ul; the dotted line represents the second section of the functional relationship, that is, the second straight line equation, which is y=10137x+212396, which is used when the aspiration volume is 80~205ul.

[0146] Assuming that the sample volume used for the hydraulic pressure curve to be tested is 55uL, substitute 55 into the first straight line equation to obtain the needle blocking threshold value NegSum_Th corresponding to the sample volume of 55uL=55*8337.4+141681=600238.

[0147] Then, the negative pulse area of ​​the hydraulic pressure curve to be tested is calculated using the above calculation method, and is set as NegSum_Test. When the negative pulse area of ​​the hydraulic pressure curve to be tested is greater than the needle blocking threshold under the preset sensitivity, the first sampling state is abnormal, which can be expressed by the formula:

[0148] ;

[0149] Among them, sens1 is the preset sensitivity.

[0150] In specific applications, if the jitter wave area of ​​the hydraulic curve to be tested is abnormal, we must first obtain a standard hydraulic curve. The hydraulic curves are different under normal suction and empty suction conditions, that is, under normal suction, there will be jitter waves at the initial position of the curve. Therefore, the empty suction calibration parameters can be calculated based on the hydraulic curve under normal suction as the standard hydraulic curve. Since the empty suction characteristics under different suction amounts are stable and do not change with the suction amount, there is no need to calculate the empty suction thresholds corresponding to different suction amounts. The jitter wave area of ​​the standard hydraulic curve can be directly compared with the hydraulic curve to be tested in terms of the jitter wave area and the number of jitter waves. The comparison result is used to determine whether the second suction state is abnormal.

[0151] In one embodiment, the second sample suction state is abnormal, including:

[0152] S11, obtaining a standard hydraulic pressure curve, and calculating the jitter wave area of ​​the standard hydraulic pressure curve as an air suction calibration parameter;

[0153] S12, calculating and obtaining the number of calibration jitter waves according to the air suction calibration parameters;

[0154] S13: When the jitter wave area of ​​the hydraulic pressure curve to be measured is smaller than the calibrated jitter wave area under a preset sensitivity, or the number of jitter waves is smaller than a preset jitter wave number threshold, the second sampling state is abnormal.

[0155] The number of calibrated jitter waves is calculated according to the air suction calibration parameters, and the formula is:

[0156] ;

[0157] Among them, PeakSum_Th is the air suction calibration parameter obtained based on the air suction hydraulic pressure curve, and K is the adjustable coefficient.

[0158] An embodiment of the present invention also provides a terminal device including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, each step of the sampling abnormality detection method described in the above embodiment is implemented.

[0159] An embodiment of the present invention further provides a storage medium, which is a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, each step in the method for detecting abnormal sampling as described in the above embodiment is implemented.

[0160] The embodiments described above are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the foregoing embodiments describe the present invention in detail, those skilled in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or replace some of the technical features therein by equivalents. Such modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included in the protection scope of the present invention.

Claims

1. A method for detecting abnormality in suction sampling, characterized in that: include: Obtain the hydraulic pressure curve to be measured when the sampling needle aspirates the sample; Calculating the negative pulse area of ​​the hydraulic pressure curve to be measured according to the hydraulic pressure curve to be measured includes: Calculating a first mean value of the hydraulic pressure curve to be measured; removing the curve portion whose hydraulic value is less than the first mean value in the hydraulic pressure curve to be measured, calculating the mean value of the remaining curve portion, and obtaining a second mean value; Calculating and obtaining the negative pulse area of ​​the hydraulic pressure curve to be measured according to the pattern area formed by the hydraulic pressure curve to be measured and the second mean value; Acquiring a first sampling state through the negative pulse area of ​​the hydraulic pressure curve to be measured; The second sampling state is obtained based on the shaking wave area and the number of shaking waves of the hydraulic pressure curve to be measured, including: Finding the minimum position of the hydraulic pressure value in the hydraulic pressure curve to be measured; Starting from the minimum position of the hydraulic pressure value and heading leftward, 200 points in the hydraulic pressure curve to be measured are selected as target data, and curve smoothing is performed on the target data to obtain a smooth curve; In the range of the minimum position of the hydraulic value and the displacement of 200 points of the minimum position of the hydraulic value, the pattern area formed by the smooth curve and the hydraulic pressure curve to be measured is calculated to obtain the jitter wave area of ​​the hydraulic pressure curve to be measured; Calculating the number of jitter waves whose peak values ​​of the pattern formed by the smooth curve and the hydraulic pressure curve to be measured exceed a set threshold, and obtaining the number of jitter waves of the hydraulic pressure curve to be measured; When the first sample suction state or the second sample suction state is abnormal, abnormal processing is performed and an alarm prompt is given; If the first sample suction state is abnormal, it means that the sampling needle is in a blocked state; if the second sample suction state is abnormal, it means that the sampling needle is in an empty suction state.

2. The method for detecting abnormality in aspiration of a sample according to claim 1, characterized in that: The first mean value of the hydraulic pressure curve to be measured is calculated using the following formula: ； Wherein, Pressure[i] is the hydraulic pressure curve to be tested, n is the data length of the hydraulic pressure curve to be tested, and average1 is the first average value of the hydraulic pressure curve to be tested; The second mean is calculated as follows: ； in, is the number of data with hydraulic pressure values ​​greater than average1, and average2 is the second average value of the hydraulic pressure curve to be tested; According to the pattern area formed by the hydraulic pressure curve to be measured and the second mean value, the negative pulse area of ​​the hydraulic pressure curve to be measured is calculated and obtained by the following formula: ; ； in, Represents an array, is the negative pulse area of ​​the hydraulic curve; n is the data length of the hydraulic curve.

3. The method for detecting abnormality in aspiration of a sample according to claim 1, characterized in that: The target data is subjected to curve smoothing to obtain a smooth curve, and the formula is as follows: ； Among them, p_min is the minimum position of the hydraulic value, and p_min-200 is the position after the minimum position of the hydraulic value is displaced by 200 points; The pattern area formed by the smooth curve and the hydraulic pressure curve to be measured is calculated to obtain the jitter wave area of ​​the hydraulic pressure curve to be measured, and the formula is as follows: ； Among them, PeakSum is the jitter wave area, Pressure[i] is the hydraulic pressure curve to be tested, and Smooth[i] is the smooth curve.

4. The method for detecting abnormality in aspiration of a sample according to claim 1, wherein: The first sample suction state is abnormal, including: Obtaining the standard functional relationship between the sample aspiration volume and the sample negative pulse area; Calculating the needle blocking calibration parameters through the standard functional relationship; Acquire a standard hydraulic pressure curve, and calculate a needle blocking threshold value according to the standard hydraulic pressure curve and the needle blocking calibration parameter; When the negative pulse area of ​​the hydraulic pressure curve to be measured is greater than the needle blocking threshold under the preset sensitivity, the first sampling state is abnormal; The second sampling state is abnormal, including: Obtain the standard hydraulic pressure curve and calculate the jitter wave area of ​​the air suction hydraulic pressure curve as the air suction calibration parameter; Calculate and obtain the number of calibration jitter waves according to the air suction calibration parameters; When the jitter wave area of ​​the hydraulic pressure curve to be measured is smaller than the calibrated jitter wave area under the preset sensitivity, or the number of jitter waves is smaller than a preset jitter wave number threshold, the second sampling state is abnormal.

5. The method for detecting abnormality in aspiration of a sample according to claim 4, characterized in that: The standard functional relationship is a piecewise linear relationship; Calculating the needle blocking threshold according to the needle blocking calibration parameter includes: Through the standard functional relationship, the needle blocking calibration parameters are calculated, and the formula is: ； ； ； ； Among them, Vol1 to Vol4, NegSum1 to NegSum4 represent the four known sample aspiration volumes and sample negative pulse areas respectively; Using the needle blocking calibration parameters to describe the standard functional relationship, obtaining a first straight line equation and a second straight line equation; The sample aspirated amount of the sampling needle during sample aspirating is substituted into the first linear equation or the second linear equation to obtain the needle blocking threshold.

6. The method for detecting abnormality in aspiration of a sample according to claim 4, characterized in that: The calibration jitter wave quantity is calculated according to the air suction calibration parameters, and the formula is: ； Among them, PeakSum_Th is the air suction calibration parameter obtained based on the standard hydraulic curve, and K is the adjustable coefficient.

7. The method for detecting abnormality in aspiration of a sample according to claim 4, characterized in that: Before obtaining the standard hydraulic curve, include: Detecting whether the currently prepared hydraulic pressure curve is a hydraulic pressure curve used under standard conditions when the same sample volume as that of the sampling needle is aspirated; If yes, use the currently prepared hydraulic pressure curve as the standard hydraulic pressure curve; If not, the currently prepared hydraulic pressure curve is replaced.

8. A terminal device, characterized in that: The method comprises a memory, a processor and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, each step of the method for detecting abnormal sampling as claimed in any one of claims 1 to 7 is implemented.

9. A storage medium, the storage medium being a computer-readable storage medium, on which a computer program is stored, characterized in that: When the computer program is executed by a processor, each step of the method for detecting abnormal sampling according to any one of claims 1 to 7 is implemented.

Images