Research article Topical Sections

An economic comparison of dedicated crops vs agricultural residues as feedstock for biogas of vehicle fuel quality

  • The vast majority of the biofuels presently used in the EU are so called first generation biofuels produced from crops. Concerns of food security, displacement of food crop production and indirect land use change (iLUC) has led to the introduction of measures to reduce the use of first generations biofuels and promote so called advanced biofuels based on feedstock that does not compete with food/feed crops, such as waste and agricultural residues. In Sweden, 60% of the biofuel consumption is already based on waste/residual feedstock, and a unique feature of the Swedish biofuel supply is the relatively large use of biogas for transport, representing 9% of the current use of biofuels. The use of waste/residues dominates the biogas production, but agricultural residues, representing a large domestic feedstock potential, are barely used at present. This could indicate that biofuels from such feedstock is non-competitive compared both to fossil fuels and to biofuels produced from crops and waste under existing policy framework. This study show that without subsidies, the production cost of biogas as biofuel from all non-food feedstocks investigated (grass, crop residues and manure) is higher than from food crops. A shift from food crops to residues, as desired according to EU directives, would thus require additional policy instruments favoring advanced biofuel feedstock. Investment or production subsidies must however be substantial in order for biogas from residues to be competitive with biogas from crops.

    Citation: Mikael Lantz, Emma Kreuger, Lovisa Björnsson. An economic comparison of dedicated crops vs agricultural residues as feedstock for biogas of vehicle fuel quality[J]. AIMS Energy, 2017, 5(5): 838-863. doi: 10.3934/energy.2017.5.838

    Related Papers:

    [1] Lixiang Zhang, Yian Zhu, Jie Ren, Wei Lu, Ye Yao . A method for detecting abnormal behavior of ships based on multi-dimensional density distance and an abnormal isolation mechanism. Mathematical Biosciences and Engineering, 2023, 20(8): 13921-13946. doi: 10.3934/mbe.2023620
    [2] Hongwei Sun, Qian Gao, Guiming Zhu, Chunlei Han, Haosen Yan, Tong Wang . Identification of influential observations in high-dimensional survival data through robust penalized Cox regression based on trimming. Mathematical Biosciences and Engineering, 2023, 20(3): 5352-5378. doi: 10.3934/mbe.2023248
    [3] Jiongxuan Zheng, Joseph D. Skufca, Erik M. Bollt . Heart rate variability as determinism with jump stochastic parameters. Mathematical Biosciences and Engineering, 2013, 10(4): 1253-1264. doi: 10.3934/mbe.2013.10.1253
    [4] Ping Zhou, Yuting Zhang, Yunlei Yu, Weijia Cai, Guangquan Zhou . 3D shape measurement based on structured light field imaging. Mathematical Biosciences and Engineering, 2020, 17(1): 654-668. doi: 10.3934/mbe.2020034
    [5] Andrew J. Whittle, Suzanne Lenhart, Louis J. Gross . Optimal control for management of an invasive plant species. Mathematical Biosciences and Engineering, 2007, 4(1): 101-112. doi: 10.3934/mbe.2007.4.101
    [6] Tian Ma, Boyang Meng, Jiayi Yang, Nana Gou, Weilu Shi . A half jaw panoramic stitching method of intraoral endoscopy images based on dental arch arrangement. Mathematical Biosciences and Engineering, 2024, 21(1): 494-522. doi: 10.3934/mbe.2024022
    [7] Sukun Tian, Ning Dai, Linlin Li, Weiwei Li, Yuchun Sun, Xiaosheng Cheng . Three-dimensional mandibular motion trajectory-tracking system based on BP neural network. Mathematical Biosciences and Engineering, 2020, 17(5): 5709-5726. doi: 10.3934/mbe.2020307
    [8] Xin Wei, Jianjun Paul Tian, Jiantao Zhao . Fairy circles and temporal periodic patterns in the delayed plant-sulfide feedback model. Mathematical Biosciences and Engineering, 2024, 21(8): 6783-6806. doi: 10.3934/mbe.2024297
    [9] Tianran Yuan, Hongsheng Zhang, Hao Liu, Juan Du, Huiming Yu, Yimin Wang, Yabin Xu . Watertight 2-manifold 3D bone surface model reconstruction from CT images based on visual hyper-spherical mapping. Mathematical Biosciences and Engineering, 2021, 18(2): 1280-1313. doi: 10.3934/mbe.2021068
    [10] Junhua Gu, Xu Deng, Ningjing Zhang, Suqi Zhang . Recommendation model based on generative adversarial network and social reconstruction. Mathematical Biosciences and Engineering, 2023, 20(6): 9670-9692. doi: 10.3934/mbe.2023424
  • The vast majority of the biofuels presently used in the EU are so called first generation biofuels produced from crops. Concerns of food security, displacement of food crop production and indirect land use change (iLUC) has led to the introduction of measures to reduce the use of first generations biofuels and promote so called advanced biofuels based on feedstock that does not compete with food/feed crops, such as waste and agricultural residues. In Sweden, 60% of the biofuel consumption is already based on waste/residual feedstock, and a unique feature of the Swedish biofuel supply is the relatively large use of biogas for transport, representing 9% of the current use of biofuels. The use of waste/residues dominates the biogas production, but agricultural residues, representing a large domestic feedstock potential, are barely used at present. This could indicate that biofuels from such feedstock is non-competitive compared both to fossil fuels and to biofuels produced from crops and waste under existing policy framework. This study show that without subsidies, the production cost of biogas as biofuel from all non-food feedstocks investigated (grass, crop residues and manure) is higher than from food crops. A shift from food crops to residues, as desired according to EU directives, would thus require additional policy instruments favoring advanced biofuel feedstock. Investment or production subsidies must however be substantial in order for biogas from residues to be competitive with biogas from crops.


    The three-dimensional reconstruction technology of micro-nano surface topography has been widely studied and applied in many fields such as industrial manufacturing, computer vision reverse engineering [1]. Especially, in the last few decades, optical methods such as coherence scanning interferometry (CSI), phase shifting interferometry (PSI) and focus variation microscopy (FVM) have become increasingly popular for three-dimensional surface reconstruction [2]. During the process of reconstruction, the point cloud contains noise by factors such as the brightness of the ambient light, measuring equipment, object materials and surface texture and man-made factors which will affect the accuracy of reconstruction [3]. Therefore, it is urgent to research an effective point cloud denoising strategy for three-dimensional reconstruction of microscopic measurement as a preprocessing step. However, automatic and effective removal of measurement outliers is challenging since the object surface model is unavailable and the estimate of the object surface shape would be inaccurate in the presence of extensive outliers.

    Traditional point cloud denoising algorithms are difficult to filter out the outliers of the point cloud reconstructed during the microscopic measurement process. There two significant reasons that some light during the process of collecting pictures for microscopic measurement and the reflective properties of the material causing different kinds of outliers, such as non-isolated and isolated outliers. Figure 1 shows the point cloud reconstructed by microscopic measurement with isolated outliers and non-isolated outliers.

    Figure 1.  Copperring with different outliers.

    Because the noise generated during the microscopic measurement is relatively random, and mixed with non-isolated outliers, point cloud filtering methods based on the new method is far from enough. This paper compares the outlier model of the point cloud with the model of the social circle, and studies a robust outlier filtering algorithm based on the social circle. The work has important innovation value and application significance for filtering outliers of ordered point cloud reconstructed by microscopy and improving the accuracy of reconstruction.

    In recent years, a large number of scholars have done in-depth research on point cloud denoising algorithms for improving the accuracy and stability of reconstruction. These filtering algorithms are mainly used to filter outliers for 3D laser scan which is also suitable for point cloud filtering of microscopic reconstruction. For example, since outliers reconstructed by 3D laser scan are usually randomly distributed, Sotoodeh [3] proposed a density-based removal method, which judges whether the point set is a noise point based on the mean value of the closest distance in the K field of the point. However, it is difficult to identify dense outliers depending on local information, so the filtering effect of dense outliers is very poor. For some isolated outliers, Yu [4] performed point cloud denoising through a three-dimensional grid which used a measure of local entropy to update mesh vertices, but this method cannot guarantee convergence, what’s more, it’s not suitable for non-isolated outliers. Some scholars proposed a point cloud outlier filtering algorithm based on statistical distribution [5] which can filtered out outliers quickly, this method performs well in isolated outlier removal, but fails in the non-isolated outlier removal. Wang et al. [6] used the majority voting method to remove non-isolated outliers, but it is not suitable for large-scale point clouds and the robustness is not good. In addition, many filtering algorithms [7,8,9] in image processing can also be applied to filter outliers in point clouds, such as mean filtering [10] and bilateral filtering [11], for example Fleishman et al. [12] applied a bilateral filtering algorithm based on image denoising to the point cloud filtering process, adjusting the position of the vertex according to its normal direction to achieve denoising processing. However, the outliers are not identified, and the original data information is modified while removing the outliers.

    Reviewing the development described in afore mentioned studies, it has poor effect on non-isolated outliers which comes from the process of microscopic measurement. Therefore the aim of this work is to develop a fast and effective method to remove the outlier of organized point cloud without modification of non-outlier, besides, the method not only filter out isolate outliers, but also filter out non-isolated outliers. Inspired by social circle [13], an effective solution is proposed removing outlier of organized point cloud. A point is defined as outlier if it lies in exterior of the social circle of all other points [14]. Firstly, a non-outlier needs to be set as the initial point and determine the radius of social circle. Then, the social circle is expanded to iteratively start with the initial point. At last, we save the point which lies in the interior of social circle as the final result.

    This study is divided into the following sections. In section 3, we not only begin with the introduction of Focus Variation Technology, but also introduce the algorithm of social circle filtering. Next, in section 4, experiments are carried out to validation of the proposed algorithm. Conclusions are finally addressed in section 5.

    Focus Variation is also called depth from focus, originally derived from the shape from focus proposed by Nayar [15]. This technology mainly utilizes the small depth of field of the microscope and precision vertical movement equipment. The core idea is to move the object to be measured relative to the microscope, at the same time take a set of sequence pictures, and record height information [16].

    As is shown in Figure 2(a), Z axis represents depth information X and Y axis represent pixel. A set of images are captured while the objective of the microscope moves upward (or downward) relative to the surface being measured. As is shown in Figure 2(b), we get sequence of FM (Foucs Measure) images after determining the sharpness of all pixels in sequence of images.

    Figure 2.  The process of FVM. (a) Images acquisition; (b) Focus measure; (c) Depth estimation; (d) Organized point cloud.

    The commonly used focus evaluation algorithm is SML which calculated as follows:

    Ml=|2(I(x,y)I(xs,y)I(x+s,y)|+|2(I(x,y)I(x,ys)I(x,y+s)| (1)
    F(i,j)=i+Nx=iNj+Ny=jNMl(x,y) (2)
    Ml(x,y)=T (3)

    s is the step size, N is the calculation window size, and T is the threshold.

    At last, a set of depth information is generated at pixel P(x, y) in sequence of images which is shown in Figure 2(c), we fit the set of depth information. The depth information of measured surface is determined after calculating all pixels in sequence of images. Finally, each pixel is fused according to the position of the corresponding image to obtain an organized point cloud.

    As we know, there are some outliers in the reconstructed point cloud due to some environmental factors, the outlier filtering method proposed in this paper is based on social circle filter. The social circle theory originated from social networks. Users in social networks have a large number of friends, organizing these friends through certain rules is our social circle. Usually, the user's social circle can be described as a self-network clustering problem. Everyone has their own social connections in the society [17]. People are divided into different groups or organizations, within which a man can easily get in touch with the other through friends or friends of friends. However, it is hard for a person to make contact with someone who does not belong to the group. When a person cannot get in touch with people in the circle, then the person can be considered as not belonging to the community, the group will consider him as a stranger. What’s more, Wang [18] proposed an improved density peaks-based clustering method which incorporates both two information of users for social circle discovery in social networks. An illustration of social circle is shown in Figure 3, in this example blue dot is the centeral point of the social circle, blue circle is the initial social circle, black circle represents the expanded social circle, and the red dot represents elements that do not belong to this social circle.

    Figure 3.  An illustration of social circle.

    Inspired by which, the concept of the social circle is also can be applied to the outlier removal. The main task of social circle filtering is to define a community and relationship between community members. Through this relationship, we will expand our social circle until the social circle includes all community members, and the members outside the community are outliers point. Our method includes two steps. The first step is to determine the initial point and social circle radius. The second step is to expand the social circle iteratively. The method is detailed as follows.

    The initial point is very important as it is the beginning point of social circle expansion. Here, the I is taken as a symbol of the initial point. If I is an outlier, the final result will be error. The final point set will lie within the social circle of outlier, which is contrary to our goal. Thus, it is important to determine a reliable non-outlier. In this paper, we find the initial point I through a “vote” strategy.

    A square area needs to be determined as the vote area. Here, users need to select a square area T(xi, yi, zi) (I = 0, 1, 2, …)from the original organized point cloud P(xi, yi, zi) (I = 0, 1, 2, …) which comes from FVM. Usually the square area T should not be less than 60 points in the x and y directions when selected. Surely, the larger the selection area, the more accurate the initial point, and try to select the smoother area in the original point cloud P. For each point p of T, it is necessary to calculate the distance between the remaining point and p. Then the vote value of p is calculated for the remaining point through the gaussian function as the following formula shown [19].

    Vp(q)=12πexp((Dp(q)/μp)22)p,qT (4)
    μp=median(Dp(q))       qT (5)

    Dp(q) is the distance between p and point q. μp is a ratio coefficient which scales Dp(q) into a smaller range. Vp(q) is the vote value of p on q. The reason why we set the μp as the median value of Dp(q) is that μp is inversely proportional to the vote accumulated value difference of adjacent point. Tiny difference may result in a mistaken choice of the initial point. We calculate the vote accumulated value of each point of T and assign it to a matrix M. The size of M is same as T.

    The final matrix M stores vote accumulated value of each point of T. Find the max value of matrix M and record its position. The point of the same position in T is the final initial point I. This method is pretty stable. It is still able to give a reliable non-outlier in the case of heavy outlier contamination. Its performance is shown in the follow Figure 4.

    Figure 4.  The determined initial point under different situation (red point). (a) The situation with little outlier; (b) the situation with many outlier.

    The picture tells that this method performs well even if the voting area is heavily polluted by outlier. All points are participated in the voting process. Thus, the final initial point is “fair” for each point. As the outlier is always the minority, their accumulated vote value is always smaller than the non-outlier. In other words, the point which has the max voting accumulated value is always a non-outlier. Therefore our method is able to give a reliable initial point.

    The social circle radius limits the rang of social circle. A larger social circle will increase the possibility that outlier is mistaken for non-outlier. So, it is important to determine a proper social circle radius. The size of search window needs to be set first, the size of search window is inversely proportional to time efficiency. Generally, it is suggested to set a 3×3 window. As the initial point I and the search window size have been determined, the window is expanded in order to get a reliable social circle radius.

    As is shown in Figure 5, we expand the window 3×3 and build a new window 5×5, whose radius is the double of original window size. Then, the distance between the initial point I and the rest point is calculated out in the new window. The distance value is assigned to an array B. We sort the array B. The final social circle radius can be obtained through the follow formula.

    B(n)=(xIxn)2+(yIyn)2+(zIzn)2n=0,1,2...15 (6)
    R=sort(B(n))8n=0,1,2,....15 (7)
    Figure 5.  The expansion of window.

    R is the social circle radius. sort(B(n))8 is the element value at index 8 of sorted array B. The purpose of expand the window is to restrain the outlier disturbance. The outlier may lead to a too large social circle radius. A large radius will weaken the outlier removal performance. We sort the array B and the large distance calculated from outlier will locate in the tail of B. In this way, we exclude the outlier disturbance.

    The expansion process of social circle is actually the overlay of social circle. Every point has its own social circle. We combine the social circle of each point and obtain a big group. In the group, each point can get in touch with other points indirectly or directly [17].

    We assume that the input point cloud P, initial point I, search window W and the search radius R is determined, the points within the search radius are used as the initial points in the next iteration, we select the point on the edge of the search window as the initial point to improve calculating effciency. By analogy, continue to iterate the initial point until the final social circle store points which is all connected to the initial point. The algorithm flow is shown in Figure 6.

    Figure 6.  Iteratively expand the social circle algorithm process.

    Operation A: Get the initial point, social circle radius R, search window, create containers V1, V2 for storing three-dimensional points, matrix B with value filled zero which is same size as point cloud P, What ’s more add the initial points to V1.

    Condition A: whether V1 is empty.

    Operation B: Find the point whose value in the corresponding matrix B is not 0 in the search window.

    Operation C: Find the point where the distance between the initial point and the current point is less than the distance from the social circle, and the height of the current point is assigned to the corresponding position in B.

    Condition B: whether the current point is at the search window boundary.

    Operation D: put the point in V2.

    Condition C: whether all points in the search window have been visited.

    Operation E: clear V1, add the elements in V2 to V1, and go back to condition A.

    Operation F: output the matrix B.

    The expansion can be divided into seven steps:

    (1) Determine the initial point I of the organized point P, the social circle radius R and search window W.

    (2) Create containers V1 and V2 to store three-dimensional points, and add the initial points to V1, what’s more create a matrix B with value filled zero which is same size as P.

    (3) Judge whether V1 is empty, if it is, output the matrix B, otherwise take each point in V1 as the initial point and execute steps (4)–(8) cyclically.

    (4) For the currently visited point in the search window, determine whether the value of the corresponding position element in B is 0, if it is, execute the next step, otherwise visit the next point.

    (5) Calculate the distance between the initial point and the access point. If it is less than the radius of the social circle, the height of the access point is assigned to the corresponding position in B, and perform the next step. Otherwise, visit the next point of the window.

    (6) Determine whether the current access point is at the window boundary, if so, put the point in V2.

    (7) Judge whether all points in the search window have been visited, if yes, clear V1, add the elements in V2 to V1, and execute step (3), otherwise visit the next point.

    Through the above social circle iteration process, matrix B is obtained. The points in A corresponding to values greater than 0 in B are non-outliers. Convert all values greater than 0 in B to three-dimensional points and output the processed points. Finally we get point cloud without outliers.

    The principle of our method has been explained in the section 3, we will perform some experiments and discuss the experiment results in this section. A point cloud of metal groove with high outlier is selected as original data. The outlier is non-uniform distribution. We compare our method with another two frequently-used outlier removal methods. They are statistical filter and radius filter which come from Point Cloud Library. They can also remove outlier fast and will not modify the non-outlier. The comparison is conducted in the aspect of time consuming and removal performance. The raw point cloud is shown in Figure 7.

    Figure 7.  The raw point cloud.

    As shown in Figure 7, the degree of outlier contamination differs in different area. The groove area is the most contaminated area. We reimplement the statistical and radius filter so that they are suitable for processing the organized point cloud. Their search window sizes are all set as 5×5(points). We use different combination of parameter to conduct experiment for each method and take the relatively good result as the final output data (Figure 8).

    Figure 8.  The outlier removal result. (a) Social circle filter; (b) Radius filter; (c) Statistical filter.

    As shown in Figure 8, they have similar performance in the low outlier contaminated area. However, their performance differs a lot in the high outlier contaminated area. The radius filter distinguishes the single outlier well but fails to detect dense outlier cluster. The statistical filter performs better than radius filter, but is hard to detect dense outlier cluster. It is due to the fact that they consider the low local point density and rely on the condition that the local neighborhood of each point contains sufficient data points. So, they cannot remove the dense outlier clusters. Compared with the radius and statistical filter, the social circle filter can not only remove single outlier but also dense outlier cluster. For the high outlier contaminated area, the social circle filter can remove outlier efficiently while the radius and statistical filter perform poorly. Their difference is shown in Table 1 which is corresponding to Figure 8.

    Table 1.  The outlier removal comparison of three methods.
    Method Number of raw 3D points Number of processed 3D points Time (ms)
    Social circle filter 1228793 1111450 148
    Radius filter 1228793 1012853 1134
    Statistical filter 1228793 1019601 474

     | Show Table
    DownLoad: CSV

    As shown in Table 1 and Figure 8, the proposed method keeps more points when outlier is removed. However, the radius and statistical filter keep fewer points while outlier still exists. Besides, the proposed method consumes less time than other two methods.

    The previous experiment intercepted part of the point cloud data for comparison, and then we filtered the entire cloud, and interpolated the filtered part to get final point cloud. The point cloud interpolation method selected in this article unified the use of cubic b-spline curve for patching. The point cloud obtained by the method is not affected by the filtering method. As before, we still compare with statistical removal and radius outlier removal, and we compare the two different occasions, one is a metal groove, the other is a composite material. The raw point cloud is shown as Figure 9.

    Figure 9.  The raw point cloud (a) Metal groove; (b) Composite material.

    As shown in Figure 9, both metal grooves and composite material have some outlier. We also use the three filtering methods, and use the same repair method at the same time, and the final output data are shown in the Figures 10 and 11.

    Figure 10.  Metal groove point cloud. (a) Social circle filter; (b) Radius filter; (c) Statistical filter.
    Figure 11.  Composite material point cloud. (a) Social circle filter; (b) Radius filter; (c) Statistical filter.

    It can be seen from Figures 10 and 11 that regardless of whether the outliers are dense or sparse, the quality of the final point cloud data reconstructed by the method of social circle filter is better than that of the traditional methods radius filter and statistical filter, especially when the outliers are dense.

    This paper proposes an outlier removal method which is applied for organized point cloud reconstructed by microscopic measurement. With comparison, we found that the method of filtering through social circles can quickly and accurately filter out some outliers without changing the non-outliers, and it is also very robust. The core idea of our method is derived from social circle. It includes three steps. First, an area is selected to determine an initial point. Then, the social circle radius needs to be calculated out. At last, the social circle is iteratively expanded and remove the outlier. Compared with the radius and statistical filter, the experiment result shows that the method proposed in the paper is able to remove outlier robustly and consumes less time.

    This work was supported by the National Key R & D program of China (2019YFB2006100, 2019YFB1707500), Fundamental Research Funds for the Central Universities (NO. NS2020030), Jiangsu Province Nature Science Fund (BK20191280).

    The authors declare there is no conflict of interest.

    [1] EU (2009) DIRECTIVE 2009/28/EC OF THE EUROPEAN PARLIAMENT AND OF THE COUNCIL of 23 April 2009 on the promotion of the use of energy from renewable sources and amending and subsequently repealing Directives 2001/77/EC and 2003/30/EC. Official Journal of the European Union.
    [2] Eurostat (2017) Energy from renewable sources. Available from: http://ec.europa.eu/eurostat/statistics-explained/index.php/Energy_from_renewable_sources.
    [3] EurObservÉR (2016) Biofuels Barometer 2016. EurObservÉR.
    [4] Fritsche UR, Wiegmann K (2011) Indirect land use change and biofuels. Committee on Environment, Public Health and Food Safety, European Parliament, Brussels.
    [5] EU (2015) DIRECTIVE (EU) 2015/1513 OF THE EUROPEAN PARLIAMENT AND OF THE COUNCIL of 9 September 2015 amending Directive 98/70/EC relating to the quality of petrol and diesel fuels and amending Directive 2009/28/EC on the promotion of the use of energy from renewable sources. Official Journal of the European Union.
    [6] SEA (2016) Transportsektorns energianvändning 2015 [Use of energy in the transportation sector 2015]. The Swedish Energy Agency. ES 2016: 03.
    [7] SEA (2016) Produktion och användning av biogas och rötrester 2015 [Production and utilization of biogas and digestate 2015]. The Swedish Energy Agency.
    [8] EurObservÉR (2016) The state of renewable energies in Europe-Edition 2016, 16th EurObservÉR Report. EurObservÉR.
    [9] FNR (2016) Bioenergy in Germany-Facts and Figures 2016. Fachagentur Nachwachsende Rohstoffe e.V.
    [10] Grahn M, Hansson J (2015) Prospects for domestic biofuels for transport in Sweden 2030 based on current production and future plans. Wires Energy Environ 4: 290–306. doi: 10.1002/wene.138
    [11] STA (2016) Styrmedel och åtgärder för att minska transportsystemets utsläpp av växthusgaser–med fokus på transportinfrastrukturen [Policy instruments and measures to reduce emissions of GHG from the transport system-with focus on transport infrastructure]. The Swedish Transport Administration.
    [12] STA (2016) Åtgärder för att minska transportsektorns utsläpp av växthusgaser–ett regeringsuppdrag [Measures to reduce GHG emissions from the transport sector-a government assignment]. The Swedish Transport Administration.
    [13] SOU (2016) En klimat-och luftvårdsstrategi för Sverige [A climate and air pollution strategi] Swedish Government Official Report.
    [14] Gissén C, Prade T, Kreuger E, et al. (2014) Comparing energy crops for biogas production–yields, energy input and costs in cultivation using digestate and mineral fertilisation. Biomass Bioenerg 64: 199–210. doi: 10.1016/j.biombioe.2014.03.061
    [15] Börjesson P, Prade T, Lantz M, et al. (2015) Energy crop-based biogas as vehicle fuel: the impact of crop selection on energy efficiency and greenhouse gas performance. Energies 8: 6033–6058. doi: 10.3390/en8066033
    [16] FNR (2010) Biogas Messprogramm II-61 Biogasanlagen im Vergleich. Gulzow: Bundesministerium fur Ernährung, Landwirtshafts und Verbraucherschutz.
    [17] Hjort-Gregersen K (2015) Udvikling og effektivisering af biogasproduktionen i Danmark - Ökonomi, nögletal og benchmark, Energistyrelsens Biogas Taskforce. AgroTech.
    [18] Möller HB, Nielsen KJ (2015) Udvikling og effektivisering af biogasproduktionen i Danmark, Faglig rapport, Biogas Taskforce. Inst. for Ingeniörvidenskab, Aarhus Universitet.
    [19] SJV (2016) Rötning av animaliska biprodukter, 2016-10-24 Jordbruksverket.
    [20] Lantz M, Ekman A, Börjesson P (2009) Systemoptimerad produktion av fordonsgas, Report 69. Lund, Sweden: Department of Technology and Society, Lund University. 69. 110 p.
    [21] Petersson A, Wellinger A (2009) Biogas upgrading technologies–development and innovations. IEA Bioenergy Task 37, European Commission Joint Research Centre, Petten, The Netherlands.
    [22] IEA (2016) Biogas upgrading plant list, data up to the end of 2015. IEA Bioenergy Task 37, European Commission Joint Research Centre, Petten, The Netherlands.
    [23] Bauer F, Hultenberg C, Persson T, et al. (2013) Biogas upgrading–Review of commercial technologies. Malmö: Svenskt Gastekniskt Center. Rapport 2013: 270.
    [24] Hoyer K, Hultenberg C, Svensson M, et al. (2016) Biogas Upgrading–Technical Review. Energiforsk.
    [25] Urban W, Girod K, Lohmann H (2008) Technologien und Kosten der Biogasaufbereitung und Einspeisung in das Erdgasnetz. Ergebnisse der Markterhebung 2007–2008. Fraunhofer-Institut für Umwelt- Sicherheits- und Energietechnik, Oberhausen, Germany.
    [26] Berglund P, Bohman M, Svensson M, et al. (2012) Teknisk och ekonomisk utvärdering av lantbruksbaserad fordonsgasproduktion. Swedish Gas Technology Centre.
    [27] Colnerud GS, Gåverud H, Glimhall A (2010) Förändrade marknadsvillkor för biogasproduktion. The Energy Market Inspectorate.
    [28] Walla C, Schneeberger W (2008) The optimal size for biogas plants. Biomass Bioenerg 32: 551–557. doi: 10.1016/j.biombioe.2007.11.009
    [29] Börjesson P, Lantz M, Andersson J, et al. (2016) METHANE AS VEHICLE FUEL–A WELL-TO-WHEEL ANALYSIS (METDRIV), report 2016:06 f3 The Swedish Knowledge Centre for Renewable Transportation Fuels, Sweden.
    [30] Berglund M, Börjesson P (2006) Assesment of energy performance in the life-cycle of biogas production. Biomass Bioenerg 30: 254–266.
    [31] Brown J, Nizami AS, Thamsiriroj T, et al. (2011) Assessing the cost of biofuel production with increasing penetration of the transport fuel market: A case study of gaseous biomethane in Ireland. Renew Sust Energ Rev 15: 4537–4547. doi: 10.1016/j.rser.2011.07.098
    [32] Jury C, Benetto E, Koster D, et al. (2010) Life cycle assessment of biogas production by monofermentation of energy crops and injection into the natural gas grid. Biomass Bioenerg 34: 54–66.
    [33] FNR (2006) Handreichung Biogasgewinnung und nutzung. Fachagentur Nachwachsende Rohstoffe e.V., Bundesministerium für Ernährung, Landwirtschaft and Verbraucherschutz, Gülzow-Prüzen, German.
    [34] Murphy J, Braun R, Weiland P, et al. (2011) Biogas from Crop digestion. Available from: https://www.nachhaltigwirtschaften.at/resources/iea_pdf/reports/iea_bioenergy_task37_biogas_from_crop_digestion.pdf.
    [35] FNR (2005) Ergebnisse des Biogas-Messprogram. Fachagentur Nachwachsende Rohstoffe e.V., Bundesministerium für Ernährung, Landwirtschaft and Verbraucherschutz, Gülzow-Prüzen, Germany.
    [36] Lundqvist P (2009) Fjärrvärme för utökad biogasproduktion–Teknisk och ekonomisk utvärdering av fjärrvärme för uppvärmning av biogasprocesser. Thermal Engineering Research Institute, Stockholm, Sweden.
    [37] Norin E (2007) Alternativa hygieniseringsmetoder. Swedish Gas Technology Centre, Malmö, Sweden.
    [38] Rosenqvist H (2017) Kalkyler för energigrödor 2017-fastbränsle, biogas, spannmål och raps. Jordbruksverket.
    [39] Björnsson L, Prade T, Lantz M (2016) Grass for biogas-Arable land as carbon sink, Report 2016:280. Energiforsk.
    [40] Valin H, Peters D, Berg Mvd, et al. (2015) The land use change impact of biofuels consumed in the EU - Quantification of area and greenhouse gas impacts. ECOFYS.
    [41] Kreuger E, Prade T, Björnsson L, et al. (2014) Biogas från Skånsk betblast - potential, teknik och ekonomi [Biogas from beet tops in Scania-potential, technology and economy]. Environmental and Energy System Studies, Lund University.
    [42] SCB (2017) Price Indices and Prices in the Food Sector-Annual and Monthly Statistics-2017:06. Statistics Sweden. JO 49 SM 1708.
    [43] Energimyndigheten (2017) Trädbränsle- och torvpriser Nr 3 2017. Swedish Statistics.
    [44] Hansen MT (2015) Standardforutsaetninger-til VE til process-ansögningsmateriale. FORCE Technology.
    [45] Odhner PB, Svensson SE, Prade T (2015) Extruder för ökad biogasproduktion [Extrusion increases the production of biogas], Report 2015: 26 Sveriges lantbruksuniversitet, Fakulteten för landskapsarkitektur, trädgårds- och växtproduktionsvetenskap.
    [46] FNR (2010) Guide for biogas-From production to use Nachwachsende Rohstoffe e.V., Bundesministerium für Ernährung, Landwirtschaft and Verbraucherschutz, Gülzow-Prüzen, Germany.
    [47] FNR (2010) Leitfaden Biogas–Von der Gewinnung zur Nutzung. Nachwachsende Rohstoffe e.V. Gülzow-Prüzen, Germany: Bundesministerium für Ernährung, Landwirtschaft and Verbraucherschutz.
    [48] FNR (2016) Leitfaden Biogas–Von der Gewinnung zur Nutzung. Nachwachsende Rohstoffe e.V. Gülzow-Prüzen, Germany: Bundesministerium für Ernährung, Landwirtschaft and Verbraucherschutz.
    [49] Lantz M, Björnsson L (2011) Biogas från gödsel och vall-analys av föreslagna styrmedel. The Federation of Swedish Farmers, Stockholm, Sweden: LRF.
    [50] Berglund P (2010) Biogödselhandbok–Biogödsel från storskaliga biogasanläggningar. Swedish Waste Management.
    [51] SEA (2017) Trädbränsle och torvpriser [Wood fuel and peat prices], Nr 1 2017. The Swedish Energy Agency.
    [52] SCB (2016) Priser på el för industrikunder [Electricity price for industries]. Available from: http://www.scb.se/hitta-statistik/statistik-efter-amne/energi/prisutvecklingen-inom-energiomradet/energipriser-pa-naturgas-och-el/pong/tabell-och-diagram/genomsnittspriser-per-halvar-2007/priser-pa-el-for-industrikunder-2007/, accessed 2017-05-23.
    [53] Smyth B, H S, Murphy J (2010) Can grass biomethane be an economically viable biofuel for the farmer and the consumer? Biofuel Bioprod Bio 4: 519–537. doi: 10.1002/bbb.238
    [54] Lantz M (2013) Biogas in Sweden-opportunities and challenges from a systems perspective. Lund, Sweden: Lund University.
    [55] SJV (2016) Rekommendationer för gödsling och kalkning 2017, Report JO16:24. The Swedish Board of Agriculture. Jönköping, Sweden.
    [56] Nutrients Fo (2011) Stallgödselkalkyl (Manure calculator) version 2011-03-31. accessed 2012-10-18.
    [57] Lantz M, Svensson M, Björnsson L, et al. (2007) The prospects for an expansion of biogas systems in Sweden–Incentives, barriers and potentials. Energy Policy 35: 1830–1843.
    [58] EC (2015) Statligt stöd-Skattebefrielser och skattenedsättningar för flytande drivmedel, C(2015) 9344. European Commission.
    [59] EC (2015) Statligt stöd-Skattebefrielser för biogas som används som motorbränsle, C(2015) 9345. European Commission.
    [60] Regeringen (2017) Reduktion av växthusgasutsläpp genom inblandning av biodrivmedel i bensin och dieselbränslen, Lagrådsremiss In: energidepartementet M-o, editor: The government of Sweden.
    [61] SJV (2017) Gödselgasstöd. Available from: http://www.jordbruksverket.se/amnesomraden/ stod/andrastod/godselgasstod.4.ac526c214a28250ac23333e.html, accessed 20170613.
    [62] SJV (2017) Investeringsstöd till biogas. Available from: http://www.jordbruksverket.se/amnesomraden/stod/stodilandsbygdsprogrammet/investeringar/biogas.4.6ae223614dda2c3dbc44f95.html, accessed.
    [63] EC (2014) COMMISSION REGULATION (EU) No 651/2014 of 17 June 2014 declaring certain categories of aid compatible with the internal market in application of Articles 107 and 108 of the Treaty In: COMMISSION TE, editor. Official Journal of the European Union.
    [64] SPBI (2017) Priser och skatter. Available from: http://spbi.se/statistik/priser/, accessed.
    [65] Karlsson S, Rodhe L (2002) Översyn av Statistiska Centralbyråns beräkning av ammoniakavgången i jordbruket–emissionsfaktorer för ammoniak vid lagring och spridning av stallgödsel.
    [66] Bengtsson B, Rasic Z (2005) Kadmium i odlingssystem med tillförsel av rötslam. Jord-skörd- och markvatten-analyser. Department of plant sciences, Swedish University of Agricultural Sciences, Alnarp, Sweden.
    [67] Roth U, Wulf S (2010) Gasausbeute in landwirtschaftlichen Biogasanlagen, KTBL-Heft 88. Kuratorium Fur Technik und Bauwesen in der Landwirtschaft e.V. (KTBL), Darmstadt, Germany.
    [68] ECN (2014) The database Phyllis2: Biomass and waste. Available from: https://www.ecn.nl/phyllis2, accessed 2017-04-03.
    [69] Tamaki Y, Mazza G (2010) Measurement of structural carbohydrates, lignins, and micro-components of straw and shives: Effects of extractives, particle size and crop species. Ind Crop Prod 31: 534–541. doi: 10.1016/j.indcrop.2010.02.004
    [70] Björnsson L, Castillo MdP, Gunnarsson C, et al. (2014) Förbehandling av lignocellulosarika råvaror för biogasproduktion-Nyckelaspekter vid jämförande utvärdering. Lund, Sweden: Environmental and Energy Systems Studies.
    [71] Moller HB, Sommer SG, Ahring B (2004) Methane productivity of manure, straw and solid fractions of manure. Biomass Bioenerg 26: 485–495. doi: 10.1016/j.biombioe.2003.08.008
    [72] Amon T, Amon B, Kryvoruchko V, et al. (2007) Biogas production from maize and dairy cattle manure-Influence of biomass composition on the methane yield. Agr Ecosyst Enviro 118: 173–182.
    [73] Mähnert P, Linke B (2009) Kinetic study of biogas production from energy crops and animal waste slurry: Effect of organic loading rate and reactor size. Environ Technol 30: 93–99.
    [74] McCarty PL (1964) Anaerobic waste treatment fundamentals, Part one, Chemistry and microbiology. Public Works 95: 107–112.
    [75] FNR (2010) Biogas-Messprogramm II, 61 Biogasanlagen im Vergleich. Fachagentur Nachwachsende Rohstoffe e.V., Bundesministerium für Ernährung, Landwirtschaft and Verbraucherschutz, Gülzow-Prüzen, Germany.
    [76] Ljung E, Palm O, Rodhe L (2013) Ökad acceptans för biogödsel inom lantbruket. Uppsala, Sweden: JTI-Institutet för jordbruks-och miljöteknik.
    [77] Procházka J, Dolejš P, Máca J, et al. (2012) Stability and inhibition of anaerobic processes caused by insufficiency or excess of ammonia nitrogen. Appl Microbiol Biot 93: 439–447. doi: 10.1007/s00253-011-3625-4
    [78] Schnürer A, Nordberg Å (2008) Ammonia, a selective agent for methane production by syntrophic acetate oxidation at mesophilic temperature. Water Sci Technol 57: 735–740.
    [79] Scherer P, Lippert H, Wolff G (1983) Composition of the major elements and trace elements of 10 methanogenic bacteria determined by inductively coupled plasma emission spectrometry Biol Trace Elem Res 5: 149–163.
    [80] Banks CJ, Zhang Y, Jiang Y, et al. (2012) Trace element required for stable food waste digestion at elevated ammonia concentrations. Bioresource Technol 104: 127–135. doi: 10.1016/j.biortech.2011.10.068
    [81] Overend R (1982) Haul distance and transportation work factors for biomass. Biomass 2: 75–79. doi: 10.1016/0144-4565(82)90008-7
    [82] Börjesson P, Gustavsson L (1996) Regional production and utilization of biomass in Sweden. Energy 21: 747–764. doi: 10.1016/0360-5442(96)00029-1
    [83] Björnsson L, Lantz M, Murto M, et al. (2011) Biogaspotential i Skåne-inventering och planeringsunderlag på översiktsnivå [The biogas potential in Scania-inventory and data for planning]. The County Administrative Board of Skåne.
  • Reader Comments
  • © 2017 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6590) PDF downloads(1326) Cited by(7)

Figures and Tables

Figures(7)  /  Tables(8)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog