Shear wave propagation velocity estimation device

Abstract

To provide a shear wave propagation velocity estimation device capable of simply and highly efficiently measuring elasticity of lung parenchyma.SOLUTION: A shear wave propagation velocity estimation device comprises: a vibration exciter for exciting shear waves in a measuring object; an ultrasonic probe for outputting an ultrasonic pulse to the measuring object and receiving an ultrasonic signal in the measuring object; an ultrasonic Doppler device for acquiring an ultrasonic Doppler signal from the received ultrasonic signal; and a velocity estimator for estimating a shear wave propagation velocity from the ultrasonic Doppler signal acquired by the ultrasonic Doppler device. The velocity estimator estimates phases of the shear waves in a plurality of or one wide B-line from the ultrasonic Doppler signal, and estimates the shear wave propagation velocity from the phases of the shear waves in the plurality of or one wide B-line.SELECTED DRAWING: Figure 1

Description

[Technical Field]

[0001] The present invention relates to a shear wave velocity estimation device that uses ultrasound to estimate the propagation velocity of shear waves inside a measurement object. The shear wave velocity estimation device of the present invention is suitable for diagnosing lung tissue. [Background technology]

[0002] Trauma and infections such as COVID-19 can cause acute, severe inflammation in lung tissue, leading to acute respiratory distress syndrome (ARDS). The primary pathology of ARDS is non-cardiogenic pulmonary edema, and as inflammation persists, lung tissue becomes fibrotic, leading to insufficient ventilation and a mortality rate of 30-40%.

[0003] ARDS is a syndrome caused by a wide variety of factors, and as there are few established treatments, patients are given symptomatic treatment such as artificial respiration while waiting for the lungs to recover. The treatment period is not fixed, and there are cases that last for more than several months, so the mortality rate is high. Therefore, the development of diagnostic techniques that can assess the disease progression and prognosis of ARDS is an urgent issue not only in peacetime but also in situations where medical resources are depleted, such as during the COVID-19 pandemic.

[0004] Chest X-ray CT is the primary diagnostic method for acute lung diseases, including ARDS. However, chest X-ray CT scans irradiate radiation and must be performed in a dedicated examination room. This requires patients with poor breathing to be moved to the examination room, and there is also a risk of radiation exposure, so the scan cannot be performed frequently. In addition, since the patient's condition changes every moment in ARDS, it is desirable to develop an evaluation method that can be performed quickly and easily at the bedside without moving the patient.

[0005] Currently, an assessment method that can be used at the bedside is an oxygen saturation meter (pulse oximeter) that measures blood oxygen saturation. However, the oximeter evaluates the respiratory condition, which is a result of lung disease, and does not evaluate the progression or prognosis of lung disease. It is also difficult to evaluate the area where dysfunction is occurring.

[0006] The functional requirements for a diagnostic method for acute lung disease are: (1) since acute lung disease can occur in anyone, anywhere, at any time, the device must be easy to use, (2) there must be no problem with radiation exposure, and the method must be highly safe for the living body; and (3) the localization of the lesion must be able to be evaluated.

[0007] On the other hand, it has been reported that in acute lung diseases, particularly ARDS, the elasticity (stiffness) of lung tissue increases as the disease progresses. Therefore, attempts have been made to quantitatively measure lung elasticity and use it to evaluate the progression and prognosis of lung diseases.

[0008] Among the methods for quantitatively measuring lung elasticity, one that can also be applied in vivo is a method in which vibrations of approximately 40 to 200 Hz are applied using a vibrator placed on the surface of the living body, and the mechanical vibration waves (shear waves) that are transmitted to the lung parenchyma are visualized using nuclear magnetic resonance imaging (MRI), and elasticity is measured from the propagation speed of the shear waves (see, for example, Non-Patent Documents 1 and 2). However, this method requires a large, non-portable MRI diagnostic device, does not satisfy the above-mentioned functional requirements, and is therefore not a desirable diagnostic method.

[0009] On the other hand, there have been attempts to apply vibrations to the surface of a living body and measure the shear waves that travel through the pleura on the surface of the lungs using an ultrasound examination device (see, for example, Non-Patent Documents 3 and 4). However, because ultrasound cannot reach the air-rich lung parenchyma, this method can only measure the sum of the elasticity of multiple biological tissues in the immediate vicinity of the pleura, and is therefore not suitable for measuring the elasticity of the lung parenchyma in clinically relevant ways.

