The subject disclosure relates to the investigation of formation fluids. More particularly, the subject disclosure relates to apparatus and methods for identifying downhole the characteristics of a formation fluid such as, e.g., the fluid compressibility and fluid bulk density of the formation fluid.
It has long been desirable to characterize fluids in a geological formation. For example, in interpreting wellbore monitoring measurements and seismic surveys for phase saturation, the thermodynamic properties of the multicomponent reservoir fluid are required, because the acoustic velocity of the fluid is determined by both the density and the isentropic compressibility, and both velocity and density are needed to decipher seismic data. Thus, it may be inferred that density and isentropic compressibility are fundamental for seismic interpretation. Correspondingly, isothermal fluid compressibility is required for well-test interpretation. In formation testing, early transients are strongly influenced by the fluid within the tool, and interpretation of data requires tool-fluid compressibility. Despite the desirability of obtaining in situ measurements of compressibility of nearly incompressible fluids, such measurements have been generally unavailable. Rather, for measuring compressibility, practice is to bring reservoir fluid samples to the surface laboratory.
When reservoir fluid samples are brought to the surface, pressure and/or temperature changes during the transfer can lead to undesirable component separations and potentially irreversible alterations of the fluid. While gas evolution may be reversed, asphaltene separation from crude oil is generally not reversible within reasonable time-scales. As a result, the results of surface measurements on the fluid can have large uncertainties, even when the fluid is reconstituted.
This summary is provided to introduce a selection of concepts that are further described below in the detailed description. This summary is not intended to identify key or essential features of the claimed subject matter, nor is it intended to be used as an aid in limiting the scope of the claimed subject matter.
Methods and apparatus are provided for in situ measurement of compressibility of reservoir fluid. The apparatus and methods derive compressibility from volumetric changes imposed by mechanical strain of a piezoelectric material induced by an applied electric field.
In some embodiments, a downhole tool is provided with a chamber that is arranged to receive formation fluids and either contains or has a wall coupled to a piezoelectric material. A pressure sensor that measures the pressure in the chamber is provided. Fluid that is to be investigated is provided to the chamber, and a voltage is applied to the piezoelectric material in order to alter the shape of the piezoelectric material and thereby change the volume of the chamber. A change in fluid pressure is measured by the pressure sensor, and the change in pressure is used in conjunction with a known change in chamber volume to derive compressibility.
In some embodiments, the piezoelectric material is located completely inside the fluid chamber. When a voltage is applied to the piezoelectric material, the piezoelectric material contracts in directions orthogonal to poling directions, but expands along its vertical poling direction, thereby changing its volume. The change in bulk volume of the piezoelectric material thereby changes the remaining volume in the chamber for the fluid. The change in pressure resulting from the change in bulk volume (or the change in fluid volume) is then related to the compressibility of the fluid.
In some embodiments, where the change in the fluid volume in the chamber is impacted by change in the volume of the piezelectric material resulting from change in the dimension of the material in multiple directions, the relationship between the change in pressure and the fluid compressibility is expressed according to
where ΔP is the change in pressure, α is the volume ratio between the fluid and the piezoelectric material, d31 and d33 are the respective tensor components of piezoelectric moduli of the piezoelectric material in the directions orthogonal and parallel to the poling direction, respectively, E3 is a measure of the applied electrical field along the poling direction, and βT is the isothermal fluid compressibility.
In another embodiment, the piezoelectric material is located either outside the fluid chamber and is coupled to a wall of the fluid chamber, or a diaphragm is located in the chamber, and the piezoelectric material is separated from the fluid in the chamber by the diaphragm. When a voltage is applied to the piezoelectric material, the material moves the diaphragm or a wall of the chamber, thereby changing the volume of the chamber (or portion of the chamber containing the fluid) and hence the volume of the fluid. The change in pressure resulting from the change in in fluid volume is then be related to the compressibility of the fluid.
In some embodiments, where the fluid volume in the chamber is changed by movement of a chamber wall or a diaphragm resulting from application of a voltage to a piezoelectric material, the relationship between the change in pressure and the fluid compressibility can be expressed according to
where A is the cross-sectional area of the chamber, lp is the height (thickness) of the piezoelectric material prior to activation, and Vf is the fluid volume.
In some embodiments, a sound speed sensor is also provided. The sound speed sensor may include a sonic transmitter and receiver that are coupled to a wall or to walls of the chamber. Using a sound speed measurement made by the sound speed sensor and the compressibility as determined resulting from the volume change in the chamber and the resulting pressure change, other characteristics of the fluid such as the bulk density or the specific heat ratio may be inferred.
Additional aspects, embodiments, objects and advantages of the disclosed methods may be understood with reference to the following detailed description taken in conjunction with the provided drawings.
The particulars shown herein are by way of example and for purposes of illustrative discussion of the examples of the subject disclosure only and are presented in the cause of providing what is believed to be the most useful and readily understood description of the principles and conceptual aspects of the subject disclosure. In this regard, no attempt is made to show details in more detail than is necessary, the description taken with the drawings making apparent to those skilled in the art how the several forms of the subject disclosure may be embodied in practice. Furthermore, like reference numbers and designations in the various drawings indicate like elements.
Turning to
As shown in
In some embodiments, valve 64 is opened, the pump 67 is activated, and valve 62 is opened in order to allow fluid from the formation (or optionally fluid stored in chamber 22 or another chamber) to enter the chamber 50, typically at or close to downhole formation (ambient) pressure and be flushed via valve 66c that is opened to the borehole. When satisfactory collection of the native formation fluid is obtained, valves 66c and 62 may be shut. With the pressure sensor 82 sensing that a desired pressure has been obtained, pump 67 is shut and valve 64 is shut, thereby establishing a closed system in chamber 52. In situations, such as in
A voltage is then applied by voltage source 80 to the piezoelectric material 70 to alter the shape of the material 70 and to cause a change in the fluid-containing volume of the chamber as described in more detail hereinafter. An indication of the pressure (i.e., the pressure or pressure change) of the fluid-containing volume of the chamber is then measured, and the pressure change is used by the electronics and processing block 18 to find an indication of the compressibility of the fluid in the chamber as described in more detail hereinafter.
In the illustrated examples, the walls of the fluid chamber (and relevant portions of pipes/tubes 60) are substantially rigid, except to the extent of walls, membranes, or other surfaces that are selectably deformable for the purpose of controllably reducing the volume of the chamber 52 as described herein. It should be understood, however, that other examples may have other walls that may not be substantially rigid, e.g., where a wall is flexible or non-rigid in a manner that may be accounted for in calculating the compressibility.
In some embodiments, the chamber 52 is also provided with a sound speed (velocity) sensor 90, which may include a sonic wave transmitter 90a and a receiver 90b that may be coupled to one or more walls of the chamber. Using a sound speed measurement made by the sound speed sensor 90 and the compressibility as determined, characteristics of the fluid such as bulk density or specific heat ratio may be inferred by the processing block 18 as described in more detail hereinafter.
In some embodiments, once an experiment has been conducted down-hole to find the compressibility of a fluid sample, valve 66 may be opened to cause the fluid in the chamber to be jettisoned from the module or stored in a desired location, or valve 66c may be opened to eject the stored fluid into the borehole so as to be replaced by a second fluid. Thus, the module 50 may then be utilized to conduct additional experiments with a different fluid samples obtained at the same or a different location in the borehole.
In one aspect, isothermal compressibility βT is defined as
where V(P,T) is the molar volume and P and T are the fluid pressure and temperature, respectively. For a finite fluid volume change ΔV, and the corresponding change in fluid pressure ΔP, isothermal compressibility of the corresponding fluid may be calculated from
Since βT varies with P and T, according to some embodiments, ΔP/P<<1; e.g., ΔP/P is less than 1/100.
In one aspect, piezoelectric materials, when deformed, develop dipoles within the solid and cause a resulting charge accumulation on surface electrodes. The developed electric displacement or charge density is proportional to the imposed mechanical stress. Conversely, in an electric field, a piezoelectric material (e.g., crystal) experiences a volumetric strain. This strain is well characterized, and as set forth hereinafter, it is this volumetric strain that is utilized to change the volume of the fluid chamber 52.
Many practical piezo-materials are polycrystalline ceramics, examples of which include lead zirconate, lead titanate, and barium titanate. Although they are apparently piezoelectric, they are polarized electrostrictive, i.e., polycrystalline materials are subjected to a poling DC potential across the material so as to exhibit piezoelectricity through grain alignment. The material undergoes a semi-permanent dimensional increase in the poling direction. A dimensional decrease occurs in the orthogonal directions.
For post-poling deformations, if the applied voltage has the same polarity as the poling voltage, the material expands parallel to the applied electric field (as seen by comparing
In some embodiments, the tendency of the piezoelectric materials to expand or shrink is utilized for the purpose of changing the volume of the fluid chamber 52. In some embodiments, the ceramic material may be well-characterized for its deformation behavior. In particular, the piezoelectric material is strained by applying an electric voltage which results in an increase in pressure in the confined fluid. The change in pressure ΔP is easily calculated from the known converse piezoelectric effect through the constitutive relationship of the ceramic. In particular, ΔP is easily calculated from the volume of the fluid (Vf), the change in Vf, (ΔVf=−ΔVp when fluid surrounds the piezoelectric material), and the temperature, provided the consequential pressure induced volume strain on the piezo-electric material is shown to be small. Since any alteration in pressure resists deformation, in some embodiments, a perturbation series is considered.
Piezoceramic (as opposed to soft piezopolymers) materials generally have bulk moduli of between 10 and 100 GPa and density of between 7000 and 8000 kg m−3. The expected piezoelectric moduli for ceramics have magnitudes of about 4×10−10 m V−1. The strain transverse to the imposed field is nearly one-half and therefore the volumetric strain is rather small, about one-tenth of the longitudinal strain. For electric fields of about 400 kV m−1, it is reasonable to expect a pressure deflection of 780 Pa or more for an arrangement having a 1:100 volume ratio of piezoceramic to aqueous fluid. The sensitivity may be improved with smaller fluid chambers.
In general, piezoelectric materials can produce strains in the range of 0.01% to 0.2% for hard materials and 10 to 100% for soft materials, thus exhibiting a bulk volume change in enclosed media.
While soft piezoelectric materials would provide a higher strain than hard ones, ceramics show no detectable degradation and are resilient to cyclic loading. In one aspect, any of many ceramic materials may be utilized, including, by way of example and not by way of limitation, a number of naturally occurring materials such as quartz, tourmaline, sodium potassium tartarate and Rochelle salts, and a number of synthetic piezoceramic materials such as PZT (Lead Zirconium Titanate) and PT (Lead Titanate). Materials such as PZT (Pb(Zr, Ti)O3, PT (Pb TiO3) and PLZT (Pb La)(Zr Ti)O3) may be manufactured with properties such that their physical, chemical and piezoelectric characteristics may be adapted for desired purposes, i.e, suitable shapes, size, with choice of axes and orientation.
The relation between the applied electric field strength and the resulting strain in a piezoceramic material is given by
εj=Sijσj+dijEi, i,j=1,2,3 (5)
where i and j are the cartesian indices, dij are tensor components of the piezoelectric moduli, Ei is the applied electrical field, and σj and Sij are the stress and compliance of the material respectively. Piezoelectric moduli tensor component values for an example of both a soft and a hard piezoelectric material are given in Table 1.
Bulk moduli of piezo-ceramic materials are quite high, in excess of (100 GPa), making the material quite incompressible compared to fluids. Therefore, eq. (5) may be rewritten as
εj≈dijEi (6)
For a vertical poling direction (z-axis) on a piezoceramic disk, displacements due to piezoelectric effect are given by
ΔR=Rd31E3 (7)
Δlp=lpd33E3 (8)
where R and lp are respectively the radius and the thickness of the piezo-ceramic disk, and E3=Vdr/lp where Vdr is the drive voltage applied to the piezo-ceramic. It should be appreciated that d31 is a negative quantity (close to d33/2 in magnitude) and, therefore, represents a radial contraction, i.e., orthogonal to the poling direction. Thus, as suggested by comparing
Equation (9) can be reduced to
ΔVf=−[2d31+d33]E3Vp (10)
where the volume of piezoceramic disk Vp=πR2lp.
In one aspect, the characteristic numbers for a specific piezoceramic crystal (ceramic APC 855(Navy VI)) subject to an electric field of 40 kV m−1 may be derived. Consider a 1 cm3 ceramic APC 855 crystal. For lp=1 cm, R=1/π0.5 cm. The piezoelectric moduli for the ceramic APC 855 are given as d33=630×10−12 m V−1, and d31=−276×10−12 m V−1. The APC 855 piezo-ceramic disk undergoes an extension of h=2.52×10−4 mm and a contraction of R=−6.23×10−5 mm under the specified 40 kV m−1 electrical field. A corresponding volumetric change in the disk (which is the same magnitude as the displaced fluid) is about 3.12×10−3 mm3. In some embodiments, it is possible to amplify this volumetric change, by having the ceramic deflect a metallic membrane or wall, as described hereinafter. Also, since the change in fluid volume ΔVf is proportional to the electrical field applied, increasing the electrical field will result in a larger volume change.
The difference in fluid pressure corresponding to fluid volume displaced in closed chamber can be calculated by using the equation
where βT is the isothermal compressibility of the fluid. For in situ conditions, the temperature change ΔT resulting from the converse piezoelectric effect is negligible. Furthermore, since this measurement can be made very rapidly and the response is essentially instantaneous, setting Vp=αVf for αϵR and substituting equation (10) into equation (11), the expression for change in pressure becomes
Thus, a volume ratio of 100:1 (α= 1/100) between the fluid and the piezo-ceramic results in a fluid volume strain of −3.12×10−8. For water, this results in a pressure increase of about 78 Pa under isothermal conditions. A more than one hundred times increase in ΔP is possible by increasing the volume ratio and the electric field (to, e.g., a maximum of about 1 kV mm−1), and amplifying the displacement by utilizing the expansion of the disk in the poling direction only, e.g., deflecting a metallic membrane through the strain of the disk in poling direction. Hence, a 78-7800 Pa pressure increase through converse piezoelectric effect is attainable, a quantity which is certainly within pressure measurement capability. For example, it is known that Crystal Quartz Gauge (CQG) sensors for downhole pressure measurements yield a maximum 6.89 kPa+0.01% of reading error in static measurements and have resolutions better than 20.7 Pa. (See, R. J. Besson, et al., “A dual-mode thickness-shear quartz pressure sensor”, IEEE Transactions on Ultrasonic, Ferroelectrics and Frequency Control, 40(5), 1993, and N. Matsumoto, et al., “Long-term stability and performance characteristics of crystal quartz gauge at high pressures and temperatures”, IEEE Transactions on Ultrasonic, Ferroelectrics and Frequency Control, 47(2), 2000)). Moreover, the 0.01% of reading error is associated with the uncertainty in static pressure and will not contribute to a pressure difference measurement. With a fluid whose compressibility is about a factor ten more than water, pressure differences are lowered by a factor of ten, a quantity whose low end is just measurable at pressures of interest, but whose high end is easily measurable. Table 2 shows attainable pressure increases in water and oil through the converse piezoelectric effect for a range of ceramic-to-fluid volume ratios (e.g., 1%, 2%, 3%, 4%, 5%, 10% and 15%) and different electrical fields (e.g., 40 Vmm−1, 50 Vmm−1, 100 Vmm−1 and 1000 Vmm−1), where compressibility βT was set to equal 4.0×10−10 Pa−1 for water and 1.0×10−9 Pa−1 for oil.
According to some embodiments, the compressibility of a fluid may be determined by providing the fluid to a chamber containing a piezoelectric material, applying a voltage to the piezoelectric material in order to alter the shape of the piezoelectric material, measuring a change in fluid pressure in the chamber, and using a processor (e.g., processing unit 18) to calculate the compressibility of the fluid according to equation (12), where the volume ratio of piezoelectric material to fluid (α), the piezoelectric modulii d31 and d33, and the drive voltage parameter E3=Vdr/lp are known. In some embodiments, the compressibility of the fluid is determined with the fluid being analyzed downhole. In some embodiments, the compressibility of the fluid is determined using an apparatus such as described above with respect to
As previously mentioned, net fluid volume displaced by the ceramic piezoelectric element due to the piezoelectric effect is affected by the shrinkage of the element in the direction orthogonal to the poling direction. This is especially true when |d31|≈d33/2, which would make the net fluid volume displaced almost zero, as is evident from equation (12).
According to one aspect, the change in net fluid volume ΔVf may be amplified by reconfiguring the system such that only the expansion of the piezoelectric element (in the poling direction) is utilized for the compression of the fluid and any contribution from d31 is avoided. For instance, rather than submerging the piezoelectric element 70 fully into the fluid F as in
membrane M aligned in the same plane. Thus, the strain due to d33 on the ceramic element is fully utilized to displace the metallic membrane M in
According to one aspect, fluid displacement may be achieved either through deflecting a membrane or pushing a plate (not shown) or a wall of the chamber to compress the fluid while making sure that fluid mass inside chamber is conserved. For example, if it is assumed that a plate or membrane sits between fluid and one or more piezoceramic disk, then ignoring the bending of the membrane or plate (this can accounted for in detailed mechanical calculation) the volume change in fluid and the corresponding pressure change will be
ΔVf=d33E3Alp (13)
where A is the cross sectional area of the chamber and lp is the thickness, or height, of the crystal. Therefore,
In this case, a substantially amplified pressure increase is obtained compared to the case of equation (12) when both d33 and d31 affect the resulting pressure change.
According to some embodiments, the compressibility of a fluid may be determined by providing the fluid to a chamber with a membrane or other element that is coupled to a piezoelectric material, applying a voltage to the piezoelectric material in order to alter the shape of the piezoelectric material, measuring a change in fluid pressure in the chamber, and using a processor to calculate the compressibility of the fluid according to equation (14), where the initial fluid volume of the chamber (Vf), the cross sectional area of the chamber (A), the piezoelectric material thickness or height (lp), the piezoelectric modulus (d33), and the drive voltage parameter E3=Vdr/lp are all known. In some embodiments, the compressibility of the fluid is determined with the fluid being analyzed downhole. In some embodiments, the compressibility of the fluid is determined using an apparatus such as described above with respect to
According to one aspect, the dimensional (and volume) changes in a piezo-ceramic element is sensitive to the make of the piezo-ceramic material. Accurate values for piezoelectric moduli d31 and d33, are useful for determining βT. While current meters have a range of 1 to 2000 pC N−1, and the error in d33 measurements has been reported to be approximately 2% for 100-2000 pC N−1, piezo-ceramic elements can be calibrated for their true piezoelectric coefficients using reference fluid measurements, and any error in compressibility contributed by d33 or d31 can be eliminated. In addition, the physical setup of the piezoelectric element and the fluid chamber may be calibrated for the variations in piezoelectric coefficients with respect to pressure and temperature in order to maintain measurement accuracy.
In some embodiments, once the compressibility of the fluid is determined (e.g., using equation (12) or (14)), other characteristics of the fluid may be determined using other measurements in conjunction with the determination of compressibility. Thus, by way of example only, the bulk density of the fluid in the chamber may be inferred from a determination of the speed of sound inside the chamber (as determined by a sound velocity sensor) and the compressibility. More particularly,
ν=√{square root over (∂P/∂ρ)}=√{square root over ((ρβs)−1)}=√{square root over (ρ−1Ks)} (15)
where the derivative is at constant entropy, and ν is the measured speed of sound, P is the pressure, ρ is the density of the fluid, βs is the adiabatic compressibility at constant entropy, and Ks is the bulk modulus. βs is the product of βT and the constant pressure to constant volume specific heat ratio commonly denoted as γ. If the density of the fluid is measured independently through any of several known mechanism, e.g a vibrating tube, βs can be inferred if the acoustic velocity is known. Thus, from the measured βT, γ may be determined from
γ=βT/βs. (16)
In one aspect, some of methods and processes described above are performed by a processor, such as determining compressibility according to equation (12) or equation (14), or determining bulk density according to equation (15), or specific heat ratio according to equation (16). The term “processor” should not be construed to limit the embodiments disclosed herein to any particular device type or system. The processor may include a computer system. The computer system may also include a computer processor (e.g., a microprocessor, microcontroller, digital signal processor, or general purpose computer) for executing any of the methods and processes described above. The computer system may further include a memory such as a semiconductor memory device (e.g., a RAM, ROM, PROM, EEPROM, or Flash-Programmable RAM), a magnetic memory device (e.g., a diskette or fixed disk), an optical memory device (e.g., a CD-ROM), a PC card (e.g., PCMCIA card), or other memory device.
The methods and processes described above may be implemented as computer program logic for use with the computer processor. The computer program logic may be embodied in various forms, including a source code form or a computer executable form. Source code may include a series of computer program instructions in a variety of programming languages (e.g., an object code, an assembly language, or a high-level language such as C, C++, or JAVA). Such computer instructions can be stored in a non-transitory computer readable medium (e.g., memory) and executed by the computer processor. The computer instructions may be distributed in any form as a removable storage medium with accompanying printed or electronic documentation (e.g., shrink wrapped software), preloaded with a computer system (e.g., on system ROM or fixed disk), or distributed from a server or electronic bulletin board over a communication system (e.g., the Internet or World Wide Web).
Alternatively or additionally, the processor may include discrete electronic components coupled to a printed circuit board, integrated circuitry (e.g., Application Specific Integrated Circuits (ASIC)), and/or programmable logic devices (e.g., a Field Programmable Gate Arrays (FPGA)). Any of the methods and processes described above can be implemented using such logic devices.
Although only a few examples have been described in detail above, those skilled in the art will readily appreciate that many modifications are possible in the examples without materially departing from this subject disclosure. Thus, by way of example only, and not by way of limitation, while various embodiments describe determining compressibility by decreasing the fluid chamber volume through the use of a piezoelectric element, it is possible to increase the fluid chamber volume using a piezoelectric element and to determine compressibility of the fluid within the fluid chamber. Also, particular arrangements of chambers and piezoelectric materials have been shown, such as a piezoelectric material located inside a container with the fluid, and a piezoelectric material located inside a container but separated from the fluid by a membrane or diaphragm, it will be appreciated that other arrangements could be provided as long as the change in one or more dimensions of the piezoelectric material causes a change in fluid volume inside the chamber. Further, while particular piezoelectric materials were described, and while particular applied voltages were described, it will be appreciated that other materials and other voltages could be used. Accordingly, all such modifications are intended to be included within the scope of this disclosure as defined in the following claims. In the claims, means-plus-function clauses are intended to cover the structures described herein as performing the recited function and not only structural equivalents, but also equivalent structures. It is the express intention of the applicant not to invoke 35 U.S.C. § 112, paragraph 6 for any limitations of any of the claims herein, except for those in which the claim expressly uses the words ‘means for’ together with an associated function. As used in the description and claims, phrases in the form of “at least one of [a] and [b]” should be construed as being disjunctive, i.e., encompassing arrangements that include [a] but not [b], arrangements that include [b] but not [a], and arrangements that include [a] and [b].
Number | Name | Date | Kind |
---|---|---|---|
4994671 | Safinya et al. | Feb 1991 | A |
5167149 | Mullins et al. | Dec 1992 | A |
5201220 | Mullins et al. | Apr 1993 | A |
5266800 | Mullins | Nov 1993 | A |
5331156 | Hines et al. | Jul 1994 | A |
5859430 | Mullins et al. | Jan 1999 | A |
6378364 | Pelletier et al. | Apr 2002 | B1 |
6474152 | Mullins et al. | Nov 2002 | B1 |
6640625 | Goodwin | Nov 2003 | B1 |
7913556 | Hsu et al. | Mar 2011 | B2 |
7921691 | DiFoggio et al. | Apr 2011 | B2 |
8321133 | Hsu et al. | Nov 2012 | B2 |
8612154 | Hsu | Dec 2013 | B2 |
20020113717 | Tang | Aug 2002 | A1 |
20120037423 | Geerits | Feb 2012 | A1 |
Entry |
---|
Besson, R. J. et al., “A Dual-Mode Thickness-Shear Quartz Pressure Sensor”, IEEE Transactions on Ultrasonic, Ferroelectrics and Frequency Control, 1993, 40(5), 8 pages. |
Cain, M. G. et al., “Degradation of Piezoelectric Materials”, National Physical Laboratory Management Ltd., 1999, 41 pages. |
Liang, L. et al., “Pressure and Electric Field Effects of Piezoelectric Responses of KNb03” Journal of Applied Physics, 2012, 112, 6 pages. |
Matsumoto, N. et al., “Long-Term Stability and Performance Characteristics of Crystal Quartz Guage at High Pressures and Temperatures”, IEEE Transactions on Ultrasonic, Ferroelectrics and Frequency Control, 2000, 47(2), 9 pages. |
Nishi, R. Y. et al., “Behavior of Plezoceramic Projector Materials Under Hydrostatic Pressure”, Journal of the Acoustical Society of America, 1964 36(7), pp. 1292-1296. |
Number | Date | Country | |
---|---|---|---|
20170074096 A1 | Mar 2017 | US |