CHARACTERIZATION OF CRUDE OIL BY SIMULATED DISTILLATION

20170363591 · 2017-12-21

    Inventors

    Cpc classification

    International classification

    Abstract

    A system and a method is provided for calculating and assigning one or more indicative properties (e.g., cetane number, pour point, cloud point and/or aniline point) of a fractions of a crude oil sample based upon gas chromatographic simulated distillation data.

    Claims

    1. A system for assigning an indicative property to a fraction of an oil sample based upon simulated distillation (SD) data, the system comprising: a non-volatile memory device that stores calculation modules and data, the data including SD data indicative of weight percent amount distilled across a range of boiling points of the oil sample; a processor coupled to the memory; and a calculation module that calculates and assigns the indicative property of a fraction of the oil sample as a function of the mid-boiling point or the weighted average boiling point determined from the SD data.

    2. A system for assigning an indicative property to a fraction of an oil sample comprising: a gas chromatograph that outputs simulated distillation (SD) data; a non-volatile memory device that stores calculation modules and data, the data including the SD data indicative of weight percent amount distilled across a range of boiling points of the oil sample; a processor coupled to the memory; and a calculation module that calculates and assigns the indicative property of a fraction of the oil sample as a function of the mid-boiling point or the weighted average boiling point determined from the SD data.

    3. The system of claim 1 wherein the calculation module calculates and assigns the indicative property of a traction of the oil sample as a function of the mid-boiling point or the weighted average boiling point of the oil sample and the density of the oil sample.

    4. A method for operating a computer to assign an indicative property to a fraction of an oil sample based upon simulated distillation (SD) data, the method comprising: entering into the computer SD data indicative of weight percent amount distilled across a range of boiling points of the oil sample; calculating and assigning the indicative property of a fraction of the oil sample as a function of the mid-boiling point or the weighted average boiling point determined from the SD data.

    5. A method for assigning an indicative property to a fraction of an oil sample comprising: operating a gas chromatograph to output simulated distillation (SD) data; entering into a computer the outputted SD data indicative of weight percent amount distilled across a range of boiling points of the oil sample; calculating and assigning the indicative property of a fraction of the oil sample as a function of the mid-boiling point or the weighted average boiling point determined from the SD data.

    6. The method of claim 4 wherein calculating and assigning the indicative property of a fraction of the oil sample is a function of the mid-boiling point or the weighted average boiling point of the oil sample and the density of the oil sample.

    7. The method of claim 4 wherein the oil sample is crude oil.

    8. The method of claim 4 wherein the oil sample is obtained from an oil well, stabilizer, extractor, or distillation tower.

    9. The method of claim 4 wherein the indicative property is a cetane number.

    10. The method of claim 4 wherein the indicative property is a pour point.

    11. The method of claim 4 wherein the indicative property is a cloud point.

    12. The method of claim 4 wherein the indicative property is an aniline point.

    13. The method of claim 4 wherein the indicative property is an octane number.

    14. The method of claim 4 wherein plural indicative properties are calculated including at least two indicative properties selected from the group consisting of cetane number, pour point, cloud point, aniline point and octane number.

    15. The method of claim 5, wherein the SD data is obtained from gas chromatography methods including ASTM D2887, ASTM D5236, ASTM D5399, ASTM D6352-04, ASTM D7213-05e1, ASTM D7398-07, ASTM D7169-05, ASTM D7096-10, ASTM D7500-10, and ASTM D5307-97.

    16. The method of claim 5, wherein the SD data is obtained from supercritical fluid chromatography methods.

    17. The method of claim 4, wherein the mid-boiling point is the mid-boiling point of whole crude oil.

    18. The method of claim 4, wherein the mid-boiling point of the crude oil is calculated at the 50 W % point of the SD data.

    19. The method of claim 4, wherein the mid-boiling point of the crude oil is calculated by taking the average of boiling points.

    20. The method of claim 4, wherein the mid-boiling point of the crude oil is calculated by taking the weighted average of boiling points.

    21. The method of claim 5, wherein correlative SD data is obtained from distillation methods selected from the group composed of ASTM D86, ASTM D1160, ASTM D2892, or any other methods based upon true boiling point distillation, supercritical fluid chromatography, and equilibrium flash.

    22. The system or method as in claim 21, wherein the true boiling point distillations were conducted in a column with the number of theoretical trays in the range 0-100.

    23. The system or method as in claim 21, wherein the true boiling point distillations were conducted in a column with the number of theoretical trays in the range 10-30.

    24. The system or method as in claim 21, wherein the true boiling point distillations were conducted in a column with the number of theoretical trays in the range 15-20.

    25. The system of claim 2 wherein the calculation module calculates and assigns the indicative property of a fraction of the oil sample as a function of the mid-boiling point or the weighted average boiling point of the oil sample and the density of the oil sample.

    26. The method of claim 5 wherein calculating and assigning the indicative property of a fraction of the oil sample is a function of the mid-boiling point or the weighted average boiling point of the oil sample and the density of the oil sample.

    Description

    BRIEF DESCRIPTION OF THE DRAWINGS

    [0020] Further advantages and features of the present invention will become apparent from the following detailed description of the invention when considered with reference to the accompanying drawings, in which:

    [0021] FIG. 1 is a graphic plot of simulated distillation data obtained from gas chromatography for three types of crude oil;

    [0022] FIG. 2 is a process flow diagram of steps carried out to establish a value for indicative properties of a gas oil fraction, using the system and method herein; and

    [0023] FIG. 3 is a block diagram of a component of a system for implementing the invention, according to one embodiment.

    DETAILED DESCRIPTION OF INVENTION

    [0024] A system and method is provided for determining one or more indicative properties of a hydrocarbon sample. Indicative properties (e.g., cetane number, pour point, cloud point and aniline point) of a gas oil fraction in crude oil samples are assigned as a function of the density and data obtained from simulated distillation data of a crude oil sample. The correlations also provide information about gas oil indicative properties without fractionation/distillation (crude oil assays) and will help producers, refiners, and marketers to benchmark the oil quality and, as a result, valuate the oils without performing the customary extensive and time-consuming crude oil assays. The currently used crude oil assay method is costly in terms of money and time. It costs about $50,000 US and takes two months to complete one assay. With the method and system herein, the crude oil can be classified as a function of SD data, and thus decisions can be made for purchasing and/or processing.

    [0025] The systems and methods are applicable for naturally occurring hydrocarbons derived from crude oils, bitumens, heavy oils, shale oils and from refinery process units including hydrotreating, hydroprocessing, fluid catalytic cracking, coking, and visbreaking or coal liquefaction.

    [0026] In the system and method herein, crude oil simulated distillation data is obtained by a suitable known or to be developed simulated distillation processes. Simulated distillation (SD) is a technique which separates individual hydrocarbon components in the order of their boiling points, and is used to simulate laboratory-scale physical distillation procedures. The separation can be accomplished with a gas chromatograph equipped with a chromatography column coated with a nonpolar (hydrocarbon-like) stationary phase, an oven and injector which can be temperature programmed. A flame ionization detector (FID) is used for detection and measurement of the hydrocarbon analytes. The SD analysis result provides a quantitative percent mass yield as a function of boiling point of the hydrocarbon components of the sample being analyzed. The chromatographic elution times of the hydrocarbon components are calibrated to the atmospheric equivalent boiling point (AEBP) of the individual n-alkane as described in a method from the ASTM by using n-alkane (n-paraffin) reference material. In the SD method ASTM D2887, the n-alkane calibration reference covers the boiling range 55-538° C. (100-1000° F.) which covers the n-alkanes with a chain length of about C5-C44.

    [0027] Alternative methods may be used, including ASTM D5236, ASTM D86, ASTM D5399, ASTM D6352-04, ASTM D7213-05e1, ASTM D7398-07, ASTM D7169-05, ASTM D7096-10, ASTM D7500-10, ASTM D5307-97, ASTM D1160, ASTM D2892, or any other methods based upon gas chromatography, true boiling point distillation, supercritical fluid chromatography, and equilibrium flash. There are well-known correlations to convert distillation numbers: for example, if the true boiling point distillation (ASTM D2892) is known, the SD data (ASTM D2887) can be estimated. While there are not yet ASTM methods using supercritical fluid chromatography, non-standardized lab methods are known.

    [0028] In the high temperature simulated distillation method (HTSD), the n-alkane calibration reference (a hydrogenated polyolefin wax, polywax 655) covers the boiling range 36-750° C. (97-1382° F.) which covers the n-alkanes with a chain length of about C5-C120. A key difference between ASTM D2887 and HTSD is the ability of HTSD to handle residue-containing samples (i.e. material boiling >538° C. or 1000° F.).

    [0029] SD and laboratory-scale physical distillation methods are routinely used for determining boiling ranges of petroleum crude oils and refined products. The boiling points with yield profile data of these materials are used by refinery engineers to make operational decisions to improve product yields and product quality. SD is valuable for, and can improve results from, computer modeling of refining processes for improvements in design and process optimization. Precise yield correlations between HTSD and crude assay distillation (a procedure which uses methods ASTM D5236 and D2892) have allowed HTSD to be successfully used in place of the physical distillation procedures. A typical simulated distillation chart obtained from a gas chromatogram of crude oil is shown in FIG. 1, where the W % of distilled fractions is plotted against the boiling temperature.

    [0030] The indicative properties (e.g., the cetane number, pour point, cloud point and aniline point) of the gas oil fraction, e.g., boiling in the range of 150-400° C. and in certain embodiments in the range of 180-370° C., can be assigned as a function of the density and the mid boiling point or weighted average boiling point of the gas oil or whole crude oil (T.sub.MBP), measured in Kelvin.

    [0031] FIG. 2 shows a process flowchart of steps in a method according to one embodiment herein that occur after gas chromatography is completed and the results are tabulated. In step 210, the cetane number is calculated. In step 220, the pour point is calculated. In step 230, the cloud point is calculated. In step 240, the aniline point is calculated. In step 250, the octane number is calculated. While FIG. 2 shows the steps performed sequentially, they can be performed in parallel or in any order. In certain embodiments, only one or more or steps 210, 220, 230, 240, 250 are carried out.

    [0032] That is,


    Indicative Property=f(density.sub.crude oil,T.sub.MBP gas oil)  (1a);


    Indicative Property=f(density.sub.crude oil,T.sub.MBP crude oil)  (1b);

    [0033] Equations (2) through (5) show that the cetane number, pour point, cloud point and aniline point can be assigned as a function of the density and simulated distillation of crude oils.


    Cetane Number (CET)=K.sub.CET+X1.sub.CET*DEN+X2.sub.CET*DEN.sup.2+X3.sub.CET*DEN.sup.3+X4.sub.CET*(T.sub.MBP/1000)+X5.sub.CET*(T.sub.MBP/1000).sup.2+X6.sub.CET*(T.sub.MBP/1000).sup.3+X7.sub.CET*DEN*(T.sub.MBP/1000)  (2);


    Pour Point (PP)=K.sub.PP+X1.sub.PP*DEN+X2.sub.PP*DEN.sup.2+X3.sub.PP*DEN.sup.3+X4.sub.PP*(T.sub.MBP/1000)+X5.sub.PP*(T.sub.MBP/1000).sup.2+X6.sub.PP*(T.sub.MBP/1000).sup.3+X7.sub.PP*DEN*(T.sub.MBP/1000)  (3);


    Cloud Point (CP)=K.sub.CP+X1.sub.CP*DEN+X2.sub.CP*DEN.sup.2+X3.sub.CP*DEN.sup.3+X4.sub.CP*(T.sub.MBP/1000)+X5.sub.CP*(T.sub.MBP/1000).sup.2+X6.sub.CP*(T.sub.MBP/1000).sup.3+X7.sub.CP*DEN*(T.sub.MBP/1000)  (4);


    Aniline Point (AP)=K.sub.AP+X1.sub.AP*DEN+X2.sub.AP*DEN.sup.2+X3.sub.AP*DEN.sup.3+X4.sub.AP*(T.sub.MBP/1000)+X5.sub.AP*(T.sub.MBP/1000).sup.2+X6.sub.AP*(T.sub.MBP/1000).sup.3+X7.sub.AP*DEN*(T.sub.MBP/1000)  (5);

    [0034] where:

    [0035] DEN=density of the crude oil sample;

    [0036] T.sub.MBP=mid boiling point of the gas oil or crude oil (derived from the simulated distillation curves of crude oils);

    [0037] and K.sub.CET, X1.sub.CET-X7.sub.CET, K.sub.PP, X1.sub.CP-X7.sub.CP, K.sub.CP, X1.sub.CP-X7.sub.CP, K.sub.AP, and X1.sub.AP-X7.sub.AP are constants.

    [0038] As an alternative to the mid boiling point a weighted average boiling point (WABP) can be used, shown in equation (6), below.

    [00001] WABP = ( T 10 * 10 ) + ( T 30 * 30 ) + ( T 50 * 50 ) + ( T 70 * 70 ) + ( T 90 * 90 ) 10 + 30 + 50 + 70 + 90 , ( 6 )

    [0039] where T.sub.10 is the boiling temperature of oil determined when 10 W % or V % of the fraction is recovered during the distillation, and where T.sub.30, T.sub.50, T.sub.70 and T.sub.90 are determined accordingly.

    [0040] An exemplary block diagram of a computer system 300 by which indicative property calculation modules can be implemented is shown in FIG. 3. Computer system 300 includes a processor 310, such as a central processing unit, an input/output interface 320 and support circuitry 330. In certain embodiments, where the computer 300 requires direct human interaction, a display 340 and an input device 350 such as a keyboard, mouse or pointer are also provided. The display 340, input device 350, processor 310, input/output interface 320 and support circuitry 330 are shown connected to a bus 360 which also connects to a memory unit 370. Memory 370 includes program storage memory 380 and data storage memory 390. Note that while computer 300 is depicted with the direct human interface components of display 340 and input device 350, programming of modules and importation and exportation of data can also be accomplished over the interface 320, for instance, where the computer 300 is connected to a network and the programming and display operations occur on another associated computer, or via a detachable input device, as are well known in the art for interfacing programmable logic controllers.

    [0041] Program storage memory 380 and data storage memory 390 can each comprise volatile (RAM) and non-volatile (ROM) memory units and can also comprise hard disk and backup storage capacity, and both program storage memory 380 and data storage memory 390 can be embodied in a single memory device or separated in plural memory devices. Program storage memory 380 stores software program modules and associated data, and in particular stores one or more indicative property calculation modules 381-385 such as cetane number calculation module 381, a pour point calculation module 382, a cloud point calculation module 383, an aniline point calculation module 384, and an octane number calculation module 385. Data storage memory 390 stores data used and/or generated by the one or more modules of the present system, including density of the crude oil sample, SD data or portions thereof used by the one or more modules of the present system, and calculated indicative properties generated by the one or more modules of the present system.

    [0042] The calculated and assigned results in accordance with the systems and methods herein are displayed, audibly outputted, printed, and/or stored to memory for use as described herein.

    [0043] It is to be appreciated that the computer system 300 can be any general or special purpose computer such as a personal computer, minicomputer, workstation, mainframe, a dedicated controller such as a programmable logic controller, or a combination thereof. While the computer system 300 is shown, for illustration purposes, as a single computer unit, the system can comprise a group/farm of computers which can be scaled depending on the processing load and database size, e.g., the total number of samples that are processed and results maintained on the system. The computer system 300 can serve as a common multi-tasking computer.

    [0044] The computing device 300 preferably supports an operating system, for example, stored in program storage memory 390 and executed by the processor 310 from volatile memory. According to the present system and method, the operating system contains instructions for interfacing the device 300 to the calculation module(s). According to an embodiment of the invention, the operating system contains instructions for interfacing computer system 300 to the Internet and/or to private networks.

    Example 1

    [0045] A set of constants K.sub.CET, X1.sub.CET-X7.sub.CET, K.sub.PP, X1.sub.PP-X7.sub.PP, K.sub.CP, X1.sub.CP-X7.sub.CP, K.sub.AP, and X1.sub.AP-X7.sub.AP was determined using linear regression. These constants are given in Table 3.

    TABLE-US-00003 TABLE 3 Cetane Number Pour Point Cloud Point Aniline Point Property (CET) (PP) (CP) (AP) K 544509.8 1344488.4 395024.0 24390.7 X1 −1932359.8 −4907366.2 −1429569.6 −49357.1 X2 2161099.3 5503008.0 1604628.0 52455.3 X3 −796440.7 −2031119.7 −592968.1 −18616.3 X4 142762.7 527938.4 136360.5 −41985.4 X5 −177339.2 −699945.0 −177392.5 65171.0 X6 90209.8 361176.8 91570.3 −33881.4 X7 −30458.6 −87436.2 −25137.2 408.8

    Example 2

    [0046] The following example is provided to demonstrate an application of equations (2) through (5). A sample of Arabian medium crude with a 15° C./4° C. density of 0.8828 Kg/l (e.g., at 15° C./4° C. using the method described in ASTM D4052) was analyzed by gas chromatography using the ASTM D2887 method. The simulated distillation data is shown in Table 4:

    TABLE-US-00004 TABLE 4 W % Temp. ° C. 0 1 2 37 3 68 4 83 5 94 6 100 7 113 8 121 9 127 10 138 11 144 12 151 13 157 14 165 15 172 16 175 17 185 18 191 19 196 20 204 21 210 22 216 23 222 24 229 25 235 26 241 27 249 28 255 29 261 30 267 31 272 32 279 33 285 34 290 35 297 36 303 37 308 38 315 39 319 40 326 41 331 42 337 43 342 44 348 45 354 46 360 47 366 48 372 49 378 50 384 51 390 52 396 53 402 54 409 55 415 56 422 57 428 58 434 59 440 60 446 61 452 62 458 63 465 64 471 65 478 66 485 67 492 68 499 69 506 70 513 71 520 72 528 73 535 74 543 75 551 76 559 77 567 78 575 79 583 80 592 81 599 82 608

    [0047] The mid boiling point of the crude oil is taken from the data at the 50 W % point, which is 384° C. (657 K). The calculations below are shown with temperature expressed in Kelvin, and using a density at 15° C./4° C. using the method described in ASTM D4052.

    [0048] Applying equation 2 and the constants from Table 3,

    [00002] Cetane .Math. .Math. Number .Math. .Math. ( CET ) = K CET + X .Math. .Math. 1 CET * DEN + X .Math. .Math. 2 CET * DEN 2 + X .Math. .Math. 3 CET * DEN 3 + X .Math. .Math. 4 CET * ( T MBP / 1000 ) + X .Math. .Math. 5 CET * ( T MBP / 1000 ) 2 + X .Math. .Math. 6 CET * T ( T MBP / 1000 ) 3 + X .Math. .Math. 7 CET * DEN * ( T MBP / 1000 ) = ( 544509.8 ) + ( - 1932359.8 ) .Math. ( 0.8828 ) + ( 2161099.3 ) .Math. ( 0.8828 ) 2 + ( - 796440.7 ) .Math. ( 0.8828 ) 3 + ( 142762.7 ) .Math. ( 657 / 1000 ) + ( - 177339.2 ) .Math. ( 657 / 1000 ) 2 + ( 90209.8 ) .Math. ( 657 / 1000 ) 3 + ( - 30458.6 ) .Math. ( 0.8828 ) .Math. ( 657 / 1000 ) = 59

    [0049] Applying equation 3 and the constants from Table 3,

    [00003] Pour .Math. .Math. Point .Math. .Math. ( PP ) = K PP + X .Math. .Math. 1 PP * DEN + X .Math. .Math. 2 PP * DEN 2 + X .Math. .Math. 3 PP * DEN 3 + X .Math. .Math. 4 PP * ( T MBP / 1000 ) + X .Math. .Math. 5 PP * ( T MBP / 1000 ) 2 + X .Math. .Math. 6 PP * ( T MBP / 1000 ) 3 + X .Math. .Math. 7 PP * DEN * ( T MBP / 1000 ) = ( 1344488.4 ) + ( - 4907366.2 ) .Math. ( 0.8828 ) + ( 5503008.06 ) .Math. ( 0.8828 ) 2 + ( - 2031119.7 ) .Math. ( 0.8828 ) 3 + ( 527938.4 ) .Math. ( 657 / 1000 ) + ( - 699945.0 ) .Math. ( 657 / 1000 ) 2 + ( 361176.8 ) .Math. ( 654 / 1000 ) 3 + ( - 87436.2 ) .Math. ( 0.8828 ) .Math. ( 657 / 1000 ) = - 10

    [0050] Applying equation 4 and the constants from Table 3,

    [00004] Cloud .Math. .Math. Point .Math. .Math. ( CP ) = K CP + X .Math. .Math. 1 CP * DEN + X .Math. .Math. 2 CP * DEN 2 + X .Math. .Math. 3 CP * DEN 3 + X .Math. .Math. 4 CP * ( T MBP / 1000 ) + X .Math. .Math. 5 CP * ( T MBP / 1000 ) 2 + X .Math. .Math. 6 CP * ( T MBP / 1000 ) 3 + X .Math. .Math. 7 CP * DEN * ( T MBP / 1000 ) = ( 395024.0 ) + ( - 1429569.6 ) .Math. ( 0.8828 ) + ( 1604628.0 ) .Math. ( 0.8828 ) 2 + ( - 592968.1 ) .Math. ( 0.8828 ) 3 + ( 136360.5 ) .Math. ( 657 / 1000 ) + ( - 177392.5 ) .Math. ( 657 / 1000 ) 2 + ( 91570.3 ) .Math. ( 657 / 1000 ) 3 + ( - 25137.2 ) .Math. ( 0.8828 ) .Math. ( 657 / 1000 ) = - 10

    [0051] Applying equation 5 and the constants from Table 3,

    [00005] Aniline .Math. .Math. Point .Math. .Math. ( AP ) = K AP + X .Math. .Math. 1 AP * DEN + X .Math. .Math. 2 AP * DEN 2 + X .Math. .Math. 3 AP * DEN 3 + X .Math. .Math. 4 AP * ( T MBP / 1000 ) + X .Math. .Math. 5 AP * ( T MBP / 1000 ) 2 + X .Math. .Math. 6 AP * ( T MBP / 1000 ) 3 + X .Math. .Math. 7 AP * DEN * ( T MBP / 1000 ) = ( 24390.7 ) + ( - 49357.1 ) .Math. ( 0.8828 ) + ( 52455.3 ) .Math. ( 0.8828 ) 2 + ( - 18616.3 ) .Math. ( 0.8828 ) 3 + ( - 41985.4 ) .Math. ( 657 / 1000 ) + ( 65171.0 ) .Math. ( 657 / 1000 ) 2 + ( - 33881.4 ) .Math. ( 657 / 1000 ) 3 + ( 408.8 ) .Math. ( 0.8828 ) .Math. ( 657 / 1000 ) = 66.

    [0052] Accordingly, as shown in the above example, indicative properties including cetane number, pour point, cloud point and aniline point can be assigned to the crude oil samples without fractionation/distillation (crude oil assays).

    Example 3

    [0053] An example calculation of WABP is presented below. When the sample is distilled, the boiling point of the sample is determined to be 149° C. when 10 W % of the sample is recovered, 230° C. when 30 W % of the sample is recovered, 282° C. when 50 W % of the sample is recovered, 325° C. when 70 W % of the sample is recovered, and 371° C. when 90 W % of the sample is recovered. Thus, T.sub.10 is 149° C., and the other values are similarly determined.

    [00006] WABT = .Math. [ 149 * 10 + 230 * 30 + 282 * 50 + 325 * 70 + 371 * 90 ] [ 10 + 30 + 50 + 70 + 90 ] = .Math. 315

    [0054] In alternate embodiments, the present invention can be implemented as a computer program product for use with a computerized computing system. Those skilled in the art will readily appreciate that programs defining the functions of the present invention can be written in any appropriate programming language and delivered to a computer in any form, including but not limited to: (a) information permanently stored on non-writeable storage media (e.g., read-only memory devices such as ROMs or CD-ROM disks); (b) information alterably stored on writeable storage media (e.g., floppy disks and hard drives); and/or (c) information conveyed to a computer through communication media, such as a local area network, a telephone network, or a public network such as the Internet. When carrying computer readable instructions that implement the present invention methods, such computer readable media represent alternate embodiments of the present invention.

    [0055] As generally illustrated herein, the system embodiments can incorporate a variety of computer readable media that comprise a computer usable medium having computer readable code means embodied therein. One skilled in the art will recognize that the software associated with the various processes described can be embodied in a wide variety of computer accessible media from which the software is loaded and activated. Pursuant to In re Beauregard, 35 USPQ2d 1383 (U.S. Pat. No. 5,710,578), the present invention contemplates and includes this type of computer readable media within the scope of the invention. In certain embodiments, pursuant to In re Nuijten, 500 F.3d 1346 (Fed. Cir. 2007) (U.S. patent application Ser. No. 09/211,928), the scope of the present claims is limited to computer readable media, wherein the media is both tangible and non-transitory.

    [0056] The system and method of the present invention have been described above and with reference to the attached figures; however, modifications will be apparent to those of ordinary skill in the art and the scope of protection for the invention is to be defined by the claims that follow.