[0010] Ultrasound examinations do not involve radiation exposure and can evaluate the lungs by placing a probe on the chest. Recently, ultrasound examination devices have become smaller, allowing results to be obtained immediately at the bedside, enabling pathological conditions to be assessed and diagnosed. Therefore, point-of-care ultrasound (POCUS) is widely used in emergency intensive care, disaster sites, home medical care, etc. (see, for example, Non-Patent Document 5). Acute lung diseases cause specific artifacts, such as "B-lines," to appear on normal ultrasound images, so ultrasound devices are also used to diagnose acute lung diseases. [Prior art documents] [Non-patent literature]

[0011] [Non-Patent Document 1] Y. Mariappan, K. Glaser, R. Hubmayer, A. Manduca, R. Ehman and K. McGee, “MR Elastography of Human Lung Parenchyma: Technical Development, Theoretical Modeling and In vivo Validation”, J. Magn Reson Imaging, 33: 1351-1361 (2011). [Non-patent document 2] J. Marineli, D. Levin, R. Vassallo, R. Carter, R. Humbmayr, R. Ehman and K. McGee, “Quantitative assessment of lung stiffness in patients with interstitial lung disease using MR elastography”, J Magn Reason Imaging, 2017, 46(2);365-374. [Non-patent document 3] X. Zhang, T. Osborn, B. Zhou, D. Meixner, R. Kinnick, B. Bartholmai, J. Greenleaf and S. Kalra, “Lung ultrasound surface wave elastography: a pilot clinical study”, IEEE Trans Ultrason Ferroelectr Freq Control. 2017 Sep; 64(9): 1298-1304. [Non-patent document 4] X. Zhang, B. Zhou, S. Karla, B. Bartholnai, J. Greenleaf, and T. Osborn, “An ultrasonic surface wave technique for assessing skin and lung diseases”, Ultrasound Med Biol. 2018, 44(2):321-331. [Non-Patent Document 5] Moore CL, Copel JA. Point-of-care ultrasonography. N Engl J Med. 2011 Feb 24;364(8):749-57. Summary of the Invention [Problem to be solved by the invention]

[0012] However, the above-mentioned methods are limited to evaluating the shape or the elasticity of the pleura on the surface of the lung or its immediate vicinity, making it difficult to measure the elasticity of lung tissue, which is expected to lead to a functional diagnosis of the lung. Therefore, it is desirable to establish as soon as possible a method for measuring the elasticity of lung parenchyma using equipment that can be used by anyone, anywhere, anytime.

[0013] In order to solve the above-mentioned problems, an object of the present invention is to provide a shear wave velocity estimation device that can measure the elasticity of lung parenchyma simply and accurately. [Means for solving the problem]

[0014] The shear wave propagation velocity estimation device of the present invention comprises an exciter for exciting shear waves in a measurement object, an ultrasonic probe for outputting ultrasonic pulses to the measurement object and receiving ultrasonic signals within the measurement object, an ultrasonic Doppler device for obtaining ultrasonic Doppler signals from the ultrasonic signals received by the ultrasonic probe, and a velocity estimation device for estimating the propagation velocity of the shear wave from the ultrasonic Doppler signals obtained by the ultrasonic Doppler device, and the velocity estimation device estimates the phase of the shear wave in multiple B-lines or a single wide B-line from the ultrasonic Doppler signals, and estimates the propagation velocity of the shear wave from the phase of the shear wave in multiple B-lines or a single wide B-line. [Effects of the Invention]

[0015] According to the shear wave velocity estimation device of the present invention described above, radiation is not used, so there is no need to provide a specialized examination room as in the case of a chest X-ray CT scan. Furthermore, since measurements can be performed by operating a vibrator and an ultrasonic probe, measurements can be performed easily. Since the propagation velocity of the shear wave in multiple B-lines or a single wide B-line is estimated, the elasticity of the lung parenchyma, which is the measurement target, can be estimated from the propagation velocity of the shear wave, thereby enabling accurate measurement of the elasticity of the lung parenchyma. Therefore, the present invention makes it possible to measure the elasticity of the lung parenchyma simply and accurately. [Brief explanation of the drawings]

[0016] [Figure 1] 1 is a diagram showing a schematic configuration of an embodiment of a shear wave propagation velocity estimation device of the present invention, together with a cross-sectional structure of the lung and surrounding body tissues. FIG. [Figure 2] 1A is a schematic diagram of a reflected ultrasonic signal at an ultrasonic scanning line where a B line appears, and FIG. 1B is a schematic diagram of a reflected ultrasonic signal at an ultrasonic scanning line where a B line does not appear. [Figure 3] FIG. 10 is a diagram showing a lung echo model in which two B-lines are formed. [Figure 4] 2 is a flowchart of signal processing in the shear wave propagation velocity estimation device of FIG. 1. [Figure 5] FIG. 1 is a diagram showing a schematic diagram of a phantom used in an experiment and an outline of an experimental device. [Figure 6] A is a B-mode image of a normal sponge phantom, and B is a B-mode image of a hardened sponge phantom. [Figure 7] This is a power Doppler image obtained when vibration was applied to a normal sponge phantom. [Figure 8] FIG. 10 shows the frame-by-frame changes in brightness values ​​of power Doppler images appearing on each B-line when a normal sponge phantom is vibrated. [Figure 9] This is a power Doppler image obtained when vibration was applied to a hardened sponge phantom. [Figure 10] FIG. 10 shows the frame-by-frame transition of the brightness value of the power Doppler image appearing on each B-line when the hardened sponge phantom is vibrated. [Figure 11] FIG. 10 shows the results of estimating the propagation velocity of shear waves from the phase difference of shear waves observed between B-lines in a normal sponge phantom and a hardened sponge phantom. DETAILED DESCRIPTION OF THE INVENTION

[0017] The shear wave propagation velocity estimation device of the present invention includes an exciter for exciting shear waves in a measurement object, an ultrasonic probe for outputting ultrasonic pulses to the measurement object and receiving ultrasonic signals within the measurement object, an ultrasonic Doppler device for obtaining ultrasonic Doppler signals from the ultrasonic signals received by the ultrasonic probe, and a velocity estimation device for estimating the propagation velocity of the shear wave from the ultrasonic Doppler signals obtained by the ultrasonic Doppler device. In the shear wave propagation velocity estimation device of the present invention, the velocity estimation device estimates the phase of the shear wave in multiple or one wide B-line from the ultrasonic Doppler signal, and estimates the propagation velocity of the shear wave from the phase of the shear wave in multiple or one wide B-line.

[0018] The shear wave velocity estimation device of the present invention focuses on B-lines, which are prominent findings that appear in ultrasound B-mode images of acute pulmonary diseases, and estimates the phase and propagation velocity of shear waves propagating through the pulmonary parenchyma from Doppler signals obtained on the B-lines. The shear wave velocity estimation device of the present invention uses ultrasound to estimate the phase and propagation velocity of shear waves propagating through the pulmonary parenchyma, making it possible to measure the elasticity of the pulmonary parenchyma from the propagation velocity of the shear waves.

[0019] B-lines are artifacts caused by multiple reflections of ultrasound waves within the alveoli, which are filled with fluid. Conversely, if a Doppler signal is detected on an ultrasound scanning line where a B-line appears, the signal obtained is the Doppler signal from the lung parenchyma itself. Therefore, for lung diseases where multiple B-lines appear, if the phase of the shear wave is measured on the B-line, the elasticity of the lung parenchyma can be measured from the phase difference between the multiple B-lines. Furthermore, even when multiple B-lines are connected to form a single wide B-line, if there is a phase difference in the shear waves within that single wide B-line, the elasticity of the lung parenchyma can be measured from the phase difference in the shear waves of that single wide B-line. In other words, this measurement method focuses on the B-line, which is an artifact of ultrasound images but contains Doppler information about the lung parenchyma due to its generation process, and is a way to measure the elasticity of the lung parenchyma, which could not be measured using ultrasound until now. This measurement method is currently the only method that can measure the elasticity of lung parenchyma, which changes due to acute lung disease, using ultrasound, which has previously been thought to be difficult to measure.

[0020] The propagation velocity of shear waves in multiple or a single wide B-line is estimated from the phase difference of shear waves in multiple or a single wide B-line. The following two methods can be used to estimate the propagation velocity of shear waves from the phase difference of shear waves in multiple B-lines or a single wide B-line.

[0021] The first method is to perform Fourier analysis on the Doppler signal. Specifically, for example, the phase of a complex quadrature detection signal (complex Doppler signal) is determined by Fourier analysis at the excitation frequency (the frequency of the shear wave from the exciter). However, the phase may also be determined from the time at which the phase reaches a predetermined phase value (for example, zero) without using Fourier analysis.

[0022] The second method is a method using the C-SWE method developed by the inventors of the present invention. C-SWE (Continuous Shear Wave Elastography) is a shear wave imaging technique that utilizes the fact that when continuous shear waves are excited inside the living body from a small vibrator placed on the surface of the living body to be measured, shear wave images appear in color Doppler images for blood imaging, which is a standard function of general-purpose ultrasound echo devices (see, for example, the specification of Japanese Patent No. 649183). Note that power Doppler images may also be used instead of color Doppler images.

[0023] In the C-SWE method, as described in the specification of Japanese Patent No. 649183, for example, vibrations having a center frequency that is n / 4 times (where n is an odd number greater than or equal to 1) the burst frequency, which corresponds to the repetition period of the ultrasonic pulse, are applied to the object to be measured. Specifically, the method for estimating the propagation velocity of shear waves using the C-SWE method involves setting the frequency of the shear waves relative to the repetition frequency of the ultrasound waves to match the frequency conditions of the C-SWE method described above (n / 4 times the repetition frequency of the ultrasound waves (where n is an odd number greater than or equal to 1)), and then performing Fourier analysis on the color Doppler image or power Doppler image obtained in this way in the frame direction.

[0024] Compared to the first method, which uses Fourier analysis of Doppler signals, the second method, which uses C-SWE, has the advantage that it can be used with conventional ultrasound diagnostic equipment without modification, and is therefore easy to incorporate into small, portable ultrasound diagnostic equipment.

[0025] According to the configuration of the shear wave propagation velocity estimation device of the present invention described above, it includes an exciter for exciting shear waves in the measurement object, an ultrasonic probe for outputting ultrasonic pulses to the measurement object and receiving ultrasonic signals from within the measurement object, an ultrasonic Doppler device for obtaining ultrasonic Doppler signals from the ultrasonic signals received by the ultrasonic probe, and a velocity estimation device for estimating the propagation velocity of the shear wave from the ultrasonic Doppler signals obtained by the ultrasonic Doppler device. This makes it possible to estimate the propagation velocity of the shear wave in the object to be measured by using the shear wave excited in the object to be measured from the vibrator and the ultrasonic pulse output to the object to be measured from the ultrasonic probe. Since it uses shear waves and ultrasonic pulses, it does not use radiation, so there is no need to set up a specialized examination room like in chest X-ray CT.In addition, measurements can be made simply by operating the vibrator and ultrasonic probe, making measurements easy.

[0026] In the shear wave propagation velocity estimation device of the present invention, the velocity estimation device estimates the phase of the shear wave in multiple or one wide B-line from the ultrasonic Doppler signal, and estimates the propagation velocity of the shear wave from the phase of the shear wave in multiple or one wide B-line. This makes it possible to estimate the elasticity of the pulmonary parenchyma, which is the measurement target, from the propagation velocity of the shear wave estimated by the velocity estimation device, thereby enabling the elasticity of the pulmonary parenchyma to be measured with high accuracy. Therefore, the shear wave propagation velocity estimation device of the present invention makes it possible to measure the elasticity of the lung parenchyma simply and accurately.

[0027] Furthermore, the shear wave propagation velocity estimation device of the present invention may be configured so that the velocity estimation device estimates the propagation velocity of the shear wave by performing Fourier analysis on the phase of a complex quadrature detection signal obtained from an ultrasonic Doppler signal at the frequency of the shear wave. That is, the velocity estimation device is configured to estimate the propagation velocity of the shear wave by the first method described above. With this configuration, the propagation velocity of the shear wave can be estimated from the ultrasonic Doppler signal.

[0028] Furthermore, in the shear wave velocity estimation device of the present invention, the frequency of the shear wave can be set to n / 4 times the repetition frequency of the ultrasound (where n is an odd number equal to or greater than 1), and the velocity estimation device can estimate the phase of the shear wave in each B-line by performing Fourier analysis in the frame direction on a color Doppler image or a power Doppler image obtained from the ultrasound Doppler signal. That is, the velocity estimation device is configured to estimate the propagation velocity of the shear wave by the second method described above. With this configuration, the propagation velocity of the shear wave can be estimated from a color Doppler image or a power Doppler image obtained from the ultrasonic Doppler signal. Furthermore, as mentioned above, since it can be used without modifying conventional ultrasonic diagnostic equipment, there is an advantage that it can be easily incorporated into small, portable ultrasonic diagnostic equipment.

[0029] Hereinafter, specific embodiments of the shear wave propagation velocity estimation device of the present invention will be described.

[0030] FIG. 1 shows a schematic configuration diagram of an embodiment of a shear wave propagation velocity estimation device according to the present invention. As shown in FIG. 1, the shear wave propagation velocity estimation device 10 of this embodiment includes an oscillator 1, a small vibrator 2, an ultrasonic probe 3, an ultrasonic Doppler device 4, a velocity estimation device 5, and an output device 6. The functions of the components 1 to 6 of the shear wave propagation velocity estimation device 10 will be described later.

[0031] FIG. 1 also shows the cross-sectional structure of the lungs and surrounding body tissues, which are the objects to be examined by the shear wave propagation velocity estimation device 10. As shown in Figure 1, the cross-sectional structure of the lung consists of, from the outside in, a body surface 11, soft tissue 12, pleura 13, and lung parenchyma 14. Although not shown, lung parenchyma 14 contains many tiny alveoli that contain air inside.

[0032] When ultrasound is transmitted from the body surface 11 to the lung structure, the cross-sectional structure of which is shown in Figure 1, and a lung echo (chest echo) is obtained, a linear echo, not shown, extending in a lateral direction (horizontal direction) with high reflection strength due to soft tissue 12, called a reflected echo or pleural echo, is observed from a position close to the body surface 11.

[0033] In acute lung diseases, if fluid accumulates in some of the alveoli in the lung parenchyma 14 and the fluid-filled alveoli (two alveoli 15 shown in Figure 1) are in contact with the pleura 13, an artifact called a "B-line" appears in the ultrasound echo image, which is a long, linear artifact extending in the depth direction. This is thought to be because the ultrasound penetrates into the alveoli 15 where liquid has accumulated, and the ultrasound causes multiple reflections inside the alveoli 15 . In FIG. 1, on the ultrasound imaging plane 7 surrounded by the dashed line, a B-line 1 and a B-line 2 are generated in the depth direction from two alveoli 15 where liquid has accumulated.

[0034] The shear wave propagation velocity estimation device 10 of this embodiment utilizes B-lines generated from alveoli 15 where fluid accumulates, obtains information on shear waves propagating through lung parenchyma 14 from Doppler signals obtained on the B-lines, and measures the elasticity of the lung parenchyma using ultrasound.

[0035] The functions of the components 1 to 6 of the shear wave propagation velocity estimation device 10 are as follows. The oscillator 1 determines the frequency of the shear wave SW that is propagated inside the lung parenchyma 14 shown in FIG. The small vibration exciter 2 excites shear waves SW inside the lung parenchyma 14. When the small vibration exciter 2 is brought into contact with the body surface 11 and shear waves SW are generated from the small vibration exciter 2, the shear waves SW can be excited inside the lung parenchyma 14 as shown in FIG. The ultrasonic probe 3 outputs ultrasonic pulses to the lung parenchyma 14 and receives the resulting ultrasonic signals SU. The ultrasonic Doppler device 4 obtains an ultrasonic Doppler signal SUD on the B line from the ultrasonic signal SU received by the ultrasonic probe 3 . The velocity estimation device 5 estimates the phase of the shear wave SW from the ultrasonic Doppler signal SUD, and estimates the propagation velocity v of the shear wave SW propagating inside the lung parenchyma 14 from the phase of the shear wave SW. The output device 6 outputs the value of the propagation velocity v of the shear wave SW estimated by the velocity estimation unit 5.

[0036] The ultrasonic probe 3 generally has a width of about 5 to 7 cm. The small vibrator 2 and the ultrasonic probe 3 are configured as small devices. Furthermore, if the ultrasonic Doppler device 4, the velocity estimation device 5, and the output device 6 are each constructed as small devices, the overall device can be made smaller. Furthermore, by making the entire device smaller and configuring it so that it can be transported on casters or the like, it will be possible to transport the device to a hospital room or to load it into an ambulance vehicle.

[0037] In the shear wave propagation velocity estimation device 10 shown in FIG. 1, the ultrasonic Doppler device 4 and the velocity estimation device 5 are separate devices, but the ultrasonic Doppler device and the velocity estimation device may be configured as an integrated device.

[0038] Using the shear wave velocity estimation device 10 of this embodiment, the propagation velocity v of the shear wave SW can be measured, for example, by the method described below.

[0039] First, a continuous shear wave SW is propagated into the lung parenchyma 14 by a small vibrator 2 placed on the body surface 11 . At this time, at the position where the alveoli 15 where the liquid is stored are present, the depth Z of the alveoli 15 is A This produces an artifact called a B-line, which appears as continuous ultrasound waves reflecting from a deeper position than the normal. The reflected ultrasonic signal y(x,z,t) generated at this time can be modeled by the following equation (1).

[0040]

number

[0041] Here, y0(x,z) is the amplitude of the reflected ultrasonic signal on the B line, and θ(x,z,t) is the phase change of the ultrasonic wave due to the shear wave SW, and is given by the following equation (2).

[0042]

number

[0043] where f0 is the center frequency of the ultrasound and c is the speed of sound. L , t) is the displacement caused by the shear wave SW and is given by the following equation (3):

[0044]

number

[0045] In equation (3), ξ0(x,z L , t) is the amplitude of the shear wave SW, and φ b (x,z L ) is the phase of the shear wave SW, and ω b is the angular frequency of the shear wave SW. In equation (3), the amplitude and phase of the shear wave SW are determined by the mechanism by which the B-line appears. A It should be noted that the amplitude and phase at the position

[0046] Here, a schematic diagram of the reflected ultrasonic signal from an ultrasonic scan line where a B line appears is shown in Fig. 2A, and a schematic diagram of the reflected ultrasonic signal from an ultrasonic scan line where a B line does not appear is shown in Fig. 2B.

[0047] In the ultrasound scan line where the B line appears, multiple reflections occur due to the fluid-filled alveoli 15, as shown in Figure 2A. In contrast, in the ultrasound scan line where the B line does not appear, multiple reflections do not occur, as shown in Figure 2B.

[0048] The received ultrasonic waves are subjected to quadrature detection using a reference signal R(t) given by the following equation (4). This quadrature detection is performed by the ultrasonic Doppler device 4 of the shear wave propagation velocity estimation device 10 in FIG.

[0049]

number

[0050] In this case, the complex Doppler signals after quadrature detection are given by the following equations (5) to (7).

[0051]

number

[0052] From the quadrature detection signal (complex Doppler signal above), the phase of the shear wave, φ b (x,z L) can be calculated in two ways: The first method is to obtain the phase of the complex quadrature detection signal (complex Doppler signal) by Fourier analysis at the excitation frequency. However, the phase may also be obtained from the time when the phase becomes a predetermined phase value (for example, zero) without using Fourier analysis. The second method involves setting the frequency of the shear waves relative to the ultrasound repetition frequency to match the frequency conditions of the C-SWE method described above (n / 4 times the ultrasound repetition frequency (where n is an odd number greater than or equal to 1)), and then performing Fourier analysis in the frame direction on the color Doppler image or power Doppler image obtained at this time.

[0053] Either of the above two methods can be used, but the second method is based on the C-SWE method and can obtain the phase of shear waves without any modification to the ultrasound echo device, making it a method that is easily applicable to small, portable ultrasound diagnostic devices.

[0054] In an actual shear wave propagation velocity estimation device, there is a delay in the transmission timing of the ultrasonic waves depending on the position of the ultrasonic scanning line, and the phase φ of the shear wave obtained from the quadrature detection signal M (x,z L ) can be written as the sum of the phase difference due to the propagation of shear waves and the phase difference due to the delay in the ultrasonic transmission timing, as shown in the following equation (8).

[0055]

number

[0056] where φ US (x) is the phase difference due to the delay in the ultrasonic transmission timing. US (x) is determined by the ultrasonic transmission and reception method of the ultrasonic diagnostic equipment, and this value is known, so φ US Using (x), the phase due to the propagation of shear waves can be calculated using the following equation (9).

[0057]

number

[0058] Next, a method for determining the propagation velocity v of the shear wave SW from the phase of the shear wave SW when a plurality of B lines are formed will be described. FIG. 3 shows, as an example, a model of a lung echo in which two B-lines are formed. In this case, the phases of the shear waves SW measured on the B-line 1 and B-line 2 can be expressed by the following equations (10) and (11), respectively.

[0059]

number

[0060] Here, if the shear wave SW is propagating in the x-axis direction (azimuth direction) as shown by arrow 8 in FIG. 3, the propagation speed v of the shear wave SW can be estimated by the following equation (12).

[0061]

number

[0062] However, f b is the frequency of the shear wave SW. From this, the propagation velocity v of the shear wave SW between the two B lines can be estimated by the following equation (13).

[0063]

number

[0064] Here, x1 and x2 are the positions of B line 1 and B line 2 in the x-axis direction, respectively, as shown in FIG. The above equation (13) is the basic equation for obtaining the propagation velocity v of the shear wave SW in the pulmonary parenchyma in the present invention.

[0065] In reality, there may be three or more B lines. This method is also applicable to cases where lung disease has progressed and multiple B-lines are connected in the lateral direction (horizontal direction) to form a single wide B-line. If there is a phase difference of the shear wave SW in the single wide B line, the propagation velocity v of the shear wave SW can be obtained from the phase difference of the shear wave SW in the single wide B line. Specifically, the phase of the shear wave SW is measured at two different points on the wide line B. From this, the propagation velocity v of the shear wave SW can be obtained from the phase difference between the two points.

[0066] Next, signal processing in the shear wave propagation velocity estimation device 10 of this embodiment will be described. FIG. 4 shows a flowchart of signal processing in the shear wave propagation velocity estimation device 10 of this embodiment.

[0067] In the flowchart shown in FIG. 4, first, in step S11, vibration is applied from the surface of the living body (body surface) 11 using the small vibrator 2, and shear waves SW are propagated to the lung parenchyma 14.

[0068] Next, in step S12, the ultrasonic Doppler device 4 performs quadrature detection on the ultrasonic signal SU received from the ultrasonic probe 3 to obtain a complex quadrature detection signal (IQ signal).

[0069] Next, in step S13, the phase of the shear wave SW is obtained from the IQ signal. As a method for obtaining the phase of the shear wave SW, the following two methods can be used, as described above. 1. Method using the phase of complex quadrature detection signals 2. Method using the C-SWE method

[0070] Thereafter, in step S14, the propagation velocity v of the shear wave SW is measured using the above equation (13) from the phase of the shear wave SW estimated from the two B lines. In this way, the shear wave velocity estimation device 10 of this embodiment can obtain the propagation velocity v and estimate the stiffness of the lung parenchyma 14.

[0071] As described above, the shear wave propagation velocity estimation device 10 of this embodiment includes a small vibrator 2 that excites shear waves SW in the lung parenchyma 14 that is the measurement target, an ultrasound probe 3 that outputs ultrasound pulses to the lung parenchyma 14 and receives ultrasound signals SU within the lung parenchyma 14, an ultrasound Doppler device 4 that obtains ultrasound Doppler signals SUD from the received ultrasound signals SU, and a velocity estimation device 5 that estimates the propagation velocity of the shear waves SW from the ultrasound Doppler signals SUD.

[0072] As a result, according to the shear wave propagation velocity estimation device 10 of this embodiment, the propagation velocity of the shear wave SW in the lung parenchyma 14 can be estimated by the velocity estimation device 5 using the shear wave SW excited from the small vibrator 2 and the ultrasonic pulse output from the ultrasonic probe 3. Furthermore, since the shear wave propagation velocity estimation device 10 of this embodiment uses shear waves SW and ultrasonic pulses, it does not use radiation, so there is no need to set up a specialized examination room like a chest X-ray CT.Furthermore, since measurements can be made by operating the small vibrator 2 and the ultrasonic probe 3, measurements can be made easily.

[0073] In the shear wave propagation velocity estimation device 10 of this embodiment, the velocity estimation device 5 estimates the phase of the shear wave SW in multiple or one wide B-line from the ultrasonic Doppler signal SUD, and estimates the propagation velocity of the shear wave SW from the phase of the shear wave SW in multiple or one wide B-line. This makes it possible to estimate the elasticity of the lung parenchyma 14, which is the measurement target, from the propagation velocity of the shear wave SW estimated by the velocity estimation device 5, and therefore the elasticity of the lung parenchyma 14 can be measured with high accuracy. Therefore, the shear wave velocity estimation device 10 of this embodiment can measure the elasticity of the lung parenchyma 14 simply and accurately.

[0074] [experiment] In order to confirm the effect of the shear wave propagation velocity estimation device of the present invention, an experiment was carried out using a phantom that mimics the structure of the lungs.

[0075] (Experimental equipment) A schematic diagram of the phantom and an outline of the experimental apparatus used in the experiment are shown in Figure 5.

[0076] As shown in Figure 5, the phantom used in the experiment has a two-layer structure. Starting from the layer closest to the ultrasound probe 3, it consists of a konjac layer 31 approximately 10 mm thick, which simulates the soft tissue interposed between the body surface and lung tissue, and a coarse sponge layer 32 located below the konjac layer 31, which simulates the air-containing lungs. Furthermore, as shown in Figure 5, three tiny glass beads 33 with a diameter of approximately 2 mm were embedded at intervals of several mm into the sponge layer 32 just below the interface between the sponge layer 32 and the konjac layer 31 to simulate alveoli filled with liquid. Furthermore, a thin layer of ultrasonic gel was applied to the boundary between the konjac layer 31 and the sponge layer 32 so that ultrasonic waves could reach the glass beads 33 .

[0077] As the shear wave propagation velocity estimation device, a shear wave propagation velocity estimation device 10 shown in FIG. 1, which is equipped with a small vibrator 2, an ultrasonic probe 3, an ultrasonic Doppler device 4, etc., was used.

[0078] The small vibrator 2 was placed approximately 10 mm away from the ultrasonic probe 3 so that shear waves would propagate from the right side to the left side in FIG. The shear wave velocity estimation device 10 was set to a power Doppler mode, and the excitation frequency was set to 72.8 Hz, which was selected so that it was within the frequency range of the C-SWE method and so that the wavefront of the shear wave could be easily visualized on the power Doppler image.

[0079] The experiments were conducted using a normal sponge (hereinafter referred to as a normal sponge phantom) as the sponge layer 32 in the phantom, and a sponge that had been thoroughly sprayed with acrylic paint and hardened to simulate hardened lungs (hereinafter referred to as a hardened sponge phantom).

[0080] (Experimental results) B-mode images, which are ultrasound echo images of the two sponge phantoms, are shown in Figures 6A and 6B. Figure 6A shows the B-mode image of the normal sponge phantom, and Figure 6B shows the B-mode image of the hardened sponge phantom.

[0081] 6A and 6B, it can be seen that in both sponge phantoms, a B line, in which a high-intensity echo trails in the vertical direction (depth direction), is formed in the portion where the glass beads 33 are embedded.

[0082] The power Doppler image obtained when vibration was applied to a normal sponge phantom is shown in Figure 7. In Figure 7, an explanation has been added to the image to make it easier to understand the areas where the Doppler signal intensity is high.

[0083] Figure 7 shows that, over time, power Doppler images of increasing intensity are visualized in the areas formed in the sponge layer, B-lines 1 to 3. This is because the shear wave front propagates from right to left. The time intervals at which Doppler signals appear on each B-line were measured and found to be approximately every 0.45 seconds.

[0084] Figure 8 shows the frame-by-frame changes in the brightness values ​​of the power Doppler images appearing on each B-line when the normal sponge phantom is vibrated.

[0085] As can be seen from Figure 8, the wavefront of a shear wave appears as abrupt changes in brightness, and in Figure 8, the time transition of the brighter areas is shown by the dashed line. The slope of the dashed line in Figure 8 also indicates the elasticity of the sponge.

[0086] On the other hand, Figure 9 shows the power Doppler image obtained when vibration was applied to the hardened sponge phantom.

[0087] As can be seen in Figure 9, similar to the observations with the normal sponge phantom, high-intensity power Doppler images appear only on the B lines in the sponge layer. The time interval at which Doppler signals appear on each B line is approximately 0.3 seconds. This time interval is shorter than that of the normal sponge, and this result indicates that the elasticity of the cured sponge is higher than that of the normal sponge.

[0088] Figure 10 shows the frame-by-frame changes in the brightness values ​​of the power Doppler images appearing on each B-line when the hardened sponge phantom is vibrated. In Figure 10, as in Figure 8, the dashed lines show the results of connecting the high brightness areas.

[0089] Figure 10 shows that the slope of the dashed line is steeper than the results for the regular sponge (Figure 8). This indicates that the elasticity of the hardened sponge phantom is higher than that of the regular sponge phantom.

[0090] The propagation velocity v of the shear wave was estimated using equation (13) from the phase difference of the shear wave observed between B lines. Specifically, the propagation velocity v was estimated at three locations: between B line 1 and B line 2, between B line 2 and B line 3, and between B line 1 and B line 3, as shown in Figures 6A and 6B. The results of estimating the propagation velocity v in the normal sponge phantom and the hardened sponge phantom are shown in Figure 11. In Figure 11, the mean value and standard deviation of the estimated propagation velocity v are shown by error bars.

[0091] As can be seen from Figure 11, the propagation velocity of the hardened sponge phantom is 1.53 times faster than that of the normal sponge phantom. This difference in propagation velocity reflects the difference in elasticity between the two sponge phantoms. This result indicates that the shear wave propagation velocity estimation device of the present invention can be used to estimate the propagation velocity of shear waves propagating within a sponge. [Explanation of symbols]

[0092] 1 Oscillator, 2 Small vibrator, 3 Ultrasound probe, 4 Ultrasound Doppler device, 5 Velocity estimator, 6 Output device, 7 Ultrasound image plane, 10 Shear wave propagation velocity estimator, 11 Body surface, 12 Soft tissue, 13 Pleura, 14 Lung parenchyma, SW Shear wave, SU Ultrasound signal, SUD Ultrasound Doppler signal, v Shear wave propagation velocity

Claims

1. a vibrator for exciting shear waves in the measurement target; an ultrasonic probe that outputs ultrasonic pulses to the measurement object and receives ultrasonic signals from within the measurement object; an ultrasonic Doppler device for obtaining an ultrasonic Doppler signal from the ultrasonic signal received by the ultrasonic probe; a velocity estimation device that estimates a propagation velocity of the shear wave from the ultrasonic Doppler signal obtained by the ultrasonic Doppler device, The velocity estimation device estimates the phase of the shear wave in a plurality of B-lines or a single wide B-line from the ultrasonic Doppler signal, and estimates the propagation velocity of the shear wave from the phase of the shear wave in the plurality of B-lines or the single wide B-line. Shear wave propagation velocity estimation device.

2. The shear wave propagation velocity estimation device according to claim 1 , wherein the velocity estimation device estimates the propagation velocity of the shear wave by Fourier analysis of the phase of a complex quadrature detection signal obtained from the ultrasonic Doppler signal at the frequency of the shear wave.

3. The frequency of the shear wave is set to n / 4 times the repetition frequency of the ultrasonic wave (where n is an odd number equal to or greater than 1), 2. The shear wave propagation velocity estimation device according to claim 1, wherein the velocity estimation device estimates the phase of the shear wave in each of the B-lines by performing Fourier analysis in a frame direction on a color Doppler image or a power Doppler image obtained from the ultrasonic Doppler signal.

Images